Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation software.”

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 289 records · Page 16

PARAMETRIC STUDY OF GROUND SOURCE HEAT PUMP SYSTEM FOR HOT AND HUMID CLMATE

The U-tube sizes and varied thermal conductivity with different grout materials are studied based on the benchmark residential building in Hot-humid Pensacola, Florida. In this study, the benchmark building is metered and the data is used to validate the simulation model. And a list of comparative simulation cases with varied parameter value are simulated to study the importance of pipe size and grout to the ground source heat pump energy consumption. The simulation software TRNSYS [1] is employed to fulfill this task. The results show the preliminary energy saving based on varied parameters. Future work needs to be conducted for the cost analysis, include the installation cost from contractor and materials cost.

15 GEOTHERMAL ENERGY↗

Simulations of Bunch Merging in a Beta Beam Decay Ring

To further study neutrino oscillation properties a Beta Beam facility has been proposed. Beta decaying ions with high kinetic energy are stored in a storage ring (''Decay Ring'') with straight sections to create pure focused (anti) electron neutrino beams. However to reach high sensitivity to neutrino oscillation parameters in the experiment the bunched beam intensity and duty cycle in the DR have to be optimized. The first CERN-based scenario, using {sup 6}He and {sup 18}Ne as neutrino sources, has been studied using a bunch merging RF scheme. Two RF cavities at different frequencies are used to capture newly injected bunches and then merge them into the stored bunches. It was shown that this scheme could satisfy the requirements on intensity and duty cycle set by the experiment. This merging scheme has now been revised with new simulation software providing new results for {sup 6}He and {sup 18}Ne. Furthermore bunch merging has been studied for the second CERN-based scenario using {sup 8}Li and {sup 8}B.

43 PARTICLE ACCELERATORS↗

2011 EnergyValue Housing Award Report

The EnergyValue Housing Award (EVHA), an integral part of the U.S. Department of Energy’s Building America (BA) program, seeks to improve and expand the energy efficiency of residential construction nationwide by honoring builders and remodelers and awarding highly efficient homes. Receiving EVHA recognition helps builders and remodelers communicate the value of energy efficiency to prospective customers and distinguish their homes from those built by more traditional residential construction companies. It also provides opportunities for EVHA winning builders and remodelers to set an example for improving the quality and performance of homes throughout the country. Each year, the project scope of the EVHA includes the application and judging processes, developing the EVHA publication, updating the program and application for the next year, hosting the EVHA banquet, and outreach and technical support. Following the award banquet, the NAHB Research Center evaluates the energy performance of that year’s EVHA winning houses. This work is primarily comprised of estimating the energy savings of each home/ or housing unit over the BA Benchmark using energy simulation software and procedures developed by the BA program. The energy evaluations serve as a quantitative measure of the energy efficiency of the units which the judges decided merited an award. As part of this work, the NAHB Research Center Industry Partnership (NAHBRC-IP) for High Performing Homes strives to make new and existing homes more energy efficient, durable, and affordable through research and development as one of the U.S. Department of Energy’s BA teams. Each year, through the EVHA, builder and remodeler partners work with industry partners in a variety of climates endeavor to reach the energy efficiency targets of the BA program. (Department of Energy 2012) This report details the simulation tool(s) and energy modeling methodology followed in making the energy efficiency estimates and documents the estimated performance of the EVHA award winning houses in comparison with the BA Benchmark and the associated House Simulation Protocols (HSPs). A summary of each building and its features is included with a brief description of the project and the judges’ comments. The purpose of this report is to assess the energy performance of the 2011 EVHA winners and to align the EVHA program with the BA program.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analysis of Factors that Influence Infiltration Rates using the HELP Model

The Hydrologic Evaluation of Landfill Performance (HELP) model is used by Savannah River National Laboratory (SRNL) in conjunction with PORFLOW groundwater flow simulation software to make longterm predictions of the fate and transport of radionuclides in the environment at radiological waste sites. The work summarized in this report supports preparation of the planned 2018 Performance Assessment for the E-Area Low-Level Waste Facility (LLWF) at the Savannah River Site (SRS). More specifically, this project focused on conducting a sensitivity analysis of infiltration (i.e., the rate at which water travels vertically in soil) through the proposed E-Area LLWF closure cap. A sensitivity analysis was completed using HELP v3.95D to identify the cap design and material property parameters that most impact infiltration rates through the proposed closure cap for a 10,000-year simulation period. The results of the sensitivity analysis indicate that saturated hydraulic conductivity (Ksat) for select cap layers, precipitation rate, surface vegetation type, and geomembrane layer defect density are dominant factors limiting infiltration rate. Interestingly, calculated infiltration rates were substantially influenced by changes in the saturated hydraulic conductivity of the Upper Foundation and Lateral Drainage layers. For example, an order-of-magnitude decrease in Ksat for the Upper Foundation layer lowered the maximum infiltration rate from a base-case 11 inches per year to only two inches per year. Conversely, an order-of-magnitude increase in Ksat led to an increase in infiltration rate from 11 to 15 inches per year. This work and its results provide a framework for quantifying uncertainty in the radionuclide transport and dose models for the planned 2018 E-Area Performance Assessment. Future work will focus on the development of a nonlinear regression model for infiltration rate using Minitab 17® to facilitate execution of probabilistic simulations in the GoldSim® overall system model for the E-Area LLWF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Verification and validation of the Alternative Nonlinear Two-phase Subchannel (ANTS) code

The Alternative Nonlinear Two-phase Subchannel solver (ANTS) code was written to provide a fast-running, steady-state, pin-resolved modeling and simulation tool for analysis of common boiling water reactor (BWR) geometry and common operating conditions. ANTS has been integrated into the Virtual Environment for Reactor Application (VERA) core simulator software, where it can be used to provide a thermal/hydraulic (T/H) subchannel solution that is then used to provide neutronic feedback as well as perform the fuel depletion and temperature solution. Herein, this paper presents the rigorous analysis performed on the ANTS code, which includes both code and solution verification testing, benchmarking with the existing two-phase subchannel capability in VERA, CTF, and validation testing using popular two-phase experiments such as PWR Sub-channel and Bundle Tests (PSBT), BWR Full-size Fine-mesh Bundle Tests (BFBT), Risø, and FRIGG. This assessment was used to qualify ANTS for its intended applications before its use for core-scale, multiphysics BWR simulations. In general, it was found that agreement with experimental data was good; errors were within the range of experimental data uncertainty. Furthermore, code and solution verification confirmed that the governing equations and the most important closure terms were correctly implemented and behaving as expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On the feasibility of polyurethane based 3D dosimeters with optical CT for dosimetric verification of low energy photon brachytherapy seeds

Purpose: To investigate the feasibility of and challenges yet to be addressed to measure dose from low energy (effective energy <50 keV) brachytherapy sources (Pd-103, Cs-131, and I-125) using polyurethane based 3D dosimeters with optical CT. Methods: The authors' evaluation used the following sources: models 200 (Pd-103), CS-1 Rev2 (Cs-131), and 6711 (I-125). The authors used the Monte Carlo radiation transport code MCNP5, simulations with the ScanSim optical tomography simulation software, and experimental measurements with PRESAGE{sup ®} dosimeters/optical CT to investigate the following: (1) the water equivalency of conventional (density = 1.065 g/cm{sup 3}) and deformable (density = 1.02 g/cm{sup 3}) formulations of polyurethane dosimeters, (2) the scatter conditions necessary to achieve accurate dosimetry for low energy photon seeds, (3) the change in photon energy spectrum within the dosimeter as a function of distance from the source in order to determine potential energy sensitivity effects, (4) the optimal delivered dose to balance optical transmission (per projection) with signal to noise ratio in the reconstructed dose distribution, and (5) the magnitude and characteristics of artifacts due to the presence of a channel in the dosimeter. Monte Carlo simulations were performed using both conventional and deformable dosimeter formulations. For verification, 2.8 Gy at 1 cm was delivered in 92 h using an I-125 source to a PRESAGE{sup ®} dosimeter with conventional formulation and a central channel with 0.0425 cm radius for source placement. The dose distribution was reconstructed with 0.02 and 0.04 cm{sup 3} voxel size using the Duke midsized optical CT scanner (DMOS). Results: While the conventional formulation overattenuates dose from all three sources compared to water, the current deformable formulation has nearly water equivalent attenuation properties for Cs-131 and I-125, while underattenuating for Pd-103. The energy spectrum of each source is relatively stable within the first 5 cm especially for I-125. The inherent assumption of radial symmetry in the TG43 geometry leads to a linear increase in sample points within the 3D dosimeter as a function of distance away from the source, which partially offsets the decreasing signal. Simulations of dose reconstruction using optical CT showed the feasibility of reconstructing dose out to a radius of 10 cm without saturating projection images using an optimal dose and high dynamic range scanning; the simulations also predicted that reconstruction artifacts at the channel surface due to a small discrepancy in refractive index should be negligible. Agreement of the measured with calculated radial dose function for I-125 was within 5% between 0.3 and 2.5 cm from the source, and the median difference of measured from calculated anisotropy function was within 5% between 0.3 and 2.0 cm from the source. Conclusions: 3D dosimetry using polyurethane dosimeters with optical CT looks to be a promising application to verify dosimetric distributions surrounding low energy brachytherapy sources.

07 ISOTOPE AND RADIATION SOURCES↗

GoldSim® Aquifer Model Calibration and Plume Interaction

The E-Area Low-Level Waste Facility (ELLWF) Probabilistic Aquifer Model (PAM) utilizes GoldSim ® Monte Carlo simulation software (GTG, 2017) to evaluate the transport of a tracer species as it travels from the water table below the disposal unit footprint, through the aquifer, to the Point of Assessment (POA) at the 100-meter boundary. This report documents the development and calibration of PAM as well as the implementation of plume interaction. PAM is a key component of the effort to include uncertainty quantification and sensitivity analysis (UQSA) in the next revision of the E-Area Performance Assessment (PA), considering recommendations from the 2015 PA strategic planning team outlined by Butcher and Phifer (2016). The Aquifer Model and associated optimized geometric parameters will be implemented in the future GoldSim® (GS) system model that will simulate subsurface flow and radionuclide transport from the ground surface to the 100-meter POA.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

High-Burnup BWR LOCA Burst Analysis Using High-Fidelity Multiphysics Simulations

The US nuclear industry is looking to improve on the operating economics of the current fleet of light-water reactors (LWRs). One way of achieving this is by operating fuel to higher burnup. In pressurized water reactors (PWRs), relaxing the current burnup limit will allow for cycle length extensions and power uprates; in boiling water reactors (BWRs) it may allow for improved fuel utilization and reduced feed assemblies, as well as more efficient power uprates and increased capacity factors that will support the Administration’s Executive Order to facilitate 5 GW of power uprates at existing nuclear facilities. However, one of the key limitations to operating fuel to higher burnup is the risk of fuel fragmentation, relocation, and dispersal (FFRD). Recognizing the high interest in extending burnup limits, the US Nuclear Regulatory Commission (NRC) has issued Draft Regulatory Guide DG-1434, which defines an approach that would be acceptable to the NRC for addressing FFRD risk. The approach defined will require better understanding of the phenomena leading to FFRD as well as best-estimate simulation methods to understand FFRD risk in high-burnup cores. The Nuclear Energy Advanced Modeling and Simulation program is supporting the FFRD industry challenge problem through development of state-of-the-art, high-fidelity modeling and simulation LWR analysis capabilities; namely, the BISON fuel performance code and the VERA core simulator software. These tools, along with the US NRC TRACE system analysis code, have been utilized for analysis of FFRD risk in both PWR and BWR cores in recent years. The work documented in this report addresses the lack of high-fidelity research for BWRs and builds on a previous activity where the framework has been applied to Cycles 16 through 18 of Limerick Unit 1, a BWR/4, with introduction of 8 high-burnup lead use assemblies (HBLUAs) that were representative of the 8 HBLUAs loaded into Limerick Unit 2 in 2021. VERA was used in this previous activity to model rod-by-rod depletion in these cycles, and its solution was used to initialize a TRACE simulation of a large-break loss-of-coolant accident (LBLOCA) at the end of Cycle 18. In the work documented in this report, the TRACE model was improved by refining the core mesh and utilizing a new feature that allows for capturing the full 3D VERA power distribution in the model. This allows for a more detailed solution for setting BISON boundary conditions. Furthermore, the solutions from VERA and TRACE were used to set up and perform BISON simulations of about 1,000 rods sampled from the core, including all burnup levels. Utilizing two cladding burst models, it was shown that no fuel rods were predicted to burst during the postulated LBLOCA transient. Additionally, a sensitivity study was performed by artificially increasing linear heat rate during the postulated LBLOCA to identify parameters that correlate with rod burst susceptibility. Burnup, fission gas release, and hoop strain were all found to be positively correlated with rod burst susceptibility. Small-break loss-of-coolant accident (SBLOCA) analyses were also performed; these analyses predicted cladding temperature increases that were bounded by the LBLOCA cladding temperatures for all small break sizes studied for this plant. However, future refinements to the plant response assumptions during the SBLOCA could impact the predicted cladding response. Finally, a benchmark study was performed between CTF and TRACE for LOCA conditions to better qualify CTF for BWR LOCA modeling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Atmospheric release model for the E-area low-level waste facility: Updates and modifications

The atmospheric release model (ARM) utilizes GoldSim® Monte Carlo simulation software (GTG, 2017) to evaluate the flux of gaseous radionuclides as they volatilize from E-Area disposal facility waste zones, diffuse into the air-filled soil pores surrounding the waste, and emanate at the land surface. This report documents the updates and modifications to the ARM for the next planned E-Area PA considering recommendations from the 2015 PA strategic planning team outlined by Butcher and Phifer.

54 ENVIRONMENTAL SCIENCES↗

Reflections of ions in electrostatic analyzers: A case study with New Horizons/Solar Wind Around Pluto

Electrostatic analyzers (ESAs), in various forms, are used to measure plasma in a range of applications. In this article, we describe how ions reflect from the interior surfaces of an ESA, the detection of which constitutes a fundamentally nonideal response of ESAs. We demonstrate this effect by comparing laboratory data from a real ESA-based space instrument, the Solar Wind Around Pluto (SWAP) instrument, aboard the NASA New Horizons spacecraft, to results from a model based on quantum mechanical simulations of particles reflected from the instrument's surfaces combined with simulations of particle trajectories through the instrument's applied electrostatic fields. Thus, we show, for the first time, how reflected ions in ESAs lead to nonideal effects that have important implications for understanding the data returned by these instruments, as well as for designing new low-background ESA-based instruments. Specifically, we show that the response of SWAP widens considerably below a level of 10{sup -3} of the peak response. Thus, a direct measurement of a plasma distribution with SWAP will have an energy-dependent background on the order of {<=}10{sup -3} of the peak of the signal due to that distribution. We predict that this order of magnitude estimate for the background applies to a large number of ESA-based instruments because ESAs operate using a common principle. However, the exact shape of the energy-dependent response will be different for different instruments. The principle of operation is that ions outside the ideal range of energy-per-charge are deflected into the walls of the ESA. Therefore, we propose that a new design paradigm is necessary to mitigate the effect of ion reflections and thus accurately and directly measure the energy spectrum of a plasma using ESAs. In this article, we build a framework for minimizing the effect of ion reflections in the design of new ESAs. Through the use of existing computer simulation software, a design team can use our method to quantify the amount of reflections in their instrument and iteratively change design parameters before fabrication, conserving resources. A possible direction for the new design paradigm is having nonsolid walls of the ESA, already used in some applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

FEM (finite element method) thermal modeling and thermal hydraulic performance of an enhanced thermal conductivity UO 2 /BeO composite fuel

An enhanced thermal conductivity UO 2 -BeO composite nuclear fuel was studied. A methodology to generate ANSYS (an engineering simulation software) FEM (Finite Element Method) thermal models of enhanced thermal conductivity oxide nuclear fuels was developed. The results showed significant increase in the fuel thermal conductivities and have good agreement with the measured ones. The reactor performance analysis showed that the decrease in centerline temperature was 250-350K for the UO 2 -BeO composite fuel, and thus we can improve nuclear reactors' performance and safety, and high-level radioactive waste generation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessing the Optical Performance Impact of Tracking Error in an Operational Concentrated Solar Power Plant Using Monte Carlo Ray-Tracing Simulation

Concentrating Solar Power (CSP) provides firm and dispatchable electricity due to its thermal storage and hybridization capabilities, which supports the decarbonization of our energy supply. Of the various CSP technologies, parabolic trough collectors are the most mature, with over 500 MW operating worldwide. The optical performance of parabolic trough systems is sensitive to tracking error, which is defined as the angular offset of a collector away from the sun position in the transversal plane. Tracking error commonly occurs due to non-continuous adjustment of the trough angle to point toward the sun, but other factors such as gravity, heating, and wind loading can also contribute to tracking error. Researchers have explored the impact of tracking error on optical performance both numerically and experimentally, but existing studies do not include measurements from operational utility-scale power plants. Tilt angle measurements of parabolic troughs at operational utility-scale power plants illustrate spatial variations in optical performance and include various sources of tracking error such as gravity, heating, and wind loading. To fully characterize wind driven loads on parabolic troughs, we are conducting a long-term field measurement campaign at the Nevada Solar One CSP plant located in Boulder City, NV, which has a nominal capacity of 72 MW and 0.5 hours of full-load storage. We record load measurements on four outer trough rows, collecting support structure bending moments, drive torque moments, dynamic accelerations of the spaceframe, mirror displacement, and tilt angles. Using the tilt measurements acquired at 20 Hz frequency, we calculate the deviation between the nominal sun position and the tracker angle. Using a Monte-Carlo ray-tracing simulation software, we assess the impact of the tracker angle deviation on optical performance throughout the diurnal cycle at various spatial locations within the CSP plant. In our view, this first-of-a-kind study will provide important guidance for future trough designs that reduce the impact of various sources of tracking error on performance.

concentrating solar power↗

Experimental evaluation of thermal-fluids performance of Helium-Cooled Flat Plate (HCFP) divertor

Various helium (He)-cooled solid-tungsten (W) divertor concepts have been proposed for long-pulse magnetic fusion energy reactors. Among these concepts, the He-cooled flat plate divertor (HCFP) modules have the largest plasma-facing surface area of ~0.2 m 2 . Simulations have shown that the most recent version of the design can withstand heat fluxes as great as 8 MW/m 2 . Earlier experimental studies of a single shortened HCFP cooling unit with a slot length of 7.6 cm used air at ambient temperature and pressures below 0.6 MPa. Here, we present initial experimental studies of a copper alloy and steel test section modeling a shortened HCFP cooling unit using He at prototypical pressure of 10 MPa, inlet temperatures T i ≤ 200 °C and steady-state incident heat fluxes q" ≤ 1.2 MW/m 2 . Results for Nusselt number Nu as a function of Reynolds number Re were obtained for Re = 1.2×10 4 –3.4×10 4 , and used to develop a Nu(Re) correlation and validate numerical models of the test section using commercial computational fluid dynamics (CFD) software. Simulations with this model are performed to evaluate the effect of the shortened slot. Furthermore, these analyses are used to estimate the thermal-fluids performance of the HCFP under prototypical conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FISPACT-II: An Advanced Simulation System for Activation, Transmutation and Material Modelling

Fispact-II is a code system and library database for modelling activation-transmutation processes, depletion-burn-up, time dependent inventory and radiation damage source terms caused by nuclear reactions and decays. The Fispact-II code, written in object-style Fortran, follows the evolution of material irradiated by neutrons, alphas, gammas, protons, or deuterons, and provides a wide range of derived radiological output quantities to satisfy most needs for nuclear applications. It can be used with any ENDF-compliant group library data for nuclear reactions, particle-induced and spontaneous fission yields, and radioactive decay (including but not limited to TENDL-2015, ENDF/B-VII.1, JEFF-3.2, JENDL-4.0u, CENDL-3.1 processed into fine-group-structure files, GEFY-5.2 and UKDD-16), as well as resolved and unresolved resonance range probability tables for self-shielding corrections and updated radiological hazard indices. The code has many novel features including: extension of the energy range up to 1 GeV; additional neutron physics including self-shielding effects, temperature dependence, thin and thick target yields; pathway analysis; and sensitivity and uncertainty quantification and propagation using full covariance data. The latest ENDF libraries such as TENDL encompass thousands of target isotopes. Nuclear data libraries for Fispact-II are prepared from these using processing codes PREPRO, NJOY and CALENDF. These data include resonance parameters, cross sections with covariances, probability tables in the resonance ranges, PKA spectra, kerma, dpa, gas and radionuclide production and energy-dependent fission yields, supplemented with all 27 decay types. All such data for the five most important incident particles are provided in evaluated data tables. The Fispact-II simulation software is described in detail in this paper, together with the nuclear data libraries. The Fispact-II system also includes several utility programs for code-use optimisation, visualisation and production of secondary radiological quantities. Included in the paper are summaries of results from the suite of verification and validation reports available with the code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multiplexed infrared photodetection using resonant radio-frequency circuits

We demonstrate a room-temperature semiconductor-based photodetector where readout is achieved using a resonant radio-frequency (RF) circuit consisting of a microstrip split-ring resonator coupled to a microstrip busline, fabricated on a semiconductor substrate. The RF resonant circuits are characterized at RF frequencies as function of resonator geometry, as well as for their response to incident IR radiation. The detectors are modeled analytically and using commercial simulation software, with good agreement to our experimental results. Though the detector sensitivity is weak, the detector architecture offers the potential for multiplexing arrays of detectors on a single read-out line, in addition to high speed response for either direct coupling of optical signals to RF circuitry, or alternatively, carrier dynamics characterization of semiconductor, or other, material systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving operational anodising process performance using simulation approach

The use of aluminium is very widespread, especially in transportation, electrical and electronics, architectural, automotive and engineering applications sectors. Therefore, the anodizing process is an important process for aluminium in order to make the aluminium durable, attractive and weather resistant. This research is focused on the anodizing process operations in manufacturing and supplying of aluminium extrusion. The data required for the development of the model is collected from the observations and interviews conducted in the study. To study the current system, the processes involved in the anodizing process are modeled by using Arena 14.5 simulation software. Those processes consist of five main processes, namely the degreasing process, the etching process, the desmut process, the anodizing process, the sealing process and 16 other processes. The results obtained were analyzed to identify the problems or bottlenecks that occurred and to propose improvement methods that can be implemented on the original model. Based on the comparisons that have been done between the improvement methods, the productivity could be increased by reallocating the workers and reducing loading time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Review of Permanent Magnet Stirring During Metal Solidification

Rather than using conventional electromagnetic stirring (EMS) with three-phase alternating current, permanent magnet stirring (PMS), based on the use of sintered NdFeB material which has excellent magnetic characteristics, can be employed to generate a magnetic field for the stirring of liquid metal during solidification. Recent experience with steel casting indicates that PMS requires less than 20 pct of the total energy compared with EMS. Despite the excellent magnetic density properties and low power consumption, this relatively new technology has received comparatively little attention by the metal casting community. This paper reviews simulation modeling, experimental studies, and industrial trials of PMS conducted during recent years. With the development of magnetic simulation software, the magnetic field and associated flow patterns generated by PMS have been evaluated. Based on the results obtained from laboratory experiments, the effects of PMS on metal solidification structures and typical defects such as surface pinholes and center cavities are summarized. The significance of findings obtained from trials of PMS within the metals processing sector, including the continuous casting of steel, are discussed with the aim of providing an overview of the relevant parameters that are of importance for further development and industrial application of this innovative technology.

36 MATERIALS SCIENCE↗

Molecular dynamics study of the mechanical loss in amorphous pure and doped silica

Gravitational wave detectors and other precision measurement devices are limited by the thermal noise in the oxide coatings on the mirrors of such devices. We have investigated the mechanical loss in amorphous oxides by calculating the internal friction using classical, atomistic molecular dynamics simulations. We have implemented the trajectory bisection method and the non-local ridge method in the DL-POLY molecular dynamics simulation software to carry out those calculations. These methods have been used to locate the local potential energy minima that a system visits during a molecular dynamics trajectory and the transition state between any two consecutive minima. Using the numerically calculated barrier height distributions, barrier asymmetry distributions, relaxation times, and deformation potentials, we have calculated the internal friction of pure amorphous silica and silica mixed with other oxides. The results for silica compare well with experiment. Finally, we use the numerical calculations to comment on the validity of previously used theoretical assumptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗