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

$\mathrm{REX}$: An analytical tool for reactor operating envelope expansion through fuel-clad thermo-mechanics

Here REX, an analytical reactor design tool for the identification and expansion of reactor operating envelopes, is introduced and detailed. Written in MATLAB, REX takes core geometries and axially discretized detector and depletion output from the Serpent neutron transport code and performs full-core, pin-specific thermal-hydraulic and thermo-mechanical calculations. REX identifies the operating envelopes of advanced, solid-fuel nuclear reactors on the basis of their fuel-cladding, fuel-gas, and coolant-cladding interfacial thermo-mechanics. It iteratively pushes coolant inlet temperatures and modifies fuel pin geometries to induce mechanical failure in their original cladding materials at the fullest extent of their fuel cycle lengths and permissible coolant inlet temperatures. It then attempts to expand their operating envelopes by determining the mechanical responses of alternative cladding materials under the same geometric conditions. The result of the REX calculation sequence is a set of 5-D variables describing the temperatures, pressures, geometries, and mechanics of the core as functions of assembly, fuel pin, axial zone, depletion step, and coolant inlet temperature for each candidate cladding material with the limiting fuel pin geometry for the original material.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Miniature DC electromagnetic pumps of molten lead and sodium to support development of Gen-IV nuclear reactors

This work developed designs of submerged, dual-regions, miniature DC-EM pumps for circulating molten lead and liquid sodium at ≤500 °C without active cooling for use in ex-pile and in-pile test loops to support of materials and fuel developments for Gen-IV sodium and molten lead fast nuclear reactor. These pumps with two pumping regions are 57 mm, 66.8 mm, 95.4 mm, and 133.5 mm in diameter and have two Alnico 5 permanent magnets with Hiperco-50 pole pieces for focusing the magnetic field lines in the flow duct. The Equivalent Circuit Model (ECM) linked to the Finite Element Method Magnetics (FEMM) software in MATLAB platform calculated the pumps characteristics for a wide range of parameters. Also determined are the pumps dimensions for the highest cumulative pumping power and peak efficiency. These dimensions are the height, width, and length of the flow duct, the thickness of the Alnico magnets, the length of current electrodes, and the separation distance between the two pumping regions. For molten lead, the pumping power increases from 368 W to 728 W, and the peak efficiency from 14.7% to 31.6% with increased pump diameter from 57 mm to 133.5 mm. In conclusion, for liquid sodium, both the pumping power and peak efficiency are higher, increasing from 392 W to 767 W and from 44.3 to 51.2%, respectively, with increased pump diameter from 57 mm to 133.5 mm.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Second-order wave excitation forces in WEC-Sim/MOST: Implementation, experimental validation, and code-to-code comparison

Accurate prediction of second-order hydrodynamic loads is essential for floating bodies, including floating offshore wind turbines, wave energy converters, and hybrid wind–wave platforms. These nonlinear effects, arising from both sum- and difference-frequency forcing, are critical for capturing key response characteristics but remain challenging to model efficiently. In this work, we extend the open-source Wave Energy Converter Simulator / MATLAB for Offshore Simulation Tool by implementing second-order wave excitation forces, supporting both the full Quadratic Transfer Function formulation and the Newman approximation. The full Quadratic Transfer Function method is used for all code-to-code comparisons and experimental validation, while the Newman approximation is provided as a computationally lighter alternative. To benchmark the new capability, we perform a code-to-code comparison with OpenFAST and OrcaFlex. We then validate the enhanced model using wave-tank measurements of a 1:96 scale DeepCwind semi-submersible, showing that second-order effects are required to reproduce platform motions. The implementation employs a computationally efficient pre-computation strategy for second-order wave excitation forces, reducing simulation cost while maintaining engineering accuracy. Overall, this work advances the tool as an open-source and versatile tool for modelling floating offshore renewable-energy systems requiring second-order hydrodynamic fidelity.

17 WIND ENERGY↗
DOE OSTI2020

Theoretical power-law relationship between permeability and formation factor

Permeability k and formation factor F are important petrophysical properties that have broad applications to reservoir engineering. Although well-log analysis has advanced reservoir characterizations, direct measurement of k is still not available. Accordingly, resistivity logs and more particularly formation factor have been frequently used to estimate k. In this study, we present a theoretical framework to relate k to F. We use power-law scaling from continuum percolation theory to develop a theoretic relationship between k and F, which is similar in form to a published power-law equation with empirical coefficients. Our result gives a physical interpretation to this empirical relationship. We also compare power-law scaling laws from percolation theory with data derived from experimental measurements and numerical simulations in the literature. Finally, using the Curve Fitting toolbox of MATLAB, the parameters of the power-law equations are determined, and the significance of the non-universal power laws are discussed.

42 ENGINEERING↗
DOE OSTI2025

On the optimal sizing of power take-off systems for wave energy converters

The power take-off (PTO) system in a wave energy converter (WEC) is the means by which the energy in an ocean wave is converted into useful energy. There has been a general lack of technology convergence in the wave energy industry and PTO design has been largely device-specific with the design process determined by the WEC developer. An optimized PTO design can significantly affect the efficiency, reliability, performance, and overall cost and viability of a device. There remains a need in the marine energy industry for best practices and efficient design processes to be identified as it could greatly benefit developers and researchers in designing an optimal WEC PTO. This research utilizes the open-source Matlab-Simulink based software, WEC-Sim, which was developed by the National Renewable Energy Laboratory, and the Eagle, a high-performance computing system, to model and optimize two WECs of different archetypes and then investigates if any trends exist that can be exploited for greater efficiency in optimal PTO design. The results from this research indicate that there may be an optimal sizing for WEC PTOs for both power and PTO force rating and that substantial reductions can be made in the cost of a device without significant loss in the energy produced. Furthermore this optimal sizing may be independent of WEC type and deployment location. With an optimal, or near optimal, approach to WEC PTO sizing demonstrated, a methodology is proposed to address the challenge of nameplate ratings in the wave energy industry.

13 HYDRO ENERGY↗
DOE OSTI2024

2-D soil zymography: Accounting for the spatial variation of pH

Soil zymography is commonly used to quantify spatial distribution of hydrolytic enzyme activities on soil and plant root surfaces. It is recommended to adjust pH in zymography substrates and calibration solutions with respect to soil/root pH. However, pH values may vary greatly within a few mm of plant rhizosphere, potentially altering the distribution of pH in zymography membranes. Despite the fact that the effect of pH on the calibration of zymography membranes is generally known, its potential impact on zymography results is unaccounted for in processing zymography images and calculations of enzyme activity. In this study we assessed the effect of pH variations on the persistency of the methylumbelliferone (MUF) calibration. The studied pH values ranged from 4.5 to 7.5. The MUF calibration curves greatly deviated from that at a reference pH of 6.5, with a marked nonlinear increase of deviation with greater membrane brightness. We suggest that the problem can be partially alleviated by reducing the membrane incubation time. However, such deviations suggest the need for a more comprehensive resolution via mapping pH and using pH-specific calibrations to process zymography images. Here, we developed a MATLAB code to implement a pixel-based correction of enzyme activity for pH in processing time-lapse zymography images.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enhanced ZIF-8-enabled colorimetric CO 2 sensing through dye-precursor synthesis

The accumulation of carbon dioxide (CO 2 ) within enclosed spaces, along with volatile organic compounds, under certain humidity, temperature, and ventilation conditions is associated with detrimental human health symptoms such as fatigue. Color-based chemical sensing is a promising approach to detect CO 2 levels relevant to indoor air quality through producing fast, quantifiable output visible to the naked eye. In a prior work, a colorimetric gas sensor was fabricated through synthesizing the metal-organic framework, ZIF-8, as the adsorbent, followed by post-synthetic mixing with a dye, phenol red (PSP), and primary amine, ethylenediamine (ED). While this sensor (termed PSP-ED/ZIF-8) maintained its structural integrity in atmospheric conditions and exhibited an increasing fuchsia-to-yellow color change with increasing CO 2 levels in dry environment, the colorimetric response greatly suffered in the presence of humid CO 2 . Here in this work, a significantly improved colorimetric CO 2 sensor (referred to as ED/PSP:ZIF-8) is accomplished through directly incorporating phenol red in the ZIF-8 metal and linker precursor solutions and then blending with ethylenediamine. MATLAB-generated color distributions and in-situ ultraviolet-visible (UV-Vis) spectroscopic studies quantitatively demonstrate an enhanced colorimetric gas response of ED/PSP:ZIF-8 compared to that of PSP-ED/ZIF-8 across an important range of CO 2 for indoor air quality monitoring (500 – 3500 ppm) and across a range of humidity. The new sensor also exhibits high selectivity to CO 2 compared to select volatile organic compounds, such as acetone and ethanol, which contribute to human health symptoms experienced indoors. The enhanced performance is attributed to the proposed incorporation of phenol red within ZIF-8, while maintaining the chemical stability of the MOF.

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

Recycling of silver from silicon solar cells by laser debonding

A laser-debonding approach to silver electrode recovery from solar cells is presented to address the critical need for efficient and eco-friendly recycling methods. The study explores the use of UV nanosecond and IR continuous-wave lasers to precisely debond silver electrodes from silicon wafers. By optimizing parameters like laser power, scan speed, and the number of passes, intact silver electrical contact lines were recovered from solar cells. Scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDS) analysis confirmed the successful separation of silver from silicon wafers. Furthermore, an increase in the number of laser passes led to the production of silver microparticles, as validated by morphology and compositional studies using EDS. The estimated silver is approximately 90 mg for the Si solar cells of 15.4 * 15.4 cm 2 used in the study. To achieve a large area recovery of silver, a MATLAB code incorporating detection algorithms and specific criteria to automate laser scanning by accurately identifying the location of silver electrodes was developed. This presented laser debonding method is not confined to silicon solar cells but can be extended to other solar cell types featuring metal electrodes, offering a versatile solution. By minimizing the use of hazardous chemicals and reducing operational costs, the developed process aligns with the imperative of environmentally responsible photovoltaic waste management.

14 SOLAR ENERGY↗
DOE OSTI2021

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

GriddingMachine, a database and software for Earth system modeling at global and regional scales

Land and Earth system modeling is moving towards more explicit biophysical representations, requiring increasing variety of datasets for initialization and benchmarking. However, researchers often have difficulties in identifying and integrating non-standardized datasets from various sources. We aim towards a standardized database and one-stop distribution method of global datasets. Here, we present the GriddingMachine as (1) a database of global-scale datasets commonly used to parameterize or benchmark the models, from plant traits to vegetation indices and geophysical information and (2) a cross-platform open source software to download and request a subset of datasets with only a few lines of code. The GriddingMachine datasets can be accessed either manually through traditional HTTP, or automatically using modern programming languages including Julia, Matlab, Octave, Python, and R. The GriddingMachine collections can be used for any land and Earth modeling framework and ecological research at the regional and global scales, and the number of datasets will continue to grow to meet the increasing needs of research communities.

58 GEOSCIENCES↗
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↗