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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments

The computational study of chemical reactions in complex, wet environments is critical for applications in many fields. It is often essential to study chemical reactions in the presence of applied electrochemical potentials, taking into account the non-trivial electrostatic screening coming from the solvent and the electrolytes. As a consequence, the electrostatic potential has to be found by solving the generalized Poisson and the Poisson-Boltzmann equations for neutral and ionic solutions, respectively. In the present work, solvers for both problems have been developed. A preconditioned conjugate gradient method has been implemented for the solution of the generalized Poisson equation and the linear regime of the Poisson-Boltzmann, allowing to solve iteratively the minimization problem with some ten iterations of the ordinary Poisson equation solver. In addition, a self-consistent procedure enables us to solve the non-linear Poisson-Boltzmann problem. Both solvers exhibit very high accuracy and parallel efficiency and allow for the treatment of periodic, free, and slab boundary conditions. The solver has been integrated into the BigDFT and Quantum-ESPRESSO electronic-structure packages and will be released as an independent program, suitable for integration in other codes.

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

SU-E-T-563: A Fast and Quantative Picket-Fence Test of a Submillimeter Patient Positioning System for Stereotactic Radiosurgery

Purpose: Picket-fence test is a qualitative TG142-recommended quality assurance (QA) test for multileaf collimators. In study, we adopted the same concept and developed a fast but quantatitive QA test for an automatic patient positioning system that requires submilleter accuracy for a radiosurgical treatment. Methods: A piece of radiochromic film was first placed inside a spherical solid water phantom and then irradiated with a sequenence of linearly placed shots of same collimator size (e.g. 4-mm) via the Leskell Gamma Knife Perfexion system (PFX). The shots were positioned with either equal or non-equal gaps of approximately 4 mm to 8 mm depending on the location of the region of interest. A pattern recognization program was developed and then applied to measure the gap spacing between two adjacent shots. The measured distance was then compared with the initial preset values for the test. Results: By introducing variable systematic and random shifts of 0.1 mm to 0.5 mm to the shot sequence, the maximum gap variation from the described test was found to be 0.35 mm or less. On average the positioning uncertainty for the PFX delivery system was found to be 0.1±0.2 mm. No significant difference in the positioning uncertainty was noted for the centrally aligned shot sequence locations versus the peripherally aligned shot sequence locations. Conclusion: A new quantitative picket-fence type test was developed and demonstrated for routine QA of the submillimeter PFX patient positioning sytem. This test also enables independent verification of any patient-specific shot positioning for a critical treatment such as a tumor in brainstem. Dr Ma is currently on the board of international society of stereotactic Radiosurgery.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Tailoring Carbide Dispersed Steels: A Path to Increased Strength and Hydrogen Tolerance

The use of transition metal carbides is reported for use as a hydrogen trapping mechanism for ferritic and austenitic steel materials. The program combined computational modeling and simulations to guide experiments towards candidate metal carbide traps, both for interfacial and interior trapping. It was found that interfacial trapping is less effective than interior trapping, with the group IVB transition metal carbides being the most effect internal traps with a loss of carbon. The sub-stoichiometric rocksalt structure accommodate the hydrogen atoms in its octahedral interstices. Using percolation theory, carbon loss of approximately 25% or more was sufficient to ensure an interconnected network of vacancies for such trapping from the surface to the internal sites within the carbide. Using this as a guide, the program developed a means to provide a uniform dispersion of ZrC nanoparticles with either Fe or 304L micron-scale powders which was then consolidated by direct current sintering. Electrolytic hydrogen diffusivity studies confirmed the reduction of hydrogen diffusion in the matrix with increasing ZrC content, which was a linear response over the sample range studied (0.01 to 1.0 wt.%). The consolidated material was micro-tensile tested in either a non-hydrogen or hydrogen charge condition and compared to a control with no carbides. Additions up to 0.05 wt.% ZrC increased the yield strength with no loss in ductility in either the non-hydrogen or hydrogen tested condition. ZrC concentrations above this amount further increased the yield strength at the expense of ductility. While these samples had a lower absolute ductility value prior to failure, the relative change in ductility between the non-hydrogen and hydrogen charge states was less for the carbides than that of the control. Metal-rich ZrC nanoparticles were fabricated through a conformal coating process yielding ZrC0.66 particles that were then incorporated into a metal matrix. Notch fatigue testing in a hydrogen environment was conducted where the number of cycles to failure was found to be less in the control than that of the carbide addition. However, the spread in experimental data and the number of samples tested limits a conclusive outcome based on defects noticed in the gauge section of all the powder processed samples. The collective outcomes of this report provide further insight into the mechanisms by which carbides act as hydrogen traps; a means to process such carbides through powder metallurgy; and their associated mechanical performance in either a non-hydrogen or hydrogen-charged condition.

08 HYDROGEN↗
DOE OSTI2025

Accelerators for the Future: R&D at the Fermilab FAST Facility

High energy physics in the U.S. has ambitious plans requiring new and improved accelerator technologies, including an upgraded complex at Fermilab for DUNE, next-generation light sources, and the potential of a future collider to be built on U.S. soil. To address these requirements, Fermilab operates the Fermilab Accelerator Science and Technology (FAST) facility, dedicated to accelerator R&D. FAST includes an electron gun and superconducting RF linac (up to 300 MeV), a storage ring, and an upcoming proton source and injector line (up to 2.5 MeV). In addition to the future of Fermilab accelerators, the broad physics program includes general R&D with potential impact across the DOE science program. This colloquium will provide a basic introduction to accelerator technologies and describe the principles of proton and electron accelerators, including the challenges associated with next-generation operations. We will discuss the exciting R&D ongoing at FAST to address the required technological advancements, focusing on Non-Linear Integrable Optics (NIO) for improved beam intensity and Optical Stochastic Cooling (OSC) for improved beam quality, and touching on a wide range of additional ongoing research. I intend to make the content interesting and accessible to those who have never taken any formal courses in accelerator physics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Structural stability of thin overhanging walls during material extrusion additive manufacturing of thermoset-based ink

Recent developments have enabled material extrusion additive manufacturing of thermoset-based composite inks on the large scale. In addition, printing out-of-plane components is of broad interest to the polymer material extrusion community. Here, we address some of the challenges associated with both large-scale and out-of-plane thermoset material extrusion additive manufacturing by studying the height at which thin overhanging walls collapse. Walls at a range of overhang angles were printed until they collapsed. An optical camera captured the profile of each wall throughout the print, allowing the collapse height to be identified and the geometric fidelity to the programmed angle to be evaluated. Using previously measured rheological properties, predictive models were generated to approximate the collapse height and profile of the deflected walls. First, an analytical model was created to predict the height at which the walls would yield. The analytical model assumes the walls exhibit a perfectly linear profile; however, experiments proved this assumption to be false. Therefore, a finite element simulation was developed to account for the elastic deflection that occurs during printing. The finite element simulation predicts both the yield height and the deflected profile after the deposition of each layer. For the properties of the thermoset ink used here, the yield height predicted by the analytical model and finite element simulation are virtually identical. These predictions match experimental data reasonably well, but minor errors are observed. Accounting for the fully plastic moment appears to explain the small mismatch between experimental data and predictions. Additionally, the finite element simulation provides an excellent prediction of the deflected profile before the wall begins to collapse. Finally, by demonstrating that the collapse height and deflected profile of thin overhanging walls can be predicted, this work illustrates how the soft viscoelastic properties of thermoset-based composite inks limit the scale of a key feature required to print some nonplanar components. It also provides a basis to tailor in-process curing systems to suppress deflection and collapse of thin overhanging walls.

36 MATERIALS SCIENCE↗
DOE OSTI2017

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM development aims for advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability for reactor transient analyses. To facilitate the code development, SAM utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage modern advanced software environments and numerical methods. SAM focuses on modeling advanced reactor concepts such as SFRs (sodium fast reactors), LFRs (lead-cooled fast reactors), and FHRs (fluoride-salt-cooled high temperature reactors) or MSRs (molten salt reactors). These advanced concepts are distinguished from light-water reactors in their use of single-phase, low-pressure, high-temperature, and low Prandtl number (sodium and lead) coolants. As a new code development, the initial effort has been focused on modeling and simulation capabilities of heat transfer and single-phase fluid dynamics responses in Sodium-cooled Fast Reactor (SFR) systems. The system-level simulation capabilities of fluid flow and heat transfer in general engineering systems and typical SFRs have been verified and validated. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As the code is still under ongoing development, this SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

An Algorithm to Retrieve Total Precipitable Water Vapor in the Atmosphere from FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) Data

The atmosphere has substantial effects on optical remote sensing imagery of the Earth’s surface from space. These effects come through the functioning of atmospheric particles on the radiometric transfer from the Earth’s surface through the atmosphere to the sensor in space. Precipitable water vapor (PWV), CO 2 , ozone, and aerosol in the atmosphere are very important among the particles through their functioning. This study presented an algorithm to retrieve total PWV from the Chinese second-generation polar-orbiting meteorological satellite FengYun 3D Medium Resolution Spectral Imager 2 (FY-3D MERSI-2) data, which have three near-infrared (NIR) water vapor absorbing channels, i.e., channel 16, 17, and 18. The algorithm was improved from the radiance ratio technique initially developed for Moderate-Resolution Imaging Spectroradiometer (MODIS) data. MODTRAN 5 was used to simulate the process of radiant transfer from the ground surfaces to the sensor at various atmospheric conditions for estimation of the coefficients of ratio technique, which was achieved through statistical regression analysis between the simulated radiance and transmittance values for FY-3D MERSI-2 NIR channels. The algorithm was then constructed as a linear combination of the three-water vapor absorbing channels of FY-3D MERSI-2. Measurements from two ground-based reference datasets were used to validate the algorithm: the sun photometer measurements of Aerosol Robotic Network (AERONET) and the microwave radiometer measurements of Energy’s Atmospheric Radiation Measurement Program (ARMP). The validation results showed that the algorithm performs very well when compared with the ground-based reference datasets. The estimated PWV values come with root mean square error (RMSE) of 0.28 g/cm 2 for the ARMP and 0.26 g/cm 2 for the AERONET datasets, with bias of 0.072 g/cm 2 and 0.096 g/cm 2 for the two reference datasets, respectively. The accuracy of the proposed algorithm revealed a better consistency with ground-based reference datasets. Thus, the proposed algorithm could be used as an alternative to retrieve PWV from FY-3D MERSI-2 data for various remote sensing applications such as agricultural monitoring, climate change, hydrologic cycle, and so on at various regional and global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Plasma-Material-Interaction Research Using PISCES Linear Plasma Devices

In this paper, recent plasma-material-interaction (PMI) studies of relevance to fusion PMI science are discussed. The results are a mix of science experiments and in situ diagnostic technology development. Scientific issues addressed include the nature of He and W atom migration in the evolution of W fuzz surfaces, understanding the nature of D retention in damaged W surfaces via modeling with the tritium migration and analysis program, and an examination of synergistic effects associated with D trapping on composite Be codeposit and W substrate surfaces. In situ diagnostic developments discussed include laser-induced desorption spectroscopy for the measurement of bulk hydrogen retention and laser-induced breakdown spectroscopy for the measurement of near-surface composition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Data-driven optimization of mixed-integer bi-level multi-follower integrated planning and scheduling problems under demand uncertainty

The coordination of interconnected elements across the different layers of the supply chain is essential for all industrial processes and the key to optimal decision-making. Yet, the modeling and optimization of such interdependent systems are still burdensome. Here we address the simultaneous modeling and optimization of medium-term planning and short-term scheduling problems under demand uncertainty using mixed-integer bi-level multi-follower programming and data-driven optimization. Bi-level multi-follower programs model the natural hierarchy between different layers of supply chain management holistically, while scenario analysis and data-driven optimization allow us to retrieve the guaranteed feasible solutions of the integrated formulation under various demand considerations. We address the data-driven optimization of this challenging class of problems using the DOMINO framework, which was initially developed to solve single-leader single-follower bi-level optimization problems to guaranteed feasibility. This framework is extended to solve single-leader multi-follower stochastic formulations and its performance is characterized by well-known single and multi-product process scheduling case studies. Through our data-driven algorithmic approach, we present guaranteed feasible solutions to linear and nonlinear mixed-integer bi-level formulations of simultaneous planning and scheduling problems and further characterize the effects of the scheduling level complexity on the solution performance, which spans over several hundred continuous and binary variables, and thousands of constraints.

42 ENGINEERING↗
DOE OSTI2023

Semidefinite programming algorithm for the quantum mechanical bootstrap

Here, we present a semidefinite programming algorithm to find eigenvalues of Schrödinger operators within the bootstrap approach to quantum mechanics. The bootstrap approach involves two ingredients: a nonlinear set of constraints on the variables (expectation values of operators in an energy eigenstate), plus positivity constraints (unitarity) that need to be satisfied. By fixing the energy we linearize all the constraints and show that the feasibility problem can be presented as an optimization problem for the variables that are not fixed by the constraints and one additional slack variable that measures the failure of positivity. To illustrate the method we are able to obtain high-precision, sharp bounds on eigenenergies for arbitrary confining polynomial potentials in one dimension.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Evaluated 182,183,184,186 W Neutron Cross Sections and Covariances in the Resolved Resonance Region

Oak Ridge National Laboratory (ORNL) has recently completed the resonance parameter evaluation of four tungsten isotopes, i.e., 182,183,184,186 W, in the neutron energy range of thermal up to several keV. This nuclear data work was performed with support from the US Nuclear Criticality Safety Program (NCSP) in an effort to provide improved tungsten cross section and covariance data for criticality safety analyses. The evaluation methodology uses the Reich-Moore approximation of the R-matrix formalism of the code SAMMY to fit high-resolution measurements performed in 2010 and 2012 at the Geel linear accelerator facility (GELINA), as well as other experimental data sets on natural tungsten available in the EXFOR library. In the analyzed energy range, this work nearly doubles the resolved resonance region (RRR) present in the latest US nuclear data library ENDF/B-VII.1. In view of the interest in tungsten for distinct types of nuclear applications and the relatively homogeneous distribution of the isotopic tungsten—namely, 182 W(26.5%), 183 W(14.31%), 184 W(30.64%), and 186 W(28.43%) - the completion of these four evaluations represents a significant contribution to the improvement of the ENDF library. This paper presents an overview of the evaluated resonance parameters and related covariances for total and capture cross sections on the four tungsten isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Evaluated {sup 182,183,184,186}W neutron cross sections and covariances in the resolved resonance region - 14617

Oak Ridge National Laboratory (ORNL) has recently completed the resonance parameter evaluation of four tungsten isotopes, i.e., {sup 182,183,184,186}W, in the neutron energy range of thermal up to several keV. This nuclear data work was performed with support from the US Nuclear Criticality Safety Program (NCSP) in an effort to provide improved tungsten cross section and covariance data for criticality safety analyses. The evaluation methodology uses the Reich-Moore approximation of the R-matrix formalism of the code SAMMY to fit high-resolution measurements performed in 2010 and 2012 at the Geel linear accelerator facility (GELINA), as well as other experimental data sets on natural tungsten available in the EXFOR library. In the analyzed energy range, this work nearly doubles the resolved resonance region (RRR) present in the latest US nuclear data library ENDF/B-VII.1. In view of the interest in tungsten for distinct types of nuclear applications and the relatively homogeneous distribution of the isotopic tungsten - namely, {sup 182}W(26.5%), {sup 183}W(14.31%), {sup 184}W(30.64%), and {sup 186}W(28.43%) - the completion of these four evaluations represents a significant contribution to the improvement of the ENDF library. This paper presents an overview of the evaluated resonance parameters and related covariances for total and capture cross sections on the four tungsten isotopes. (authors)

42 ENGINEERING↗
DOE OSTI2015

TU-EF-304-09: Quantifying the Biological Effects of Therapeutic Protons by LET Spectrum Analysis

Purpose: To correlate in vitro cell kill with linear energy transfer (LET) spectra using Monte Carlo simulations and knowledge obtained from previous high-throughput in vitro proton relative biological effectiveness (RBE) measurements. Methods: The Monte Carlo simulation toolkit Geant4 was used to design the experimental setups and perform the dose, dose-averaged LET, and LET spectra calculations. The clonogenic assay was performed using the H460 lung cancer cell line in standard 6-well plates. Using two different experimental setups, the same dose and dose-averaged LET (12.6 keV/µm) was delivered to the cell layer; however, each respective energy or LET spectrum was different. We quantified the dose contributions from high-LET (≥10 keV/µm, threshold determined by previous RBE measurements) events in the LET spectra separately for these two setups as 39% and 53%. 8 dose levels with 1 Gy increments were delivered. The photon reference irradiation was performed using 6 MV x-rays from a LINAC. Results: The survival curves showed that both proton irradiations demonstrated an increased RBE compared to the reference photon irradiation. Within the proton-irradiated cells, the setup with 53% dose contribution from high-LET events exhibited the higher biological effectiveness. Conclusion: The experimental results indicate that the dose-averaged LET may not be an appropriate indicator to quantify the biological effects of protons when the LET spectrum is broad enough to contain both low- and high-LET events. Incorporating the LET spectrum distribution into robust intensity-modulated proton therapy optimization planning may provide more accurate biological dose distribution than using the dose-averaged LET. NIH Program Project Grant 2U19CA021239-35.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

A mathematical programming approach for the optimal placement of flame detectors in petrochemical facilities

Flame detectors provide an important layer of protection for personnel in petrochemical plants, but effective placement can be challenging. A mixed-integer nonlinear programming formulation is proposed for optimal placement of flame detectors while considering non-uniform probabilities of detection failure. We show that this approach allows for the placement of fire detectors using a fixed sensor budget and outperforms models that do not account for imperfect detection. In this study, we develop a linear relaxation to the formulation and an efficient solution algorithm that achieves global optimality with reasonable computational effort. We integrate this problem formulation into the Python package, Chama, and demonstrate the effectiveness of this formulation on a small test case and on two real-world case studies using the fire and gas mapping software, Kenexis Effigy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2013

Singular linear-quadratic control problem for systems with linear delay

A singular linear-quadratic optimization problem on the trajectories of non-autonomous linear differential equations with linear delay is considered. The peculiarity of this problem is the fact that this problem has no solution in the class of integrable controls. To ensure the existence of solutions is required to expand the class of controls including controls with impulse components. Dynamical systems with linear delay are used to describe the motion of pantograph from the current collector with electric traction, biology, etc. It should be noted that for practical problems fact singularity criterion of quality is quite commonly occurring, and therefore the study of these problems is surely important. For the problem under discussion optimal programming control contained impulse components at the initial and final moments of time is constructed under certain assumptions on the functional and the right side of the control system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Characterization of unconventional electron fields for the treatment of mycosis fungoides using the total skin irradiation technique

Mycosis fungoides is a cutaneous lymphoma that accounts for 2–3% of all lymphomas. Several clinical studies have demonstrated the effectiveness of TSEBT (Total Skin Electron Beam Therapy) in patients with mycosis fungoides. It is important to develop this technique and make it available to a larger number of patients in Mexico. Because large fields for electron TSEBT are required in order to cover the entire body of the patient, beam characterization at conventional treatment distances is not sufficient and a calibration distance of 500cm or higher is required. Materials and methods: Calibration of radiochromic Gafchromic® EBT2 film (RCF) for electrons was performed in a solid water phantom (Scanditronix Wellhöfer) at a depth of 1.4cm and a Source Axis Distance (SAD) of 100cm. A polynomial fit was applied to the calibration curve, in order to obtain the equation relating dose response with optical density. The spatial distribution is obtained in terms of percentage of the dose, placing 3×3cm samples of RCF on the acrylic screen, which is placed in front of the patient in order to obtain maximum absorbed dose on the skin, covering an area of 200×100cm{sup 2}. The Percentage Depth Dose (PDD) curve was obtained placing RCF samples at depths of 0, 1, 1.2, 1.4, 1.5, 2, 3, 4, 5, 6, 7, 8 and 9cm in the solid water phantom, irradiated with an ELEKTA SINERGY Linear Accelerator electron beam, with an energy of 6 MeV, at a Source Skin Distance (SSD) of 500cm, with 1000MU = 100Gy, with a cone of 40×40cm and gantry angle of 90°. The RCFs were scanned on a flatbed scanner (EPSON EXPRESSION 10000 XL) and the images were processed with the ImageJ program using a region of interest (ROI) of 1×1cm{sup 2}. Results: The relative spatial dose distribution and the percentage depth dose for a SSD of 500±0.5cm, over an area of 200×100cm{sup 2} was obtained, resulting to an effective maximum dose depth (Z{sub ref}) for electrons of 1.4±0.05cm. Using the same experimental data, horizontal and vertical beam profiles were also graphed, showing a horizontal symmetry of ±035%, horizontal flatness of ±3.62%, vertical symmetry of ±2.1% and vertical flatness of ±14.2%. Conclusions: The electron beam was characterized and the data obtained were useful to determine the spatial dose distribution to a SSD of 500±0.5cm, in an area of 200×100cm{sup 2}. Dose profiles were obtained both horizontally and vertically, thus allowing to assess electron beam symmetry and flatness. PDD analysis up to a depth of 9±0.05cm, has made possible to establish the depth of electron penetration, assuring an only skin irradiation treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

PETSc/TAO Users Manual V.3.21

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations (PDEs) and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for implementing large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and optimizers that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started ). The library is organized hierarchically, enabling users to employ the abstraction level most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users.

97 MATHEMATICS AND COMPUTING↗