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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
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 OSTI2016

Coupled modeling of a directly heated tubular solar receiver for supercritical carbon dioxide Brayton cycle: Optical and thermal-fluid evaluation

In single phase performance and appealing thermo-physical properties supercritical carbon dioxide (s-CO 2 ) make a good heat transfer fluid candidate for concentrating solar power (CSP) technologies. The development of a solar receiver capable of delivering s-CO 2 at outlet temperatures ~973 K is required in order to merge CSP and s-CO 2 Brayton cycle technologies. A coupled optical and thermal-fluid modeling effort for a tubular receiver is undertaken to evaluate the direct tubular s-CO 2 receiver’s thermal performance when exposed to a concentrated solar power input of ~0.3–0.5 MW. Ray tracing, using SolTrace, is performed to determine the heat flux profiles on the receiver and computational fluid dynamics (CFD) determines the thermal performance of the receiver under the specified heating conditions. Moreover, an in-house MATLAB code is developed to couple SolTrace and ANSYS Fluent. CFD modeling is performed using ANSYS Fluent to predict the thermal performance of the receiver by evaluating radiation and convection heat loss mechanisms. Understanding the effects of variation in heliostat aiming strategy and flow configurations on the thermal performance of the receiver was achieved through parametric analyses. Finally, a receiver thermal efficiency ~85% was predicted and the surface temperatures were observed to be within the allowable limit for the materials under consideration.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2020

Spectral Mesh-Free Quadrature for Planar Regions Bounded by Rational Parametric Curves

This article presents spectral, mesh-free, Green’s theorem-based numerical quadrature schemes for integrating functions over planar regions bounded by rational parametric curves. Our algorithm proceeds in two steps: (1) We first find intermediate quadrature rules for line integrals along the region’s boundary curves corresponding to Green’s theorem. (2) We then use a high-order quadrature rule to compute the numerical antiderivative of the integrand along a coordinate axis, which is used to evaluate the Green’s theorem line integral. We present two methods to compute the intermediate quadrature rule. The first is spectrally accurate (it converges faster than any algebraic order with respect to number of quadrature points) and is relatively easy to implement, but has no guarantee of polynomial exactness. The second guarantees exactness for polynomial integrands up to a pre-specified degree with an a priori-known number of quadrature points and retains the convergence properties of the first, but is slightly more complicated. The quadrature schemes have applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order meshes, and initialization of multi-material simulations with rational geometry. We compare the quadrature schemes produced using our method to other methods in the literature and show that they are much more efficient both in terms of number of quadrature points and computational time. We provide an open-source implementation of the algorithm in MATLAB.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2017

The introspective may achieve more: Enhancing existing Geoscientific models with native-language emulated structural reflection

Geoscientific models manage myriad and increasingly complex data structures as trans -disciplinary models are integrated. They often incur significant redundancy with cross-cutting tasks. Reflection, the ability of a program to inspect and modify its structure and behavior at runtime, is known as a powerful tool to improve code reusability, abstraction, and separation of concerns. Reflection is rarely adopted in high-performance Geoscientific models, especially with Fortran, where it was previously deemed implausible. Practical constraints of language and legacy often limit us to feather-weight, native-language solutions. We demonstrate the usefulness of a structural-reflection-emulating, dynamically-linked metaObjects, gd . We show real-world examples including data structure self-assembly, effortless save/restart and upgrade to parallel I/O, recursive actions and batch operations. We share gd and a derived module that reproduces MATLAB-like structure in Fortran and C++. We suggest that both a gd representation and a Fortran-native representation are maintained to access the data, each for separate purposes. In conclusion, embracing emulated reflection allows generically-written codes that are highly re-usable across projects.

58 GEOSCIENCES↗
DOE OSTI2018

Alternating Anderson–Richardson method: An efficient alternative to preconditioned Krylov methods for large, sparse linear systems

Here, we present the Alternating Anderson–Richardson (AAR) method: an efficient and scalable alternative to preconditioned Krylov solvers for the solution of large, sparse linear systems on high performance computing platforms. Specifically, we generalize the recently proposed Alternating Anderson–Jacobi (AAJ) method (Pratapa et al., 2016) to include preconditioning, discuss efficient parallel implementation, and provide serial MATLAB and parallel C/C++ implementations. In serial applications to nonsymmetric systems, we find that AAR is comparably robust to GMRES, using the same preconditioning, while often outperforming it in time to solution; and find AAR to be more robust than Bi-CGSTAB for the problems considered. In parallel applications to the Helmholtz and Poisson equations, we find that AAR shows superior strong and weak scaling to GMRES, Bi-CGSTAB, and Conjugate Gradient (CG) methods, using the same preconditioning, with consistently shorter times to solution at larger processor counts. Finally, in massively parallel applications to the Poisson equation, on up to 110,592 processors, we find that AAR shows superior strong and weak scaling to CG, with shorter minimum time to solution. We thus find that AAR offers a robust and efficient alternative to current state-of-the-art solvers, with increasing advantages as the number of processors grows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Simulated impedance of diffusion in porous media

This paper describes the use of a frequency domain, finite-difference scheme to simulate the impedance spectra of diffusion in porous microstructures. We investigate both open and closed systems for a range of ideal geometries, as well as some randomly generated synthetic volumes and tomographically derived microstructural data. In many cases, the spectra deviate significantly from the conventional Warburg-type elements typically used to represent diffusion in equivalent circuit analysis. Furthermore, a key finding is that certain microstructures show multiple peaks in the complex plane, which may be misinterpreted as separate electrochemical processes in real impedance data. This is relevant to battery electrode design as the techniques for nano-scale fabrication become more widespread. This simulation tool is provided as an open-source MatLab application and is freely available online as part of the TauFactor platform.

36 MATERIALS SCIENCE↗
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 OSTI2020

A novel triple power cycle featuring a gas turbine cycle with supercritical carbon dioxide and organic Rankine cycles: Thermoeconomic analysis and optimization

In this study, a novel triple power cycle is proposed where waste heat from a gas turbine cycle is utilized to drive a supercritical carbon dioxide (s-CO 2 ) recompression cycle and a recuperative organic Rankine cycle (ORC) in sequence. A detailed thermoeconomic model is developed and implemented in MATLAB to evaluate the performance of the proposed cycle under different operating conditions. Optimization using a particle swarm optimization (PSO) algorithm is performed to minimize the levelized cost of electricity (LCOE) and determine the optimum design conditions of the cycle. The optimization results show that for a 100 MW cycle, the overall thermal efficiency and LCOE are 0.521 and $\$52.819$/MWh, respectively. The turbine inlet temperature of the gas turbine and s-CO 2 cycles are found as the most influential parameters on the thermoeconomic performance of the triple cycle. The proposed triple cycle shows an excellent waste energy recovery potential and superiority over a very thermodynamically efficient cycle from the literature, which included a gas turbine topping cycle and a complex cascade s-CO 2 power cycle used a bottoming cycle. As a result, the proposed triple cycle shows up to 0.9% points higher efficiency while it has fewer heat exchangers and turbomachinery than the cycle from the literature.

30 DIRECT ENERGY CONVERSION↗
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 OSTI2019

A Generic Modeling and Development Approach for WECC Composite Load Model

The composite load model (CLM) proposed by Western Electricity Coordinating Council (WECC) gains increasing traction in the industry, particularly in North America. At the same time, it is recognized that further improvements in the structure, initialization and aggregation methods are needed. However, lack of an open implementation of WECC CLM becomes a roadblock for many researchers to work on improving it. This paper presents the first open reference implementation of WECC CLM. Individual load components and the CLM are first developed and tested in Matlab, then translated to the high performance computing (HPC), parallel simulation platform GridPACK. Main contributions of the paper include: 1) presenting some important yet undocumented details of modeling and initializing the CLM, particularly for a parallel simulation platform like GridPACK; 2) implementation details of the load components such as the single-phase air conditioner motor; 3) implementing the CLM in a modular and extensible manner. The implementation has been tested at both the component as well as system levels and bench-marked against commercial simulation programs.

Huang, Qiuhua↗
DOE OSTI2018

Techno-economic evaluation of a process for direct conversion of methane to aromatics

Natural gas is widely abundant and relatively inexpensive. As such, large quantities are flared every year. This presents an opportunity to use natural gas as a chemical feedstock. Non-oxidative dehydroaromatization of methane, the major component of natural gas is proposed as an option to monetize natural gas. This study evaluates the economics of the methane to benzene process. A catalytic membrane reactor model written in MATLAB has been directly implemented in an Aspen Plus V9 process model using CAPE-OPEN interoperability. Minimum utility requirements and heat exchanger designs were identified in Aspen Energy Analyzer, and process economics for a variety of feed and product prices were assessed using Aspen Process Economic Analyzer. The results indicate that the venture is profitable (profitability index of 1.17, IRR of 35.5%, NRR of 18.2%). The capital cost was found to be $35,500 per daily standard barrel of benzene. Additionally, the process remains economically attractive so long as either the price of benzene remains above $470/ton, or the price of hydrogen remains above $0.8/kg.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2018

Computational examination of two-phase microchannel heat transfer correlations with conjugate heat spreading

Microchannel flow boiling is an attractive thermal management strategy for ever-growing volumetric heat dissipation demands associated with electronic systems. Due to difficulties related to measurement at the microscales, the majority of researchers have chosen to study relatively simple situations with uniform heat flux. However, many applications involve local hotspots which give rise to highly non-uniform heat flux and temperature gradients due to heat spreading. This necessitates the consideration of conjugate heat transfer for accurate analysis. The current work is aimed at investigating conjugate heat transfer in a two-phase microchannel array. Experimental data was collected on R134a flow boiling heat transfer for very small hydraulic diameter (<100 µm) silicon microchannels with very high heat fluxes (>1 kW cm -2 ) applied via platinum strip heaters with a footprint size of 1 cm × 1 mm. The collected experimental data was then combined with detailed computational modeling utilizing finite element modeling with COMSOL Multiphysics and MATLAB to examine the applicability of five published heat transfer correlations for use in determining local heat transfer coefficients. A two-phase correlation from Agostini and Bontemps, developed for markedly different test parameters, provided the best computational agreement with an RMS temperature difference from experiment of 3.3 °C and a predicted peak perimeter-averaged heat transfer coefficient of 116 MW m -2 K -1 . Modeling confirms the presence of highly non-uniform local heat flux and correspondingly non-uniform local heat transfer coefficient. Finally, the results of this study make clear the need for better micro-scale two-phase correlations developed to predict local heat transfer coefficients at these small scales and high local heat fluxes.

42 ENGINEERING↗
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↗