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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Joint Optimization for Transport and Bucket Loading Phases of Automated Wheel Loaders

This article investigates optimization of fuel-efficiency and productivity for automated wheel loaders. A control-oriented model for both the transport phase and bucket loading phase is proposed. Here, the vehicle model includes an automatic gear shift schedule that can be incorporated into the optimization problem. Based on the model, the multistage optimization problem is formulated to simultaneously consider all phases of a short cycle with physical constraints. Cycle time and fuel efficiency are used as the weighted performance indexes in a multiobjective cost function. Bucket fill factor is included as a constraint during the bucket loading phase. A nonlinear programming problem is created with collocation using MATLAB and CasADi. The optimization solver IPOPT solves the problem to obtain the optimal state and control trajectories, which can be used as a reference for automated wheel loaders or even as a driver advisory for human-driven wheel loaders.

42 ENGINEERING↗
DOE OSTI2020

Multi-objective Optimization Paradigm for Toroidal Inductors with Spatially Tuned Permeability

Spatially tuning core permeability of an electromagnetic device enables superior performance. A permeability profile can be heuristically selected to improve the flux distribution in a device with a given geometry, but in order to fully leverage the capacity of spatial dependent permeability engineering, the geometry and the permeability should be optimized simultaneously. The work herein presented sets forth a multi-physics design optimization paradigm that includes the permeability profile tuning in the context of an inductor design. This approach enables the determination of Pareto optimal fronts consisting of a set of optimal solutions against competing objectives (e.g. mass and loss) under imposed constraints. To this end, analytical solutions of the heat transfer and electromagnetic formulations are derived for toroidal inductors. The software implemented in Matlab 2018b is available online as an attachment to this paper.

42 ENGINEERING↗
DOE OSTI2017

An Introduction to Goodness of Fit for PMU Parameter Estimation

New results of measurements of phasor-like signals are presented based on our previous work on the topic. In this document an improved estimation method is described. The algorithm (which is realized in MATLAB software) is discussed. We examine the effect of noisy and distorted signals on the Goodness of Fit metric. The estimation method is shown to be performing very well with clean data and with a measurement window as short as a half a cycle and as few as 5 samples per cycle. The Goodness of Fit decreases predictably with added phase noise, and seems to be acceptable even with visible distortion in the signal. While the exact results we obtain are specific to our method of estimation, the Goodness of Fit method could be implemented in any phasor measurement unit.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Neural Networks-Based Inverter Control: Modeling and Adaptive Optimization for Smart Distribution Networks

The optimal voltage control of inverter-based resources, especially under the high penetration of solar photovoltaics, is critical to the stability of the distribution power system. However, the computational complexity as well as the coordinated operation performance of the voltage control optimization in the distribution power system limits the real-time applications. To mitigate this issue, a model-free based adaptive optimal control scheme for the smart inverter is proposed to maximize the active power generation, minimize the power loss, and maintain the bus voltages in smart distribution networks. An inverter-based optimization model for coordinated operation is first established, considering the uncertainties of renewable power generation. Subsequently, by collecting the data and control strategies, the neural networks (NNs) based algorithm is proposed to efficiently predict the best possible control strategy. The main objective of this scheme is to accurately predict candidate optimal solutions with near-negligible feasibility and optimization gaps, with the advantage of avoiding complicated iteration-based numerical algorithms. Thereafter, the co-simulation among OpenDSS, MATLAB, and Python is set up to fully take advantage of the three individual software. Experiments are conducted based on different control parameter characteristics and structures of NNs. Finally, the results reveal that an average mean squared error of 0.013 and 1 ms response time are achieved, which is lower than some state-of-the-art methods.

42 ENGINEERING↗
DOE OSTI2019

Optimization of Air-Gap Profile in Interior Permanent-Magnet Synchronous Motors for Torque Ripple Mitigation

Interior permanent magnet synchronous motors (IPMSMs) have been widely used due to their high efficiency and high power density. Minimization of torque pulsation resulting in vibration and acoustic noise is one of the important design considerations for IPMSMs. In this paper, a grid on/off search method for the rotor profile is proposed to mitigate torque pulsation. Selection of the rotor profile is due to the fact that air-gap is the most sensitive parameter in electric machines wherein changes in flux densities can cause substantial differences in the distribution of forces. A layer comprised of 20 partitions with a 0.1 mm thickness and 3 degrees wide grids have been introduced to the rotor surface for each pole and the possible geometries have been analyzed using finite element method (FEM) in ANSYS Maxwell. An optimal design was found that has the lowest torque ripple with a higher average torque compared to the original design. Genetic algorithm (GA) has also been applied to the method to automate the coupling between Maxwell and Matlab thereby saving the simulation time. Finally, complete structural analysis has been done for both of the original and optimal designs to verify the superiority and feasibility of the proposed design.

42 ENGINEERING↗
DOE OSTI2015

A Method for Determining Bulk Density, Material Density, and Porosity of Melter Feed During Nuclear Waste Vitrification

Glassmelting efficiency largely depends on heat transfer to reacting glass batch (melter feed), which in turn is influenced by the bulk density (ρ b ) and porosity (ϕ) of the reacting feed as functions of temperature ( T ). Neither ρ b ( T ) nor ϕ( T ) functions are readily accessible from direct measurements. For the determination of ρ b , we monitored the profile area of heated feed pellets and calculated the pellet volume using numerical integration. For the determination of ϕ, we measured the material density of feeds quenched at various stages of conversion via pycnometry and then computed the feed density at heat‐treatment temperature using thermal expansion values of basic feed constituents.

Hilliard, Zachary↗
DOE OSTI2012

New Optical Evaluation Approach for Parabolic Trough Collectors: First-Principle OPTical Intercept Calculation

A new analytical method -- First-principle OPTical Intercept Calculation (FirstOPTIC) -- is presented here for optical evaluation of trough collectors. It employs first-principle optical treatment of collector optical error sources and derives analytical mathematical formulae to calculate the intercept factor of a trough collector. A suite of MATLAB code is developed for FirstOPTIC and validated against theoretical/numerical solutions and ray-tracing results. It is shown that FirstOPTIC can provide fast and accurate calculation of intercept factors of trough collectors. The method makes it possible to carry out fast evaluation of trough collectors for design purposes. The FirstOPTIC techniques and analysis may be naturally extended to other types of CSP technologies such as linear-Fresnel collectors and central-receiver towers.

14 SOLAR ENERGY↗
DOE OSTI2020

Metamodeling Time-Temperature Creep Parameters

Abstract There exist many time-temperature parameter (TTP) models for creep rupture prediction of components including the Larson–Miller (LM), Manson–Haferd (MH), Manson–Brown (MB), Orr–Sherby–Dorn (OSD), Manson–Succop (MS), Graham–Walles (GW), Chitty–Duval (CD), Goldhoff–Sherby (GS) models. It remains a challenge to determine which model is “best”, capable of accurate interpolation and physically realistic extrapolation of creep rupture data for a given material. In this study, metamodeling is applied to create a unified TTP metamodel that combines and regresses into twelve TTP models (eight existing and four newly derived). An analysis of the mathematical problems that exist in TTP models is provided. A matlab code is written that can: (1) calibrate the material constants of any of the twelve TTP models (using the metamodel); (2) determine the most suitable stress-parameter function; (3) and report the normalized mean square error (NMSE) of rupture predictions for a given material database. Using the metamodel, and code, a design engineer can make an intelligent selection of the “best” TTP model for creep resistant design. This process is demonstrated using four isotherms of alloy P91 creep rupture data. To assess the influence of material, further validation is performed on alloys Hastelloy X, 304SS, and 316SS. It is determined that the “best” model is dependent on material type and the quality and quantity of available data.

Engineering↗
DOE OSTI2016

Organizing Cells Within Non-Periodic Microarchitectured Materials That Achieve Graded Thermal Expansions

The aim of this paper is to introduce an approach for optimally organizing a variety of different unit cell designs within a large lattice such that the bulk behavior of the lattice exhibits a desired Young’s modulus with a graded change in thermal expansion over its geometry. This lattice, called a graded microarchitectured material, can be sandwiched between two other materials with different thermal expansion coefficients to accommodate their different expansions or contractions caused by changing temperature while achieving a desired uniform stiffness. First, this paper provides the theory necessary to calculate the thermal expansion and Young’s modulus of large multi-material lattices that consist of periodic (i.e., repeating) unit cells of the same design. Then it introduces the theory for calculating the graded thermal expansions of a large multimaterial lattice that consists of non-periodic unit cells of different designs. An approach is then provided for optimally designing and organizing different unit cells within a lattice such that both of its ends achieve the same thermal expansion as the two materials between which the lattice is sandwiched. A MATLAB tool is used to generate images of the undeformed and deformed lattices to verify their behavior and various examples are provided as case studies. The theory provided is also verified and validated using finite element analysis and experimentation.

Young’s modulus↗
DOE OSTI2017

MECHANICAL MOTION RECTIFIER BASED EFFICIENT POWER TAKEOFF FOR OCEAN WAVE ENERGY HARVESTING

This paper proposes a novel mechanical-motion-rectifier (MMR) based power-takeoff (PTO) for ocean wave energy harvesting. The proposed PTO directly converts irregular oscillatory wave motion into regular unidirectional rotation of the generator. It is mainly composed of two ball screws, three bevel gears, two one-way clutches, and a generator. The two one-way clutches and the bevel gears change the bi-directional rotation of the two ball screws into unidirectional ration of the generator. The MMR rectifies the irregular reciprocating motion into unidirectional rotation; similar to the way the electric voltage rectifier regulates an AC voltage. The proposed PTO can be integrated into a heaving point wave energy converter (WEC). The dynamics and modelling of the PTO are presented. The frequency-domain dynamics of the WEC are then formulated for operating condition and control. The power generation capability of the proposed WEC has been evaluated in MATLAB and WAMIT. The simulation results demonstrate that the power generation capability can be improved by using the MMR method.

Yin, Xiuxing↗
DOE OSTI2017

Surface recogniton for cars: A comprehensive approach for neural networks

This paper explores the viability of neural-networkbased classification of ground surface for vehicles. By classifying road surface in near realtime, improvements in vehicle performance (e.g. braking and cornering) may be possible. Classification performance for many combinations of feature encoding and neural network types are compared. The vehicle used here was an An Audi “S3” with a magnetic suspension system on the sport mode. An NI CompactRIO (or cDAQ) module was used to record from a lowing the cDAQ to communicate with the PCB 352C03 one-axis accelerometer. The accelerometer was firmly attached to the windshield of the car. This work focuses on the classification of four road surfaces (asphalt, dirt, concrete, and sand), though larger target sets were also considered. The most accurate method involved a MATLAB feature extraction package with a back-propagation neural network, yielding an overall accuracy of 97%. Lessons learned from this wide exploration of options may extend to other related classification problems.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

MODELING AND OPTIMIZATION OF A PCM-BASED OCEAN THERMAL ENERGY HARVESTER FOR POWERING UNMANNED UNDERWATER VEHICLES

As oceans cover over 70% of the planet’s surface, they represent a large reservoir of resources that remained vastly untapped. Uncrewed Underwater Vehicles (UUVs) are becoming key technology for ocean exploration. Ocean thermal gradient is a permanent and reliable energy source that can be used to power UUVs using phase change material (PCM)-based thermal engines. However, the main drawback of this technology is its low conversion efficiency. When using PCM-based thermal engines to power UUVs, there are different energy conversion stages, thermal-hydraulic, hydraulic-kinetic, kinetic-electrical, dependent of a wide variety of parameters. Thus, optimization of the overall energy conversion is still a challenge for powering the increasing energy demanding UUVs for long missions. The goal of this study is to propose a PCM-based ocean thermal energy harvesting system for powering the Solo II float. Using the ocean thermal energy reduces the cost for battery replacement and expands the float’s mission time. For this purpose, we developed a multiphysics model for a system consisting of thermal material, hydraulic, and electrical systems. This system is designed to provide the electrical power needed for the instruments aboard the underwater vehicle. The models of the hydraulic and the electric systems are implemented in MATLAB-Simulink/Simscape environment. The final model developed can provide 13.66 kJ of electrical energy, which is more than 1.5 times the energy requirement per cycle for the SOLO II float.

Ouro-Koura, Habilou↗
DOE OSTI2022

A crystal plasticity finite element method modeling of Zircaloy with hydride phases based on scanning electron micrographs

This work developed a microstructure-based finite element model to predict the stress state of alloys with second phase inclusions. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. Generation of digital microstructure was achieved through two steps of tessellations using software Neper and Matlab. The process is demonstrated using an example of Zircaloy material with secondary phases of hydrides dispersed within and stress-strain response of Zircaloy containing hydrides was predicted. The constitutive material model for Zircaloy in this study was based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as isotropic elasto-plastic material. A parametric study had been conducted to understand the effect of volume fraction, and lamellae thickness of the hydride phase on the mechanical properties of the overall material. Results can help designers to alter the manufacturing process to obtain the enhanced mechanical properties for components used in nuclear applications made by Zircaloy material.

Kulkarni, Shank S.↗
DOE OSTI2023

New Developments and Capabilities Within WEC-Sim

WEC-Sim is an open-source software for simulating wave energy converters and has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory and Sandia National Laboratories within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with the latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1, which spans a 2-year period from June 2020 to September 2022. New features described here include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗
DOE OSTI2022

Large-scale parametric modeling of spent nuclear fuel dynamics in the 30 cm package drop scenario

Packages used to transport spent nuclear fuel (SNF) are required by the U.S. Code of Federal Regulations 10 CFR 71.71 to demonstrate satisfactory performance during a drop scenario. While the CFR is meant to ensure safe package function, it does not evaluate survival of the SNF within. The U.S. Department of Energy Spent Fuel and Waste Science and Technology program is working on closing the knowledge gap related to the response of SNF to external mechanical loads, including the hypothetical 30 cm package drop scenario in the CFR. In support of this effort, LS-DYNA finite element simulations were developed by Pacific Northwest National Laboratory (PNNL) to model generic drop scenarios at both the package and fuel assembly level. The models were validated against one-third scale package and full scale fuel assembly drop test data and were exercised to predict fuel cladding strains in a narrow range of model configurations. This work describes a large-scale parametric study conducted by PNNL using the previously developed and validated PWR finite element model, with the addition of a new generic BWR assembly model. The motivation for the parametric study was to characterize the broad range of SNF responses in the 30 cm package drop scenario. This was accomplished by varying the drop orientation, fuel assembly type (17x17 PWR and 10x10 BWR), burnup, cladding temperature, spacer grid buckling load, package mass, impact limiter stiffness, and mechanical gap conditions within the basket. A MATLAB framework was developed to automate LS-DYNA model generation and execution on PNNL institutional computing resources. In total, over 2000 simulations were performed. For each simulation, the SNF response was quantified in terms of permanent grid deformation, fuel rod contact pressure, and strains within the fuel rods, guide tubes, and water rods. The results provide valuable insight into the range of responses that could be reasonably expected from SNF in the hypothetical drop scenario, as well as the sensitivity to each input parameter. The results of this parametric study are a key component of the testing and modeling strategy the Spent Fuel and Waste Science and Technology program is using to close the external loads knowledge gap.

Kadooka, Kevin↗
DOE OSTI2010

Fabrication of ZnO photonic crystals by nanosphere lithography using inductively coupled-plasma reactive ion etching with CH{sub 4}/H{sub 2}/Ar plasma on the ZnO/GaN heterojunction light emitting diodes

This article reports fabrication of n-ZnO photonic crystal/p-GaN light emitting diode (LED) by nanosphere lithography to further booster the light efficiency. In this article, the fabrication of ZnO photonic crystals is carried out by nanosphere lithography using inductively coupled plasma reactive ion etching with CH{sub 4}/H{sub 2}/Ar plasma on the n-ZnO/p-GaN heterojunction LEDs. The CH{sub 4}/H{sub 2}/Ar mixed gas gives high etching rate of n-ZnO film, which yields a better surface morphology and results less plasma-induced damages of the n-ZnO film. Optimal ZnO lattice parameters of 200 nm and air fill factor from 0.35 to 0.65 were obtained from fitting the spectrum of n-ZnO/p-GaN LED using a MATLAB code. In this article, we will show our recent result that a ZnO photonic crystal cylinder has been fabricated using polystyrene nanosphere mask with lattice parameter of 200 nm and radius of hole around 70 nm. Surface morphology of ZnO photonic crystal was examined by scanning electron microscope.

36 MATERIALS SCIENCE↗