Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MATLAB”

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 649 records · Page 36

Performance evaluation of photovoltaic/thermal–HDH desalination system

The efficiency of photovoltaic (PV) panel drops with increase in cell temperature. The temperature of the PV panel can be controlled with various cooling techniques. In the proposed work the PV panel is cooled by circulating water and the recovered heat energy is used to run a humidification and dehumidification desalination to produce distilled water from sea water (or) brackish water. This work deals with a detailed analysis of performance of combined power and desalination (Photovoltaic/Thermal–Humidification and Dehumidification) system. A mathematical model of PV/thermal–humidification dehumidification plant was developed and simulations were carried out in MATLAB environment. The performance of photovoltaic/ thermal desalination (Photovoltaic/Thermal–Humidification and Dehumidification) system was investigated under various solar radiation levels (800–1000 W/m{sup 2}). For each solar radiation level the effect of mass flow rate of coolant water (30–110 kg/h) on water outlet temperature, PV efficiency, PVT thermal efficiency, distilled water production, and plant efficiency was studied. Results show that under each solar radiation level increasing coolant flow rate increases efficiency of PV panel and reduces the plant efficiency. The highest PV efficiency (16.598%) was reached under 800 W/m{sup 2} at mass flow rate of 110 kg/h and the highest plant efficiency (43.15%) was reached under 800 W/m{sup 2} at a mass flow rate of 30 kg/h. The maximum amount of distilled water production rate (0.82 L/h) was reached under 1000 W/m{sup 2} at water mass flow rate of 30 kg/h.

14 SOLAR ENERGY↗

Fault detection in PV systems

A practical fault detection approach for PV systems intended for online implementation is developed. The fault detection model here is built using artificial neural network. initially the photovoltaic system is simulated using MATLAB software and output power is collected for various combinations of irradiance and temperature. Data is first collected for normal operating condition and then four types of faults are simulated and data are collected for faulty conditions. Four faults are considered here and they are: Line to Line faults with a small voltage difference, Line to line faults with a large voltage difference, degradation fault and open-circuit fault. This data is then used to train the neural network and to develop the fault detection model. The fault detection model takes irradiance, temperature and power as the input and accurately gives the type of fault in the PV system as the output. This system is a generalized one as any PV module datasheet can be used to simulate the Photovoltaic system and also this fault detection system can be implemented online with the use of data acquisition system.

14 SOLAR ENERGY↗

Time Delay of Wide Area Damping Control in Urban Power Grid: Model-Based Analysis and Data-Driven Compensation

Due to the rapid development of economies, large urban cities consume an increasing amount of energy and have a higher requirement for power quality. Voltage source converter based high voltage direct current (VSC-HVDC) is a promising device to transmit clean power from remote regions to urban power systems, while also providing wide area damping control (WADC) for frequency stabilization. However, the time-delay naturally existing in the VSC-HVDC system may degrade the performance of WADC and even result in instability. To address this issue, this paper develops a time-delay correction control strategy for VSC-HVDC damping control in urban power grids. First, a small signal model of WADC is built to analyze the negative impacts of time delay. Then, a data-driven approach is proposed to compensate for the inherent time delay in VSC-HVDC damping control. The extensive training data will be generated under various disturbances. After offline training, the long short-term memory network (LSTM) can be implemented online to predict the actual frequency deviation based on real-time measurements. Finally, the proposed method is validated through MATLAB-Simulink in a two-area four-machine system. The results indicate that the data-driven compensation has a strong generalization ability for random delay time constants and can improve the performance of WADC significantly.

24 POWER TRANSMISSION AND DISTRIBUTION↗

MetaboTools: A comprehensive toolbox for analysis of genome-scale metabolic models

Metabolomic data sets provide a direct read-out of cellular phenotypes and are increasingly generated to study biological questions. Previous work, by us and others, revealed the potential of analyzing extracellular metabolomic data in the context of the metabolic model using constraint-based modeling. With the MetaboTools, we make our methods available to the broader scientific community. The MetaboTools consist of a protocol, a toolbox, and tutorials of two use cases. The protocol describes, in a step-wise manner, the workflow of data integration, and computational analysis. The MetaboTools comprise the Matlab code required to complete the workflow described in the protocol. Tutorials explain the computational steps for integration of two different data sets and demonstrate a comprehensive set of methods for the computational analysis of metabolic models and stratification thereof into different phenotypes. The presented workflow supports integrative analysis of multiple omics data sets. Importantly, all analysis tools can be applied to metabolic models without performing the entire workflow. Taken together, the MetaboTools constitute a comprehensive guide to the intra-model analysis of extracellular metabolomic data from microbial, plant, or human cells. In conclusion, this computational modeling resource offers a broad set of computational analysis tools for a wide biomedical and non-biomedical research community.

59 BASIC BIOLOGICAL SCIENCES↗

A One-Dimensional Stack Model for Redox Flow Battery Analysis and Operation

Current redox flow battery (RFB) stack models are not particularly conducive to accurate yet high-throughput studies of stack operation and design. To facilitate system-level analysis, we have developed a one-dimensional RFB stack model through the combination of a one-dimensional Newman-type cell model and a resistor-network to evaluate contributions from shunt currents within the stack. Inclusion of hydraulic losses and membrane crossover enables constrained optimization of system performance and allows users to make recommendations for operating flow rate, current densities, and cell design given a subset of electrolyte and electrode properties. Over the range of experimental conditions explored, shunt current losses remain small, but mass-transfer losses quickly become prohibitive at high current densities. Attempting to offset mass-transfer losses with high flow rates reduces system efficiency due to the increase in pressure drop through the porous electrode. The development of this stack model application, along with the availability of the source MATLAB code, allows for facile approximation of the upper limits of performance with limited empiricism. This work primarily presents a readily adaptable tool to enable researchers to perform either front-end performance estimates based on fundamental material properties or to benchmark their experimental results.

Barton, John L. (ORCID:0000000225729127)↗

An Optimized H5 Hysteresis Current Control with Clamped Diodes in Transformer-Less Grid-PV Inverter

With the rise of renewable energy penetration in the grid, photovoltaic (PV) panels are connected to the grid via inverters to supply solar energy. Transformer-less grid-tied PV inverters are gaining popularity because of their improved efficiency, reduced size, and lower costs. However, they can induce a path for leakage currents between the PV and the grid due to the absence of galvanic isolation. This leads to serious electromagnetic interference, loss in efficiency, and safety concerns. The leakage current is primarily influenced by the nature of the common mode voltage (CMV), which is determined by the switching techniques of the inverter. In this paper, a novel inverter topology of Hysteresis Controlled H5 with Two Clamping Diodes (HCH5-D2) is derived. The HCH5-D2 topology helps decouple the AC part (Grid) and DC part (PV) during the freewheeling period to make the CMV constant, thereby reducing the leakage current. Additionally, the extra diodes help reduce voltage spikes generated during the freewheeling period and maintain the CMV at a constant value. Finally, a 2.2 kW grid-connected single-phase HCH5-D2 PV inverter system’s MATLAB simulation is presented, showing better results compared to a traditional H4 inverter.

42 ENGINEERING↗

Per-Phase and 3-Phase Optimal Coordination of Directional Overcurrent Relays Using Genetic Algorithm

Penetration of the power grid by renewable energy sources, distributed storage, and distributed generators is becoming increasingly common. Increased utilization of these distributed energy resources (DERs) has given rise to additional protection coordination concerns, particularly where they are utilized in an unbalanced manner or where loading among phases is unbalanced. Digital relays such as the SEL-751 (produced by Schweitzer Engineering Laboratories, Pullman, WA, USA) series have the capability of being set on a per-phase basis. This capability is underutilized in common practice. Additionally, in optimization algorithms for determining relay settings, the time-overcurrent characteristics (TOCs) of relays are generally not treated as variables and are assigned before running the optimization algorithm. In this paper, TOC options themselves are treated as discrete variables to be considered in the optimization algorithm. A mixed integer nonlinear programming problem (MINLP) is set up where the goal is to minimize relay operating times. A genetic algorithm (GA) approach is implemented in MATLAB where two cases are considered. In the first case, the TOC and Time dial setting (TDS) of each relay is set on a three-phase basis. In the second case, per-phase settings are considered. Relay TDSs and TOCs are both considered as simultaneous discrete control variables. Despite the stochastic nature of using per-phase settings for unbalanced systems is found to generally allow for shorter operating times. However, for relatively balanced systems, it is best to use three-phase settings if computation time is of importance.

Matthews, Ronald C.↗

Proactive Frequency Stability Scheme: A Distributed Framework Based on Particle Filters and Synchrophasors

The reactive nature of traditional under-frequency load shedding schemes can lead to delayed response and unnecessary loss of load. This work presents a proactive framework for power system frequency stability. Bayesian filters and synchrophasors are leveraged to produce predictions after disturbances are detected. By being able to estimate the future state of frequency corrective actions can be taken before the system reaches a critical condition. This proactive approach makes it possible to optimize the response to a disturbance, which results in a decrease in the amount of compensation utilized. The framework is tested via Matlab simulations based on Kundur’s Two-Area System, and the IEEE 14-Bus System. Performance metrics are provided and evaluated against other contemporary solutions found in literature. During testing this framework outperformed other solutions by drastically reducing the amount of load dropped during compensation.

42 ENGINEERING↗

GA-Based Voltage Optimization of Distribution Feeder with High-Penetration of DERs Using Megawatt-Scale Units

In this paper, genetic algorithm (GA)-based voltage optimization of a modified IEEE-34 node distribution feeder with high penetration of distributed energy resources (DERs) is proposed using two megawatt-scale reactive power sources. Traditional voltage support units present in distribution grids are not suitable for DER-rich feeders, while voltage support using small-scale DERs present in the feeder requires considerable communication effort to reach a global solution. In this work, two megawatt-scale units are placed to improve the voltage profile across the IEEE 34-node feeder, which has been modified to include several PV units and an energy storage unit. The megawatt-scale units are optimized using GA for fast and accurate operation. The performance of the proposed scheme is verified using simulation results with a multi-platform setup where the modified IEEE-34 node feeder is modeled in OpenDSS while the GA optimization scheme is programmed in MATLAB.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Method for Spatiotemporal Solar Power Profile Estimation for a Proposed U.S.–Caribbean–South America Super Grid under Hurricanes

Solar photovoltaic (PV) generation technology stands out as a scalable and cost-effective solution to enable the transition toward decarbonization. However, PV solar output, beyond the daily solar irradiance variability and unavailability during nights, is very sensitive to weather events like hurricanes. Hurricanes nucleate massive amounts of clouds around their centers, shading hundreds of kilometers in their path, reducing PV power output. This research proposes a spatiotemporal method, implemented in MATLAB R2023b coding, to estimate the shading effect of hurricanes over a wide distribution of PV solar plants connected to a high-voltage power infrastructure called the U.S.–Caribbean–South America super grid. The complete interconnection of the U.S., the Caribbean, and South America results in the lowest power valley levels, i.e., an overall percentual reduction in PV power output caused by hurricane shading. The simulations assess the impact of hurricanes in 10 synthetic trajectories spanning from Texas to Florida. The Caribbean would also experience lower power valleys with expanded interconnectivity schemes. The U.S.–Caribbean–South America super grid reduces Caribbean variability from 37.8% to 8.9% in the case study. The proposed spatiotemporal method for PV power profile estimation is a valuable tool for future solar power generation expansion, transmission planning, and system design considering the impact of hurricanes.

14 SOLAR ENERGY↗

Multibody for Everybody (M4E): A Symbolic Dynamics Modeling Tool with Applications in Simulation, Control, and Optimization

Developing the analytical model of a multibody system is often the initial step in control and optimization. The analytical model (equations of motion) describes a system’s time evolution under specified forcing conditions. Although developing these equations is easy for simple systems, this process becomes more complex for systems composed of multiple bodies. Deriving equations of motion for complex multibody systems requires specialized expertise in multibody dynamics, is time-consuming, and is susceptible to error. To address this issue, this paper presents an open-source, easy-to-use, systematic framework to derive symbolic equations of motion in both Python and MATLAB using the joint coordinate formulation. This formulation results in a set of ordinary differential equations that use the minimum set of coordinates needed to model a system. The symbolic representation provides better insight into the influence of design parameters on system performance, facilitates sensitivity analysis and parameter studies, and supports direct implementation of control and optimization routines. The tool enables numerical simulation for specified parameter sets, is modular for straightforward integration with other tools and libraries, and allows incorporation of hydrodynamics, mooring, and other external forces. The result is a reproducible, extensible pipeline for modeling, simulation, and design of complex multibody systems. The proposed tool is versatile and can be applied to domains such as robotics, control, and design. In addition, we integrated external libraries that provide capabilities for modeling offshore systems such as underwater robots and marine energy converters.

16 TIDAL AND WAVE POWER↗

Picometer-Precision Atomic Position Tracking through Electron Microscopy

The modern aberration-corrected scanning electron microscopes have successfully achieved direct visualization of atomic columns with sub-angstrom resolution. With this significant progress, advanced image quantification and analysis is still at its early stages. In this work, we present the complete pathway for the metrology of atomic resolution STEM images. This includes: 1) tips for acquiring high-quality STEM images; 2) denoising and drift-correction for enhancing measurement accuracy; 3) obtaining initial atom positions; 4) indexing the atoms based on unit cell vectors; 4) quantifying the atom column positions with either 2D-gaussian single peak fitting or 5) multi-peak fitting routines for slightly overlapping atomic columns; 6) quantification of lattice distortion/strain within the crystal structures or at the defects/interfaces where the lattice periodicity is disrupted, and 7) some common methods to visualize and present the analysis. Furthermore, a simple self-developed free MATLAB app (EASY-STEM) with a graphical user interface (GUI) will be introduced that can help with the analysis of STEM images without the need for writing dedicated analysis code or software. The advanced data analysis methods presented here can be applied for the local quantification of defect relaxations, local structural distortions, local phase transformations, and non-centrosymmetry in a wide range of materials.

47 OTHER INSTRUMENTATION↗

Optimal Matched Filter in the Low-number Count Poisson Noise Regime and Implications for X-Ray Source Detection

Detection of templates (e.g., sources) embedded in low-number count Poisson noise is a common problem in astrophysics. Examples include source detection in X-ray images, γ-rays, UV, neutrinos, and search for clusters of galaxies and stellar streams. However, the solutions in the X-ray-related literature are sub-optimal in some cases by considerable factors. Using the lemma of Neyman–Pearson, we derive the optimal statistics for template detection in the presence of Poisson noise. We demonstrate that, for known template shape (e.g., point sources), this method provides higher completeness, for a fixed false-alarm probability value, compared with filtering the image with the point-spread function (PSF). In turn, we find that filtering by the PSF is better than filtering the image using the Mexican-hat wavelet (used by wavdetect). For some background levels, our method improves the sensitivity of source detection by more than a factor of two over the popular Mexican-hat wavelet filtering. This filtering technique can also be used for fast PSF photometry and flare detection; it is efficient and straightforward to implement. We provide an implementation in MATLAB. The development of a complete code that works on real data, including the complexities of background subtraction and PSF variations, is deferred for future publication.

79 ASTRONOMY AND ASTROPHYSICS↗

An Easy Algorithm to Generate Colored Noise Sequences

The need for accurate and efficient simulation of the noise background arises in statistical significance tests for periodic signals buried in colored noise. This paper discusses techniques for generating colored-noise sequences which simulate processes with a given spectral density. Matlab routines based on the fast fractional difference algorithm are presented. These routines can create various stochastic models (e.g., first-order autoregressive (AR(1)), power law (PL), autoregressive fractionally integrated moving average ((1, d, 0)), and generalized Gauss Markov (GGM)) that serve as possible candidate null hypotheses to test against in various scenarios. Allan variance and power spectral density (PSD) show that our algorithms are accurate and efficient, and can be easily implemented for stationary noise models and non-stationary PLs with spectral indices up to 2. Our algorithms can also be extended to produce a non-stationary PL with a range-limited steeper PSD (spectral indices up to about 4) using a GGM approximation with a proper break frequency. The red-noise leakage effect on the periodogram is further discussed. The result shows that our GGM approximation has a potential for alleviating the red-noise leakage in the PSD estimates. The Timmer and Koenig procedure is also included for comparison.

79 ASTRONOMY AND ASTROPHYSICS↗

An accurate and efficient algorithm for detection of radio bursts with an unknown dispersion measure, for single-dish telescopes and interferometers

Astronomical radio signals are subjected to phase dispersion while traveling through the interstellar medium. To optimally detect a short-duration signal within a frequency band, we have to precisely compensate for the unknown pulse dispersion, which is a computationally demanding task. We present the “fast dispersion measure transform” algorithm for optimal detection of such signals. Our algorithm has a low theoretical complexity of 2N{sub f}N{sub t}+N{sub t}N{sub Δ}log{sub 2}(N{sub f}), where N{sub f}, N{sub t}, and N{sub Δ} are the numbers of frequency bins, time bins, and dispersion measure bins, respectively. Unlike previously suggested fast algorithms, our algorithm conserves the sensitivity of brute-force dedispersion. Our tests indicate that this algorithm, running on a standard desktop computer and implemented in a high-level programming language, is already faster than the state-of-the-art dedispersion codes running on graphical processing units (GPUs). We also present a variant of the algorithm that can be efficiently implemented on GPUs. The latter algorithm’s computation and data-transport requirements are similar to those of a two-dimensional fast Fourier transform, indicating that incoherent dedispersion can now be considered a nonissue while planning future surveys. We further present a fast algorithm for sensitive detection of pulses shorter than the dispersive smearing limits of incoherent dedispersion. In typical cases, this algorithm is orders of magnitude faster than enumerating dispersion measures and coherently dedispersing by convolution. We analyze the computational complexity of pulsed signal searches by radio interferometers. We conclude that, using our suggested algorithms, maximally sensitive blind searches for dispersed pulses are feasible using existing facilities. We provide an implementation of these algorithms in Python and MATLAB.

79 ASTRONOMY AND ASTROPHYSICS↗

How to COAAD Images. II. A Coaddition Image that is Optimal for Any Purpose in the Background-dominated Noise Limit

Image coaddition is one of the most basic operations that astronomers perform. In Paper I, we presented the optimal ways to coadd images in order to detect faint sources and to perform flux measurements under the assumption that the noise is approximately Gaussian. Here, we build on these results and derive from first principles a coaddition technique that is optimal for any hypothesis testing and measurement (e.g., source detection, flux or shape measurements, and star/galaxy separation), in the background-noise-dominated case. This method has several important properties. The pixels of the resulting coadded image are uncorrelated. This image preserves all the information (from the original individual images) on all spatial frequencies. Any hypothesis testing or measurement that can be done on all the individual images simultaneously, can be done on the coadded image without any loss of information. The PSF of this image is typically as narrow, or narrower than the PSF of the best image in the ensemble. Moreover, this image is practically indistinguishable from a regular single image, meaning that any code that measures any property on a regular astronomical image can be applied to it unchanged. In particular, the optimal source detection statistic derived in Paper I is reproduced by matched filtering this image with its own PSF. This coaddition process, which we call proper coaddition, can be understood as the maximum signal-to-noise ratio measurement of the Fourier transform of the image, weighted in such a way that the noise in the entire Fourier domain is of equal variance. This method has important implications for multi-epoch seeing-limited deep surveys, weak lensing galaxy shape measurements, and diffraction-limited imaging via speckle observations. The last topic will be covered in depth in future papers. We provide an implementation of this algorithm in MATLAB.

79 ASTRONOMY AND ASTROPHYSICS↗