Search NASA⌕ Search

SEARCH · Search NASA

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A high-precision instrument for analyzing nonlinear dynamic behavior of bearing cage

The high-precision ball bearing is fundamental to the performance of complex mechanical systems. As the speed increases, the cage behavior becomes a key factor in influencing the bearing performance, especially life and reliability. This paper develops a high-precision instrument for analyzing nonlinear dynamic behavior of the bearing cage. The trajectory of the rotational center and non-repetitive run-out (NRRO) of the cage are used to evaluate the instability of cage motion. This instrument applied an aerostatic spindle to support and spin test the bearing to decrease the influence of system error. Then, a high-speed camera is used to capture images when the bearing works at high speeds. A 3D trajectory tracking software TEMA Motion is used to track the spot which marked the cage surface. Finally, by developing the MATLAB program, a Lissajous’ figure was used to evaluate the nonlinear dynamic behavior of the cage with different speeds. The trajectory of rotational center and NRRO of the cage with various speeds are analyzed. The results can be used to predict the initial failure and optimize cage structural parameters. In addition, the repeatability precision of instrument is also validated. In the future, the motorized spindle will be applied to increase testing speed and image processing algorithms will be developed to analyze the trajectory of the cage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

MHD effects and heat transfer for the UCM fluid along with Joule heating and thermal radiation using Cattaneo-Christov heat flux model

Present study examines the numerical analysis of MHD flow of Maxwell fluid with thermal radiation and Joule heating by considering the recently developed Cattaneo-Christov heat flux model which explains the time relaxation characteristics for the heat flux. The objective is to analyze the governing parameters such as viscoelastic fluid parameter, Magnetic parameter, Eckert and Prandtl number’s impact on the velocity and temperature profiles through graphs and tables. Suitable similarity transformations have been used to reduce the formulated PDEs into a system of coupled non-linear ODEs. Shooting technique has been invoked for finding the numerical solutions of the dimensionless velocity and temperature profiles. Additionally, the MATLAB built-in routine bvp4c has also been used to verify and strengthen the results obtained by shooting method. From some special cases of the present work, a comparison with the previously published results has been presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Quantitative analysis of total reflection X-ray fluorescence from finely layered structures using XeRay

Total reflection x-ray fluorescence (TXRF) is a widely applicable experimental technique for studying chemical element distributions across finely layered structures at extremely high sensitivity. To promote and facilitate scientific discovery using TXRF, we developed a MATLAB-based software package with a graphical user interface, named XeRay, for quick, accurate, and intuitive data analysis. XeRay lets the user model any layered system, each layer with its independent chemical composition and thickness, and enables fine-tuned data fitting. The accuracy of XeRay has been tested in the analysis of TXRF data from both air/liquid interface and liquid/liquid interfacial studies and has been compared to literature results. In an air/liquid interface study, Ca 2+ sequestration was measured at a Langmuir monolayer of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidic acid (SOPA) on a buffer solution of 1 mM CaCl2 at pH 7.5. Data analysis with XeRay reveals that each 1 nm2 of interfacial area contains 2.38 ± 0.06 Ca 2+ ions, which corresponds to a 1:1 ratio between SOPA headgroups and Ca 2+ ions, consistent with several earlier reports. For the liquid/liquid interface study of Sr 2+ enrichment at the dodecane/surfactant/water interface, analysis using XeRay gives a surface enrichment of Sr 2+ at $68^{+6}_{-5}$ Å 2 per ion, consistent with the result published for the same dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Accurate, rapid identification of dislocation lines in coherent diffractive imaging via a min-max optimization formulation

Defects such as dislocations impact materials properties and their response during external stimuli. Imaging these defects in their native operating conditions to establish the structure-function relationship and, ultimately, to improve performance via defect engineering has remained a considerable challenge for both electron-based and x-ray-based imaging techniques. While Bragg coherent x-ray diffractive imaging (BCDI) is successful in many cases, nuances in identifying the dislocations has left manual identification as the preferred method. Derivative-based methods are also used, but they can be inaccurate and are computationally inefficient. Here we demonstrate a derivative-free method that is both more accurate and more computationally efficient than either derivative-or human-based methods for identifying 3D dislocation lines in nanocrystal images produced by BCDI. We formulate the problem as a min-max optimization problem and show exceptional accuracy for experimental images. We demonstrate a 227x speedup for a typical experimental dataset with higher accuracy over current methods. We discuss the possibility of using this algorithm as part of a sparsity-based phase retrieval process. Here, we also provide MATLAB code for use by other researchers.

47 OTHER INSTRUMENTATION↗
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 OSTI2014

Bayesian Proteoform Modeling Improves Protein Quantification of Global Proteomic Measurements

As the capability of mass spectrometry-based proteomics has matured, tens of thousands of peptides can be measured simultaneously, which has the benefit of offering a systems view of protein expression. However, a major challenge is that with an increase in throughput, protein quantification estimation from the native measured peptides has become a computational task. A limitation to existing computationally-driven protein quantification methods is that most ignore protein variation, such as alternate splicing of the RNA transcript and post-translational modifications or other possible proteoforms, which will affect a significant fraction of the proteome. The consequence of this assumption is that statistical inference at the protein level, and consequently downstream analyses, such as network and pathway modeling, have only limited power for biomarker discovery. Here, we describe a Bayesian model (BP-Quant) that uses statistically derived peptides signatures to identify peptides that are outside the dominant pattern, or the existence of multiple over-expressed patterns to improve relative protein abundance estimates. It is a research-driven approach that utilizes the objectives of the experiment, defined in the context of a standard statistical hypothesis, to identify a set of peptides exhibiting similar statistical behavior relating to a protein. This approach infers that changes in relative protein abundance can be used as a surrogate for changes in function, without necessarily taking into account the effect of differential post-translational modifications, processing, or splicing in altering protein function. We verify the approach using a dilution study from mouse plasma samples and demonstrate that BP-Quant achieves similar accuracy as the current state-of-the-art methods at proteoform identification with significantly better specificity. BP-Quant is available as a MatLab ® and R packages at https://github.com/PNNL-Comp-Mass-Spec/BP-Quant.

Protein Quantitation, proteomics↗
DOE OSTI2017

A New Modeling Technique to Analyze Safeguards Measurements in Large Systems

SafeGuards Analysis (SGA) is a computational toolbox able to simulate different safeguards scenarios across a number of different fuel cycles and at many different scales within the Matlab Simulink framework. SGA functions by simulating Material Balance Areas (MBAs) under safeguards materials control and accountability and allows the user to define the uncertainty parameters of the associated flow and inventory measurements. The simulated safeguard system uses the uncertain measurement estimates to calculate a mass-balance across the MBA. This mass balance is then evaluated by one of a number of different statistical tests to determine if a significant amount of material has been removed from the MBA. This work describes the design of SGA, the results of testing each element of the toolbox, and a number of single MBA example scenarios. In all of the test cases SGA performed as expected, and produced acceptable results from the single MBA scenario.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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 OSTI2017

A Fault Location Algorithm for Two-End Series-Compensated Double-Circuit Transmission Lines Using the Distributed Parameter Line Model

A new fault location algorithm for two-end series-compensated double-circuit transmission lines utilizing unsynchronized two-terminal current phasors and local voltage phasors is presented in this paper. The distributed parameter line model is adopted to take into account the shunt capacitance of the lines. The mutual coupling between the parallel lines in the zero-sequence network is also considered. The boundary conditions under different fault types are used to derive the fault location formulation. The developed algorithm directly uses the local voltage phasors on the line side of series compensation (SC) and metal oxide varistor (MOV). However, when potential transformers are not installed on the line side of SC and MOVs for the local terminal, these measurements can be calculated from the local terminal bus voltage and currents by estimating the voltages across the SC and MOVs. MATLAB SimPowerSystems is used to generate cases under diverse fault conditions to evaluating accuracy. The simulation results show that the proposed algorithm is qualified for practical implementation.

current phasors, distributed parameter line model,↗
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 OSTI2012

Constructing realistic Szekeres models from initial and final data

The Szekeres family of inhomogeneous solutions, which are defined by six arbitrary metric functions, offers a wide range of possibilities for modelling cosmic structure. Here we present a model construction procedure for the quasispherical case using given data at initial and final times. Of the six arbitrary metric functions, the three which are common to both Szekeres and Lemaître-Tolman models are determined by the model construction procedure of Krasinski and Hellaby. For the remaining three functions, which are unique to Szekeres models, we derive exact analytic expressions in terms of more physically intuitive quantities — density profiles and dipole orientation angles. Using MATLAB, we implement the model construction procedure and simulate the time evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

catsHTM: A Tool for Fast Accessing and Cross-matching Large Astronomical Catalogs

Fast access to large catalogs is required for some astronomical applications. Here we introduce the catsHTM tool, consisting of several large catalogs reformatted into HDF5-based file format, which can be downloaded and used locally. To allow fast access, the catalogs are partitioned into hierarchical triangular meshes and stored in HDF5 files. Several tools are provided to perform efficient cone searches at resolutions spanning from a few arc-seconds to degrees, within a few milliseconds time. The first released version includes the following catalogs (by alphabetical order): 2MASS, 2MASS extended sources, AKARI, APASS, Cosmos, DECaLS/DR5, FIRST, GAIA/DR1, GAIA/DR2, GALEX/DR6Plus7, HSC/v2, IPHAS/DR2, NED redshifts, NVSS, Pan-STARRS1/DR1, PTF photometric catalog, ROSAT faint source, SDSS sources, SDSS/DR14 spectroscopy, SkyMapper, Spitzer/SAGE, Spitzer/IRAC galactic center, UCAC4, UKIDSS/DR10, VST/ATLAS/DR3, VST/KiDS/DR3, WISE and XMM. We provide Python code that allows to perform cone searches, as well as MATLAB code for performing cone searches, catalog cross-matching, general searches, as well as load and create these catalogs.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2017

Preliminary Mechanical Design Study of the Hollow Electron Lens for HL-LHC

A Hollow Electron Lens (HEL) has been proposed in order to improve performance of halo control and collimation in the Large Hadron Collider in view of its High Luminosity upgrade (HL-LHC). The concept is based on a hollow beam of electrons that travels around the protons for a few meters. The electron beam is produced by a cathode and then guided by a strong magnetic field. The first step of the design is the definition of the magnetic field that drives the electron trajectories. The estimation of such trajectories by means of a dedicated MATLAB tool is presented. The influence of the main geometrical and electrical parameters is analyzed and discussed. Then, the main mechanical design choices for the solenoids, cryostats gun and collector are described. The aim of this paper is to provide an overview of the feasibility study of the Electron Lens for LHC. The methods used in this study also serve as examples for future mechanical and integration designs of similar devices.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Online helium inventory monitoring of JLab cryogenic systems

There are five cryogenic plants at Jefferson Lab which support the LINAC, experiment hall end-stations and test facility. The majority of JLab's helium inventory, which is around 15 tons, is allocated in the LINAC cryo-modules, with the majority of the balance of helium distributed at the cryogenic-plant level mainly as stored gas and liquid for stable operation. Due to the organic evolution of the five plants and independent actions within the experiment halls, the traditional inventory management strategy suffers from rapid identification of potential leaks. This can easily result in losses many times higher than the normally accepted (average) loss rate. A real-time program to quickly identify potential excessive leakage was developed and tested. This program was written in MATLAB© for portability, easy diagnostics and modification. It interfaces directly with EPICS to access the cryogenic system state, and with and NIST REFPROP© for real fluid properties. This program was validated against the actual helium offloaded into the system. The present paper outlines the details of the inventory monitoring program, its validation and a sample of the achieved results.

42 ENGINEERING↗
DOE OSTI2019

Neutron imaging analysis using jupyter Python notebook

Abstract Independently of the image modality (x-rays, neutrons, etc), image data analysis requires normalization, a preprocessing step. While the normalization can sometimes easily be generalized, the analysis is, in most cases, specific to an experiment and a sample. Although many tools (MATLAB, ImageJ, VG Studio…) offer a large collection of pre-programmed image analysis tools, they usually require a learning step that can be lengthy depending on the skills of the end user. We have implemented Jupyter Python notebooks to allow easy and straightforward data analysis, along with live interaction with the data. Jupyter notebooks require little programming knowledge and the steep learning curve is bypassed. Most importantly, each notebook can be tailored to a specific experiment and sample with minimized effort. Here, we present the pros and cons of the main methods to analyse data and show the reason why we have found that Jupyter Python notebooks are well suited for imaging data processing, visualization and analysis.

Bilheux, Jean-Christophe (ORCID:0000000321726487)↗
DOE OSTI2014

Causal network inference using biochemical kinetics

Motivation: Networks are widely used as structural summaries of biochemical systems. Statistical estimation of networks is usually based on linear or discrete models. However, the dynamics of biochemical systems are generally non-linear, suggesting that suitable non-linear formulations may offer gains with respect to causal network inference and aid in associated prediction problems. Results: We present a general framework for network inference and dynamical prediction using time course data that is rooted in nonlinear biochemical kinetics. This is achieved by considering a dynamical system based on a chemical reaction graph with associated kinetic parameters. Both the graph and kinetic parameters are treated as unknown; inference is carried out within a Bayesian framework. This allows prediction of dynamical behavior even when the underlying reaction graph itself is unknown or uncertain. Results, based on (i) data simulated from a mechanistic model of mitogen-activated protein kinase signaling and (ii) phosphoproteomic data from cancer cell lines, demonstrate that non-linear formulations can yield gains in causal network inference and permit dynamical prediction and uncertainty quantification in the challenging setting where the reaction graph is unknown. Availability and implementation: MATLAB R2014a software is available to download from warwick.ac.uk/chrisoates.

60 APPLIED LIFE SCIENCES↗