Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear model”

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.

777 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-T-295: Simultaneous Beam Sampling and Aperture Shape Optimization for Station Parameter Optimized Radiation Therapy (SPORT)

Purpose: Station Parameter Optimized Radiation Therapy (SPORT) was recently proposed to fully utilize the technical capability of emerging digital LINACs, in which the station parameters of a delivery system, (such as aperture shape and weight, couch position/angle, gantry/collimator angle) are optimized altogether. SPORT promises to deliver unprecedented radiation dose distributions efficiently, yet there does not exist any optimization algorithm to implement it. The purpose of this work is to propose an optimization algorithm to simultaneously optimize the beam sampling and aperture shapes. Methods: We build a mathematical model whose variables are beam angles (including non-coplanar and/or even nonisocentric beams) and aperture shapes. To solve the resulting large scale optimization problem, we devise an exact, convergent and fast optimization algorithm by integrating three advanced optimization techniques named column generation, gradient method, and pattern search. Column generation is used to find a good set of aperture shapes as an initial solution by adding apertures sequentially. Then we apply the gradient method to iteratively improve the current solution by reshaping the aperture shapes and updating the beam angles toward the gradient. Algorithm continues by pattern search method to explore the part of the search space that cannot be reached by the gradient method. Results: The proposed technique is applied to a series of patient cases and significantly improves the plan quality. In a head-and-neck case, for example, the left parotid gland mean-dose, brainstem max-dose, spinal cord max-dose, and mandible mean-dose are reduced by 10%, 7%, 24% and 12% respectively, compared to the conventional VMAT plan while maintaining the same PTV coverage. Conclusion: Combined use of column generation, gradient search and pattern search algorithms provide an effective way to optimize simultaneously the large collection of station parameters and significantly improves quality of resultant treatment plans as compared with conventional VMAT or IMRT treatments.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Physics and technology considerations for the deuterium–tritium fuel cycle and conditions for tritium fuel self sufficiency

The tritium aspects of the DT fuel cycle embody some of the most challenging feasibility and attractiveness issues in the development of fusion systems. The review and analyses in this paper provide important information to understand and quantify these challenges and to define the phase space of plasma physics and fusion technology parameters and features that must guide a serious R&D in the world fusion program. We focus in particular on components, issues and R&D necessary to satisfy three 'principal requirements': (1) achieving tritium self-sufficiency within the fusion system, (2) providing a tritium inventory for the initial start-up of a fusion facility, and (3) managing the safety and biological hazards of tritium. A primary conclusion is that the physics and technology state-of-the-art will not enable DEMO and future power plants to satisfy these principal requirements. We quantify goals and define specific areas and ideas for physics and technology R&D to meet these requirements. A powerful fuel cycle dynamics model was developed to calculate time-dependent tritium inventories and flow rates in all parts and components of the fuel cycle for different ranges of parameters and physics and technology conditions. Dynamics modeling analyses show that the key parameters affecting tritium inventories, tritium start-up inventory, and tritium self-sufficiency are the tritium burn fraction in the plasma (f b ), fueling efficiency (η f ), processing time of plasma exhaust in the inner fuel cycle (t p ), reactor availability factor (AF), reserve time (tr) which determines the reserve tritium inventory needed in the storage system in order to keep the plant operational for time t r in case of any malfunction of any part of the tritium processing system, and the doubling time (t d ). Results show that η f f b > 2% and processing time of 1–4 h are required to achieve tritium self-sufficiency with reasonable confidence. For η f f b = 2% and processing time of 4 h, the tritium start-up inventory required for a 3 GW fusion reactor is ~11 kg, while it is <5 kg if η f f b = 5% and the processing time is 1 h. To achieve these stringent requirements, a serious R&D program in physics and technology is necessary. The EU-DEMO direct internal recycling concept that carries fuel directly from the plasma exhaust gas to the fueling systems without going through the isotope separation system reduces the overall processing time and tritium inventories and has positive effects on the required tritium breeding ratio (TBR R ). A significant finding is the strong dependence of tritium self-sufficiency on the reactor availability factor. Simulations show that tritium self-sufficiency is: impossible if AF < 10% for any η f f b , possible if AF > 30% and 1% ≤ η f f b ≤ 2%, and achievable with reasonable confidence if AF > 50% and η f f b > 2%. These results are of particular concern in light of the low availability factor predicted for the near-term plasma-based experimental facilities (e.g. FNSF, VNS, CTF), and can have repercussions on tritium economy in DEMO reactors as well, unless significant advancements in RAMI are made. There is a linear dependency between the tritium start-up inventory and the fusion power. The required tritium start-up inventory for a fusion facility of 100 MW fusion power is as small as 1 kg. Since fusion power plants will have large powers for better economics, it is important to maintain a 'reserve' tritium inventory in the tritium storage system to continue to fuel the plasma and avoid plant shutdown in case of malfunctions of some parts of the tritium processing lines. But our results show that a reserve time as short as 24 h leads to unacceptable reserve and start-up inventory requirements. Therefore, high reliability and fast maintainability of all components in the fuel cycle are necessary in order to avoid the need for storing reserve tritium inventory sufficient for continued fusion facility operation for more than a few hours. The physics aspects of plasma fueling, tritium burn fraction, and particle and power exhaust are highly interrelated and complex, and predictions for DEMO and power reactors are highly uncertain because of lack of experiments with burning plasma. Fueling by pellet injection on the high field side of tokamak has evolved to be the preferred method to fuel a burning plasma. Extrapolation from the DIII-D penetration scaling shows fueling efficiency expected in DEMO to be <25%, but such extrapolations are highly uncertain. The fueling efficiency of gas in a reactor relevant regime is expected to be extremely poor and not very useful for getting tritium into the core plasma efficiently. Gas fueling will nonetheless be useful for feedback control of the divertor operating parameters. Extensive modeling has been carried out to predict burn fraction, fueling requirements, and fueling efficiency for ITER, DEMO, and beyond. The fueling rate required to operate Q = 10 ITER plasmas in order to provide the required core fueling, helium exhaust and radiative divertor plasma conditions for acceptable divertor power loads was calculated. If this fueling is performed with a 50–50 DT mix, the tritium burn fraction in ITER would be ~0.36%, which is too low to satisfy the self-sufficiency conditions derived from the dynamics modeling for fusion reactors. Extrapolation to DEMO using this approach would also yield similarly low burn fraction. Extensive analysis presented shows that specific features of edge neutral dynamics in ITER and fusion reactors, which are different from present experiments, open possibilities for optimization of tritium fueling and thus to improve the burn fraction. Using only tritium in pellet fueling of the plasma core, and only deuterium for edge density, divertor power load and ELM control results in significant increase of the burn fraction to 1.8–3.6%. These estimates are performed with physics models whose results cannot be fully validated for ITER and DEMO plasma conditions since these cannot be achieved in present tokamak experiments. Thus, several uncertainties remain regarding particle transport and scenario requirements in ITER and DEMO. The safety standard requirements for protection of the public and release guidelines for tritium have been reviewed. General safety approaches including minimizing tritium inventories, reducing tritium permeation through materials, and decontaminating material for waste disposal have been suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗