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

Prototype HB650 Cryomodule Heat Loads Simulations

During the design stages of the PIP-II cryomodules, many analytical calculations and FEA have been performed on simpler geometry in order to estimate the heat loads and also to optimize the design. To better analyze the cryomodule cold tests, simulations have been performed with MATLAB to determine the temperature of the main components during cool down and to determine the heat loads of the cryomodule. These simulations have been applied to the High Beta 650 MHz prototype cryomodule design and compared to the cold tests performed on it.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Automatic Calibration of a Geomechanical Model from Sparse Data for Estimating Stress in Deep Geological Formations

Summary In this study, we demonstrate geomechanical modeling with fully automatic parameter calibration to estimate the full geomechanical stress fields of a prospective US carbon dioxide (CO2) storage site, based on sparse measurement data. The goal is to compute full stress tensor field estimates (principal stresses and orientations) that are maximally compatible with observations within the constraints of the model assumptions, thereby extending pointwise, incomplete partial stress measurement to a simulated full formation stress field, as well as a rough assessment of the associated error. We use the Perch site, located in Otsego County, Michigan, USA, as our case study. The input data consist of partial stress tensor information inferred from in-situ borehole tests, geophysical well logs, and processing of seismic data. A static earth model (SEM) of the site was developed, and geomechanical simulation functionality of the open-source MATLAB Reservoir Simulation Toolbox (MRST) was used to model the stress field. Adjoint-based nonlinear optimization was used to adjust boundary conditions and material properties to calibrate simulated results of observations. Results were interpreted through a Bayesian framework. The focus of this paper is to demonstrate how the fully automatic calibration procedure works and discuss the results obtained; it does not attempt a detailed analysis of the stress field in the context of the proposed CO2 storage initiatives. Our work is part of a larger effort to noninvasively determine in-situ stresses in deep formations considered for CO2 storage. Guided by previously published research on geomechanical model calibration, our work presents a novel calibration approach supporting a potentially large number of linear or nonlinear calibration parameters to produce results optimally agreeing with available measurements and thus extend partial pointwise estimates to full tensor fields compatible with the physics of the site.

Engineering↗
DOE OSTI2024

QSpace - An open-source tensor library for Abelian and non-Abelian symmetries

This is the documentation for the tensor library QSpace (v4.0), a toolbox to exploit ‘quan tum symmetry spaces’ in tensor network states in the quantum many-body context. QSpace permits arbitrary combinations of symmetries including the abelian symmetries $\mathbb{Z}_n$ and U(1), as well as all non-abelian symmetries based on the semisimple classical Lie algebras: A n , B n , C n , and D n , or respectively, the special unitary group SU(n), the odd orthogonal group SO(2n+1), the symplectic group Sp(2n), and the even orthogonal group SO(2n). The code (C++ embedded via the MEX interface into Matlab) is available open source as of QSpace v4.0 on bitbucket under the Apache 2.0 license. QSpace is designed as a bottom-up approach for non-abelian symmetries. It starts from the defining representation and the respective Lie algebra. By explicitly comput ing and tabulating generalized Clebsch-Gordan coefficient tensors, QSpace is versatile in the type of operations that it can perform across all symmetries. At the level of an ap plication, much of the symmetry-related details are hidden within the QSpace C++ core libraries. Hence when developing tensor network algorithms with QSpace, these can be coded (nearly) as if there are no symmetries at all, despite being able to fully exploit general non-abelian symmetries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2010

TORCH Computational Reference Kernels - A Testbed for Computer Science Research

For decades, computer scientists have sought guidance on how to evolve architectures, languages, and programming models in order to improve application performance, efficiency, and productivity. Unfortunately, without overarching advice about future directions in these areas, individual guidance is inferred from the existing software/hardware ecosystem, and each discipline often conducts their research independently assuming all other technologies remain fixed. In today's rapidly evolving world of on-chip parallelism, isolated and iterative improvements to performance may miss superior solutions in the same way gradient descent optimization techniques may get stuck in local minima. To combat this, we present TORCH: A Testbed for Optimization ResearCH. These computational reference kernels define the core problems of interest in scientific computing without mandating a specific language, algorithm, programming model, or implementation. To compliment the kernel (problem) definitions, we provide a set of algorithmically-expressed verification tests that can be used to verify a hardware/software co-designed solution produces an acceptable answer. Finally, to provide some illumination as to how researchers have implemented solutions to these problems in the past, we provide a set of reference implementations in C and MATLAB.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2010

Project Report on Development of a Safeguards Approach for Pyroprocessing

The Idaho National Laboratory has undertaken an effort to develop a standard safeguards approach for international commercial pyroprocessing facilities. This report details progress for the fiscal year 2010 effort. A component by component diversion pathway analysis has been performed, and has led to insight on the mitigation needs and equipment development needed for a valid safeguards approach. The effort to develop an in-hot cell detection capability led to the digital cloud chamber, and more importantly, the significant potential scientific breakthrough of the inverse spectroscopy algorithm, including the ability to identify energy and spatial location of gamma ray emitting sources with a single, non-complex, stationary radiation detector system. Curium measurements were performed on historical and current samples at the FCF to attempt to determine the utility of using gross neutron counting for accountancy measurements. A solid cost estimate of equipment installation at FCF has been developed to guide proposals and cost allocations to use FCF as a test bed for safeguards measurement demonstrations. A combined MATLAB and MCNPX model has been developed to perform detector placement calculations around the electrorefiner. Early harvesting has occurred wherein the project team has been requested to provide pyroprocessing technology and safeguards short courses.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2010

The generation of shared cryptographic keys through channel impulse response estimation at 60 GHz.

Methods to generate private keys based on wireless channel characteristics have been proposed as an alternative to standard key-management schemes. In this work, we discuss past work in the field and offer a generalized scheme for the generation of private keys using uncorrelated channels in multiple domains. Proposed cognitive enhancements measure channel characteristics, to dynamically change transmission and reception parameters as well as estimate private key randomness and expiration times. Finally, results are presented on the implementation of a system for the generation of private keys for cryptographic communications using channel impulse-response estimation at 60 GHz. The testbed is composed of commercial millimeter-wave VubIQ transceivers, laboratory equipment, and software implemented in MATLAB. Novel cognitive enhancements are demonstrated, using channel estimation to dynamically change system parameters and estimate cryptographic key strength. We show for a complex channel that secret key generation can be accomplished on the order of 100 kb/s.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2011

Velocimetry signal synthesis with fringen.

An important part of velocimetry analysis is the recovery of a known velocity history from simulated data signals. The fringen program synthesizes VISAR and PDV signals, given a specified velocity history, using exact formulations for the optical signal. Time-dependent light conditions, non-ideal measurement conditions, and various diagnostic limitations (noise, etc.) may be incorporated into the simulated signals. This report describes the fringen program, which performs forward VISAR (Velocity Interferometer System for Any Reflector) and PDV (Photonic Doppler Velocimetry, also known as heterodyne velocimetry) analysis. Nearly all effects that might occur in VISAR/PDV measurement of a single velocity can be modeled by fringen. The program operates in MATLAB, either within a graphical interface or as a user-callable function. The current stable version of fringen is 0.3, which was released in October 2010. The following sections describe the operation and use of fringen. Section 2 gives a brief overview of VISAR and PDV synthesis. Section 3 illustrates the graphical and console interface of fringen. Section 4 presents several example uses of the program. Section 5 summarizes program capabilities and discusses potential future work.

47 OTHER INSTRUMENTATION↗
DOE OSTI2011

Data Acquisition in a High Harmonic Generation Lab and at LCLS

In this paper, we examine data acquisition in a high harmonic generation (HHG) lab and preliminary data analysis with the Cyclohexadiene Collaboration at the Linac Coherent Lightsource (LCLS) at SLAC National Accelerator Laboratory. HHG experiments have a large number of parameters that need to be monitored constantly. In particular, the pressure of the target is critical to HHG yield. However, this pressure can fluctuate wildly and without a tool to monitor it, it is difficult to analyze the correlation between HHG yield and the pressure. I used the Arduino microcontroller board and created a complementary MATLAB graphical user interface (GUI), thereby enhancing the ease with which users can acquire time-stamped parameter data. Using the Arduino, it is much easier to match the pressure to the corresponding HHG yield. Collecting data by using the Arduino and the GUI is flexible, user-friendly, and cost-effective. In the future, we hope to be able to control and monitor parts of the lab with the Arduino alone. While more parameter information is needed in the HHG lab, we needed to reduce the amount of data during the cyclohexadiene collaboration. This was achieved by sorting the data into bins and filtering out unnecessary details. This method was highly effective in that it minimized the amount of data without losing any valuable information. This effective preliminary data analysis technique will continue to be used to decrease the size of the collected data.

43 PARTICLE ACCELERATORS↗
DOE OSTI2010

Model Predictive Control of Integrated Gasification Combined Cycle Power Plants

The primary project objectives were to understand how the process design of an integrated gasification combined cycle (IGCC) power plant affects the dynamic operability and controllability of the process. Steady-state and dynamic simulation models were developed to predict the process behavior during typical transients that occur in plant operation. Advanced control strategies were developed to improve the ability of the process to follow changes in the power load demand, and to improve performance during transitions between power levels. Another objective of the proposed work was to educate graduate and undergraduate students in the application of process systems and control to coal technology. Educational materials were developed for use in engineering courses to further broaden this exposure to many students. ASPENTECH software was used to perform steady-state and dynamic simulations of an IGCC power plant. Linear systems analysis techniques were used to assess the steady-state and dynamic operability of the power plant under various plant operating conditions. Model predictive control (MPC) strategies were developed to improve the dynamic operation of the power plants. MATLAB and SIMULINK software were used for systems analysis and control system design, and the SIMULINK functionality in ASPEN DYNAMICS was used to test the control strategies on the simulated process. Project funds were used to support a Ph.D. student to receive education and training in coal technology and the application of modeling and simulation techniques.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2011

Discretization error estimation and exact solution generation using the method of nearby problems.

The Method of Nearby Problems (MNP), a form of defect correction, is examined as a method for generating exact solutions to partial differential equations and as a discretization error estimator. For generating exact solutions, four-dimensional spline fitting procedures were developed and implemented into a MATLAB code for generating spline fits on structured domains with arbitrary levels of continuity between spline zones. For discretization error estimation, MNP/defect correction only requires a single additional numerical solution on the same grid (as compared to Richardson extrapolation which requires additional numerical solutions on systematically-refined grids). When used for error estimation, it was found that continuity between spline zones was not required. A number of cases were examined including 1D and 2D Burgers equation, the 2D compressible Euler equations, and the 2D incompressible Navier-Stokes equations. The discretization error estimation results compared favorably to Richardson extrapolation and had the advantage of only requiring a single grid to be generated.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2011

Investigations on Repository Near-Field Thermal Modeling - Repository Science/Thermal Load Management & Design Concepts (M41UF033302) (Rev.1)

The modeling tools generated in MathCAD, Excel, and MatLab were used to calculate the temperature evolution for combinations of waste form and geologic medium, assuming a particular emplacement layout for each concept (Sutton et al. 2011). Two types of spent nuclear fuel assemblies (SNFA) were considered, namely UOX from an open fuel cycle and MOX from a modified-­‐open fuel cycle. Four types of high level waste (HLW) canisters were considered, namely “co-­‐extraction” glass from a modified-­‐open fuel cycle and “new extraction” glass, electrochemical ceramic (EC-­‐C) and electrochemical metal (EC-­‐M) from a closed fuel cycle. While the “coextraction” process is similar in function to the industrial Co-­‐Extraction™ (COEX) process deployed by AREVA, the two processes assume different processing methods and steps and so the product and waste streams cannot be directly compared. Similar is true for the “new extraction” process documented in this report and the NUEX industrial process proposed by Energy Solutions, which also cannot be directly compared.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2012

C++ Tensor Toolbox User Manual (V.1.0)

The C++ Tensor Toolbox is a software package for computing tensor decompositions. It is based on the Matlab Tensor Toolbox, and is particularly optimized for sparse data sets. This user manual briefly overviews tensor decomposition mathematics, software capabilities, and installation of the package. Tensors (also known as multidimensional arrays or N-way arrays) are used in a variety of applications ranging from chemometrics to network analysis. The Tensor Toolbox provides classes for manipulating dense, sparse, and structured tensors in C++. The Toolbox compiles into libraries and is intended for use with custom applications written by users.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2012

Final Progress Report: Internship at Los Alamos National Laboratory

Originally I was tasked fluidized bed modeling, however, I changed projects. While still working with ANSYS Fluent, I performed a study of particle tracks in glove boxes. This is useful from a Health-Physics perspective, dealing respirable particles that can be hazardous to the human body. I iteratively tested different amounts of turbulent particles in a steady-state flow. The goal of this testing was to discover how Fluent handles built-in Rosin-Rammler distributions for particle injections. I worked on the health physics flow problems and distribution analysis under the direction of two mentors, Bruce Letellier and Dave Decroix. I set up and ran particle injection calculations using Fluent. I tried different combinations of input parameters to produce sets of 500,000, 1 million, and 1.5 million particles to determine what a good test case would be for future experiments. I performed a variety of tasks in my work as an Undergraduate Student Intern at LANL this summer, and learned how to use a powerful CFD application in addition to expanding my skills in MATLAB. I enjoyed my work at LANL and hope to be able to use the experience here to further my career in the future working in a security-conscious environment. My mentors provided guidance and help with all of my projects and I am grateful for the opportunity to work at Los Alamos National Laboratory.

36 MATERIALS SCIENCE↗
DOE OSTI2012

A Comparison of Image Quality Evaluation Techniques for Transmission X-Ray Microscopy

Beamline 6-2c at Stanford Synchrotron Radiation Lightsource (SSRL) is capable of Transmission X-ray Microscopy (TXM) at 30 nm resolution. Raw images from the microscope must undergo extensive image processing before publication. Since typical data sets normally contain thousands of images, it is necessary to automate the image processing workflow as much as possible, particularly for the aligning and averaging of similar images. Currently we align images using the 'phase correlation' algorithm, which calculates the relative offset of two images by multiplying them in the frequency domain. For images containing high frequency noise, this algorithm will align noise with noise, resulting in a blurry average. To remedy this we multiply the images by a Gaussian function in the frequency domain, so that the algorithm ignores the high frequency noise while properly aligning the features of interest (FOI). The shape of the Gaussian is manually tuned by the user until the resulting average image is sharpest. To automatically optimize this process, it is necessary for the computer to evaluate the quality of the average image by quantifying its sharpness. In our research we explored two image sharpness metrics, the variance method and the frequency threshold method. The variance method uses the variance of the image as an indicator of sharpness while the frequency threshold method sums up the power in a specific frequency band. These metrics were tested on a variety of test images, containing both real and artificial noise. To apply these sharpness metrics, we designed and built a MATLAB graphical user interface (GUI) called 'Blur Master.' We found that it is possible for blurry images to have a large variance if they contain high amounts of noise. On the other hand, we found the frequency method to be quite reliable, although it is necessary to manually choose suitable limits for the frequency band. Further research must be performed to design an algorithm which automatically selects these parameters.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2011

Space Charge Correction on Emittance Measurement of Low Energy Electron Beams

The goal of any particle accelerator is to optimize the transport of a charged particle beam along a set path by confining the beam to a small region close to the design trajectory and directing it accurately along the beamline. To do so in the simplest fashion, accelerators use a system of magnets that exert approximately linear electromagnetic forces on the charged beam. These electromagnets bend the beam along the desired path, in the case of bending magnets, and constrain the beam to the desired area through alternating focusing and defocusing effects, in the case of quadrupole magnets. We can model the transport of such a beam through transfer matrices representing the actions of the various beamline elements. However, space charge effects, produced from self electric fields within the beam, defocus the beam and must be accounted for in the calculation of beam emittance. We present below the preliminary results of a MATLAB code built to model the transport of a charged particle beam through an accelerator and measure the emittance under the influence of space charge effects. We demonstrate the method of correctly calculating the emittance of a beam under space charge effects using a least square fit to determine the initial properties of the beam given the beam size measured at a specific point after transport.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Graphite Oxidation Simulation in HTR Accident Conditions (3 rd Year and Final Technical Report)

Massive air and water ingress, following a pipe break or leak in steam-generator tubes, is a design-basis accident for high-temperature reactors (HTRs). Analysis of these accidents in both prismatic and pebble bed HTRs requires state-of-the-art capability for predictions of: 1) oxidation kinetics, 2) air helium gas mixture stratification and diffusion into the core following the depressurization, 3) transport of multi-species gas mixture, and 4) graphite corrosion. This project will develop a multi-dimensional, comprehensive oxidation kinetics model of graphite in HTRs, with diverse capabilities for handling different flow regimes. The chemical kinetics/multi-species transport model for graphite burning and oxidation will account for temperature-related changes in the properties of graphite, oxidants (O 2 , H 2 O, CO), reaction products (CO, CO 2 , H 2 , CH 4 ) and other gases in the mixture (He and N2). The model will treat the oxidation and corrosion of graphite in geometries representative of HTR core component at temperatures of 900°C or higher. The developed chemical reaction kinetics model will be user-friendly for coupling to full core analysis codes such as MELCOR and RELAP, as well as computational fluid dynamics (CFD) codes such as CD-adapco. The research team will solve governing equations for the multi-dimensional flow and the chemical reactions and kinetics using Simulink, an extension of the MATLAB solver, and will validate and benchmark the model's predictions using reported experimental data. Researchers will develop an interface to couple the validated model to a commercially available CFD fluid flow and thermal-hydraulic model of the reactor , and will perform a simulation of a pipe break in a prismatic core HTR, with the potential for future application to a pebble-bed type HTR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2012

Unique Methodologies for Nano/Micro Manufacturing Job Training Via Desktop Supercomputer Modeling and Simulation

This project establishes an initiative in high speed (Teraflop)/large-memory desktop supercomputing for modeling and simulation of dynamic processes important for energy and industrial applications. It provides a training ground for employment of current students in an emerging field with skills necessary to access the large supercomputing systems now present at DOE laboratories. It also provides a foundation for NIU faculty to quantum leap beyond their current small cluster facilities. The funding extends faculty and student capability to a new level of analytic skills with concomitant publication avenues. The components of the Hewlett Packard computer obtained by the DOE funds create a hybrid combination of a Graphics Processing System (12 GPU/Teraflops) and a Beowulf CPU system (144 CPU), the first expandable via the NIU GAEA system to ~60 Teraflops integrated with a 720 CPU Beowulf system. The software is based on access to the NVIDIA/CUDA library and the ability through MATLAB multiple licenses to create additional local programs. A number of existing programs are being transferred to the CPU Beowulf Cluster. Since the expertise necessary to create the parallel processing applications has recently been obtained at NIU, this effort for software development is in an early stage. The educational program has been initiated via formal tutorials and classroom curricula designed for the coming year. Specifically, the cost focus was on hardware acquisitions and appointment of graduate students for a wide range of applications in engineering, physics and computer science.

Kimball, Clyde↗
DOE OSTI2012

Reprocessing plant scale model integration

The national laboratory complex contains a number of niche modeling capabilities for fuel cycle facilities, at different scales and utilizing a range of codes or modeling platforms. The goal of this work has been to integrate modeling capabilities for reprocessing plants to develop a plant scale model. The Separations and Safeguards Performance Model (SSPM), developed in Matlab Simulink, has been used as a base to integrate more detailed codes. In the past, both the AMUSE and SEPHIS solvent extraction codes were integrated into the SSPM to provide more realistic modeling of the separations in aqueous reprocessing plants. For this work, the AMUSE integration was expanded to include time dependence and more fully integrated for both UREX and TRUEX contactor banks. The Nitron code is a new continuous dissolver model which has been examined for integration in the SSPM as well. Finally, in order to improve the model run time, a parallel processing capability was added. The integration progress along with initial results from testing are presented.

Cipiti, Benjamin B.↗