Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix multiplication”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Post-irradiation examination of the Sirius-1 nuclear thermal propulsion fuel test

Nuclear Thermal Propulsion (NTP) systems hold promise in reducing transit times for exploration of the solar system by both crewed and uncrewed missions. NTP systems currently under investigation include a once-through high-temperature gas-cooled fission reactor to provide thermal energy to heat the coolant, which also serves as the propellant. The fuel system consists of angular UN fuel particles dispersed in a matrix of W/Re, creating a ceramic and metallic composite or cermet, which has been irradiated in Idaho National Laboratory's (INL) Transient Reactor Test Facility (TREAT). This enabled evaluation of these materials under representative nuclear heating rates (~95 K/s) and peak surface temperatures (~2527 K). Further, these surface conditions were indicative of achievement of the target peak internal temperature of approximately 2850 K. These tests named Sirius-1 (UN–W/Re), have been irradiated and this paper will present post-irradiation examination results. The Sirius-1 test produced minor, stable cracks in the fuel specimen and spalling of surface material. Volatilized uranium ‘soot’ was found deposited on the inner wall of the irradiation capsule indicating loss of some fissile material from the fuel specimen. Spalling from the surfaces was also noted upon visual inspection. Uranium diffusion from the fuel particles resulted in the formation of U/Re phases and the production of a laminar microstructure on the edge of the fuel system. Overall, the specimen was stable and remained within a coolable geometry upon conclusion of multiple thermal cycles to prototypical operating temperatures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2002

Development of Three-Dimensional Multicellular Tissue-Like Constructs for Mutational Analysis Using Macroporous Microcarriers

A three-dimensional (3-D), tissue-like model was developed for the genotoxic assessment of space environment. In previous experiments, we found that culturing mammalian cells in a NASA-designed bioreactor, using Cytodex-3 beads as a scaffold, generated 3-D multicellular spheroids. In an effort to generate scaffold-free spheroids, we developed a new 3-D tissue-like model by coculturing fibroblast and epithelial cell in a NASA bioreactor using macroporous Cultispher-S(TradeMark) microcarriers. Big Blue(Registered Trademark) Rat 2(Lambda) fibroblasts, genetically engineered to contain multiple copies (>60 copies/cell) of the Lac I target gene, were cocultured with radio-sensitive human epithelial cells, H184F5. Over an 8-day period, samples were periodically examined by microscopy and histology to confirm cell attachment, growth, and viability. Immunohistochemistry and western analysis were used to evaluate the expression of specific cytoskeletal and adhesion proteins. Key cell culture parameters (glucose, pH, and lactate concentrations) were monitored daily. Controls were two-dimensional mono layers of fibroblast or epithelial cells cultured in T-flasks. Analysis of 3-D spheroids from the bioreactor suggests fibroblast cells attached to and completely covered the bead surface and inner channels by day 3 in the bioreactor. Treatment of the 3-day spheroids with dispase II dissolved the Cultisphers(TradeMark) and produced multicellular, bead-less constructs. Immunohistochemistry confirmed the presence of vi.mentin, cytokeratin and E-cadherin in treated spheroids. Examination of the dispase II treated spheroids with transmission electron microscopy (TEM) also showed the presence of desmosomes. These results suggest that the controlled enzymatic degradation of an artificial matrix in the low shear environment of the NASA-designed bioreactor can produce 3-D tissue-like spheroids. 2

Jordan, Jacqueline A.↗
DOE OSTI2023

Investigation of desorption of hydrogen gas from polymer matrix using thermal desorption analysis and finite element modeling

Polymer materials are used widely as sealing materials and liners in hydrogen storage and transportation infrastructure. During operating conditions, polymers are exposed to high pressure hydrogen gas periodically. When in contact, the hydrogen gas diffuses through the polymer material with ease due to its smaller sized molecules compared to the polymer molecules. In the event of rapid decompression of outside hydrogen gas, the hydrogen molecules within the polymer diffuse out. This process is known as the desorption of polymer and can cause blisters or cracking causing permanent damage to the polymer. Here, we used a combined experimental and modeling approach to investigate the desorption of the hydrogen gas in the polymer. Thermal desorption analysis (TDA) experiments were performed on the polymer sample exposed to the high-pressure hydrogen gas to measure the diffusion coefficient and equilibrium hydrogen concentration. Experiments also provided the quantity of hydrogen within the polymer during the desorption process which followed Fick’s law. However, for higher pressurization values, some anomaly was observed in the quantity of hydrogen remaining in the polymer sample. A continuum mechanics based coupled diffusion–deformation-damage model was deployed to simulate the gas desorption process in the material using the finite element method to understand its effect on the material. It was found that the anomaly in TDA results were due to the blister formation on the polymer surface which resulted in longer hydrogen retention. Simulations also showed that blister formation can be affected by multiple parameters such as the addition of filler particles, and pre-existing cavity size.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

A New Concurrent Engineering Tool for the New Mission Design Center at NASA Ames Research Center

The NASA Ames Mission Design Center (MDC) has undergone a significant transformation process in recent years, culminating in becoming a true Concurrent Engineering Center. The underlying goal of this transformation is to better serve the Ames PI community in a cost-effective and rapidly-responsive way, providing quality products pivotal to the decision-making proposal strategy of our center. The four pillars of the change focus on four areas: personnel, physical places, tools, and training. Personnel: the MDC switched from using a dedicated pool of engineers—mostly entry-level career—for a matrixed approach with more experienced Subject Matter Experts, usually involved in spaceflight projects. These SMEs get detailed to the MDC for the duration of the study on a part-time or full-time basis. Physical places: The physical aspect pertains mostly to the layout of the new Engineering building and the dedicated concurrent engineering layout to foster collaboration. The current plan includes a dedicated space for concurrent engineering sessions as well as dedicated project rooms for concept studies. Tool: the MDC team is developing a new in-house concurrent engineering tool to facilitate concurrent engineering sessions. The team leveraged the many lessons learned and insight gained from using the current concurrent engineering tool, in use since 2007. SMEs have been involved throughout its design and development, ensuring accuracy and providing validity to the implementation approach. This new tool is database-driven and allows a team to work concurrently on the same model of a mission concept. “Workspaces” are provided for each of the traditional spacecraft mission design disciplines, including cost accounting and systems engineering. There are multiple tiers of fidelity available within each workspace, which can be completed parametrically or independently from the other disciplines as needed or synchronously with the other disciplines. A “commit” step synchronizes a workspace with the mission database and identifies any potential conflicts, along with the user who entered the conflicting data. A Mission Summary workspace enables a Study Lead to run an effective CE session, displaying key graphics, workspace status information, and technical and monetary budget information. This new CE tool will be undergoing Beta Testing with its primary user base and a limited Parts database in the summer of 2020. The team plans to release the first fully operational version in December 2020. Training: working in early concept development and maturation for a space mission requires a unique mindset and being comfortable with uncertainties in an environment where requirements are not yet fully developed, and where changes and trade-offs happen at a high tempo. The MDC started a weekly training program to ensure a common knowledge base on how do develop and mature early concept studies, as well as how to operate in a concurrent engineering environment. In conclusion, the MDC transformation will provide the overall PI community focused on cost-effective small satellite scientific missions with a powerful tool to address and mature early concepts more effectively and efficiently.

David Mauro↗
NASA NTRS2016

Analysis and Characterization of Damage and Failure Utilizing a Generalized Composite Material Model Suitable for Use in Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in state-of-the art composite impact models is under development. In particular, a next generation composite impact material model, jointly developed by the FAA and NASA, is being implemented into the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage, and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters (such as modulus and strength). The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in the various coordinate directions. Due to the fact that the plasticity and damage models are uncoupled, test procedures and methods to both characterize the damage model and to covert the material stress-strain curves from the true (damaged) stress space to the effective (undamaged) stress space have been developed. A methodology has been developed to input the experimentally determined composite failure surface in a tabulated manner. An analytical approach is then utilized to track how close the current stress state is to the failure surface.

Finite Element Method↗
DOE OSTI2024

Application Guide for the Cyber-Resilient Design Framework for Hybrid Systems

As the energy landscape evolves, hybrid power plants—integrating multiple renewable energy sources (e.g., solar, wind, and battery storage) with the bulk power electric system—play a crucial role in meeting growing energy demands. However, with a hybrid power plant’s increased number of components, complex network connectivity and digitization of controls, these systems face growing cybersecurity risks. To help mitigate these risks, Idaho National Laboratory (INL) presents this application guide specifically designed for operators and systems engineers to evaluate the cybersecurity resilience of their hybrid power plant’s design.

14 SOLAR ENERGY↗
DOE OSTI2022

Sulphur variations in annually layered stalagmites using benchtop micro-XRF

Variation of sulphur in annually laminated stalagmites can be used to infer the impact of past volcanic activities, anthropogenic pollution, and climate change due to increased bushfire activity. The synchrotron radiation micro-X-Ray fluorescence (SR-XRF) microprobe is a powerful tool to analyse and image sulphur recorded in stalagmites with micrometre resolution. However, access to SR-XRF beamlines can be limited, so researchers must select the most promising stalagmites for imaging. Benchtop micro-XRF is an effective tool for trace elemental analysis of speleothem samples and is a candidate for routine laboratory measurement of sulphur along stalagmite laminae and screening for SR-XRF. Here, this study describes a protocol using matrix-matched standards to measure annual variations of sulphur at trace to percent level along the laminae of two Western Australian stalagmites, one of which already having been analysed using SR-XRF. Parameters that affected quantitation include X-ray tube voltage and current, spot size of the X-ray beam and stalagmite surface roughness and porosity. The use of a 20 μm X-ray spot size provides sub-annual spatial resolution that can be completed in an overnight scan. The features in a 1000 point micro-XRF analysis of sulphur along a 20 mm transect show good consistency with SR-XRF microprobe data. Micro-XRF mapping was also performed to produce chemical images on the stalagmite and compared with Raman and X-ray diffraction (XRD) to confirm that the stalagmite is exclusively calcite, with no aragonite, and that the source of sulphur in the samples was gypsum and anhydrite. Regions of very high sulphur in the micro-XRF maps were found to be artefacts due to diffraction of the incident beam but these could be efficiently removed by using a multiple point statistics approach to produce a clean image suitable for analysis of the laminae. This work shows the potential of micro-XRF for routine analysis of sulphur in stalagmites, and to streamline sample characterisation before SR-XRF imaging.

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

High-performance strategies for the recent MRSF-TDDFT in GAMESS

Multiple ERI (Electron Repulsion Integral) tensor contractions (METC) with several matrices are ubiquitous in quantum chemistry. In response theories, the contraction operation, rather than ERI computations, can be the major bottleneck, as its computational demands are proportional to the multiplicatively combined contributions of the number of excited states and the kernel pre-factors. Here, this paper presents several high-performance strategies for METC. Optimal approaches involve either the data layout reformations of interim density and Fock matrices, the introduction of intermediate ERI quartet buffer, and loop-reordering optimization for a higher cache hit rate. The combined strategies remarkably improve the performance of the MRSF (mixed reference spin flip)-TDDFT (time-dependent density functional theory) by nearly 300%. The results of this study are not limited to the MRSF-TDDFT method and can be applied to other METC scenarios.

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

Assessing sensitivities of climate model weighting to multiple methods, variables, and domains in the south-central United States

Given the increasing use of climate projections and multi-model ensemble weighting for a diverse array of applications, this project assesses the sensitivities of climate model weighting strategies and their resulting ensemble means to multiple components, such as the weighting schemes, climate variables, or spatial domains of interest. The purpose of this study is to assess the sensitivities associated with multi-model weighting strategies. The analysis makes use of global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and their statistically downscaled counterparts created with the localized constructed analogs (LOCA) method. This work focuses on historical and projected future mean precipitation and daily high temperatures of the south-central United States. Results suggest that the model weights and the corresponding weighted model means can be sensitive to the weighting strategy that is applied. For instance, when estimating model weights based on Louisiana precipitation, the weighted projections show a wetter and cooler south-central domain in the future compared to other weighting strategies. Alternatively, for example, when estimating model weights based on New Mexico temperature, the weighted projections show a drier and warmer south-central domain in the future. However, when considering the entire south-central domain in estimating the model weights, the weighted future projections show a compromise in the precipitation and temperature estimates. As for uncertainty, our matrix of results provided a more certain picture of future climate compared to the spread in the original model ensemble. If future impact assessments utilize weighting strategies, then our findings suggest that how the specific weighting strategy is used with climate projections may depend on the needs of an impact assessment or adaptation plan.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Effect of multiple strain-anneal cycles on the 1000 C creep behaviour of gamma/gamma prime-alpha

Various multiple strain-anneal cycles (1000 C) were imposed on specimens of the directionally solidified eutectic (DSE) alloy gamma/gamma prime-alpha to identify thermomechanical processing methods (TMP) which would improve the creep behavior. Specimens of the Ni-32.3Mo-6.3Al wt pct alloy were grown with a modified Bridgeman technique. Some of the cylindrical specimens were alternately heat-treated at 900 C, then strained, or heat-treated only, while other specimens were annealed at 900 C after swaging and then worked at ambient temperature. The specimens were all examined microstructurally using transmission electron microscopy, some before and after being exposed to constant-load compression tests at 1000 C. The creep strain increased for all TMP specimens for strain rates of at least 2 millionths per sec. Strain rates of about 2 ten millionths per sec were only improved with strain annealing with 13 percent work at ambient temperature. A slight improvement, compared to as-grown materials, was observed in the 1000 C creep behavior of materials annealed at 900 C. Strain-annealing was found to introduce three-dimensional dislocation networks into the gamma-prime matrix.

Whittenberger, J. D.↗
NASA NTRS1997

Dynamic Reconstruction and Multivariable Control for Force-Actuated, Thin Facesheet Adaptive Optics

The Multiple Mirror Telescope (MMT) under development at the University of Arizona takes a new approach in adaptive optics placing a large (0.65 m) force-actuated, thin facesheet deformable mirror at the secondary of an astronomical telescope, thus reducing the effects of emissivity which are important in IR astronomy. However, The large size of the mirror and low stiffness actuators used drive the natural frequencies of the mirror down into the bandwidth of the atmospheric distortion. Conventional adaptive optics takes a quasi-static approach to controlling the, deformable mirror. However, flexibility within the control bandwidth calls for a new approach to adaptive optics. Dynamic influence functions are used to characterize the influence of each actuator on the surface of the deformable mirror. A linearized model of atmospheric distortion is combined with dynamic influence functions to produce a dynamic reconstructor. This dynamic reconstructor is recognized as an optimal control problem. Solving the optimal control problem for a system with hundreds of actuators and sensors is formidable. Exploiting the circularly symmetric geometry of the mirror, and a suitable model of atmospheric distortion, the control problem is divided into a number of smaller decoupled control problems using circulant matrix theory. A hierarchic control scheme which seeks to emulate the quasi-static control approach that is generally used in adaptive optics is compared to the proposed dynamic reconstruction technique. Although dynamic reconstruction requires somewhat more computational power to implement, it achieves better performance with less power usage, and is less sensitive than the hierarchic technique.

Grocott, Simon C. O.↗
NASA NTRS2020

Analysis of the Honeywell Uncertified Research Engine with Ice Crystal Cloud Ingestion at Simulated Altitudes

The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. To develop the test matrix of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The goal of the test matrix was to provide operating conditions such that ice would accrete in either the fan-stator through the inlet guide vane region of the compression system or within the first stator of the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDESMELT) which includes an ice particle model. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method (COMDES-MELT) was enhanced by computing key parameters through the fan-stator at multiple span wise locations, in order to increase the fidelity with the current mean-line method. The Icing Wedge static wet-bulb temperature thresholds were applied for determining the risk of ice accretion in the fan-stator, which is thought to be an adiabatic region. At some operating conditions near the splitter-lip region, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion, and the Icing Wedge was not applied to predict accretion at that location. A simple order-of-magnitude heat transfer model was implemented into the COMDES-MELT code to estimate the wall temperature minimum and maximum thresholds that support ice accretion, as observed by video confirmation. The results from this model spanned the range of wall temperatures measured on a previous engine that experienced ice accretion at certain operating conditions. The goal of this study is to show that the computational process developed on earlier engine icing tests can be used to provide an icing risk assessment in adiabatic regions for other engines.

Turbofan engine↗
DOE OSTI2020

Measurement of the 239 Pu$(n, f$) prompt fission neutron spectrum from 10 keV to 10 MeV induced by neutrons of energy 1–20 MeV

Although the prompt fission neutron spectrum (PFNS) is an essential component of neutron-driven systems that has been measured for decades, there are still multiple glaring unknowns regarding the PFNS of major actinides in the fission neutron incident energy range, specifically with regard to multichance fission and pre-equilibrium neutron emission processes. The only impactful experimental 239 Pu PFNS measurements included in recent nuclear data evaluations were measured over a limited outgoing neutron energy range at thermal and 1.5-MeV average incident neutron energy, while other potentially impactful measurements have been shown to contain errors that resulted in either large uncertainty increases or in complete exclusion from nuclear data evaluation. Here, we report a measurement of the 239 Pu PFNS over a wide range of incident neutron energy (1–20 MeV) and three orders of magnitude in outgoing neutron energy (0.01–10 MeV) resulting from the Chi-Nu experiment at the Los Alamos Neutron Science Center. These results are the combination of separate PFNS measurements in the same experimental area, one using a Li-glass and the other a liquid scintillator detector array. Covariances between all PFNS data points from each detector and within each incident energy range were generated between all other data in both detector arrays and within all other incident neutron energy bins, yielding a single covariance matrix for all 1300 PFNS data points reported here. These covariances are based on a thorough assessment of systematic bias and uncertainties associated with the measurement, PFNS extraction technique, combination of data from each detector type, and other aspects of the analysis. The existence of covariances between PFNS data points in different incident neutron energy ranges yielded covariances between average PFNS energy values at each incident energy to be reported here as well, which allowed for firm statements to be made regarding a shape of a purely experimental mean PFNS energy trend for the first time. Although minor PFNS shape differences exist between the results reported here and recent nuclear data evaluations, the ENDF/B-VIII.0 and JEFF-3.3 PFNS evaluations agree reasonably well with the present results from 1-to 10-MeV incident neutron energy, which spans the well-measured 1.5-MeV incident neutron energy PFNS from Lestone and Shores as well as the onset of second-chance fission. However, while the pre-equilibrium component of the PFNS above 12-MeV incident neutron energy roughly agrees in position and magnitude with ENDF/B-VIII.0 and JEFF-3.3, clear differences relating to the relative magnitude of third-chance fission PFNS features are present in the PFNS shape and in the mean PFNS energy trends.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2004

Error Modeling of Multi-baseline Optical Truss: Application to SIM Metrology Truss Field Dependent Error - Part II

The current design of the Space Interferometry Mission (SIM) employs a 19 laser-metrology-beam system (also called L19 external metrology truss) to monitor changes of distances between the fiducials of the flight system's multiple baselines. The function of the external metrology truss is to aid in the determination of the time-variations of the interferometer baseline. The largest contributor to truss error occurs in SIM wide-angle observations when the articulation of the siderostat mirrors (in order to gather starlight from different sky coordinates) brings to light systematic errors due to offsets at levels of instrument components (which include comer cube retro-reflectors, etc.). This error is labeled external metrology wide-angle field-dependent error. Physics-based model of field-dependent error at single metrology gauge level is developed and linearly propagated to errors in interferometer delay. In this manner delay error sensitivity to various error parameters or their combination can be studied using eigenvalue/eigenvector analysis. Also validation of physics-based field-dependent model on SIM testbed lends support to the present approach. As a first example, dihedral error model is developed for the comer cubes (CC) attached to the siderostat mirrors. Then the delay errors due to this effect can be characterized using the eigenvectors of composite CC dihedral error. The essence of the linear error model is contained in an error-mapping matrix. A corresponding Zernike component matrix approach is developed in parallel, first for convenience of describing the RMS of errors across the field-of-regard (FOR), and second for convenience of combining with additional models. Average and worst case residual errors are computed when various orders of field-dependent terms are removed from the delay error. Results of the residual errors are important in arriving at external metrology system component requirements. Double CCs with ideally co-incident vertices reside with the siderostat. The non-common vertex error (NCVE) is treated as a second example. Finally combination of models, and various other errors are discussed.

comer cube retro-reflector↗
DOE OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
NASA NTRS2010

Fast Solution in Sparse LDA for Binary Classification

An algorithm that performs sparse linear discriminant analysis (Sparse-LDA) finds near-optimal solutions in far less time than the prior art when specialized to binary classification (of 2 classes). Sparse-LDA is a type of feature- or variable- selection problem with numerous applications in statistics, machine learning, computer vision, computational finance, operations research, and bio-informatics. Because of its combinatorial nature, feature- or variable-selection problems are NP-hard or computationally intractable in cases involving more than 30 variables or features. Therefore, one typically seeks approximate solutions by means of greedy search algorithms. The prior Sparse-LDA algorithm was a greedy algorithm that considered the best variable or feature to add/ delete to/ from its subsets in order to maximally discriminate between multiple classes of data. The present algorithm is designed for the special but prevalent case of 2-class or binary classification (e.g. 1 vs. 0, functioning vs. malfunctioning, or change versus no change). The present algorithm provides near-optimal solutions on large real-world datasets having hundreds or even thousands of variables or features (e.g. selecting the fewest wavelength bands in a hyperspectral sensor to do terrain classification) and does so in typical computation times of minutes as compared to days or weeks as taken by the prior art. Sparse LDA requires solving generalized eigenvalue problems for a large number of variable subsets (represented by the submatrices of the input within-class and between-class covariance matrices). In the general (fullrank) case, the amount of computation scales at least cubically with the number of variables and thus the size of the problems that can be solved is limited accordingly. However, in binary classification, the principal eigenvalues can be found using a special analytic formula, without resorting to costly iterative techniques. The present algorithm exploits this analytic form along with the inherent sequential nature of greedy search itself. Together this enables the use of highly-efficient partitioned-matrix-inverse techniques that result in large speedups of computation in both the forward-selection and backward-elimination stages of greedy algorithms in general.

Moghaddam, Baback↗