Search NASA⌕ Search

SEARCH · Search NASA

Results for “Design Methods Modeling and Validation”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling of deposit formation in mesoporous substrates via atomic layer deposition: Insights from pore‐scale simulation

Atomic layer deposition (ALD) has been a promising technique in fabricating membranes and tuning their properties with a precision at the atomic level. Fabrication of zeolitic imidazolate framework (ZIF) membranes using the ligand-induced permselectivation (LIPS) method starts with the formation of an oxide in a mesoporous substrate by ALD and is followed by the transformation of this oxide to ZIF using imidazolate vapor treatment. The objective of the ALD step is to block the mesopores with a thin deposit, that is, one with small penetration depth and small thickness on the top surface of the substrate. Unlike typical ALD on nonporous substrates, where all available sites react per ALD cycle, thin deposit formation in a mesoporous substrate requires that only a small fraction of the available deposition sites (i.e., close to the substrate surface) is subjected to ALD. Consequently, reactant dosing and duration of pulses are important process variables which, together with diffusion and reaction kinetics determine the deposit structure. Quantitative understanding of the interplay of these variables and phenomena can enable the rational design of ALD within mesoporous substrates. Here, we extend our earlier modeling effort considering the coexistence of ALD both inside the pores and on the external surface of the substrate. Finite-volume based models were developed and validated to simulate the two distinct modes of deposition cycle by cycle. The total mass uptake of the substrate with ALD cycles can be predicted using the combined surface deposition and pore reaction–diffusion models as affirmed by in situ quartz crystal microbalance experimental data. The ALD reactor model combined with the deposition model can accurately capture the number of ALD cycles needed to block the pores of the substrate. Based on the model, we designed a modified ALD process and examined the performance of the corresponding LIPS membranes. Furthermore, the present modeling work provides a new understanding of the deposit formation via ALD within mesoporous substrates for a variety of membrane applications.

atomic layer deposition↗
DOE OSTI2022

Flow Measurements in Pressurized Helium Gas Cooling Channels

Production of metastable Technetium-99 (Tc-99m) is vital to the medical imaging community. It is a radioactive tracer that emits gamma rays as it decays and it is these gamma rays that are then detected through imaging. Tc-99m is extracted from the decay of Molybdenum-99 (Mo-99) which has a half-life of about 2-3 days. One method of generating Mo-99 is using accelerator based technology, through the irradiation of Mo-100 using an electron beam. This is part of the NNSA’s mission to produce Mo-99 commercially, within the US, without the use of highly enriched uranium (HEU) in support of nonproliferation and global security. Los Alamos National Laboratory (LANL) is working with NorthStar medical Radioisotopes (NMR) on their efforts to produce Mo-99 from the irradiation of Mo-100 targets. The NMR target consists of an Inconel window that allows the electron beam to penetrate and irradiate a stack of Mo-100 discs. The irradiation process generates large amounts of heat and requires extensive cooling of the Mo-100 targets and target housing using pressurized helium gas. Before beam on target, the design of the target stack and target housing is tested at LANL using laser based diagnostics experiments. We quantify the velocity field of the helium gas around the target housing to estimate whether the cooling provided is sufficient for the design and to validate computational fluid dynamics (CFD) models. Fiber optic imaging is also developed for use during cold testing (without beam on target) as well as during irradiation, to obtain real time information on disc performance.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Expanding the design space of stratospheric aerosol geoengineering to include precipitation-based objectives and explore trade-offs

Previous climate modeling studies demonstrate the ability of feedback-regulated, stratospheric aerosol geoengineering with injection at multiple independent latitudes to meet multiple simultaneous temperature-based objectives in the presence of anthropogenic climate change. However, the impacts of climate change are not limited to rising temperatures but also include changes in precipitation, loss of sea ice, and many more; knowing how a given geoengineering strategy will affect each of these climate metrics is vital to understanding the limits and trade-offs of geoengineering. Here, we first introduce a new method of visualizing the design space in which desired climate outcomes are represented by 2-D surfaces on a 3-D graph. Surface orientations represent how different injection choices influence that objective, and intersecting surfaces represent objectives which can be met simultaneously. Using this representation as a guide, we present simulations of two new strategies for feedback-regulated aerosol injection, using the Community Earth System Model with the Whole Atmosphere Community Climate Model – CESM1(WACCM). The first simultaneously manages global mean temperature, tropical precipitation centroid, and Arctic sea ice extent, while the second manages global mean precipitation, tropical precipitation centroid, and Arctic sea ice extent. Both simulations control the tropical precipitation centroid to within 5 % of the goal, and the latter controls global mean precipitation to within 1% of the goal. Additionally, the first simulation overcompensates sea ice, while the second undercompensates sea ice; all of these results are consistent with the expectations of our design space model. In addition to showing that precipitation-based climate metrics can be managed using feedback alongside other goals, our simulations validate the utility of our design space visualization in predicting our climate model behavior under a given geoengineering strategy, and together they help illustrate the fundamental limits and trade-offs of stratospheric aerosol geoengineering.

58 GEOSCIENCES↗
NASA NTRS1988

Methods for heat transfer and temperature field analysis of the insulated diesel, phase 3

Work during Phase 3 of a program aimed at developing a comprehensive heat transfer and thermal analysis methodology for design analysis of insulated diesel engines is described. The overall program addresses all the key heat transfer issues: (1) spatially and time-resolved convective and radiative in-cylinder heat transfer, (2) steady-state conduction in the overall structure, and (3) cyclical and load/speed temperature transients in the engine structure. These are all accounted for in a coupled way together with cycle thermodynamics. This methodology was developed during Phases 1 and 2. During Phase 3, an experimental program was carried out to obtain data on heat transfer under cooled and insulated engine conditions and also to generate a database to validate the developed methodology. A single cylinder Cummins diesel engine was instrumented for instantaneous total heat flux and heat radiation measurements. Data were acquired over a wide range of operating conditions in two engine configurations. One was a cooled baseline. The other included ceramic coated components (0.050 inches plasma sprayed zirconia)-piston, head and valves. The experiments showed that the insulated engine has a smaller heat flux than the cooled one. The model predictions were found to be in very good agreement with the data.

Morel, Thomas↗
DOE OSTI2019

An improved laboratory-based x-ray absorption fine structure and x-ray emission spectrometer for analytical applications in materials chemistry research

X-ray absorption fine structure (XAFS) and x-ray emission spectroscopy (XES) are advanced x-ray spectroscopies that impact a wide range of disciplines. Furthermore, unlike the majority of other spectroscopic methods, XAFS and XES are accompanied by an unusual access model, wherein the dominant use of the technique is for premier research studies at world-class facilities, i.e., synchrotron x-ray light sources. In this paper, we report the design and performance of an improved XAFS and XES spectrometer based on the general conceptual design of Seidler et al. [Rev. Sci. Instrum. 85, 113906 (2014)]. New developments include reduced mechanical degrees of freedom, much-increased flux, and a wider Bragg angle range to enable extended x-ray absorption fine structure (EXAFS) measurement and analysis for the first time with this type of modern laboratory XAFS configuration. This instrument enables a new class of routine applications that are incompatible with the mission and access model of the synchrotron light sources. To illustrate this, we provide numerous examples of x-ray absorption near edge structure (XANES), EXAFS, and XES results for a variety of problems and energy ranges. Highlights include XAFS and XES measurements of battery electrode materials, EXAFS of Ni with full modeling of results to validate monochromator performance, valence-to-core XES for 3d transition metal compounds, and uranium XANES and XES for different oxidation states. Taken en masse , these results further support the growing perspective that modern laboratory-based XAFS and XES have the potential to develop a new branch of analytical chemistry.

47 OTHER INSTRUMENTATION↗
NASA NTRS1993

On the Accuracy of CFD-Based Pressure Drop Predictions for Right-Angle Ducts

The predictive capability of computational fluid dynamics (CFD) codes for turbulent flow through curved ducts is of significant importance to the design and performance analysis of modern rocket engine flowpaths. Code calibration and validation studies for this class of flow are desireable to estimate the performance margin and operating range of components designed using Navier-Stokes methods. Parametric experimental studies such as that of Weske (NACA ARR W-39) provided a wealth of performance data for the design of single- and compound elbow configurations with various cross-sections, curvature and aspect ratios at varying Reynolds numbers. In that work, the majority of data is presented in the form of loss coefficients, characterizing pressure losses due to duct curvature, and including losses due to wall friction. Using measured friction coefficients, losses of equivalent straight lengths of duct are subtracted, resulting in performance curves useful for design computations. These data are currently used in a CFD-based parametric study covering a broad range of operating conditions. Of particular interest for the accuracy of CFD predictions are the effects on pressure loss due to inlet boundary layer thickness (dependent on upstream development length), and the wall treatment for the turbulence equations (conventional wall functions vs. wall integration using a two-layer model). The experimental data are reassessed in the form of an error analysis, and are compared with CFD predictions for 18 computational cases. Grid-independence, grid spacing, and convergence requirements of the cases are discussed. Conclusions regarding the relative importance of the parametric variables will be presented.

Brankovic, Andreja↗
NASA NTRS1998

Flight Acoustics Measurement Techniques and Applications

Careful consideration must be given to data acquisition and analysis techniques in the design of experiments for the measurement of noise generated by flight vehicles. Although noise measurement locations and data reduction procedures are specified for aircraft certification by FAA and ICAO directives, for example, there are virtually no established procedures for aircraft noise measurement for other purposes. To optimize the quality and quantity of information obtained in a flight acoustics experiment, microphone layout, data acquisition, and analysis must be tailored to the specific test objective. This paper will review flight acoustics technology at NASA Langley Research Center developed over the past decade. In particular, the paper will focus on flight experiments performed for three diverse objectives: (1) research applications, such as noise prediction code validation, (2) noise impact modeling, and (3) noise abatement flight procedures. To best achieve these diverse objectives, different deployments of microphone systems on the ground are required, and different data analysis techniques are needed. In all cases, accurate positioning of the aircraft synchronized in time with the data recording is necessary. However, there are some restrictions on flight operations unique to each case for the methods to properly work.

Preisser, J. S.↗
DOE OSTI2018

Conceptual design of extended magnetic probe set to improve 3D field detection in NSTX-U

Adding toroidal arrays of magnetic probes at the top and bottom of NSTX-U would improve both the detection of the multimodal plasma response to applied magnetic perturbations and the identification of the poloidal structure of unstable plasma modes, as well as contribute to the validation of MHD models, improve the understanding of the plasma response to external fields, and improve the error field correction. In this paper, the linear MHD code MARS-F/K has been used to identify poloidal locations that would improve the capability to measure stationary or near-stationary 3D fields that may result from the plasma response to external sources of non-axisymmetric fields. The study highlighted 6 poloidal positions where new arrays of both poloidal and radial magnetic field sensors would improve the poloidal resolution. Here, the proposed set of new arrays combined with the present ones is shown to be capable of measuring the poloidal structure of perturbations with n ≤ 6 and of detecting the multimodal plasma response. Assessment of the trade-off in the poloidal length of the probes leads to an ideal length between 10 cm and 30 cm. A method to configure the probes of a toroidal array based on the singular value decomposition condition number is proposed, and an ideal solution and a low-cost one are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Inadequacy of current approaches for characterizing membrane transport properties at high salinities

Cost optimal design of osmotic membrane processes requires an accurate estimate of membrane transport parameters across their full operational range. However, standard approaches for estimating these parameters rely on empirical methods, the accuracy of which remains unquantified as a function of temperature, salinity, and measurement error. Herein, we present a systematic accuracy analysis of previously developed methods for estimation of membrane transport properties in reverse osmosis, high-pressure reverse osmosis, forward osmosis, pressure retarded reverse osmosis, and osmotically assisted reverse osmosis. We use a Monte Carlo approach to sample the full range of feasible membrane water permeabilities, salt permeabilities, structural parameters, and operating conditions for these processes. These material and process parameters are then incorporated into a physical transport model for each process. Our analysis shows that the statistical uncertainty of current empirical methods for estimating membrane parameters increases by 5 times from low-salinity to high-salinity conditions. Here, the result of this work demonstrates that empirical methods are inadequate for precisely estimating membrane transport properties at high salinity and highlight a critical need for the development of statistically validated higher accuracy methods.

42 ENGINEERING↗
DOE OSTI2019

Development of a Generalized Residual Stress Inversion Technique

Residual stress is a common result of manufacturing processes, but it is one that is often overlooked in design and qualification activities. There are many reasons for this oversight, such as lack of observable indicators and difficulty in measurement. Traditional relaxation-based measurement methods use some type of material removal to cause surface displacements, which can then be used to solve for the residual stresses relieved by the removal. While widely used, these methods may offer only individual stress components or may be limited by part or cut geometry requirements. Diffraction-based methods, such as X-ray or neutron, offer non-destructive results but require access to a radiation source. With the goal of producing a more flexible solution, this LDRD developed a generalized residual stress inversion technique that can recover residual stresses released by all traction components on a cut surface, with much greater freedom in part geometry and cut location. The developed method has been successfully demonstrated on both synthetic and experimental data. The project also investigated dislocation density quantification using nonlinear ultrasound, residual stress measurement using Electronic Speckle Pattern Interferometry Hole Drilling, and validation of residual stress predictions in Additive Manufacturing process models.

42 ENGINEERING↗
NASA NTRS2015

Investigation of Spiral Bevel Gear Condition Indicator Validation via AC-29-2C Combining Test Rig Damage Progression Data with Fielded Rotorcraft Data

This is the final of three reports published on the results of this project. In the first report, results were presented on nineteen tests performed in the NASA Glenn Spiral Bevel Gear Fatigue Test Rig on spiral bevel gear sets designed to simulate helicopter fielded failures. In the second report, fielded helicopter HUMS data from forty helicopters were processed with the same techniques that were applied to spiral bevel rig test data. Twenty of the forty helicopters experienced damage to the spiral bevel gears, while the other twenty helicopters had no known anomalies within the time frame of the datasets. In this report, results from the rig and helicopter data analysis will be compared for differences and similarities in condition indicator (CI) response. Observations and findings using sub-scale rig failure progression tests to validate helicopter gear condition indicators will be presented. In the helicopter, gear health monitoring data was measured when damage occurred and after the gear sets were replaced at two helicopter regimes. For the helicopters or tails, data was taken in the flat pitch ground 101 rotor speed (FPG101) regime. For nine tails, data was also taken at 120 knots true airspeed (120KTA) regime. In the test rig, gear sets were tested until damage initiated and progressed while gear health monitoring data and operational parameters were measured and tooth damage progression documented. For the rig tests, the gear speed was maintained at 3500RPM, a one hour run-in was performed at 4000 in-lb gear torque, than the torque was increased to 8000 in-lbs. The HUMS gear condition indicator data evaluated included Figure of Merit 4 (FM4), Root Mean Square (RMS) or Diagnostic Algorithm 1(DA1), + 3 Sideband Index (SI3) and + 1 Sideband Index (SI1). These were selected based on their sensitivity in detecting contact fatigue damage modes from analytical, experimental and historical helicopter data. For this report, the helicopter dataset was reduced to fourteen tails and the test rig data set was reduced to eight tested gear sets. The damage modes compared were separated into three cases. For case one, both the gear and pinion showed signs of contact fatigue or scuffing damage. For case two, only the pinion showed signs of contact fatigue damage or scuffing. Case three was limited to the gear tests when scuffing occurred immediately after the gear run-in. Results of this investigation highlighted the importance of understanding the complete monitored systems, for both the helicopter and test rig, before interpreting health monitoring data. Further work is required to better define these two systems that include better state awareness of the fielded systems, new sensing technologies, new experimental methods or models that quantify the effect of system design on CI response and new methods for setting thresholds that take into consideration the variance of each system.

spiral bevel gears↗
DOE OSTI2019

Turbulent nonpremixed cool flames: Experimental measurements, Direct Numerical Simulation, and manifold-based combustion modeling

Turbulence, low-temperature chemistry, and their interactions in the form of turbulent cool flames are critical to understanding and improving advanced engines. Design of such engines requires tractable simulations which in turn necessitate turbulent combustion models that can account for cool flames. While manifold-based turbulent combustion models are an attractive option for hot flames, their applicability to cool flames is not yet fully understood. This is partially due to the lack of turbulent cool flame experiments, which has made model validation difficult. This work addresses these points with a combined experimental and computational investigation of turbulent nonpremixed dimethyl ether cool flames. First, a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner is developed and tested. Turbulent cool flames are studied using the formaldehyde planar laser-induced fluorescence (PLIF), acetone PLIF, and planar Rayleigh scattering techniques. The acetone PLIF signals are converted into mixture fraction values, and measurements of time-averaged temperature are derived from the Rayleigh scattering signals by taking advantage of the similarities between cool flames and unburned mixtures. Second, a Direct Numerical Simulation (DNS) of the CARAT burner is conducted. The flame is shown to be sensitive to thermal boundary conditions such that the stabilization method is dependent on a 10 K change in inlet temperature. Comparisons of first- and second-order statistics of temperature, mixture fraction, and formaldehyde between the DNS and experiment show good agreement. The validity of manifold-based turbulent combustion models in turbulent cool flames is then explored by first studying the effective Lewis numbers in the flame using a local differential diffusion parameter analysis. Turbulent cool flames are shown to have effective unity Lewis number transport even at a much lower Reynolds number than typically required for hot flames. Good agreement is shown for temperature and formaldehyde between DNS conditional means and one-dimensional nonpremixed flame solutions. Additionally, the motion of the cool flame in mixture fraction space is shown to be well described by one-dimensional nonpremixed flame solutions. These points indicate that the fundamental physics of turbulent cool flames are analogous to turbulent hot flames, implying that the same reduced-order manifold-based modeling approaches utilized for turbulent hot flames can be utilized for turbulent cool flames.

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

New capabilities of the MORET 6 Monte Carlo neutron transport code

The MORET code is a simulation tool that solves the transport equation for neutrons using the Monte Carlo method. It allows users to model complex three-dimensional geometrical configurations in a user-friendly way. New features have been introduced to extend the application field of MORET beyond the usual criticality calculations for which it has been initially designed. The most important change is the addition of an analog fixed source mode which allows studies of systems of any reactivity combined with very flexible outputs. Other useful improvements have been added concerning the geometric part, the fission matrix, the multigroup sensitivity coefficients and the outputs. This paper presents an overview of these new features. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2017

Water Hammer Simulations of MMH Propellant - New Capability Demonstration of the Generalized Fluid Flow Simulation Program

Fluid Transient analysis is important for the design of spacecraft propulsion system to ensure structural stability of the system in the event of sudden closing or opening of the valve. Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/MSFC is capable of simulating pressure surge due to sudden opening or closing of valve when thermodynamic properties of real fluid are available for the entire range of simulation. Specifically GFSSP needs an accurate representation of pressure-density relationship in order to predict pressure surge during a fluid transient. Unfortunately, the available thermodynamic property programs such as REFPROP, GASP or GASPAK does not provide the thermodynamic properties of Monomethylhydrazine (MMH). This paper will illustrate the process used for building a customized table of properties of state variables from available properties and speed of sound that is required by GFSSP for simulation. Good agreement was found between the simulations and measured data. This method can be adopted for modeling flow networks and systems with other fluids whose properties are not known in detail in order to obtain general technical insight. Rigorous code validation of this approach will be done and reported at a future date.

Burkhardt, Z.↗
DOE OSTI2018

Reconstruction of phonon relaxation times from systems featuring interfaces with unknown properties

We present a method for reconstructing the phonon relaxation-time function $τ_ω = τ(ω)$ (including polarization) and associated phonon free-path distribution from thermal spectroscopy data for systems featuring interfaces with unknown properties. Our method does not rely on the effective thermal-conductivity approximation or a particular physical model of the interface behavior. The reconstruction is formulated as an optimization problem in which the relaxation times are determined as functions of frequency by minimizing the discrepancy between the experimentally measured temperature profiles and solutions of the Boltzmann transport equation for the same system. Interface properties such as transmissivities are included as unknowns in the optimization; however, because for the thermal spectroscopy problems considered here the reconstruction is not very sensitive to the interface properties, the transmissivities are only approximately reconstructed and can be considered as byproducts of the calculation whose primary objective is the accurate determination of the relaxation times. The proposed method is validated using synthetic experimental data obtained from Monte Carlo solutions of the Boltzmann transport equation. The method is shown to remain robust in the presence of uncertainty (noise) in the measurement.

42 ENGINEERING↗
DOE OSTI2017

Multiresolution molecular mechanics: Surface effects in nanoscale materials

Surface effects have been observed to contribute significantly to the mechanical response of nanoscale structures. The newly proposed energy-based coarse-grained atomistic method Multiresolution Molecular Mechanics (MMM) (Yang, To (2015), ) is applied to capture surface effect for nanosized structures by designing a surface summation rule SR{sup S} within the framework of MMM. Combined with previously proposed bulk summation rule SR{sup B}, the MMM summation rule SR{sup MMM} is completed. SR{sup S} and SR{sup B} are consistently formed within SR{sup MMM} for general finite element shape functions. Analogous to quadrature rules in finite element method (FEM), the key idea to the good performance of SR{sup MMM} lies in that the order or distribution of energy for coarse-grained atomistic model is mathematically derived such that the number, position and weight of quadrature-type (sampling) atoms can be determined. Mathematically, the derived energy distribution of surface area is different from that of bulk region. Physically, the difference is due to the fact that surface atoms lack neighboring bonding. As such, SR{sup S} and SR{sup B} are employed for surface and bulk domains, respectively. Two- and three-dimensional numerical examples using the respective 4-node bilinear quadrilateral, 8-node quadratic quadrilateral and 8-node hexahedral meshes are employed to verify and validate the proposed approach. It is shown that MMM with SR{sup MMM} accurately captures corner, edge and surface effects with less 0.3% degrees of freedom of the original atomistic system, compared against full atomistic simulation. The effectiveness of SR{sup MMM} with respect to high order element is also demonstrated by employing the 8-node quadratic quadrilateral to solve a beam bending problem considering surface effect. In addition, the introduced sampling error with SR{sup MMM} that is analogous to numerical integration error with quadrature rule in FEM is very small. - Highlights: • Surface effect captured by Multiresolution Molecular Mechanics (MMM) is presented. • A novel surface summation rule within the framework of MMM is proposed. • Surface, corner and edges effects are accuterly captured in two and three dimension. • MMM with less 0.3% degrees of freedom of atomistics reproduces atomistic results.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Flexibility Estimation and Control of Thermostatically Controlled Loads with Lock Time for Regulation Service

A virtual battery model is a simple and general method to quantify aggregate flexibility from thermostatically controlled loads (TCLs), enabling grid operators to effectively coordinate a large number of flexible building loads with supplyside resources in power systems. Lockout controls are designed to avoid wear and tear resulting from short-cycling of hardware. The lock on/off time could significantly affect aggregate flexibility from TCLs to provide ancillary services and may even fail control algorithms designed without considering the lock time constraints. This paper focuses on flexibility estimation and control design for TCLs with lock time constraints to provide frequency regulation service. We first investigate the potential impacts of lock time on TCLs’ aggregate flexibility and control performance. Both control-dependent and control-independent power bounds are derived, based on either previous TCL switching operations or regulation signals. While the control-dependent method provides aggregate flexibility for a given control method, the control-independent method calculates the theoretical maximum of power bounds. Two control algorithms are proposed to better distribute flexibility over time and thereby improve signal tracking performance. The proposed methods are illustrated and validated through simulations.

Wang, Peng↗