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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Perovskite- and Dye-Sensitized Solar-Cell Device Databases Auto-generated Using ChemDataExtractor

The number of scientific publications reporting cutting-edge third-generation photovoltaic devices is increasing rapidly, owing to the pressing need to develop renewable-energy technologies that address the climate-change crisis. Consequently, the field could benefit from a central repository where photovoltaic-performance metrics, such as the power-conversion efficiency (η) are recorded. We present two automatically generated databases that contain photovoltaic properties and device material data for dye-sensitized solar cells (DSCs) and perovskite solar cells (PSCs), totalling 660,881 data entries representing 57,678 photovoltaic devices. The databases were generated by applying the text-mining toolkit ChemDataExtractor on a corpus of 25,720 articles. A multi-faceted evaluation, incorporating manual and automatic methods, was applied to ensure that the data contained therein were of the highest quality, with precision metrics ranging from 73.1% to 95.8%. The DSC database contains 475,045 entries representing 41,680 devices, and the PSC database contains 185,836 entries representing 15,818 devices. The databases are available in MongoDB and JSON formats, which can be queried in Python, R, Java and MATLAB for data-driven photovoltaic materials discovery.

14 SOLAR ENERGY↗
DOE OSTI2023

Accurate real space iterative reconstruction (RESIRE) algorithm for tomography

Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Real space iterative reconstruction for vector tomography (RESIRE-V)

Tomography has had an important impact on the physical, biological, and medical sciences. To date, most tomographic applications have been focused on 3D scalar reconstructions. However, in some crucial applications, vector tomography is required to reconstruct 3D vector fields such as the electric and magnetic fields. Over the years, several vector tomography methods have been developed. Here, we present the mathematical foundation and algorithmic implementation of REal Space Iterative REconstruction for Vector tomography, termed RESIRE-V. RESIRE-V uses multiple tilt series of projections and iterates between the projections and a 3D reconstruction. Each iteration consists of a forward step using the Radon transform and a backward step using its transpose, then updates the object via gradient descent. Incorporating with a 3D support constraint, the algorithm iteratively minimizes an error metric, defined as the difference between the measured and calculated projections. The algorithm can also be used to refine the tilt angles and further improve the 3D reconstruction. To validate RESIRE-V, we first apply it to a simulated data set of the 3D magnetization vector field, consisting of two orthogonal tilt series, each with a missing wedge. Our quantitative analysis shows that the three components of the reconstructed magnetization vector field agree well with the ground-truth counterparts. We then use RESIRE-V to reconstruct the 3D magnetization vector field of a ferromagnetic meta-lattice consisting of three tilt series. Our 3D vector reconstruction reveals the existence of topological magnetic defects with positive and negative charges. We expect that RESIRE-V can be incorporated into different imaging modalities as a general vector tomography method. To make the algorithm accessible to a broad user community, we have made our RESIRE-V MATLAB source codes and the data freely available at https://github.com/minhpham0309/RESIRE-V.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Kinetic modelling of an environmentally friendly carbamazepine synthesis via urea and iminostilbene in batch and continuous processes

Accurate kinetic models for reaction systems allow for improved process understanding and greater quality control, which is particularly beneficial as the pharmaceutical industry shifts from batch to continuous manufacturing (CM). In this work, a first principles kinetic model has been developed for the synthesis of carbamazepine (CBZ) from iminostilbene and urea, starting in a batch reactor and subsequently in a continuous flow reactor. An eco-friendly reaction pathway using urea was selected to avoid the toxic reagents that are typically used for synthesis of CBZ. The kinetic parameters determined from batch reactions were utilized in a MATLAB based kinetic model to simulate the yield for the continuous process. Overall, good agreement between the model prediction and corresponding experimental values was observed for the batch and continuous reaction systems within the full factorial design space. However, the model slightly overpredicted the yield of the continuous reaction system for higher conversion values (>60%) since it did not account for the reverse reaction that can occur at the studied reaction conditions. Here, the use of broken order kinetics was compared with whole number orders, and it was determined that the whole number orders resulted in better agreement at all conversion values for the continuous system.

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

Modular horizontal pendulum wave energy converter: Exploring feasibility to power ocean observation applications in the U.S. pacific northwest

Abstract Marine renewable energy as a power source for ocean observation applications has the potential to allow longer deployment operations due to the consistent, higher, and denser energy available from this resource. This additionally could encourage deployments in remote locations where maintenance is costly or resource availability is low if dependent on solar power. More importantly, gaps in spatial data could be filled. This paper examines the feasibility of a modular horizontal pendulum wave energy converter to power National Data Buoy Center's Self‐Contained Ocean Observations Payload (SCOOP) off the coast of Washington State, U.S. The effect on power output was studied when the pendulum's radius arm, mass, and power take‐off damping were varied. Results using Matlab toolbox WEC‐Sim revealed positive correlation between radius arm length and mass to power output, where power maximised for optimal damping values. Seasonal trends in power were not significant, where a 20 kg pendulum mass was needed to meet the SCOOP base power requirement of 5 W throughout the year.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Evaluating Connected and Automated Vehicles in Co-Simulation Environment of Traffic Microsimulation and Vehicle Dynamics

Connected and automated vehicles (CAVs) have the potential to improve many aspects of the current transportation systems such as safety, mobility, and energy efficiency. In order to evaluate the benefits and impacts of a CAV, the CAV control algorithm is typically implemented on vehicles simulated in a traffic microsimulation environment. However, traffic microsimulation usually lacks detailed vehicle and powertrain dynamics, making it challenging to fully understand how a CAV control algorithm will perform and respond on an actual vehicle. Whether the same benefits measured in the simulation will also be observed in real-world remains an open question. One potential approach to fill in this gap is to conduct a co-simulation of traffic microsimulation with detailed vehicle and powertrain dynamics models, often developed in MATLAB Simulink. However, current microsimulation tools such as VISSIM and SUMO do not have a ready-to-use interface for co-simulation with vehicle dynamics and Simulink. Also, even if such an interface exists, it will be tool-specific, making it challenging to shift from one tool to another or test CAV controls in different tools. There are needs for tool-agnostic co-simulation as different microsimulation tools have their pros and cons, and researchers often need to use different tools based on the purposes of the simulation, project needs, and applications. In this work, Flexible Interface for X-in-the-loop Simulation (FIXS) is developed that can support the co-simulation of microsimulation, CAV control algorithm, and vehicle dynamics model in Simulink. Enabled by the FIXS, the benefit and performance of a CAV control algorithm can be better understood with the consideration of vehicle responses and dynamics. The connection to VISSIM and SUMO is handled internally by the interface, and users can easily switch tools by changing a configuration file. The co-simulation capability is demonstrated for a VISSIM eco-approach and departure CAV scenario and a SUMO cooperative merging scenario for both a passenger CAV and a class 8 heavy-duty connected and automated trucks.

Shao, Yunli↗
DOE OSTI2021

Versatile semi-batch apparatus for manometric measurement of gas–solid reaction rates

In this study, we describe the design, construction, operation, and performance of a simple and versatile semi-batch reactor that is especially useful for measurement of gas/solid reaction rates at pressures in the range of 1 mTorr to 1500 Torr. The reactor operates by repeatedly imposing small AC modulations of reactant gas pressure on top of a much larger DC pressure background. Based on the rate of pressure relaxation following each AC pulse, the reaction rate is determined. Our design is characterized by modular construction from off-the-shelf, ultra-high-vacuum-compatible components, which facilitate easy retrofitting and adaptation to a range of experimental conditions. Automated experiment control and data acquisition is accomplished via a custom National Instruments © LabView virtual instrument. Data analysis is automated using a custom series of Mathworks © Matlab scripts. We demonstrate reactor performance through measurements of hydrogenation kinetics for a composite H 2 getter material consisting of 1,4-bis(phenylethynyl)benzene mixed with a palladium/carbon catalyst.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Highly Available Nuclear Power for Mission-Critical Applications

Some energy consumers require power on an anytime, all-year-round basis with a high level of certainty, including defense installations, isolated communities, and some industrial processes. For these customers, interruptions in electricity or heat can mean substantial financial loss or even loss of life. In the absence of grid-scale energy storage, a high level of power availability can be accomplished only through the robustness and redundancy of power generators. The NuScale small modular reactor design is well suited to provide highly available power because of several features related to both the nuclear steam supply system and the overall plant design. In analogy to Redundant Array of Independent Disks (RAID) systems used to provide highly reliable data storage, a NuScale plant can assure sustained power generation by virtue of its Redundant Array of Integral Reactors (RAIR).This paper describes the NuScale RAIR plant features and summarizes the results of a rigorous analysis of RAIR availability as a function of power or, conversely, the RAIR plant output power as a function of power availability. The analysis utilized the Matrix Laboratory code (MATLAB) and included probability distributions for the frequency and duration of module outages due to planned and unplanned events. Further, the study evaluated the impact of implementing turbine bypass rather than module shutdown and using one or more modules to supply house loads in the case of loss of off-site power (LOOP). Availability results are presented for a 12-module RAIR plant with and without turbine bypass enabled during a LOOP and for different possible connections to the off-site power distribution grid and dedicated service loads. Results indicate that a very high level of availability can be achieved at relatively high power output levels, regardless of turbine bypass and dedicated load connection, compared to the operating fleet.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗
DOE OSTI2024

A fast, matrix-based method to perform omnidirectional pressure integration

Abstract Experimentally-measured pressure fields play an important role in understanding many fluid dynamics problems. Unfortunately, pressure fields are difficult to measure directly with non-invasive, spatially resolved diagnostics, and calculations of pressure from velocity have proven sensitive to error in the data. Omnidirectional line integration methods are usually more accurate and robust to these effects as compared to implicit Poisson equations, but have seen slower uptake due to the higher computational and memory costs, particularly in 3D domains. This paper demonstrates how omnidirectional line integration approaches can be converted to a matrix inversion problem. This novel formulation uses an iterative approach so that the boundary conditions are updated each step, preserving the convergence behavior of omnidirectional schemes while also keeping the computational efficiency of Poisson solvers. This method is implemented in Matlab and also as a GPU-accelerated code in CUDA-C++. The behavior of the new method is demonstrated on 2D and 3D synthetic and experimental data. Three-dimensional grid sizes of up to 125 million grid points are tractable with this method, opening exciting opportunities to perform volumetric pressure field estimation from 3D PIV measurements.

42 ENGINEERING↗
DOE OSTI2023

Modeling and analysis of the tritium fuel cycle for ARC- and STEP-class D-T fusion power plants

Abstract The limited tritium resources available for the first fusion power plants (FPPs) make fuel self-sufficiency and tritium inventory minimization leading issues in FPP design. This work builds on the model proposed by Abdou et al (2020 Nucl. Fusion 61 013001), which analyzed the fuel cycle (FC) of a DEMOnstration nuclear FPP-class FPP with a time-dependent system-level model. Here, we use a modified version of their model to analyze the FC of an Affordable, Robust, Compact (ARC)-class tokamak and two versions of a Spherical Tokamak for Energy Production (STEP)-class tokamak. The ARC-class tokamak breeds tritium in a 2LiF + BeF 2 liquid immersion blanket, while the STEP-class tokamak breeds tritium utilizing either a liquid-lithium blanket design or an encapsulated breeding blanket. A time-dependent system-level model is developed in Matlab Simulink ® to simulate the evolution of tritium flows and tritium inventories in the FC. The main goals of this work are to assess tritium self-sufficiency of the ARC- and STEP-class designs and to determine quantitative design requirements that can be used to analyze the adequacy of a proposed FC system. These design requirements are aimed at achieving a low tritium inventory doubling time ( t d ) and a low start-up inventory ( I s t a r t u p ) while keeping the required tritium breeding ratio (TBR r ) as low as possible. We also consider how improvements in FC technology and POs affect TBR r and I s t a r t u p . The model results show that TBR r for ARC- and STEP-class FPPs should be achievable if the tritium burn efficiency (TBE) reaches 0.5%–1% (TBR r < 1.2). This assumes significant, but attainable, improvements over current abilities. However, the model results indicate that an FPP must achieve ambitious performance targets, including FPP availability > 70%, tritium processing time < 4 h, and the implementation of direct internal recycling (DIR). If future research yields major improvements to achievable TBE, it may be possible to achieve tritium self-sufficiency while operating at lower availability and without implementing DIR.

Physics↗
DOE OSTI2020

GIWAXS-SIIRkit: scattering intensity, indexing and refraction calculation toolkit for grazing-incidence wide-angle X-ray scattering of organic materials

Grazing-incidence wide-angle X-ray scattering (GIWAXS) has become an increasingly popular technique for quantitative structural characterization and comparison of thin films. For this purpose, accurate intensity normalization and peak position determination are crucial. At present, few tools exist to estimate the uncertainties of these measurements. Here, a simulation package is introduced called GIWAXS-SIIRkit , where SIIR stands for scattering intensity, indexing and refraction. Furthermore, the package contains several tools that are freely available for download and can be executed in MATLAB. The package includes three functionalities: estimation of the relative scattering intensity and the corresponding uncertainty based on experimental setup and sample dimensions; extraction and indexing of peak positions to approximate the crystal structure of organic materials starting from calibrated GIWAXS patterns; and analysis of the effects of refraction on peak positions. Additionally, each tool is based on a graphical user interface and designed to have a short learning curve. A user guide is provided with detailed usage instruction, tips for adding functionality and customization, and exemplary files.

36 MATERIALS SCIENCE↗
DOE OSTI2020

XtalCAMP: a comprehensive program for the analysis and visualization of scanning Laue X-ray micro-/nanodiffraction data

XtalCAMP is a software package based on the MATLAB platform, which is suitable for, but not limited to, the analysis and visualization of scanning Laue X-ray micro-/nanodiffraction data. The main objective of the software is to provide complementary functionalities to the Laue indexing software packages used at several synchrotron beamlines. Here, the graphical user interfaces allow the easy analysis of characteristic microstructure features, including real-time intensity mapping for a quick examination of phase, grain and defect distribution, 2D color-coded mapping of microstructural properties from the output of other Laue indexing software, crystal orientation visualization, grain boundary characterization based on orientation/misorientation calculation, principal strain/stress analysis, and strain ellipsoid representation, as well as a series of additional toolkits. As an example, XtalCAMP is applied to the microstructural investigation of a solution-heat-treated Ni-based superalloy manufactured using a laser 3D-printing technique, and a deformed natural quartzite from Val Bregaglia in the Central Alps.

36 MATERIALS SCIENCE↗
DOE OSTI2022

AXEAP : a software package for X-ray emission data analysis using unsupervised machine learning

The Argonne X-ray Emission Analysis Package ( AXEAP ) has been developed to calibrate and process X-ray emission spectroscopy (XES) data collected with a two-dimensional (2D) position-sensitive detector. AXEAP is designed to convert a 2D XES image into an XES spectrum in real time using both calculations and unsupervised machine learning. AXEAP is capable of making this transformation at a rate similar to data collection, allowing real-time comparisons during data collection, reducing the amount of data stored from gigabyte-sized image files to kilobyte-sized text files. With a user-friendly interface, AXEAP includes data processing for non-resonant and resonant XES images from multiple edges and elements. AXEAP is written in MATLAB and can run on common operating systems, including Linux, Windows, and MacOS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Framework for X-ray mirror surface shape fitting

For accurate characterization of grazing-incidence X-ray mirrors, we present a comprehensive framework to fit measured surface shapes (either slope or height) of X-ray mirrors used in synchrotron radiation and free-electron laser facilities. We summarize the closed-form expressions of some typical surface shapes of X-ray mirrors including elliptic cylinders, hyperbolic cylinders, ellipsoids, hyperboloids, and diaboloids. This framework is composed of four layers: definition of standard shapes with closed-form expressions, generation of theoretical surface with pose parameters (six degrees of freedom defining an object's position and orientation relative to a coordinate system), parameter optimization with the ability to select which parameters are fit and which are held constant, and the development of user-friendly fitting function wrappers for particular fitting tasks. A few practical fitting examples are demonstrated to verify the effectiveness of the proposed fitting framework. We discuss the physical meanings of the fitting parameters, and provide several examples using the elliptic cylinder and ellipsoid shapes to highlight some features of the framework. Moreover, we provide the presented framework as open-source codes (MATLAB and Python codes available at https://github.com/nsls2omf/xmf) to the community to encourage academic collaboration and further improvements.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Genetic Algorithm-Based Commutation Angle Control for Torque Ripple Mitigation in Switched Reluctance Motor Drives

This work addresses the application of the Genetic Algorithm (GA) technique to optimize the commutation angles of a 2 kW 8/6 switched reluctance machine (SRM). The primary goal is to reduce the well-known drawback of SRMs: the torque ripple. Firstly, the machine was modeled in Matlab /Simulink ® using lookup tables obtained via finite element method (FEM) simulations. Subsequently, the model was used to perform the GA routine aiming to find the optimal phase commutation angles that minimize the torque ripple factor. Notably, the torque performance of the SRM was significantly affected by the commutation angles during the search for the optimal solution. Afterwards, the GA results for four different operation points were verified experimentally through a developed drive platform with digital signal processor-based (DSP) control and an asymmetric bridge converter. As showed by the experiments, the proposed approach was suitable to reduce the torque ripple by more than 50% for one of the evaluated operating points. Furthermore, it was confirmed that the torque ripple mitigation led to acoustic noise improvement.

42 ENGINEERING↗
DOE OSTI2022

A Comprehensive Approach towards Multi-Objective EMI Filter Design Optimization in High-Frequency SiC-Based Motor Drives

With rapid penetration of high-frequency, highdensity power electronics into the industry, controlling the system generated EMI becomes one of the major design challenges. Designing practical optimized EMI filters require simultaneous consideration of multiple aspects - including electromagnetic couplings, magnetic materials, practical component parasitics and their impacts at high-frequencies, and also the overall manufacturing costs. The multi-dimensional nature of the problem often results in overcompensated designs that hurt the system efficiency and power density. This paper focuses on a comprehensive design platform (developed in MATLAB) towards designing high-density, optimized, highly efficient EMI filters by analyzing the generated noise spectrum and evaluating a vast array of potential solutions before outputting an optimized filter solution. The proposed tool includes database-based component selection, optimal topology selection, multistage filter design, overall optimization for volume, mass, cost and total loss. The proposed filter has been thoroughly evaluated in a PSIM simulation environment emulating the hardware in CRD300DA12EXM3, a 300kW three-phase inverter from Wolfspeed with an RL load representing a motor drive application. Furthermore, a 20kHz switching frequency is chosen and impacts of the motor high-frequency characteristics on the generated EMI noise spectrum of the overall system have been analyzed.

42 ENGINEERING↗
DOE OSTI2025

A Grid-Forming Split-Phase Three-Leg Inverter with Unbalanced Loading and Active Power Decoupling

This paper presents a split-phase three-leg inverter that regulates 180-degree phase-shifted AC voltages to support residential house loads under unbalanced loading conditions and decouples double line frequency power pulsation by utilizing the third half-bridge switching leg. The active power decoupling function is added without extra active switching components. The control feasibility and effectiveness are verified through MATLAB Simulink simulation and controller hardware-in-the-loop test.

Kim, Namwon [ORNL] (ORCID:0000000200438489)↗