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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Hysteresis Effects on a Non-uniform Transmission Line with Induced Quantum Mechanical Atomic Transitions

This study describes the effect of hysteresis on the generation of the higher harmonics in the line voltage and current in a non-uniform transmission line when we include a quadratic non-linear term having memory to account for the hysteresis effect in the inductance part of line. The simulational analysis using MATLAB is presented which works with a non-uniform line having hysteresis. These simulations are based on a theoretical spatial Fourier analysis combined with perturbation theory which is used without bothering about source and load. The line differential equations are solved approximately in the space- time-frequency domain using perturbation theory and the numerical computations have been carried out for the voltage and current amplitude characteristics. In Appendix 1, some theoretical explanations that support our non-linear memory-based hysteresis model based on Landau’s equation of precession of a magnetic moment in a magnetic field are included. This paper also includes a brief derivation of atomic transition probability produced by e-m radiation from the transmission line falling on an atom in the far field zone. This calculation enables us to derive maximum likelihood estimation(MLE) of the hysteresis parameters by measuring which atomic transitions have occurred.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Error performance enhancement of satellite to earth FSO downlink employing hybrid subcarrier intensity modulation scheme with aperture averaging technique

In this work, we have evaluated the bit error rate (BER) of free space optical downlink using hybrid-subcarrier-intensity-modulation (hybrid-SIM) scheme. Also, we have employed the aperture averaging technique for the improvement of its error performance. The required BER analysis has been performed over the gamma–gamma channel model with the scintillation effects. The BER analysis is simplified using Meijer-G function and it is simulated in MATLAB. The simulation results show that the improvement in error performance is observed by hybrid-SIM technique compared to other modulation schemes like PPM and BPSK-SIM and also the BER improves significantly with the increase of aperture size due to the aperture averaging effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Design and performance analysis of all-optical 1:4 and 1:8 high speed demultiplexer using InGaAsP–InP optical microring resonator in Z-domain

All-optical 1:4 and 1:8 demultiplexers are designed and demonstrated using the process of photon absorption that is induced using an external optical pumping of wavelength 532 nm in InGaAsP–InP add-drop microring resonators. The 1:4 demultiplexer is designed with two cascaded optical microring resonators (OMRRs) and 1:8 demultiplexer is designed with four OMRRs. The mathematical modeling is pursued in Z-domain. The delay has been reduced by implementing a three level parallel circuit for 1:8 demultiplexer. The nonlinearity that occurs as a resonance frequency shift in OMRR due to photon energized electron–hole pairs is used to realize the optical switching operation of OMRR. The simulation is carried using MATLAB. The circuit has less propagation delay, fast switching time and can be integrated on CMOS photonic chips.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A Gibbs Energy Minimization Approach for Modeling of Chemical Reactions in a Basic Oxygen Furnace

In modern steelmaking, the decarburization of hot metal is converted into steel primarily in converter processes, such as the basic oxygen furnace. The objective of this work was to develop a new mathematical model for top blown steel converter, which accounts for the complex reaction equilibria in the impact zone, also known as the hot spot, as well as the associated mass and heat transport. An in-house computer code of the model has been developed in Matlab. The main assumption of the model is that all reactions take place in a specified reaction zone. The mass transfer between the reaction volume, bulk slag, and metal determine the reaction rates for the species. The thermodynamic equilibrium is calculated using the partitioning of Gibbs energy (PGE) method. The activity model for the liquid metal is the unified interaction parameter model and for the liquid slag the modified quasichemical model (MQM). The MQM was validated by calculating iso-activity lines for the liquid slag components. The PGE method together with the MQM was validated by calculating liquidus lines for solid components. The results were compared with measurements from literature. The full chemical reaction model was validated by comparing the metal and slag compositions to measurements from industrial scale converter. The predictions were found to be in good agreement with the measured values. Furthermore, the accuracy of the model was found to compare favorably with the models proposed in the literature. The real-time capability of the proposed model was confirmed in test calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Binary Alkali Nitrate and Chloride Phase Diagrams: NaNO{sub 3}-KNO{sub 3}, LiNO{sub 3}-NaNO{sub 3}, LiNO{sub 3}-KNO{sub 3}, and NaCl-KCl

LiNO{sub 3}-NaNO{sub 3}, LiNO{sub 3}-KNO{sub 3,} NaNO{sub 3}-KNO{sub 3}, and NaCl-KCl phase relations were formulated using the regular solution model to obtain desired compositions for their potential application as a heat transfer fluid. The simultaneous nonlinear equations formed were solved numerically using a contour plot technique in Matlab. An algorithm was developed to determine the real roots of the coupled nonlinear equations and subsequent phase diagrams were constructed. The phase diagrams were compared with those predicted using FactSage thermochemical package and literature data from experimental studies. The calculated NaNO{sub 3}-KNO{sub 3} phase diagram shows that this system has a eutectic point at 0.5 mole fraction of KNO{sub 3} with a eutectic temperature of 220.85 °C, while the LiNO{sub 3}-KNO{sub 3} system was found to have a eutectic point at 0.44 mole fraction LiNO{sub 3} at 137 °C. The LiNO{sub 3}-NaNO{sub 3} system shows a eutectic at 0.524 mole fraction LiNO{sub 3} at 192.8 °C, and finally the NaCl-KCl system has a eutectic point at 0.48 mole fraction NaCl at 634 °C. The liquidus and solidus curves obtained fitted close against values obtained from the quasi-chemical model used in FactSage and experimental values from the literature. This study showed that the regular solution model combined with the contour plots approach can sufficiently describe the liquidus–solidus curves of the alkali binary nitrate and chloride systems and could be a useful method for predicting the phase diagrams for higher order alkali nitrates and chlorides.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Determination of Material Constitutive Laws for Inconel 718 Superalloy Under Different Strain Rates and Working Temperatures

In this paper, a special procedure for the prediction of parameters of the Johnson–Cook constitutive material models is proposed based on the experimental data and specially developed MATLAB scripts which allow advanced modeling of complex 3D response surfaces. Experimental investigations concern two various strain rates of 10{sup −3} and 10{sup 1} 1/s and the testing temperature ranging from the ambient up to 700 °C. As a result, a set of mathematical equations which fit the experimental data is determined. The applicability of the experimentally derived constitutive models to the FEM modeling of real machining processes of Inconel 718 alloy is verified.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effect of Grain Boundaries on the Performance of Thin-Film-Based Polycrystalline Silicon Solar Cells: A Numerical Modeling

Solar cells/photovoltaic, a renewable energy source, is appraised to be the most effective alternative to the conventional electrical energy generator. A cost-effective alternative of crystalline wafer-based solar cell is thin-film polycrystalline-based solar cell. This paper reports the numerical analysis of dependency of the solar cell parameters (i.e., efficiency, fill factor, open-circuit voltage and short-circuit current density) on grain size for thin-film-based polycrystalline silicon (Si) solar cells. A minority carrier lifetime model is proposed to do a correlation between the grains, grain boundaries and lifetime for thin-film-based polycrystalline Si solar cells in MATLAB environment. As observed, the increment in the grain size diameter results in increase in minority carrier lifetime in polycrystalline Si thin film. A non-equivalent series resistance double-diode model is used to find the dark as well as light (AM1.5) current–voltage (I-V) characteristics for thin-film-based polycrystalline Si solar cells. To optimize the effectiveness of the proposed model, a successive approximation method is used and the corresponding fitting parameters are obtained. The model is validated with the experimentally obtained results reported elsewhere. The experimentally reported solar cell parameters can be found using the proposed model described here.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Experimental Study and Fractal Analysis on the Anisotropic Performance of Explosively Welded Interfaces of 304 Stainless Steel/245 Carbon Steel

The mechanical anisotropy of an explosive welding composite plate made of 304 stainless steel/245 steel was studied through shear experiments performed on explosively welded wavy interfaces along several orientation angles. The results indicated that the strength and the fracture energy of samples significantly varied with the orientation angles. The fracture surfaces of all samples were observed using a scanning electron microscope and through three-dimensional structure microscopy. The periodic features of all the fracture surfaces were clearly shown in different fracture modes. The fractal dimension of the fracture surfaces was calculated based on the fractal geometry by the box-counting method in MATLAB. The cohesive element model was used to analyze the fracture energy according to the physical dependence of the fractal dimension on thermodynamic entropy and interface separation energy. The fracture energy was an exponential function of the fractal dimension value, which was in good agreement with the experimental results. All results were validated for effective use in the application of anisotropy analysis to the welded interface and structural optimization of explosively welded composite plates.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Prediction and Computation of Corrosion Rates of A36 Mild Steel in Oilfield Seawater

The parameters which primarily control the corrosion rate and life of steel structures are several and they vary across the different ocean and seawater as well as along the depth. While the effect of single parameter on corrosion behavior is known, the conjoint effects of multiple parameters and the interrelationship among the variables are complex. Millions sets of experiments are required to understand the mechanism of corrosion failure. Statistical modeling such as ANN is one solution that can reduce the number of experimentation. ANN model was developed using 170 sets of experimental data of A35 mild steel in simulated seawater, varying the corrosion influencing parameters SO{sub 4}{sup 2−}, Cl{sup −}, HCO{sub 3}{sup −},CO{sub 3}{sup 2−}, CO{sub 2}, O{sub 2}, pH and temperature as input and the corrosion current as output. About 60% of experimental data were used to train the model, 20% for testing and 20% for validation. The model was developed by programming in Matlab. 80% of the validated data could predict the corrosion rate correctly. Corrosion rates predicted by the ANN model are displayed in 3D graphics which show many interesting phenomenon of the conjoint effects of multiple variables that might throw new ideas of mitigation of corrosion by simply modifying the chemistry of the constituents. The model could predict the corrosion rates of some real systems.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Microstructural Analysis of Powder Metallurgy Tool Steels in the Context of Abrasive Wear Behavior: A New Computerized Approach to Stereology

The present work describes a new methodology designed to characterize the microstructures of tool steels containing carbide hard phases, with the focus set on their abrasive wear resistance. A series of algorithms were designed and implemented in MATLAB® to (i) recognize each of the features of interest, (ii) measure relevant quantities and (iii) characterize each of the phases and the alloy in function of attributes usually neglected in wear description applications: size distribution, shape and contiguity of the hard phases. The new framework incorporates new parameters to describe each one of these attributes, as observed in SEM micrographs. All three aforementioned stages contain novel contributions that can be potentially beneficial to the field of materials design in general and to the field of alloy design for severely abrasive environments in particular. Models of known geometry and micrographs of different powder metallurgy steels were analyzed, and the obtained results were compared with the obtained by the linear intercept method. The relation between the new parameters and the ones available in the scientific literature is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Hybrid Inverse Parameter Identification of Fully Coupled Ductile Damage Model for Steel Sheet DP600 with Two Different Algorithms: Trust Region and Genetic Algorithms

Ductile damage model has been widely regarded as a valuable method to predict the failure of sheet metal. Based on the thermodynamic theory and continuum damage mechanics, the fully coupled ductile damage model can be developed, which also can better predict the initiation and growth of the fracture. But the identifications of model parameters with theoretical methodology are difficult due to the complex coupling relationships existing among all state variables. The inverse methodology is regarded as a good method to resolve the problem. In this paper, the recently proposed fully coupled ductile damage model is chosen to investigate the deformation behavior of DP600 steel, in which the mixed saturation isotropic and kinematic hardenings are taken into account and fully coupled with the ductile damage. During the identification process, the least square formula of the error between experimental and numerical results is selected as the target function. The trust region algorithm and genetic algorithm are used with the help of MATLAB for the identification of three damage parameters. Finally, by comparing the experimental and numerical results, the validations of two algorithms are proved. The efficiency of the optimization process with trust region algorithm is higher, but with lower accuracy. Meanwhile, the optimization process is greatly affected by the chosen initial values of the ductile damage parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An Automated, High-Throughput Method for Interpreting the Tandem Mass Spectra of Glycosaminoglycans

The biological interactions between glycosaminoglycans (GAGs) and other biomolecules are heavily influenced by structural features of the glycan. The structure of GAGs can be assigned using tandem mass spectrometry (MS{sup 2}), but analysis of these data, to date, requires manually interpretation, a slow process that presents a bottleneck to the broader deployment of this approach to solving biologically relevant problems. Automated interpretation remains a challenge, as GAG biosynthesis is not template-driven, and therefore, one cannot predict structures from genomic data, as is done with proteins. The lack of a structure database, a consequence of the non-template biosynthesis, requires a de novo approach to interpretation of the mass spectral data. We propose a model for rapid, high-throughput GAG analysis by using an approach in which candidate structures are scored for the likelihood that they would produce the features observed in the mass spectrum. To make this approach tractable, a genetic algorithm is used to greatly reduce the search-space of isomeric structures that are considered. The time required for analysis is significantly reduced compared to an approach in which every possible isomer is considered and scored. The model is coded in a software package using the MATLAB environment. This approach was tested on tandem mass spectrometry data for long-chain, moderately sulfated chondroitin sulfate oligomers that were derived from the proteoglycan bikunin. The bikunin data was previously interpreted manually. Our approach examines glycosidic fragments to localize SO{sub 3} modifications to specific residues and yields the same structures reported in literature, only much more quickly. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Evaluation of Digital Image Recognition Methods for Mass Spectrometry Imaging Data Analysis

Analyzing mass spectrometry imaging data can be laborious and time consuming, and as the size and complexity of datasets grow, so does the need for robust automated processing methods. We here present a method for comprehensive, semi-targeted discovery of molecular distributions of interest from mass spectrometry imaging data, using widely available image similarity scoring algorithms to rank images by spatial correlation. A fast and powerful batch search method using a MATLAB implementation of structural similarity (SSIM) index scoring with a pre-selected reference distribution is demonstrated for two sample imaging datasets, a plant metabolite study using Artemisia annua leaf, and a drug distribution study using maraviroc-dosed macaque tissue. .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A modified homotopy method for solving the principal-agent bilevel programming problem

In this paper, the principal-agent bilevel programming problem with integral operator is considered, in which the upper-level object is the agent that maximizes its expected utility with respect to an agreed compensation contract. The constraints are the principal’s participation and the agent’s incentive compatibility. The latter is a lower-level optimization problem with respect to its private action. To solve an equivalent single-level nonconvex programming problem with integral operator, a modified homotopy method for solving the Karush–Kuhn–Tucker system is proposed. This method requires only an interior point and, not necessarily, a feasible initial approximation for the constraint shifting set. Global convergence is proven under some mild conditions. Numerical experiments were performed by our homotopy method as well as by fmincon in Matlab, LOQO and MINOS. The experiments showed that: designing a piecewise linear contract is much better than designing a piecewise constant contract and only needs to solve a much lower-dimensional optimization problem and hence needs much less computation time; the optimal value of the principal-agent model with designing piecewise linear contract tends to a limitation, while the discrete segments gradually increase; and finally, the proposed modified homotopy method is feasible and effective.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

J-variational inequalities and zeroes of a family of maximal monotone operators by sunny generalized nonexpansive retraction

In this paper, using generalized resolvents and retractions onto shrinking closed subsets, we present new extragradient algorithms for finding the solution of a J-variational-like inequality, the zero points of a family of maximal monotone operators, and the solution set of mixed equilibrium problems for a J-α-inverse-strongly monotone-like operator in a 2-uniformly convex and 2-uniformly smooth Banach space. By introducing the new definitions, we prove strong convergence of generated iterates in the extragradient methods. Using FMINCON optimization toolbox in MATLAB, we give some numerical examples and compare them with several existence results in the literature to illustrate the usability of our results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Initial conceptual demonstration of control co-design for WEC optimization

Abstract While some engineering fields have benefited from systematic design optimization studies, wave energy converters have yet to successfully incorporate such analyses into practical engineering workflows. The current iterative approach to wave energy converter design leads to sub-optimal solutions. This short paper presents an open-source MATLAB toolbox for performing design optimization studies on wave energy converters where power take-off behavior and realistic constraints can be easily included. This tool incorporates an adaptable control co-design approach, in that a constrained optimal controller is used to simulate device dynamics and populate an arbitrary objective function of the user’s choosing. A brief explanation of the tool’s structure and underlying theory is presented. To demonstrate the capabilities of the tool, verify its functionality, and begin to explore some basic wave energy converter design relationships, three conceptual case studies are presented. In particular, the importance of considering (and constraining) the magnitudes of device motion and forces in design optimization is shown.

16 TIDAL AND WAVE POWER↗
DOE OSTI2018

Analysis of N-dimensional Klein–Gordon equation for hydrogen molecule in the non-central potential field

We investigate a complete normalized relativistic solutions for N-dimensional Hartmann potential and relativistic effects of non-central potential field for hydrogen molecule within the framework of asymptotic iteration approach. We also discuss in detail the influence of dimensions on the angular probability distribution of hydrogen molecule. N-dimensional bound state wavefunction solutions including angular and radial solutions are expressed in terms of the Jacobi and Laguerre polynomials, respectively. N-dimensional angular wavefunctions of hydrogen molecule were visualized with Matlab software. In hyperspherical coordinates, the recurrence relations between diagonal matrix elements have been derived for N-dimensional Hartmann potential in relativistic theory. We also obtain analytical expressions of diagonal matrix elements and the recurrence formula between off-diagonal matrix elements in N-dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗