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 469 records · Page 26

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

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

Interpretable machine learning↗

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

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

36 MATERIALS SCIENCE↗

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

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

SOEC↗

Application of the Hartmann–Tran profile to precise experimental data sets of 12 C 2 H 2

Self- and nitrogen-broadened line shape data for the P e (11) line of the ν₁ + ν₃ band of acetylene, recorded using a frequency comb-stabilized laser spectrometer, have been analyzed using the Hartmann–Tran profile (HTP) line shape model in a multispectrum fitting. In total, the data included measurements recorded at temperatures between 125 K and 296 K and at pressures between 4 and 760 Torr. New, sub-Doppler, frequency comb-referenced measurements of the positions of multiple underlying hot band lines have also been made. These underlying lines significantly affect the P e (11) line profile at temperatures above 240 K and poorly known frequencies previously introduced errors into the line shape analyses. Thus, the behavior of the HTP model was compared to the quadratic speed dependent Voigt profile (QSDVP) expressed in the frequency and time domains. A parameter uncertainty analysis was carried out using a Monte Carlo method based on the estimated pressure, transmittance and frequency measurement errors. From the analyses, the P e (11) line strength was estimated to be 1.2014(50) × 10 -20 in cm.molecules⁻¹ units at 296 K with the standard deviation in parenthesis. For analyzing these data, we found that a reduced form of the HTP, equivalent to the QSDVP, was most appropriate because the additional parameters included in the full HTP were not well determined. As a supplement to this work, expressions for analytic derivatives and a lineshape fitting code written in Matlab for the HTP are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An alternative to the crystallographic reconstruction of austenite in steels

An alternative crystallographic austenite reconstruction programme written in Matlab is developed by combining the best features of the existing models: the orientation relationship refinement, the local pixel-by-pixel analysis and the nuclei identification and spreading strategy. This programme can be directly applied to experimental electron backscatter diffraction mappings. Its applicability is demonstrated on both quenching and partitioning and as-quenched lath-martensite steels. - Highlights: • An alternative crystallographic austenite reconstruction program is developed. • The method combines a local analysis and a nuclei identification/spreading strategy. • The validity of the calculated orientation relationship is verified on a Q and P steel. • The accuracy of the reconstructed microtexture is investigated on a martensite steel.

36 MATERIALS SCIENCE↗

Tutorial: Crystal orientations and EBSD — Or which way is up?

Electron backscatter diffraction (EBSD) is an automated technique that can measure the orientation of crystals in a sample very rapidly. There are many sophisticated software packages that present measured data. Unfortunately, due to crystal symmetry and differences in the set-up of microscope and EBSD software, there may be accuracy issues when linking the crystal orientation to a particular microstructural feature. In this paper we outline a series of conventions used to describe crystal orientations and coordinate systems. These conventions have been used to successfully demonstrate that a consistent frame of reference is used in the sample, unit cell, pole figure and diffraction pattern frames of reference. We establish a coordinate system rooted in measurement of the diffraction pattern and subsequently link this to all other coordinate systems. A fundamental outcome of this analysis is to note that the beamshift coordinate system needs to be precisely defined for consistent 3D microstructure analysis. This is supported through a series of case studies examining particular features of the microscope settings and/or unambiguous crystallographic features. These case studies can be generated easily in most laboratories and represent an opportunity to demonstrate confidence in use of recorded orientation data. Finally, we include a simple software tool, written in both MATLAB® and Python, which the reader can use to compare consistency with their own microscope set-up and which may act as a springboard for further offline analysis. - Highlights: • Presentation of conventions used to describe crystal orientations • Three case studies that outline how conventions are consistent • Demonstrates a pathway for calibration and validation of EBSD based orientation measurements • EBSD computer code supplied for validation by the reader.

36 MATERIALS SCIENCE↗

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

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

36 MATERIALS SCIENCE↗

Tissue Heterogeneity in IMRT Dose Calculation for Lung Cancer

The aim of this study was to evaluate the differences in accuracy of dose calculation between 3 commonly used algorithms, the Pencil Beam algorithm (PB), the Anisotropic Analytical Algorithm (AAA), and the Collapsed Cone Convolution Superposition (CCCS) for intensity-modulated radiation therapy (IMRT). The 2D dose distributions obtained with the 3 algorithms were compared on each CT slice pixel by pixel, using the MATLAB code (The MathWorks, Natick, MA) and the agreement was assessed with the {gamma} function. The effect of the differences on dose-volume histograms (DVHs), tumor control, and normal tissue complication probability (TCP and NTCP) were also evaluated, and its significance was quantified by using a nonparametric test. In general PB generates regions of over-dosage both in the lung and in the tumor area. These differences are not always in DVH of the lung, although the Wilcoxon test indicated significant differences in 2 of 4 patients. Disagreement in the lung region was also found when the {Gamma} analysis was performed. The effect on TCP is less important than for NTCP because of the slope of the curve at the level of the dose of interest. The effect of dose calculation inaccuracy is patient-dependent and strongly related to beam geometry and to the localization of the tumor. When multiple intensity-modulated beams are used, the effect of the presence of the heterogeneity on dose distribution may not always be easily predictable.

61 RADIATION PROTECTION AND DOSIMETRY↗

Comp Plan: A computer program to generate dose and radiobiological metrics from dose-volume histogram files

Treatment planning studies often require the calculation of a large number of dose and radiobiological metrics. To streamline these calculations, a computer program called Comp Plan was developed using MATLAB. Comp Plan calculates common metrics, including equivalent uniform dose, tumor control probability, and normal tissue complication probability from dose-volume histogram data. The dose and radiobiological metrics can be calculated for the original data or for an adjusted fraction size using the linear quadratic model. A homogeneous boost dose can be added to a given structure if desired. The final output is written to an Excel file in a format convenient for further statistical analysis. Comp Plan was verified by independent calculations. A lung treatment planning study comparing 45 plans for 7 structures using up to 6 metrics for each structure was successfully analyzed within approximately 5 minutes with Comp Plan. The code is freely available from the authors on request.

61 RADIATION PROTECTION AND DOSIMETRY↗

Comparison of 3D dose distributions for HDR {sup 192}Ir brachytherapy sources with normoxic polymer gel dosimetry and treatment planning system

Radiation fluence changes caused by the dosimeter itself and poor spatial resolution may lead to lack of 3-dimensional (3D) information depending on the features of the dosimeter and quality assurance of dose distributions for high–dose rate (HDR) iridium-192 ({sup 192}Ir) brachytherapy sources is challenging and experimental dosimetry methods used for brachytherapy sources are limited. In this study, we investigated 3D dose distributions of {sup 192}Ir brachytherapy sources for irradiation with single and multiple dwell positions using a normoxic gel dosimeter and compared them with treatment planning system (TPS) calculations. For dose calibration purposes, 100-mL gel-containing vials were irradiated at predefined doses and then scanned in an magnetic resonance (MR) imaging unit. Gel phantoms prepared in 2 spherical glasses were irradiated with {sup 192}Ir for the calculated dwell positions, and MR scans of the phantoms were obtained. The images were analyzed with MATLAB software. Dose distributions and profiles derived with 1-mm resolution were compared with TPS calculations. Linearity was observed between the delivered dose and the reciprocal of the T2 relaxation time constant of the gel. The x-, y-, and z-axes were defined as the sagittal, coronal, and axial planes, respectively, the sagittal and axial planes were defined parallel to the long axis of the source while the coronal plane was defined horizontally to the long axis of the source. The differences between measured and calculated profile widths of 3-cm source length and point source for 70%, 50%, and 30% isodose lines were evaluated at 3 dose levels using 18 profiles of comparison. The calculations for 3-cm source length revealed a difference of > 3 mm in 1 coordinate at 50% profile width on the sagittal plane and 3 coordinates at 70% profile width and 2 coordinates at 50% and 30% profile widths on the axial plane. Calculations on the coronal plane for 3-cm source length showed > 3-mm difference in 1 coordinate at 50% and 70% and 2 coordinates at 30% profile widths. The point source measurements and calculations for 50% profile widths revealed a difference > 3 mm in 1 coordinate on the sagittal plane and 2 coordinates on the axial plane. The doses of 3 coordinates on the sagittal plane and 4 coordinates on the axial plane could not be evaluated in 30% profile width because of low doses. There was good agreement between the gel dosimetry and TPS results. Gel dosimetry provides dose distributions in all 3 planes at the same time, which enables us to define the dose distributions in any plane with high resolution. It can be used to obtain 3D dose distributions for HDR {sup 192}Ir brachytherapy sources and 3D dose verification of TPS.

61 RADIATION PROTECTION AND DOSIMETRY↗

Clinical analysis of the approximate, 3-dimensional, biological effective dose equation in multiphase treatment plans

A multiphase, approximate biological effective dose (BED{sub A}) equation was introduced because most treatment planning systems (TPS) are incapable of calculating the true BED (BED{sub T}). This work investigates the accuracy and precision of the multiphase BED{sub A} relative to the BED{sub T} in clinical cases. Ten patients with head and neck cancer and 10 patients with prostate cancer were studied using their treatment plans from Pinnacle{sup 3} 9.2 (Philips Medical, Fitchburg, WI). The organs at risk (OARs) that were studied are the normal brain, left and right optic nerves, optic chiasm, spinal cord, brainstem, bladder, and rectum. BED{sub A} and BED{sub T} distributions were calculated using MATLAB 2010b (MathWorks, Natick, MA) and analyzed on a voxel basis for percent error, percent error volume histograms (PEVHs), Pearson correlation coefficient, and Bland-Altman analysis. The maximum BED values that were calculated using the BED{sub A} and BED{sub T} methods were also analyzed. BED{sub A} was found to always underestimate BED{sub T}. The accuracy and precision of BED{sub A} distributions varied between the organs: for optic chiasm and brainstem, 50% of the patients had an overall BED{sub A} percent error of A percent error of A percent error of A distributions had maximum errors ranging from 2% to 11%, with the 11% error occurring for bladder. BED{sub A} produced much more accurate maximum BED values with adjacent organs such as normal brain, bladder, and rectum. This study has shown that BED{sub A} can calculate BED distributions with acceptable accuracy under certain circumstances. However, its consistency and accuracy strongly depend on the dose distributions of the different treatment phases. One should be cautious when using BED{sub A}.

61 RADIATION PROTECTION AND DOSIMETRY↗

Exudate-based diabetic macular edema detection in fundus images using publicly available datasets

Diabetic macular edema (DME) is a common vision threatening complication of diabetic retinopathy. In a large scale screening environment DME can be assessed by detecting exudates (a type of bright lesions) in fundus images. In this work, we introduce a new methodology for diagnosis of DME using a novel set of features based on colour, wavelet decomposition and automatic lesion segmentation. These features are employed to train a classifier able to automatically diagnose DME through the presence of exudation. We present a new publicly available dataset with ground-truth data containing 169 patients from various ethnic groups and levels of DME. This and other two publicly available datasets are employed to evaluate our algorithm. We are able to achieve diagnosis performance comparable to retina experts on the MESSIDOR (an independently labelled dataset with 1200 images) with cross-dataset testing (e.g., the classifier was trained on an independent dataset and tested on MESSIDOR). Our algorithm obtained an AUC between 0.88 and 0.94 depending on the dataset/features used. Additionally, it does not need ground truth at lesion level to reject false positives and is computationally efficient, as it generates a diagnosis on an average of 4.4 s (9.3 s, considering the optic nerve localization) per image on an 2.6 GHz platform with an unoptimized Matlab implementation.

59 BASIC BIOLOGICAL SCIENCES↗

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↗

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 ↗

$\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↗

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↗

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↗