Search NASA⌕ Search

SEARCH · Search NASA

Results for “LINEAR PROGRAMMING”

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 649 records · Page 36

High Gradient Accelerator Research

The goal of the MIT program of research on high gradient acceleration is the development of advanced acceleration concepts that lead to a practical and affordable next generation linear collider at the TeV energy level. Other applications, which are more near-term, include accelerators for materials processing; medicine; defense; mining; security; and inspection. The specific goals of the MIT program are: • Pioneering theoretical research on advanced structures for high gradient acceleration, including photonic structures and metamaterial structures; evaluation of the wakefields in these advanced structures • Experimental research to demonstrate the properties of advanced structures both in low-power microwave cold test and high-power, high-gradient test at megawatt power levels • Experimental research on microwave breakdown at high gradient including studies of breakdown phenomena induced by RF electric fields and RF magnetic fields; development of new diagnostics of the breakdown process • Theoretical research on the physics and engineering features of RF vacuum breakdown • Maintaining and improving the Haimson / MIT 17 GHz accelerator, the highest frequency operational accelerator in the world, a unique facility for accelerator research • Providing the Haimson / MIT 17 GHz accelerator facility as a facility for outside users • Active participation in the US DOE program of High Gradient Collaboration, including joint work with SLAC and with Los Alamos National Laboratory; participation of MIT students in research at the national laboratories • Training the next generation of Ph. D. students in the field of accelerator physics.

43 PARTICLE ACCELERATORS↗

Optimal design and dispatch of a system of diesel generators, photovoltaics and batteries for remote locations

Renewable energy technologies, specifically, solar photovoltaic cells, combined with battery storage and diesel generators, form a hybrid system capable of independently powering remote locations, i.e., those isolated from larger grids. If sized correctly, hybrid systems reduce fuel consumption compared to diesel generator-only alternatives. We present an optimization model for establishing a hybrid power design and dispatch strategy for remote locations, such as a military forward operating base, that models the acquisition of different power technologies as integer variables and their operation using nonlinear expressions. Our cost-minimizing, nonconvex, mixed-integer, nonlinear program contains a detailed battery model. Due to its complexities, we present linearizations, which include exact and convex under-estimation techniques, and a heuristic, which determines an initial feasible solution to serve as a “warm start” for the solver. We determine, in a few hours at most, solutions within 5% of optimality for a candidate set of technologies; these solutions closely resemble those from the nonlinear model. Lastly, our instances contain real data spanning a yearly horizon at hour fidelity and demonstrate that a hybrid system could reduce fuel consumption by as much as 50% compared to a generator-only solution.

14 SOLAR ENERGY↗

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An inertia-free filter line-search algorithm for large-scale nonlinear programming

We present a filter line-search algorithm that does not require inertia information of the linear system. This feature enables the use of a wide range of linear algebra strategies and libraries, which is essential to tackle large-scale problems on modern computing architectures. The proposed approach performs curvature tests along the search step to detect negative curvature and to trigger convexification. We prove that the approach is globally convergent and we implement the approach within a parallel interior-point framework to solve large-scale and highly nonlinear problems. Our numerical tests demonstrate that the inertia-free approach is as efficient as inertia detection via symmetric indefinite factorizations. We also demonstrate that the inertia-free approach can lead to reductions in solution time because it reduces the amount of convexification needed.

Chiang, Nai-Yuan↗

Performance estimation of resonance-enhanced dual-buoy wave energy converter using coupled time-domain simulation

This article presents the modeling methodology and performance evaluation of the resonance-enhanced dual-buoy WEC (Wave Energy Converter) by HEM (hydrodynamic & electro-magnetic) fully-coupled-dynamics time-domain-simulation program. The numerical results are systematically compared with the authors’ 1/6-scale experiment. With a direct-drive linear generator, the WEC consists of dual floating cylinders and a moon-pool between the cylinders, which can utilize three resonance phenomena from moon-pool dynamics as well as heave motions of inner and outer buoys. The contact and friction between the two buoys observed in the experiment are also properly modeled in the time-domain simulation by the Coulomb-friction model. Moon-pool resonance peaks significantly exaggerated in linear potential theory are empirically adjusted through comparisons with measured values. A systematic comparative study between the simulations and experiments with and without PTO (power-take-off) is conducted, and the relative heave displacements/velocities and power outputs are well matched. Then, parametric studies are carried out with the simulation program to determine optimum generator parameters. The performance with various wave conditions is also assessed.

16 TIDAL AND WAVE POWER↗

Discrete Time McKean–Vlasov Control Problem: A Dynamic Programming Approach

We consider the stochastic optimal control problem of nonlinear mean-field systems in discrete time. We reformulate the problem into a deterministic control problem with marginal distribution as controlled state variable, and prove that dynamic programming principle holds in its general form. We apply our method for solving explicitly the mean-variance portfolio selection and the multivariate linear-quadratic McKean–Vlasov control problem.

97 MATHEMATICS AND COMPUTING↗

Development of a clostridia-based cell-free system for prototyping genetic parts and metabolic pathways

Gas fermentation by autotrophic bacteria, such as clostridia, offers a sustainable path to numerous bioproducts from a range of local, highly abundant, waste and low-cost feedstocks, such as industrial flue gases or syngas generated from biomass or municipal waste. Unfortunately, designing and engineering clostridia remains laborious and slow. The ability to prototype individual genetic part function, gene expression patterns, and biosynthetic pathway performance in vitro before implementing designs in cells could help address these bottlenecks by speeding up design. Unfortunately, a high-yielding cell-free gene expression (CFE) system from clostridia has yet to be developed. Here in this paper, we report the development and optimization of a high-yielding (236 ± 24 μg/mL) batch CFE platform from the industrially relevant anaerobe, Clostridium autoethanogenum. A key feature of the platform is that both circular and linear DNA templates can be applied directly to the CFE reaction to program protein synthesis. We demonstrate the ability to prototype gene expression, and quantitatively map aerobic cell-free metabolism in lysates from this system. We anticipate that the C. autoethanogenum CFE platform will not only expand the protein synthesis toolkit for synthetic biology, but also serve as a platform in expediting the screening and prototyping of gene regulatory elements in non-model, industrially relevant microbes.

59 BASIC BIOLOGICAL SCIENCES↗

Experimental and numerical analysis of pre-compressed masonry walls in two-way-bending with second order effects

The buckling behavior of slender unreinforced masonry (URM) walls subjected to axial compression and out-of-plane lateral loads is investigated through a combined experimental and numerical homogenizedapproach. After a preliminary analysis performed on a unit cell meshed by means of elastic FEs and non-linear interfaces, macroscopic moment-curvature diagrams so obtained are implemented at a structural level, discretizing masonry by means of rigid triangular elements and non-linear interfaces. The non-linear incremental response of the structure is accounted for a specific quadratic programming routine. In parallel, a wide experimental campaign is conducted on walls in two way bending, with the double aim of both validating the numerical model and investigating the behavior of walls that may not be reduced to simple cantilevers or simply supported beams. Panels investigated are dry-joint in scale square walls simply supported at the base and on a vertical edge, exhibiting the classical Rondelet’s mechanism. The results obtained are compared with those provided by the numerical model.

42 ENGINEERING↗

The linear and nonlinear response of infinite periodic systems to static and/or dynamic electric fields. Implementation in CRYSTAL code

An implementation of the vector potential approach (VPA) for treating the response of infinite periodic systems to static and dynamic electric fields has been initiated within the CRYSTAL code. The VPA method is based on the solution of a time-dependent Hartree-Fock or Kohn-Sham equation for the crystal orbitals wherein the usual scalar potential, that describes interaction with the field, is replaced by the vector potential. This equation may be solved either by perturbation theory or by finite field methods. With some modification all the computational procedures of molecular ab initio quantum chemistry can be adapted for periodic systems. Accessible properties include the linear and nonlinear responses of both the nuclei and the electrons. The programming of static field pure electronic (hyper)polarizabilities has been successfully tested. Dynamic electronic (hyper)polarizabilities, as well as infrared and Raman intensities, are in progress while the addition of finite fields for calculation of vibrational (hyper)polarizabilities, through nuclear relaxation procedures, will begin shortly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

qsp4pde v1.0

The software program is a collection of Python implementations of quantum circuits for solving linear partial differential equations using quantum signal processing. The circuits are implemented using the qiskit SDK.

Kim, Hyeongjin [Lawrence Berkeley National Laborat↗

Analysis of Installed Measures and Energy Savings for Single-Family Residential Better Buildings Projects

This report presents an analysis of data for residential single-family projects reported by 37 organizations that were awarded federal financial assistance (cooperative agreements or grants) by the U.S. Department of Energy’s Better Buildings Neighborhood Program.1 The report characterizes the energy-efficiency measures installed for single-family residential projects and analyzes energy savings and savings prediction accuracy for measures installed in a subset of those projects.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NCSP Related Nuclear Data Research at RPI [Slides]

This lecture is on NCSP related nuclear data research at the Rensselaer Polytechnic Institute (RPI). This presentation includes an overview of fiscal year 2022 activity. Additionally, it has a talk on the RPI Nuclear Data (ND) Group Research, an update on Linear Accelerators (LINAC) refurbishment, and the Enhanced Thermal Target & Cold Moderator (ETTC) Target System.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Using Whisper-1.1 to Guide Improvements to Nuclear Data Evaluations

Along with the latest release of the general purpose Monte Carlo radiation transport code, MCNP6.2, came the first officially released version of Whisper, version 1.1, a statistical analysis code used for nuclear criticality safety (NCS) validation. The Whisper software package is primarily used to calculate the upper subcritical limit (USL) of a particular application case along with components of the margin of subcriticality (MOS). In the calculation of the MOS due to nuclear data uncertainties, Whisper uses the generalized linear-least-squares (GLLS) approach to obtain adjusted nuclear data uncertainties which are then propagated to estimate the nuclear data uncertainties in the k eff of the application. Due to the fact that Whisper contains the simulated and measured k eff values for over 1100 benchmarks, nuclear data sensitivity profiles for each benchmark along with the nuclear data covariance matrices, all of which are used in the GLLS method, it is possible to easily utilize these data available to help guide improvements to nuclear data evaluations. In the present work, the Whisper data and the GLLS method are used together to suggest an adjusted nuclear data cross section evaluation specifically for 233 U such that the agreement between the calculated and experimental k eff measurements of the relevant benchmarks are improved. With this adjusted nuclear data available, feedback to the theoretical and experimental nuclear data community is provided.

97 MATHEMATICS AND COMPUTING↗

Low Precision and Efficient Programming Languages for Sustainable AI: Final Report for the Summer Project of 2024

This document contains all relevant material generated during the authors' summer internship at NREL in 2024. This report shows how to improve energy efficiency of a few code samples by using low-precision data types combined with mixed-precision algorithms. The main applications considered here are (i) linear system solvers using mixed precision, and (ii) neural networks using mixed precision. This report also discusses how programming languages affect energy consumption of algorithms, energy metrics for a code and tools, and the available current software and hardware infrastructure.

97 MATHEMATICS AND COMPUTING↗

SU-E-T-479: IMRT Plan Recalculation in Patient Based On Dynalog Data and the Effect of a Single Failing MLC Motor

Purpose: Using Linac dynamic logs (Dynalogs) we evaluate the impact of a single failing MLC motor on the deliverability of an IMRT plan by assessing the recalculated dose volume histograms (DVHs) taking the delivered MLC positions and beam hold-offs into consideration. Methods: This is a retrospective study based on a deteriorating MLC motor (leaf 36B) which was observed to be failing via Dynalog analysis. To investigate further, Eclipse-importable MLC files were generated from Dynalogs to recalculate the actual delivered dose and to assess the clinical impact through DVHs. All deliveries were performed on a Varian 21EX linear accelerator equipped with Millennium-120 MLC. The analysis of Dynalog files and subsequent conversion to Eclipse-importable MLC files were all performed by in-house programming in Python. Effects on plan DVH are presented in the following section on a particular brain-IMRT plan which was delivered with a failing MLC motor which was then replaced. Results: Global max dose increased by 13.5%, max dose to the brainstem PRV increased by 8.2%, max dose to the optic chiasm increased by 7.6%, max dose to optic nerve increased by 8.8% and the mean dose to the PTV increased by 7.9% when comparing the original plan to the fraction with the failing MLC motor. The reason the dose increased was due to the failure being on the B-bank which is the lagging side on a sliding window delivery, therefore any failures on this side will cause an over-irradiation as the B-bank leaves struggles to keep the window from growing. Conclusion: Our findings suggest that a single failing MLC motor may jeopardize the entire delivery. This may be due to the bad MLC motor drawing too much current causing all MLCs on the same bank to underperform. This hypothesis will be investigated in a future study.

60 APPLIED LIFE SCIENCES↗

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗

SU-E-T-465: Implementation of An Automated Collision Detection Program Using Open Source Software for the Pinnacle Treatment Planning System

Purpose: Potential collisions between the gantry head and the patient or table assembly are difficult to detect in most treatment planning systems. We have developed and implemented a novel software package for the representation of potential gantry collisions with the couch assembly at the time of treatment planning. Methods: Physical dimensions of the Varian Edge linear accelerator treatment head were measured and reproduced using the Visual Python display package. A script was developed for the Pinnacle treatment planning system to generate a file with the relevant couch, gantry, and isocenter positions for each beam in a planning trial. A python program was developed to parse the information from the TPS and produce a representative model of the couch/gantry system. Using the model and the Visual Python libraries, a rendering window is generated for each beam that allows the planner to evaluate the possibility of a collision. Results: Comparison against heuristic methods and direct verification on the machine validated the collision model generated by the software. Encounters of <1 cm between the gantry treatment head and table were visualized as collisions in our virtual model. Visual windows were created depicting the angle of collision for each beam, including the anticipated table coordinates. Visual rendering of a 6 arc trial with multiple couch positions was completed in under 1 minute, with network bandwidth being the primary bottleneck. Conclusion: The developed software allows for quick examination of possible collisions during the treatment planning process and helps to prevent major collisions prior to plan approval. The software can easily be implemented on future planning systems due to the versatility and platform independence of the Python programming language. Further integration of the software with the treatment planning system will allow the possibility of patient-gantry collision detection for a range of treatment machines.

60 APPLIED LIFE SCIENCES↗

A deterministic solution of the first order linear Boltzmann transport equation in the presence of external magnetic fields

Purpose: Accurate radiotherapy dose calculation algorithms are essential to any successful radiotherapy program, considering the high level of dose conformity and modulation in many of today’s treatment plans. As technology continues to progress, such as is the case with novel MRI-guided radiotherapy systems, the necessity for dose calculation algorithms to accurately predict delivered dose in increasingly challenging scenarios is vital. To this end, a novel deterministic solution has been developed to the first order linear Boltzmann transport equation which accurately calculates x-ray based radiotherapy doses in the presence of magnetic fields. Methods: The deterministic formalism discussed here with the inclusion of magnetic fields is outlined mathematically using a discrete ordinates angular discretization in an attempt to leverage existing deterministic codes. It is compared against the EGSnrc Monte Carlo code, utilizing the emf-macros addition which calculates the effects of electromagnetic fields. This comparison is performed in an inhomogeneous phantom that was designed to present a challenging calculation for deterministic calculations in 0, 0.6, and 3 T magnetic fields oriented parallel and perpendicular to the radiation beam. The accuracy of the formalism discussed here against Monte Carlo was evaluated with a gamma comparison using a standard 2%/2 mm and a more stringent 1%/1 mm criterion for a standard reference 10 × 10 cm{sup 2} field as well as a smaller 2 × 2 cm{sup 2} field. Results: Greater than 99.8% (94.8%) of all points analyzed passed a 2%/2 mm (1%/1 mm) gamma criterion for all magnetic field strengths and orientations investigated. All dosimetric changes resulting from the inclusion of magnetic fields were accurately calculated using the deterministic formalism. However, despite the algorithm’s high degree of accuracy, it is noticed that this formalism was not unconditionally stable using a discrete ordinate angular discretization. Conclusions: The feasibility of including magnetic field effects in a deterministic solution to the first order linear Boltzmann transport equation is shown. The results show a high degree of accuracy when compared against Monte Carlo calculations in all magnetic field strengths and orientations tested.

60 APPLIED LIFE SCIENCES↗