Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 253 records · Page 14

Combining model-based and model-free methods for stochastic control of distributed energy resources

Modern distribution systems are experiencing a fast transformation with the growing penetration of distributed energy resources (DERs). Along with the economic and environmental benefits of DERs, challenges arise to address the uncertainties caused by their inherent volatility. If properly coordinated, however, DERs have the potential to provide the controllability that grid operators need. Here, we propose a hierarchical control framework that combines the model-based and model-free methods for stochastic DER control in distribution systems. The upper-level scheduler considers a chance-constrained optimal power flow problem (model-based) that schedules DER setpoints to minimize the operational cost and maintain the operating reserve. The lower-level distributed DER controllers absorb real-time disturbances and uncertainties using the extremum seeking control (model-free) to achieve grid objectives. The combination of model-based and model-free methods allows us to take the advantages of both methods to effectively manage the uncertainty in distribution systems. The proposed work is demonstrated on the IEEE 13-node feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Approximate Dynamic Programming With Enhanced Off-Policy Learning for Coordinating Distributed Energy Resources

Herein this paper proposes an innovative approximate dynamic programming (ADP) method for distributed energy resource coordination with the loss of life of battery energy storage system (BESS) explicitly modeled. The dispatch policy is designed to account for both calendrical and cyclical aging effects on BESS, explicitly modeling the impacts of ambient temperature on BESS lifespan. The proposed ADP employs an adaptive critic method and enhanced off-policy deterministic policy gradient (DPG) strategy, addressing the limitations of the on-policy gradient-based ADP approaches, including inadequate exploration, low data usage, and computational complexity. In particular, a customized policy is proposed to guide the algorithm to explore some promising decisions and thereby improve exploration capability and learning efficiency compared to conventional DPG-based learning approaches, which may struggle to find a global optimum due to random noisy action-based exploration or require expert demonstration with extra effort. The proposed method is illustrated using the IEEE 123-node system and compared with the existing ADP methods to prove solution accuracy and demonstrate the effects of incorporating degradation models into control design. Case studies showed that the proposed ADP effectively coordinates DERs with a 10 times smaller optimization gap compared to existing methods, and the incorporation of the BESS life loss model ensures the expected lifespan and results in significant cost savings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Device and methods for determination of molecular weight distributions of polymers and distributions of other polymer properties without physical separation

Devices and methods for determining the cumulative distribution of a polymer property in a reactor without physical separation of reaction subcomponents. The device includes a means of measuring an instantaneous property of the polymers being produced in a reaction vessel a plurality of times during a polymerization reaction as well as a means of determining the corresponding change in polymer concentration in the reaction vessel between measurements of the instantaneous polymer property The device also includes a means of computing a statistical distribution appropriate to the polymer characteristic and applying the statistical distribution to a recently measured instantaneous value of the polymer property so as to have an instantaneous distribution of the polymer property and a means of adding together the instantaneous distributions of the polymer property in order to obtain the cumulative distribution of the polymer property in the reactor.

Reed, Wayne Frederick↗

Device and methods for determination of molecular weight distributions of polymers and distributions of other polymer properties without physical separation

Devices and methods for determining the cumulative distribution of a polymer property in a reactor without physical separation of reaction subcomponents. The device includes a means of measuring an instantaneous property of the polymers being produced in a reaction vessel a plurality of times during a polymerization reaction as well as a means of determining the corresponding change in polymer concentration in the reaction vessel between measurements of the instantaneous polymer property. The device also includes a means of computing a statistical distribution appropriate to the polymer characteristic and applying the statistical distribution to a recently measured instantaneous value of the polymer property so as to have an instantaneous distribution of the polymer property and a means of adding together the instantaneous distributions of the polymer property in order to obtain the cumulative distribution of the polymer property in the reactor.

Reed, Wayne Frederick↗

EV Forecasting-Based Model Predictive Control for Distribution System Congestion Mitigation

The uncoordinated charging of electric vehicles (EVs) in time and space brings congestion issues to the distribution network. This paper proposes an EV charging demand forecasting-based model predictive control (MPC) method for distribution system congestion management. To effectively forecast the time-series EV station charging demand, a hybrid forecasting model that integrates the long short-term memory network (LSTM) and Transformer is proposed. The Transformer-LSTM model is trained using a one-year real historical charging dataset of EV stations to forecast future charging demand in 15-minute intervals. This informs the MPC for distribution network congestion management and minimization of PV curtailment. Numerical results carried out on the modified IEEE 123-bus distribution system demonstrate that the proposed method can effectively resolve line congestion issues through EV smart charging and PV curtailment while outperforming other benchmarks.

ADVANCED PROPULSION SYSTEMS,SOLAR ENERGY↗

A Mixed integer linear programming‐based distributed energy management for networked microgrids considering network operational objectives and constraints

Abstract Mixed integer linear programming (MILP)–based distributed energy management for networked microgrids embedded modern distribution systems is proposed. Considering the diverse ownership of microgrids, distributed energy resources (DERs) that interface directly with utilities and responsive loads, an alternating direction method of multipliers–based distributed framework was formulated for the scheduling of networked microgrids embedded modern distribution systems by adjusting nodal price signals iteratively. In addition, to make the formulated optimization problems resolvable through more accessible and popular MILP solvers, different linearisation techniques were employed to transform the nonlinear terms into linear or mixed integer linear formats. The proposed MILP‐based distributed method preserves all participants' autonomy (e.g., microgrids, DERs that interface directly with utilities and responsive loads), while incentivising them to actively participate in the distribution system operation with price signals. The proposed method is validated with results of numerical simulation using a modern distribution system consisting of multiple networked microgrids, DERs that interface directly with utilities, as well as responsive loads.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Numerical methods and hypoexponential approximations for gamma distributed delay differential equations

Abstract Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs). In this work, we address the lack of appropriate numerical tools for gamma distributed DDEs in two ways. First, we develop a functional continuous Runge–Kutta (FCRK) method to numerically integrate the gamma distributed DDE without resorting to Erlang approximation. We prove the fourth-order convergence of the FCRK method and perform numerical tests to demonstrate the accuracy of the new numerical method. Nevertheless, FCRK methods for infinite delay DDEs are not widely available in existing scientific software packages. As an alternative approach to solving gamma distributed DDEs, we also derive a hypoexponential approximation of the gamma distributed DDE. This hypoexponential approach is a more accurate approximation of the true gamma distributed DDE than the common Erlang approximation but, like the Erlang approximation, can be formulated as a system of ODEs and solved numerically using standard ODE software. Using our FCRK method to provide reference solutions, we show that the common Erlang approximation may produce solutions that are qualitatively different from the underlying gamma distributed DDE. However, the proposed hypoexponential approximations do not have this limitation. Finally, we apply our hypoexponential approximations to perform statistical inference on synthetic epidemiological data to illustrate the utility of the hypoexponential approximation.

97 MATHEMATICS AND COMPUTING↗

Three-dimensional finite-element analysis of chevron-notched fracture specimens

Stress-intensity factors and load-line displacements were calculated for chevron-notched bar and rod fracture specimens using a three-dimensional finite-element analysis. Both specimens were subjected to simulated wedge loading (either uniform applied displacement or uniform applied load). The chevron-notch sides and crack front were assumed to be straight. Crack-length-to-specimen width ratios (a/w) ranged from 0.4 to 0.7. The width-to-thickness ratio (w/B) was 1.45 or 2. The bar specimens had a height-to-width ratio of 0.435 or 0.5. Finite-element models were composed of singularity elements around the crack front and 8-noded isoparametric elements elsewhere. The models had about 11,000 degrees of freedom. Stress-intensity factors were calculated by using a nodal-force method for distribution along the crack front and by using a compliance method for average values. The stress intensity factors and load-line displacements are presented and compared with experimental solutions from the literature. The stress intensity factors and load-line displacements were about 2.5 and 5 percent lower than the reported experimental values, respectively.

Raju, I. S.↗

System and Method for Providing a Climate Data Analytic Services Application Programming Interface Distribution Package

A system, method and computer-readable storage devices for providing a climate data analytic services application programming interface distribution package. The example system can provide various components. The system provides a climate data analytic services application programming interface library that enables software applications running on a client device to invoke the capabilities of a climate data analytic service. The system provides a command-line interface that provides a means of interacting with a climate data analytic service by issuing commands directly to the system's server interface. The system provides sample programs that call on the capabilities of the application programming interface library and can be used as templates for the construction of new client applications. The system can also provide test utilities, build utilities, service integration utilities, and documentation.

Schnase, John L.↗

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗

ARDP Natrium Neutronic Methodology: Argonne Neutronic Assessment of ABR-1000

The Natrium Sodium-cooled Fast Reactor (SFR) concept developed by TerraPower, LLC was selected as one of the advanced reactor designs for demonstration under the Advanced Reactor Demonstration Program Natrium Demo project. In collaboration with TerraPower, the Argonne National Laboratory (ANL) team provide independent reviews of the neutronic methodology being used for the Natrium design. As part of this activity, both the TerraPower and ANL teams agreed to perform independent neutronics analyses of the ABR-1000 reactor problem as a demonstration on the usage of the neutronics methodologies. Comparison of the results from TerraPower and ANL would expose any methodological differences in the modeling approaches. This report describes and presents the Argonne analysis work on the ABR-1000, and will be provided to TerraPower, along with data files, to complete the comparison study. The comparison will cover the methodology to evaluate the reactivity coefficients for a safety analysis, the shutdown margin, and the reconstructed pin power distribution. The methods and models used to calculate the kinetic parameters and the reactivity coefficients (density, Doppler, axial/radial expansion, etc.) are described in detail with recommendations of specific method options and modelling techniques. The control system reactivity worth and the excess reactivity due to temperature defect are calculated for the shutdown margin analysis. The method to evaluate the excess reactivity is discussed in detail. The pin power reconstruction methodology used in SE2ANL, SE2RCT, and DASSH are discussed in this report. The pin power distribution evaluated by different methods are compared and discussed in the report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Energy conservation, using remote thermal scanning

Airborne thermal infrared scans and thermal maps utilized in NASA's energy conservation program have proven to be efficient cost-effective method for identifying heat losses from building roofs and heating system distribution lines. Method employs commercially available equipment in highly developed way.

Bowman, R. L.↗

Use of In Situ Cloud Condensation Nuclei, Extinction, and Aerosol Size Distribution Measurements to Test a Method for Retrieving Cloud Condensation Nuclei Profiles From Surface Measurements

If the aerosol composition and size distribution below cloud are uniform, the vertical profile of cloud condensation nuclei (CCN) concentration can be retrieved entirely from surface measurements of CCN concentration and particle humidification function and surface-based retrievals of relative humidity and aerosol extinction or backscatter. This provides the potential for long-term measurements of CCN concentrations near cloud base. We have used a combination of aircraft, surface in situ, and surface remote sensing measurements to test various aspects of the retrieval scheme. Our analysis leads us to the following conclusions. The retrieval works better for supersaturations of 0.1% than for 1% because CCN concentrations at 0.1% are controlled by the same particles that control extinction and backscatter. If in situ measurements of extinction are used, the retrieval explains a majority of the CCN variance at high supersaturation for at least two and perhaps five of the eight flights examined. The retrieval of the vertical profile of the humidification factor is not the major limitation of the CCN retrieval scheme. Vertical structure in the aerosol size distribution and composition is the dominant source of error in the CCN retrieval, but this vertical structure is difficult to measure from remote sensing at visible wavelengths.

Ghan, Stephen J.↗

The Effect of General Statistical Fiber Misalignment on Predicted Damage Initiation in Composites

A micromechanical method is employed for the prediction of unidirectional composites in which the fiber orientation can possess various statistical misalignment distributions. The method relies on the probability-weighted averaging of the appropriate concentration tensor, which is established by the micromechanical procedure. This approach provides access to the local field quantities throughout the constituents, from which initiation of damage in the composite can be predicted. In contrast, a typical macromechanical procedure can determine the effective composite elastic properties in the presence of statistical fiber misalignment, but cannot provide the local fields. Fully random fiber distribution is presented as a special case using the proposed micromechanical method. Results are given that illustrate the effects of various amounts of fiber misalignment in terms of the standard deviations of in-plane and out-of-plane misalignment angles, where normal distributions have been employed. Damage initiation envelopes, local fields, effective moduli, and strengths are predicted for polymer and ceramic matrix composites with given normal distributions of misalignment angles, as well as fully random fiber orientation.

Micromechanical↗

Method of measuring elastic strain distribution in specimens used for high frequency fatigue testing

A technique using a point contact piezoelectric detector has been developed which permits the displacement amplitude distributions over the entire length of high frequency fatigue specimens to be measured. Elastic strain distributions can be derived with a computer analysis that uses curve smoothing and numerical derivative methods. The equipment and techniques for measuring and processing the data are described, and displacement amplitudes, strain distributions, and magnification factors for straight rods and dumbbell-shaped specimens are obtained and shown to agree with theoretically predicted values. It is shown that the technique presented here can be used to design specimens with unusual displacement amplitudes and strain distributions.

Sirian, C. R.↗

Retrieval of aerosol size distribution moments from multiwavelength particulate extinction data

An iterative method for deriving the particle size distribution from multiwavelength satellite remote sensing data is presented and its use illustrated with an application to SAGE II data. The compression size distribution method (CSD) features interpolation and look-up features for relating the coefficients of a least-squares fit as a function of wavelength to parameters of the compression function that have been obtained by Mie calculations. Results are reported of comparisons between particulate radii predicted for the El Chichon volcano aerosol based on SAGE II data and 14 radii measured with NASA flights through the aerosol with aircraft carrying wire impactors.

Livingston, J. M.↗