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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulation and power quality analysis of a Loose-Coupled bipolar DC microgrid in an office building

With distributed generation and battery storage technologies thriving in microgrids, the use of direct current (DC) microgrids in the building sector offers multiple advantages in energy efficiency and power quality compared with alternating current (AC) systems. This study developed a new concept of a loose-coupled bipolar DC building power system. In this work, the concept was used to design a real-world office building in Shenzhen, China. A power system model was developed to study the stability and control of the DC power system and to verify DC power quality. The design and modeling of the DC power control system is discussed in detail. The study developed a few common fault scenarios in DC building microgrids that were simulated in the MATLAB-Simulink environment to validate the design of a loose-coupled bipolar DC system. The results indicate that the proposed loose-coupled bipolar DC system schema, when implemented with proper control algorithms, can achieve good fault-tolerant performance with reliable power quality, even during disruptive system events.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Distributed Wind-Hybrid Microgrids with Autonomous Controls and Forecasting

Distributed wind-hybrid microgrids have the potential to provide key resilience and economic benefits to both the customers they serve and the utility grids they are connected to. Such microgrids will likely be a key part of the grid of the future, whether connected to large utility grids or linked together in multi-microgrid systems. Through the hybridization of distributed wind and solar photovoltaics, autonomous device-level and system-level controls, battery energy storage systems with smart inverters, and forecasting, these microgrids could maintain local stability and provide grid services - all with renewable power. In the literature, these elements have been considered individually. However, they have not been combined and demonstrated at a high fidelity, which is essential to prove the concept's operation before moving to hardware-in-the-loop and physical demonstrations. In this work, we develop a high-fidelity MATLAB-Simulink model of a real distributed wind-hybrid microgrid that includes all these elements. We demonstrate the microgrid maintaining stability and production in a variety of islanded, grid-connected, and transition scenarios. This includes riding through faults and grid transitions, handling resource variability, and providing grid services. The results demonstrate, at a high fidelity, how distributed wind-hybrid microgrids can operate in an economic and resilient fashion. Finally, we provide recommendations for future research to move advanced distributed wind-hybrid microgrids toward deployment.

ancillary services↗
DOE OSTI2021

Empirical validation and comparison of methodologies to simulate micro and macro-encapsulated PCMs in the building envelope

Thermal Energy Storage (TES) has the potential to shift peak electricity demand. Passive TES is usually implemented in building envelope as micro and macro encapsulated phase change materials (PCM) to shift electric energy demand and therefore requires careful heat transfer analysis. Whole building energy modelling with simplified heat transfer analysis has become extremely important for designers, architects, engineers, and researchers to predict energy performance of buildings. It is important to validate PCM modelling algorithms used in building energy programs to quantify their error and prove their capacity to model different PCM encapsulation types. This study uses data from a microencapsulated PCM and two macroencapsulated PCMs (Bio based PCM and hydrate salts) tested in full-scale using the Advanced Multiscale Building Energy Research (AMBER) Lab located at the Colorado School of Mines and is used to validate a numerical algorithm written in MATLAB language. To approximate the heat transfer through a wall assembly with macroencapsulated PCM pouches, several modelling techniques that can reduce 3D heat transfer characteristics to 1D are explored in this research. A parallel path heat transfer modelling approach is found to give the closest agreement with the experimental data for the pouched PCMs in building envelope applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. Additionally, the scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Coupling thermal energy storage with a thermally anisotropic building envelope for building demand-side management across various US climate conditions

Here, the thermally anisotropic building envelope (TABE) is a novel active building envelope that enhances energy efficiency and thermal comfort in buildings by transferring heat and cold between building envelopes and hydronic loops. When coupled with thermal energy storage (TES) units, the TABE + TES enables the storage of both heat and cold energy captured by the TABE roof or exterior walls. This stored energy can be later released by the TABE floor for indoor heating and cooling, benefiting both the grid and the end user. This paper evaluates the merits of TABE + TES for building demand-side management across various US climate conditions, focusing on peak load shaving, annual energy savings, and cost savings under time-of-use (TOU) electric rate schedules. Simulations were conducted by integrating time-of-day–informed, rule-based control strategies in MATLAB, TABE components and TES units in COMSOL Multiphysics, and whole-building energy analysis in EnergyPlus. A case study using the US Department of Energy’s prototype single-family detached house model in Birmingham, Alabama; Los Angeles, California; Oak Ridge, Tennessee; and Denver, Colorado, showed that the TABE + TES system achieved (1) 70 % peak load shaving in Los Angeles and Denver and 20 % in Birmingham and Oak Ridge; (2) significant peak electricity savings of 351–497 kWh, reducing peak energy consumption by 38 %–78 %; and (3) annual heating cost savings of 0.79 $\$$/m2–1.17 $\$$/m 2 and cooling cost savings of 0.60 $\$$/m 2 –1.17 $\$$/m 2 using a normal utility rate or low-TOU rate. The benefits of employing the TABE + TES system are even more significant under high TOU rates.

25 ENERGY STORAGE↗
DOE OSTI2025

Dynamic modelling and control strategy of a temperature-driven metal hydride cooling system for buildings

A temperature-driven coupled metal hydride (MH) based thermal energy storage (TES) system can allow to shave and shift the peak energy demand in buildings. The high energy density and long-term (seasonal) energy storage capability are its major advantages over other energy storage methods. The dynamic nature of the MH operation, however, requires controlled hydrogen transfer between the coupled MHs at a rate needed to meet the building's transient load. While temperature-driven MH systems are studied in the literature, their application in buildings and control are scarcely reported. Here, this paper presents a control-based dynamic modeling of the temperature-driven coupled MH-TES system for building cooling applications. The dynamic model is developed in MATLAB(R) Simulink environment, considering the thermodynamic and kinetic behaviors of the MH systems. Based on a preliminary analysis of a property database of over 337 hydrides, we select around 1600 MH pairs suitable for building cooling applications. Each of these MH pairs is studied for their performance using the dynamic model, and among all, Zr 0.76 Ti 0.24 Ni 1.16 Mn 0.63 V 0.14 Fe 0.18 -Ti 0.85 Zr 0.15 Cr 1.2 Mn 0.8 MH pair showed fast dynamics along with high coefficient of performance (COP) of 0.71. A parametric investigation is performed on this MH pair to understand the effect of operating temperatures. Finally, three proportional-integral (PI) feedback controllers are investigated to regulate the temperature, pressure and mass exchange between the coupled MH pairs. The developed PI controller is sufficiently capable of rejecting the signal noise from the hydrogen flow and internal heat exchange processes with root mean square error of 5.78 W between reference and actual cooling load.

08 HYDROGEN↗
DOE OSTI2024

Predictive model using artificial neural network to design phase change material-based ocean thermal energy harvesting systems for powering uncrewed underwater vehicles

Uncrewed Underwater Vehicles (UUVs) are a major beneficiary of the phase change material (PCM)-based ocean thermal energy harvesting technology for their mission needs. However, this technology relies on different parameters and energy conversion steps that could be critical to the general energy generation efficiency. Sea trials showed that the design performed lower than their laboratory design specifications. This underperformance results from different factors, mainly the UUV’s trajectory, travel time, underwater ocean currents, temperature fluctuations, and biofouling on the heat exchanger due to long term underwater operations. Therefore, there exists a need to continuously monitor the ambient energy harvesting system and predict system performance, for mission planning purposes. Two major parameters influencing the energy harvesting system include the final pressure inside the hydraulic energy storage vessel or accumulator, and the electrical load value. Here, this work focuses on the hydraulic to electric energy conversion system. Therefore, a combination of numerical model and experimental testing is used to develop a predictive model using artificial neural network using MATLAB. After validation with experimental testing, 1000 data samples obtained from the numerical model are used to train the ANN. Compared to the experimental results, the developed ANN model can predict in less than a second the designed benchtop system’s total efficiency with less than 15 percent maximum error range. This predictive model development represents a cost-effective way for optimization and a computational energy efficient mode aboard UUVs for mission planning for deployed UUVs using PCM-based ocean thermal energy harvesting technology.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

Enhancing operational safety: The NSTX-U Shorted Turn Protection (STP) system

The NSTX-U Shorted Turn Protection (STP) system is a vital safety feature designed to protect the tokamak’s coils during operations. It employs a Kalman filter-based algorithm to monitor coil currents and detect faults in real-time. If an anomaly is detected, the system swiftly terminates the tokamak pulse to prevent damage. Developed with Matlab and Simulink, the STP system has been validated through real-time simulations, ensuring its effectiveness for the upcoming upgrades of the NSTX-U facility. In conclusion, this article presents the theory for implementation of the STP algorithm and the results from simulation and real-time implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Relationship between the contact force strength and numerical inaccuracies in piecewise-smooth systems

This work studies the different types of behavior and inaccuracies that can occur when contact is not adequately accounted for in a dynamical system with freeplay, as the strength of the contact stiffness increases. The MATLAB® ode45 time integration solver, with the built-in Event Location capability, is first validated using past experimental data from a forced Duffing oscillator with freeplay. Next, numerical results utilizing event location are compared to results neglecting event location in order to highlight possible numerical errors and effects on multistable dynamical responses. Inaccuracies tend to occur in two different ways. First, neglecting event location can affect the boundaries between basins of attraction. Second, neglecting event location has little effect on the behaviors of the attractor solutions themselves besides merely resembling poorly converged solutions. Errors are less pronounced at the limits of soft or hard contact stiffness. Finally, this study shows the importance of accurately solving piecewise-smooth systems and the existing correlation between the strength of the contact force and possible numerical inaccuracies.

42 ENGINEERING↗
DOE OSTI2022

Screening and evaluation of biomass upgrading strategies for sustainable transportation fuel production with biomass-derived volatile fatty acids

Biomass conversion to fuels and chemicals is crucial to decarbonization, but choosing an advantageous upgrading pathway out of many options is challenging. Rigorously evaluating all candidate pathways (process simulation, product property testing) requires a prohibitive amount of research effort; even simple upgrading schemes have hundreds of possible permutations. We present a method enabling high-throughput screening by approximating upgrading unit operations and drop-in compatibility of products (e.g., fuel properties) and apply it to volatile fatty acid (VFA) conversion to liquid transportation fuels via a MATLAB script, VFA Upgrading to Liquid Transportation fUels Refinery Estimation (VULTURE). VULTURE selects upgrading configurations that maximize fuel blend bio-derived content. We validate VULTURE's approximations through surrogate fuel property testing and process simulation. Techno-economic and life cycle analyses suggest that VFA upgrading processes down-selected by VULTURE are profitable and have low carbon intensities, demonstrating the potential for the strategy to accelerate process development timelines at decreased costs.

09 BIOMASS FUELS↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

RandONets: Shallow networks with random projections for learning linear and nonlinear operators

Deep neural networks have been extensively used for the solution of both the forward and the inverse problem for dynamical systems. However, their implementation necessitates optimizing a high-dimensional space of parameters and hyperparameters. This fact, along with the requirement of substantial computational resources, pose a barrier to achieving high numerical accuracy, but also interpretability. Here, to address the above challenges, we present Random Projection-based Operator Networks (RandONets): shallow networks with random projections and tailor-made numerical analysis methods that learn accurately and fast linear and nonlinear operators. Building on previous works, we prove that RandOnets are universal approximators of linear and nonlinear operators. Due to their simplicity, RandONets provide a one-step transformation of the input space, facilitating interpretability. For the evaluation of their performance, we focus on operators of PDEs. We show, that RandONets outperform by several orders of magnitude, both in terms of numerical approximation accuracy and computational cost, the “vanilla” DeepONets. Hence, we believe that our method will trigger further developments in the field of scientific machine learning, for the development of new ‘’light”schemes that will provide high accuracy while reducing dramatically the computational cost. A MATLAB toolbox for RandONets, including demos, is available on GitHub at https://github.com/GianlucaFabiani/RandONets.

Interpretable machine learning↗
DOE OSTI2021

Quantifying the effect of hydride microstructure on zirconium alloys embrittlement using image analysis

Since hydride microstructure can significantly influence Zr alloy nuclear fuel cladding’s ductility, a new metric has been developed to quantify hydride microstructure in 2D micrographs and relate it to crack propagation. As cladding failure usually results from a hoop stress, this new metric, called the Radial Hydride Continuous Path (RHCP), is based on quantifying the continuity of brittle hydride particles along the radial direction of the cladding tube. Compared to previous metrics, this approach more closely relates to the propensity of a crack to propagate radially through the cladding tube thickness. The RHCP takes into account hydride length, orientation, and connectivity to choose the optimal path for crack propagation through the cladding thickness. The RHCP can therefore be more closely linked to hydride embrittlement of the Zr alloy material, thus creating a relationship between material structure, properties, and performance. The new definition, along with previously proposed metrics such as the Radial Hydride Fraction (RHF), the Hydride Continuity Coefficient (HCC), and the Radial Hydride Continuity Factor (RHCF), have been implemented and automated in MATLAB. These codes are made available with this publication as supplementary materials. These metrics were verified by comparing their predictions of hydride morphology against expected values in simple cases, and the implementation of the new metric was validated by comparing its predictions with manual measurements of hydride microstructure performed on ImageJ. The RHCP was also validated against experimental measurements of fracture behavior and it was shown to correlate with cladding failure better than previous metrics. The information provided by these metrics will help accurately assess cladding integrity during operation, transportation, and storage.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Current interrupt method for calculating the electrochemical impedance in a solid oxide electrolysis stack

Here, in this work the time domain response of Solid Oxide Electrolysis Cells (SOEC) to a current interruption was transformed into the frequency domain using a carrier function Laplace transform, which is fit to the experimental data using a MATLAB Complex Nonlinear Least Squares (CNLS) solver. The hardware implementation, consisting principally of a high-speed switch and a fast-logging Analog to Digital Converter (ADC), was assembled and tested using a calibration module to assess the accuracy, repeatability, and speed of acquisition of the prototype device as compared against a calibrated commercial impedance spectrometer. Additionally, the current interrupt device and commercial FRA were used to acquire the impedance spectra of a four cell SOEC stack with a large, 300 cm 2 , active cell area.

SOEC↗
DOE OSTI2021

An automated procedure built on MTEX for reconstructing deformation twin hierarchies from electron backscattered diffraction datasets of heavily twinned microstructures

Here we present a set of algorithms built on the MTEX and MATLAB graph toolboxes for automatic reconstruction of deformation twin hierarchies from Electron Backscatter Diffraction (EBSD) datasets with a focus on developing methods for heavily twinned microstructures (twin fractions >0.5). The algorithms address key issues arising at large strains, mainly: missing twin relationships, grouping of heavily deformed grain fragments into families of similar orientation originating from a single initial grain, identification of parent fragments for large twin volume fractions, and classification of families having twin relationships with multiple families. To facilitate the development of these algorithms, large-grained ultra-high purity α-Ti deformed in compression along two directions is investigated. Graphs are utilized to handle non-local geometric merging and to represent relationships throughout the reconstruction process. When determining if a grain fragment is from the undeformed microstructure, the combined metrics of the fragment's orientation volume fraction in the initial texture and the directed graph centrality measure of out-closeness (the number of nodes reached in a graph from a given node) are essential. To address automation in reconstructing the sequence of twinning and relating fragments originating from a single grain in the initial microstructure, the twin family tree is formulated as a minimum spanning tree emanating from the initial grain family. A scheme constructing the distances associated with twin relationship comprising the spanning tree is developed, and a novel quasi-directional Prim spanning tree algorithm is used to determine the twin family tree. The procedure is demonstrated to significantly improve the level of automation in reconstructing twin hierarchies in heavily twinned microstructure compared to other methodologies in literature. The procedure can readily be applied to analyses of twinning in metals, as well as provide an approach for routinely extracting twin statistics at larger deformation levels than previously possible. Significantly, the procedure is demonstrated to be capable of identifying third generation twinning in α-Ti microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dissolved gas recovery from water using a sidestream hollow-fiber membrane module: First principles model synthesis and steady-state validation

This paper presents a first-principles model for the recovery of dissolved gases from liquids using a sidestream hollow-fiber membrane module. The model avoids the use of new empirical coefficients, thus providing a parametric understanding of the process behavior for future design and optimization of membrane modules. This type of first-principles model could be particularly useful when gas recovery is beneficial to biological or chemical reactions of interest, such as the acetogenesis reactions in two-stage anaerobic digesters. The steady-state behavior of the model was validated against both new experimental data for the recovery of H 2 , CH 4 and H 2 –CH 4 mixtures from pure water, as well as existing published data. The modeled gas recovery predictions agreed with experimental data to an absolute average error of 13%, and an average R value of 0.98. Parametric analysis of mixed-gas recovery suggests possible key transition points in the composition of the recovered gases. For example, at 40 °C, increasing trans-membrane pressure while keeping hydraulic residence time (HRT) under 0.5 s will result in an increase in the ratio of H 2 to CH 4 recovered. Otherwise, increasing trans-membrane pressure will instead decrease the ratio of H 2 to CH 4 recovered. The model has potential to be extended to transient analysis, but has yet to be validated with transient experimental data. Furthermore, this model was successfully implemented in both Python and MATLAB, and provides valuable insights for future net-energy optimization for anaerobic digestion systems with in-situ gas recovery.

Anaerobic Digestion↗
DOE OSTI2021

SHarD: A beam dynamics simulation code for dielectric laser accelerators based on spatial harmonic field expansion

In order to demonstrate acceleration of electrons to relativistic scales by an on chip dielectric laser accelerator (DLA), a ponderomotive focusing scheme capable of capturing and transporting electrons through nanometer-scale apertures over extended interaction lengths has been proposed. Here we present a Matlab-based numerical code (SHarD) utilizing a spatial harmonic expansion of the fields within the dielectric structure to simulate the evolution of the beam phase space distribution in this scheme. The code can be used to optimize key-parameters for the accelerator performance such as the final energy, transverse spot size evolution and total number of electrons accelerated through currently fabricated structures. Eventually, the simulation model will be applied to inform the phase mask profile to be added to a pulse front tilt drive laser pulse using a liquid crystal mask in the experimental setup being assembled at UCLA Pegasus Laboratory.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗