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 451 records · Page 25

Engineering the mechanical properties of ultrabarrier films grown by atomic layer deposition for the encapsulation of printed electronics

Direct deposition of barrier films by atomic layer deposition (ALD) onto printed electronics presents a promising method for packaging devices. Films made by ALD have been shown to possess desired ultrabarrier properties, but face challenges when directly grown onto surfaces with varying composition and topography. Challenges include differing nucleation and growth rates across the surface, stress concentrations from topography and coefficient of thermal expansion mismatch, elastic constant mismatch, and particle contamination that may impact the performance of the ALD barrier. In such cases, a polymer smoothing layer may be needed to coat the surface prior to ALD barrier film deposition. We present the impact of architecture on the performance of aluminum oxide (Al{sub 2}O{sub 3})/hafnium oxide (HfO{sub 2}) ALD nanolaminate barrier films deposited on fluorinated polymer layer using an optical calcium (Ca) test under damp heat. It is found that with increasing polymer thickness, the barrier films with residual tensile stress are prone to cracking resulting in rapid failure of the Ca sensor at 50 °C/85% relative humidity. Inserting a SiN{sub x} layer with residual compressive stress between the polymer and ALD layers is found to prevent cracking over a range of polymer thicknesses with more than 95% of the Ca sensor remaining after 500 h of testing. These results suggest that controlling mechanical properties and film architecture play an important role in the performance of direct deposited ALD barriers.

36 MATERIALS SCIENCE↗

Evaluating Methods for Measuring Grid Frequency in Low-Inertia Power Systems: Preprint

Accurate measurement of grid frequency is a critical component of reliable grid control. Traditionally, methods such as phase locked loops (PLLs) and discrete Fourier transforms (DFTs) have been used in inverters and phasor measurement units (PMUs) to measure frequency. However, as the percentage of inverter-based resources (IBRs) such as solar and wind has increased, these conventional frequency measurement methods are proving unable to guarantee reliable control in some cases. One challenge is measuring frequency during transient events, where there is a disruption in the steady state sinusoidal voltage. During these events, the underlying frequency of the grid may barely change, but measurement methods report a large spike in frequency due to the disrupted waveform. New methods must balance between suppressing spikes in frequency during faults, and providing fast, accurate, measurements in all other grid operation conditions, especially during events with high rate-of-change-of frequency (ROCOF), which are more prevalent in high-IBR power systems. This paper first surveys frequency measurement methods that have been proposed to reduce measurement errors during transient events. Then, both conventional and more novel frequency measurement methods are tested against an IEEE standard and industry recommendations, and their performance is evaluated for events simulated in PSCAD. Results quantify the trade-offs in performance during different grid conditions and lead to suggestions for the most appropriate frequency and ROCOF measurement methods for low inertia grids.

frequency↗

A measured energy use, solar production, and building air leakage dataset for a zero energy commercial building

This paper provides an open dataset of measured energy use, solar energy production, and building air leakage data from a 328 m 2 (3,531 ft 2 ) all-electric, zero energy commercial building in Virginia, USA. Over two years of energy use data were collected at 1-hour intervals using circuit-level energy monitors. Over six years of solar energy production data were measured at 1-hour resolution by 56 microinverters (presented as daily and monthly data in this dataset). The building air leakage data was measured post-construction per ASTM-E779 Standard Test Method for Determining Air Leakage Rate by Fan Pressurization and the United States Army Corps (USACE) Building Enclosure Testing procedure; both pressurization and depressurization results are provided. The architectural and engineering (AE) documents are provided to aid researchers and practitioners in reliable modeling of building performance. The paper describes the data collection methods, cleaning, and convergence with weather data. This dataset can be employed to predict, benchmark, and calibrate operational outcomes in zero energy commercial buildings.

14 SOLAR ENERGY↗

Clinical Outcomes of Biological Effective Dose-Based Fractionated Stereotactic Radiation Therapy for Metastatic Brain Tumors From Non-Small Cell Lung Cancer

Purpose: To evaluate the efficacy and toxicity of fractionated stereotactic radiation therapy (FSRT) based on biological effective dose (BED), a novel approach to deliver a fixed BED irrespective of dose fractionation, for brain metastases from non-small cell lung cancer (NSCLC). Methods and Materials: Between March 2005 and March 2009 we treated 299 patients with 1 to 5 lesions from NSCLC (573 total brain metastases) with FSRT using Novalis. The dose fractionation schedules were individually determined to deliver a peripheral BED10 (α/β ratio = 10) of approximately 80 Gy{sub 10}. The median number of fractions was 3 (range, 2-10), the median peripheral BED10 was 83.2 Gy (range, 19.1-89.6 Gy). Patients were followed up with magnetic resonance imaging (MRI) studies performed at 1- to 2-month intervals. The local tumor control rate and overall local progression-free and intracranial relapse-free survival were calculated by the Kaplan-Meier method. Results: Local control rates for all 573 lesions at 6 and 12 months were 96.3% and 94.5%, respectively. By multivariate analysis the tumor diameter was the only factor predictive of the local control rate (P=.001). The median overall survival, local progression-free survival, and intracranial relapse-free survival were 17.1, 14.9, and 4.4 months, respectively. The overall survival, local progression-free survival, and intracranial relapse-free survival rates at 6 and 12 months were 78.5% and 63.3%, 74.3% and 57.8%, and 41.0% and 21.8%, respectively. Six patients (2%) manifested progressive radiation injury to the brain even during therapy with corticosteroids; they underwent hyperbaric oxygen therapy, and follow-up MRI showed improvement. Conclusions: This study showed that BED-based FSRT for brain metastases from NSCLC is a promising strategy that may yield excellent outcomes with acceptable toxicity. Criteria must be established to determine the optimal dose fractionation for individual patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Energy levels, lifetimes and radiative rates for transitions in the bromine isoelectronic sequence La XXIII-Dy XXXII, W XL

Using the multiconfiguration Dirac–Fock method, calculations for the lowest 62 levels of the ([Ar] 3d{sup 10})4s{sup 2}4p{sup 5}, ([Ar] 3d{sup 10})4s{sup 2}4p{sup 3}4d{sup 2}, ([Ar] 3d{sup 10})4s{sup 2}4p{sup 4}4d, ([Ar] 3d{sup 10})4s4p{sup 6}, and ([Ar] 3d{sup 10})4s4p{sup 5}4d configurations are performed for the bromine isoelectronic sequence La XXIII-Dy XXXII, W XL. Results of energy levels, lifetimes, wavelengths, and electric dipole, magnetic dipole, electric quadrupole, and magnetic quadrupole radiative rates are presented. In order to assess the accuracy of results, independent calculations for W XL have been carried out using the many-body perturbation theory (MBPT) method. Comparisons are made with available theoretical results from other calculations and the observed values of the Atomic Spectra Database of the National Institute of Standards and Technology. Energy levels are estimated to be accurate to better than 1%, and radiative rates (and lifetimes) are accurate to better than 20% for a majority of strong transitions. These results should be useful in many applications of lanthanide ions related to broad area of research such as applied physics, laser physics and fusion science.

74 ATOMIC AND MOLECULAR PHYSICS↗

Prediction of liver remnant regeneration after living donor liver transplantation using preoperative CT texture analysis

Purpose: To predict the rate of liver regeneration after living donor liver transplantation (LDLT) using pre-operative computed tomography (CT) texture analysis. Materials and methods: 112 living donors who performed right hepatectomy for LDLT were included retrospectively. We measured the volume of future remnant liver (FLR) on pre-operative CT and the volume of remnant liver (LR) on follow-up CT, taken at a median of 123 days after transplantation. The regeneration index (RI) was calculated using the following equation: [(V{sub LR}−V{sub FLR})/V{sub FLR}] × 100. Computerized texture analysis of the semi-automatically segmented FLR was performed. We used a stepwise, multivariable linear regression to assess associations of clinical features and texture parameters in relation to RI and to make the best-fit predictive model. Results: The mean RI was 110.7 ± 37.8%, highly variable ranging from 22.4% to 247.0%. Among texture parameters, volume of FLR, standard deviation, variance, and gray level co-occurrence matrices (GLCM) contrast were found to have significant correlations between RI. In multivariable analysis, smaller volume of FLR (ß − 0.17, 95% CI − 0.22 to − 0.13) and lower GLCM contrast (ß − 1.87, 95% CI − 3.64 to − 0.10) were associated with higher RI. The regression equation predicting RI was following: RI = 203.82 + 10.42 × pre-operative serum total bilirubin (mg/dL) − 0.17 × V{sub FLR} (cm{sup 3}) − 1.87 × GLCM contrast (× 100). Conclusion: Volume of FLR and GLCM contrast were independent predictors of RI, showing significant negative correlations. Pre-operative CT with texture analysis can be useful for predicting the rate of liver regeneration in living donor of liver transplantation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Simulations and Experimental Verifications of an Algorithm for Radiation Source Mapping and Navigational Path Generation

Accurate and efficient mapping and localization of both ionizing and non-ionizing radiation sources are important across many different fields. As such, a versatile mapping and navigational path generation algorithm, which can be applied to any point source measurements that follow an inverse-square characteristic, was developed using non-linear least squares methods. Forty thousand simulations were performed on the algorithm, which located sources successfully in a 10 m × 10 m × 10 m three-dimensional space with a success rate of over 80% across different noise functions, given a proportional constant of 10 to 1,000. The algorithm was also verified experimentally with small-scale radioactive decontamination of a 70 cm × 70 cm surface and localization of a lost Wi-Fi router in a 70 m × 70 m open field. One hundred twenty-one measurements were taken from each experiment, which were then fed into the algorithm for navigation. For the radioactive 137Cs source, the estimated locations were within 7 cm × 7 cm of the answer in 79.3% of the scenarios, while the Wi-Fi router was located to within 7 m × 7 m in 57.9% of the tests. In general, the method requires much less information and data than a geographically comprehensive survey and thus shows a lot of potential for practical applications, such as lost source retrieval with unmanned aerial vehicles, small-scale decontamination, mapping undocumented Wi-Fi routers or radio towers, and radiation simulation with radio signals. Here, the different failure modes, desirable features, and potential improvements were also identified but remain as future work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparison of One- and Two-Variable Linear Regression Models and Classic Energy Intensity for Energy Performance Tracking of Two Manufacturing Sectors

We report manufacturing facilities consumed about 32% of total domestic energy in the United States in 2016. To evaluate the energy savings achieved through the implementation of energy conservation projects and to establish powerful arguments for future projects, accurate energy performance tracking methods are necessary. The classic energy intensity method (i.e., the ratio of annual total energy over annual total production) is the means most widely used to measure savings because it can be understood and calculated easily. This method considers the variation in production rates to some extent; however, it fundamentally assumes facilities' base energy consumption (energy consumption at zero production) to be zero, which rarely holds true. Furthermore, this method does not consider variations in other relevant parameters, such as weather conditions. Therefore, the regression models approach is commonly recommended to track energy performance improvements. However, because it requires more data and statistical expertise, the regression models approach has been adopted by only a few facilities, and many others suspect it is not worth the effort for their specific cases. For this reason, the improved accuracy this approach offers needs to be demonstrated. By analyzing 477 monthly energy (electricity and natural gas) data sets, this study quantitatively compared the accuracy of classic energy intensity, one-variable (production only) linear regression models, and two-variable (i.e., production and weather) linear regression models for two manufacturing sectors (primary metal and transportation equipment manufacturing). Results showed that significant improvements in accuracy were achieved with one-variable regression models when compared with the classic energy intensity method and with two-variable regression models when compared with one-variable regression models. These improvements were demonstrated using the p values of intercept, cooling degree days, and heating degree days. Based on these results, to achieve a good balance between accuracy improvements and resources requirements, one-variable (production only) linear regression models for electricity consumption and two-variable (production and heating degree days) linear regression models for natural gas consumption are recommended for both sectors in the case of limited resources. Lastly, facilities can also use the results to decide which approaches fit better in their cases.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Development of a High Throughput Platform for Screening Glycoside Hydrolases Based on Oxime-NIMS

Cost-effective hydrolysis of biomass into sugars for biofuel production requires high-performance low-cost glycoside hydrolase (GH) cocktails that are active under demanding process conditions. Improving the performance of GH cocktails depends on knowledge of many critical parameters, including individual enzyme stabilities, optimal reaction conditions, kinetics, and specificity of reaction. With this information, rate- and/or yield-limiting reactions can be potentially improved through substitution, synergistic complementation, or protein engineering. Given the wide range of substrates and methods used for GH characterization, it is difficult to compare results across a myriad of approaches to identify high performance and synergistic combinations of enzymes. Here, we describe a platform for systematic screening of GH activities using automatic biomass handling, bioconjugate chemistry, robotic liquid handling, and nanostructure-initiator mass spectrometry (NIMS). Twelve well-characterized substrates spanning the types of glycosidic linkages found in plant cell walls are included in the experimental workflow. To test the application of this platform and substrate panel, we studied the reactivity of three engineered cellulases and their synergy of combination across a range of reaction conditions and enzyme concentrations. We anticipate that large-scale screening using the standardized platform and substrates will generate critical datasets to enable direct comparison of enzyme activities for cocktail design.

09 BIOMASS FUELS↗

Comparative performance of modern digital mammography systems in a large breast screening program

Purpose: To compare physical measures pertaining to image quality among digital mammography systems utilized in a large breast screening program. To examine qualitatively differences in these measures and differences in clinical cancer detection rates between CR and DR among sites within that program. Methods: As part of the routine quality assurance program for screening, field measurements are made of several variables considered to correlate with the diagnostic quality of medical images including: modulation transfer function, noise equivalent quanta, d′ (an index of lesion detectability) and air kerma to allow estimation of mean glandular dose. In addition, images of the mammography accreditation phantom are evaluated. Results: It was found that overall there were marked differences between the performance measures of DR and CR mammography systems. In particular, the modulation transfer functions obtained with the DR systems were found to be higher, even for larger detector element sizes. Similarly, the noise equivalent quanta, d′, and the phantom scores were higher, while the failure rates associated with low signal-to-noise ratio and high dose were lower with DR. These results were consistent with previous findings in the authors’ program that the breast cancer detection rates at sites employing CR technology were, on average, 30.6% lower than those that used DR mammography. Conclusions: While the clinical study was not large enough to allow a statistically powered system-by-system assessment of cancer detection accuracy, the physical measures expressing spatial resolution, and signal-to-noise ratio are consistent with the published finding that sites employing CR systems had lower cancer detection rates than those using DR systems for screening mammography.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Solid-State Transformer and Hybrid Transformer With Integrated Energy Storage in Active Distribution Grids: Technical and Economic Comparison, Dispatch, and Control

Solid-state transformer (SST) and hybrid transformer (HT) are promising alternatives to the line-frequency transformer (LFT) in smart grids. The SST features medium-frequency isolation, full controllability for voltage regulation, reactive power compensation, and the capability of battery energy storage system (BESS) integration with multiport configuration. The HT has a partially-rated converter for fractional controllability and can integrate a small BESS. Fast grid-edge voltage fluctuations from increased solar photovoltaic (PV) and electric vehicle (EV) penetration are difficult to manage for mechanical load tap changers. Hence, along with the trend towards more BESS in the grid, the controllability and the storage integration capability of the SST and HT are of strong interest. However, a review of literature shows existing SST and HT research is mostly at converter level, while system-level assessments are scarce. Assessing technical and economic impacts is critical to understanding the benefits and role of the SST and HT to guide future research, which is presented for the first time in this article. Experimental results from medium-voltage (MV) SST and MV HT prototypes are shown to confirm equipment-level feasibility, where the voltage controllability waveforms of a MV HT prototype are reported for the first time. Comparative simulations are performed on a modified IEEE 34-bus system. Here, a grid-model-less decentralized grid-edge voltage control method and a day-ahead BESS dispatch method are proposed for the SST and HT. The simulations show that the SST and HT with integrated storage can host more PV, achieve peak shaving, mitigate voltage fluctuation and reverse power flow, and support energy arbitrage for operational cost reduction, as compared to the LFT. Moreover, comprehensive analyses of net present value (NPV) and internal rate of return (IRR) are performed under different installed PV capacities, HT’s partial converter ratings, and BESS capacities. Sensitivities to future cost reductions of the PV and BESS are studied. Although the NPV and IRR are currently negative, 60% capital cost reduction or 150% revenue increase will make the SST and HT economically viable in the use case studied.

14 SOLAR ENERGY↗

Energy Efficient Thermoplastic Composite Manufacturing

The objective of the project was to establish an effective and affordable method to lay-up and consolidate/join large thermoplastic composite aerospace structure with cycle times measured in minutes rather than hours. Composite airplane designs have proven efficient and effective however future potential product production rates are challenging what the current systems can efficiently produce due to material lay-down constraints and extended thermal cycle times. The ability to lay-up then rapidly heat, consolidate, and cool large complex composite structures plus very accurately tool them (i.e. matching CTE of composite materials) along with very precise thermal control is a difficult challenge. Current systems such as autoclave processing of thermoset materials require long cycle times due to method of heating and the large associated thermal masses. These extended cycle times inhibit the ability to meet higher rate production scenarios due to the need for multiple sets of equipment and tools. Thermoplastic composite materials were used to facilitate more rapid cycle times via the elimination of a need for a cure dwell at temperature. Also, utilized induction heating along with smart susceptors to enable the quick cycle times needed while providing precise intrinsic thermal control. The intrinsic control of the induction heating process via the smart susceptors along with the use of laser assisted fiber placement were the key innovations areas developed. The ability to rapidly heat and then precisely control the temperature of the consolidation component was paramount to the success of this project. It is this unique processing attribute coupled with the thermoplastic material characteristics along with rapid lay-up methods that are novel and advantageous. While many requirements and influencing factors decide the materials and processes utilized for future airplane construction, the forecast of accelerated production rates and the recent performance successes of composites in airplane construction provide an opportunity for this processing technology to have significant influence. During the course of this project a number of key accomplishments were completed. They are as follows: The needed laser assisted fiber placement capabilities have been developed and validated on a component of meaningful scale and complexity. Subsequently, a thermoplastic preform of the part scale-up part design was laser assisted fiber placed and made availalbe for consolidation. An induction consolidation system was sized, designed, and fabricated resulting in the full capability to consolidate the scale-up component. The tool for induction consolidating of the scale-up component was designed and fabricated. An aluminum pressure bladder was welded and formed for use as the pressurization membrane to perform the needed application of pressure on the part at temperature to achieve consolidation. The induction consolidation system consisting of the restraint, induction power supply, induction tool with integrated induction coils and smart susceptor liners, along with the aluminum bladder were successfully used to consolidate a large thermoplastic skin. The capability of accomplishing rapid heat-up rates and cool down rates for consolidation of large thermoplastic skins with precise thermal control and even pressure application at the consolidation temperature was validated. The general conclusion is that the use of induction consolidation with smart susceptors is an effective method to rapidly and efficiently consolidate large thermoplastic composite components and enable affordable high rate manufacturing of composites. More specific conclusions generated by this project are as follows: A restraint type induction consolidation system provides an affordable option over the standard press when using the induction consolidation processing method. Aluminum bladders can be formed to the complex shapes needed and supply uniform pressure for consolidation and co-consolidation of thermoplastic composite materials. Straightforward modifications can be made to existing fiber placement robotic systems to enable their use for thermoplastic lay-up applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Maximum Inner Container Dose Rate That Will Ensure Meeting Regulatory Limit at Package Surface - 19206

Often, the shipper of radioactive material does not know the exact radioactive isotopic composition of the materials that need to be shipped. This is especially important for legacy materials that need to be removed from a site that is in the process of being decommissioned. A method of determining the acceptability of such material for shipment before it is placed in into the shipping package will allow more flexibility in planning for decommissioning activities. A series of Monte-Carlo calculations were performed to calculate the dose rate at the surface of a convenience container or a shielding pig separately and then at the surface of a new shipping package being designed by the Savannah River National Laboratory (SRNL). The purpose of these calculations is to determine the maximum dose rate at the surface of the convenience container or the shielded pig that will ensure the dose rate at the outer surface the shipping package is less that the regulatory limit of 2 mSv/h. The calculations also considered the effect of shape of the radioactive material and its possible movement inside the convenience container along with movement of the convenience container inside the confinement vessel of the shipping package. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Enabling high rate charge and discharge capability, low internal resistance, and excellent cycleability for Li-ion batteries utilizing graphene additives

A liquid-phase mixing method is adopted to uniformly disperse the graphene nanosheets onto LiNi 1/3 Co 1/3 Mn 1/3 O 2 cathode for high-performance Li-ion batteries (LIBs). The electrochemical performance was characterized using a full pouch cells with state-of-the-art electrode areal loading (compared to half coin cells). The addition of graphene sheets (i.e., only 1 wt%) significantly improves the high rate capability for charging and discharging operation. For example, 6 times improvement in 5 C charging was achieved providing further insights in enabling extreme fast charging for LIBs. Other benefits include longer cycleability, lower internal resistance, and higher lithium ion diffusion coefficient, demonstrated by charge-discharge cycling tests and electrochemical impedance spectroscopy. Higher capacity retention of 88.2% and decreased internal resistance of ~0.9 Ω are observed after 400 cycles. The diffusion coefficient of Li ions is 6.49 × 10 -8 cm 2 s -1 when charged to 4.2 V, which is approximately 1.37 times higher compared to the configuration with no graphene sheet (4.74 × 10 -8 cm 2 s -1 ). To conclude, the improved performance is ascribed to a robust network among the active materials formed by graphene sheets, which serves as an extended current conductor and facilitates charge transfer, ionic reversibility, and ionic transportation.

25 ENERGY STORAGE↗

SU-E-T-795: Validations of Dose Calculation Accuracy of Acuros BV in High-Dose-Rate (HDR) Brachytherapy with a Shielded Cylinder Applicator Using Monte Carlo Simulation

Purpose: Acuros BV has become available to perform accurate dose calculations in high-dose-rate (HDR) brachytherapy with phantom heterogeneity considered by solving the Boltzmann transport equation. In this work, we performed validation studies regarding the dose calculation accuracy of Acuros BV in cases with a shielded cylinder applicator using Monte Carlo (MC) simulations. Methods: Fifteen cases were considered in our studies, covering five different diameters of the applicator and three different shielding degrees. For each case, a digital phantom was created in Varian BrachyVision with the cylinder applicator inserted in the middle of a large water phantom. A treatment plan with eight dwell positions was generated for these fifteen cases. Dose calculations were performed with Acuros BV. We then generated a voxelized phantom of the same geometry, and the materials were modeled according to the vendor’s specifications. MC dose calculations were then performed using our in-house developed fast MC dose engine for HDR brachytherapy (gBMC) on a GPU platform, which is able to simulate both photon transport and electron transport in a voxelized geometry. A phase-space file for the Ir-192 HDR source was used as a source model for MC simulations. Results: Satisfactory agreements between the dose distributions calculated by Acuros BV and those calculated by gBMC were observed in all cases. Quantitatively, we computed point-wise dose difference within the region that receives a dose higher than 10% of the reference dose, defined to be the dose at 5mm outward away from the applicator surface. The mean dose difference was ∼0.45%–0.51% and the 95-percentile maximum difference was ∼1.24%–1.47%. Conclusion: Acuros BV is able to accurately perform dose calculations in HDR brachytherapy with a shielded cylinder applicator.

60 APPLIED LIFE SCIENCES↗

Combining shallow-water and analytical wake models for tidal array micro-siting

For tidal-stream energy to become a competitive renewable energy source, clustering multiple turbines into arrays is paramount. Array optimisation is thus critical for achieving maximum power performance and reducing cost of energy. However, ascertaining an optimal array layout is a complex problem, subject to specific site hydrodynamics and multiple inter-disciplinary constraints. In this work, we present a novel optimisation approach that combines an analytical-based wake model, FLORIS, with an ocean model, Thetis. The approach is demonstrated through applications of increasing complexity. By utilising the method of analytical wake superposition, the addition or alteration of turbine position does not require re-calculation of the entire flow field, thus allowing the use of simple heuristic techniques to perform optimisation at a fraction of the computational cost of more sophisticated methods. Using a custom condition-based placement algorithm, this methodology is applied to the Pentland Firth for arrays with turbines of 3.05 m/s rated speed, demonstrating practical implications whilst considering the temporal variability of the tide. For a 24-turbine array case, micro-siting using this technique delivered an array 15.8% more productive on average than a staggered layout, despite flow speeds regularly exceeding the rated value. Performance was evaluated through assessment of the optimised layout within the ocean model that treats turbines through a discrete turbine representation. Used iteratively, this methodology could deliver improved array configurations in a manner that accounts for local hydrodynamic effects.

16 TIDAL AND WAVE POWER↗

Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface

A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model—the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4-fold for a two-level implementation, and can be increased up to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of r{sub s}=0.912.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preservation of physical properties of stochastic Maxwell equations with additive noise via stochastic multi-symplectic methods

Stochastic Maxwell equations with additive noise are a system of stochastic Hamiltonian partial differential equations intrinsically, possessing the stochastic multi-symplectic conservation law. It is shown that the averaged energy increases linearly with respect to the evolution of time and the flow of stochastic Maxwell equations with additive noise preserves the divergence in the sense of expectation. Moreover, we propose three novel stochastic multi-symplectic methods to discretize stochastic Maxwell equations in order to investigate the preservation of these properties numerically. We make theoretical discussions and comparisons on all of the three methods to observe that all of them preserve the corresponding discrete version of the averaged divergence. Meanwhile, we obtain the corresponding dissipative property of the discrete averaged energy satisfied by each method. Especially, the evolution rates of the averaged energies for all of the three methods are derived which are in accordance with the continuous case. Numerical experiments are performed to verify our theoretical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗