SEARCH · Engineering Papers

Results for “parameterization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Predicting X-ray diffuse scattering from translation–libration–screw structural ensembles

A method of simulating X-ray diffuse scattering from multi-model PDB files is presented. Despite similar agreement with Bragg data, different translation–libration–screw refinement strategies produce unique diffuse intensity patterns. Identifying the intramolecular motions of proteins and nucleic acids is a major challenge in macromolecular X-ray crystallography. Because Bragg diffraction describes the average positional distribution of crystalline atoms with imperfect precision, the resulting electron density can be compatible with multiple models of motion. Diffuse X-ray scattering can reduce this degeneracy by reporting on correlated atomic displacements. Although recent technological advances are increasing the potential to accurately measure diffuse scattering, computational modeling and validation tools are still needed to quantify the agreement between experimental data and different parameterizations of crystalline disorder. A new tool, phenix.diffuse, addresses this need by employing Guinier’s equation to calculate diffuse scattering from Protein Data Bank (PDB)-formatted structural ensembles. As an example case, phenix.diffuse is applied to translation–libration–screw (TLS) refinement, which models rigid-body displacement for segments of the macromolecule. To enable the calculation of diffuse scattering from TLS-refined structures, phenix.tls-as-xyz builds multi-model PDB files that sample the underlying T, L and S tensors. In the glycerophosphodiesterase GpdQ, alternative TLS-group partitioning and different motional correlations between groups yield markedly dissimilar diffuse scattering maps with distinct implications for molecular mechanism and allostery. These methods demonstrate how, in principle, X-ray diffuse scattering could extend macromolecular structural refinement, validation and analysis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

A Rietveld refinement method for angular- and wavelength-dispersive neutron time-of-flight powder diffraction data

This paper introduces a two-dimensional extension of the well established Rietveld refinement method for modeling neutron time-of-flight powder diffraction data. The novel approach takes into account the variation of two parameters, diffraction angle 2θ and wavelength λ, to optimally adapt to the varying resolution function in diffraction experiments. By doing so, the refinement against angular- and wavelength-dispersive data gets rid of common data-reduction steps and also avoids the loss of high-resolution information typically introduced by integration. In a case study using a numerically simulated diffraction pattern of Rh 0.81 Fe 3.19 N taking into account the layout of the future POWTEX instrument, the profile function as parameterized in 2θ and λ is extracted. As a proof-of-concept, the resulting instrument parameterization is then utilized to perform a typical refinement of the angular- and wavelength-dispersive diffraction pattern of CuNCN, yielding excellent residuals within feasible computational efforts. Another proof-of-concept is carried out by applying the same approach to a real neutron diffraction data set of CuNCN obtained from the POWGEN instrument at the Spallation Neutron Source in Oak Ridge. The paper highlights the general importance of the novel approach for data analysis at neutron time-of-flight diffractometers and its possible inclusion within existing Rietveld software packages.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Generalized skew-symmetric interfacial probability distribution in reflectivity and small-angle scattering analysis

Generalized skew-symmetric probability density functions are proposed to model asymmetric interfacial density distributions for the parameterization of any arbitrary density profiles in the `effective-density model'. The penetration of the densities into adjacent layers can be selectively controlled and parameterized. A continuous density profile is generated and discretized into many independent slices of very thin thickness with constant density values and sharp interfaces. The discretized profile can be used to calculate reflectivities via Parratt's recursive formula, or small-angle scattering via the concentric onion model that is also developed in this work.

36 MATERIALS SCIENCE↗

DOE OSTI2023

An Assessment of Additively Manufactured Bonded Permanent Magnets for a Distributed Wind Generator

In this paper, we examine and compare the performance of a generator design optimized using additively manufactured NdFeB-SmFeN in nylon-polymer-bonded permanent magnets (PMs) against a generator design with conventional NdFeB sintered PMs. To realize this, a commercially available 15-kW wind generator's rotor is re-optimized using both additively manufactured and sintered NdFeB magnets using simple geometric parameterization that allowed for two specific magnet shapes, namely, arc-shaped and crown-shaped designs. Results showed that for a similar generator performance, the designs with additively manufactured bonded PMs are more cost-competitive in terms of the estimated PM material cost and also have negligible eddy current magnet losses.

additive manufacturing↗

DOE OSTI2017

Modeling DNP3 Traffic Characteristics of Field Devices in SCADA Systems of the Smart Grid

In the generation, transmission, and distribution sectors of the smart grid, intelligence of field devices is realized by programmable logic controllers (PLCs). Many smart-grid subsystems are essentially cyber-physical energy systems (CPES): For instance, the power system process (i.e., the physical part) within a substation is monitored and controlled by a SCADA network with hosts running miscellaneous applications (i.e., the cyber part). To study the interactions between the cyber and physical components of a CPES, several co-simulation platforms have been proposed. However, the network simulators/emulators of these platforms do not include a detailed traffic model that takes into account the impacts of the execution model of PLCs on traffic characteristics. As a result, network traces generated by co-simulation only reveal the impacts of the physical process on the contents of the traffic generated by SCADA hosts, whereas the distinction between PLCs and computing nodes (e.g., a hardened computer running a process visualization application) has been overlooked. To generate realistic network traces using co-simulation for the design and evaluation of applications relying on accurate traffic profiles, it is necessary to establish a traffic model for PLCs. In this work, we propose a parameterized model for PLCs that can be incorporated into existing co-simulation platforms. We focus on the DNP3 subsystem of slave PLCs, which automates the processing of packets from the DNP3 master. To validate our approach, we extract model parameters from both the configuration and network traces of real PLCs. Simulated network traces are generated and compared against those from PLCs. Our evaluation shows that our proposed model captures the essential traffic characteristics of DNP3 slave PLCs, which can be used to extend existing co-simulation platforms and gain further insights into the behaviors of CPES.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2020

Simulating the Autonomous Future: A Look at Virtual Vehicle Environments and How to Validate Simulation Using Public Data Sets

The rapid evolution of autonomous vehicles (AVs) has exposed the need for fast-paced development and testing processes of a variety of perception, planning, and control algorithms. To expedite development, the AV industry and researchers leverage virtual vehicle environments to simulate a range of test scenarios that may otherwise be costly or difficult to conduct on a real test track. However, the various virtual environments may have different results depending on the fidelity of various simulation features, such as vehicle dynamics, sensor simulation, and environment recreation. Herein, this tutorial article examines a proposed framework for constructing, parameterizing, and validating a virtual vehicle environment using an existing AV data set. First, an overview of several open source and commercially available simulation tools, including their associated workflows, for scene and scenario creation is presented. Next, various open AV data sets are examined to inform the data set selection for the validation framework. Then, an example workflow of recreating a real-world scene from the selected data set in a simulation tool with various emulated sensors parameterized to match the data set is demonstrated. Finally, an example AV-perception algorithm is subjected to data streams from virtual and real-world environments and suggested metrics for analyzing the results are discussed.

42 ENGINEERING↗

DOE OSTI2016

A Scalable Method for Extracting Soiling Rates from PV Production Data

We present a method for analyzing time series production data from photovoltaic systems to extract the rate at which energy yield is affected by the accumulation of dust, dirt, and other forms of soiling. We describe an approach that is based on prevailing methods which consider the change in energy production during dry periods. The method described here builds upon these methods by considering a statistical sample of soiling intervals from each site under consideration. The method enables straightforward application to a large number of sites with minimal parameterization of data-filtering requirements. Furthermore, it enables statistical confidence intervals and comparisons between sites.

14 SOLAR ENERGY↗

DOE OSTI2020

Digital zero noise extrapolation for quantum error mitigation

Zero-noise extrapolation (ZNE) is an increasingly popular technique for mitigating errors in noisy quantum computations without using additional quantum resources. We review the fundamentals of ZNE and propose several improvements to noise scaling and extrapolation, the two key components in the technique. We introduce unitary folding and parameterized noise scaling. Furthermore, these are digital noise scaling frameworks, i.e. one can apply them using only gate-level access common to most quantum instruction sets. We also study different extrapolation methods, including a new adaptive protocol that uses a statistical inference framework. Benchmarks of our techniques show error reductions of 18X to 24X over non-mitigated circuits and demonstrate ZNE effectiveness at larger qubit numbers than have been tested previously. In addition to presenting new results, this work is a self-contained introduction to the practical use of ZNE by quantum programmers.

quantum computing↗

DOE OSTI2025

Optimal Power Management for Large-Scale Battery Energy Storage Systems via Bayesian Inference

Large-scale battery energy storage systems (BESS) have found ever-increasing use across industry and society to accelerate clean energy transition and improve energy supply reliability and resilience. However, their optimal power management poses significant challenges: the underlying high-dimensional nonlinear nonconvex optimization lacks computational tractability in real-world implementation, and the uncertainty of the exogenous power demand makes exact optimization difficult. This paper presents a new solution framework to address these bottlenecks. The solution pivots on introducing power-sharing ratios to specify each cell’s power quota from the output power demand. To find the optimal power-sharing ratios, we formulate a nonlinear model predictive control (NMPC) problem to achieve power-loss-minimizing BESS operation while complying with safety, cell balancing, and power supply-demand constraints. We then propose a parameterized control policy for the power-sharing ratios, which utilizes only three parameters, to reduce the computational demand in solving the NMPC problem. This policy parameterization allows us to translate the NMPC problem into a Bayesian inference problem for the sake of 1) computational tractability, and 2) overcoming the nonconvexity of the optimization problem. We leverage the ensemble Kalman inversion technique to solve the parameter estimation problem. Concurrently, a low-level control loop is developed to seamlessly integrate our proposed approach with the BESS to ensure practical implementation. This low-level controller receives the optimal power-sharing ratios, generates output power references for the cells, and maintains a balance between power supply and demand despite uncertainty in output power. We conduct extensive simulations and experiments on a 20-cell prototype to validate the proposed approach.

Battery energy storage systems (BESSs)↗

DOE OSTI2021

Parameter Reduction of Composite Load Model Using Active Subspace Method

Over the past decades, the increasing penetration of distributed energy resources (DERs) has dramatically changed the power load composition in the distribution networks. The traditional static and dynamic load models can hardly capture the dynamic behavior of modern loads especially for fault-induced delayed voltage recovery (FIDVR) events. Thus, a more comprehensive composite load model with combination of static load, different types of induction motors, single-phase A/C motor, electronic load and DERs has been proposed by Western Electricity Coordinating Council (WECC). However, due to the large number of parameters and model complexity, the WECC composite load model (WECC CMLD) raises new challenges to power system studies. To overcome these challenges, in this paper, a cutting-edge parameter reduction (PR) approach for WECC CMLD based on active subspace method (ASM) is proposed. Firstly, the WECC CMLD is parameterized in a discrete-time manner for the application of the proposed method. Then, parameter sensitivities are calculated by discovering the active subspace, which is a lower-dimensional linear subspace of the parameter space of WECC CMLD in which the dynamic response is most sensitive. The interdependency among parameters can be taken into consideration by our approach. Finally, the numerical experiments validate the effectiveness and advantages of the proposed approach for WECC CMLD model.

active subspace↗

DOE OSTI2020

The response of stomatal conductance to seasonal drought in tropical forests

Stomata regulate CO2 uptake for photosynthesis and water loss through transpiration. The approaches used to represent stomatal conductance (gs) in models vary. In particular, current understanding of drivers of the variation in a key parameter in those models, the slope parameter (i.e. a measure of intrinsic plant water-use-efficiency), is still limited, particularly in the tropics. Here we collected diurnal measurements of leaf gas exchange and water potential (?leaf), and a suite of plant traits from the upper canopy of 15 tropical trees in two contrasting Panamanian forests throughout the dry season of the 2016 El Niño. The plant traits included wood density, leaf-mass-per-area (LMA), leaf carboxylation capacity (Vc,max25), leaf water content, the degree of isohydry, and predawn ?leaf. We first investigated how the choice of four commonly used leaf-level gs models with and without the inclusion of ?leaf as an additional predictor variable influence the ability to predict gs, and then explored the abiotic (i.e. month, site-month interaction) and biotic (i.e. tree-species-specific characteristics) drivers of slope parameter variation. Our results show that the inclusion of ?leaf did not improve model performance and that the models that represent the response of gs to vapor pressure deficit performed better than corresponding models that respond to relative humidity. Within each gs model, we found large variation in the slope parameter, and this variation was attributable to the biotic driver, rather than abiotic drivers. We further investigated potential relationships between the slope parameter and the six available plant traits mentioned above, and found that only one trait, LMA, had a significant correlation with the slope parameter (R2=0.66, n=15), highlighting a potential path towards improved model parameterization. This study advances understanding of gs dynamics over seasonal drought, and identifies a practical, trait-based approach to improve modeling of carbon and water exchange in tropical forests.

carbon and water exchange, stomatal conductance mo↗

DOE OSTI2016

A frequentist approach to computer model calibration

The paper considers the computer model calibration problem and provides a general frequentist solution. Under the framework proposed, the data model is semiparametric with a non-parametric discrepancy function which accounts for any discrepancy between physical reality and the computer model. In an attempt to solve a fundamentally important (but often ignored) identifiability issue between the computer model parameters and the discrepancy function, the paper proposes a new and identifiable parameterization of the calibration problem. It also develops a two-step procedure for estimating all the relevant quantities under the new parameterization. This estimation procedure is shown to enjoy excellent rates of convergence and can be straightforwardly implemented with existing software. For uncertainty quantification, bootstrapping is adopted to construct confidence regions for the quantities of interest. As a result, the practical performance of the methodology is illustrated through simulation examples and an application to a computational fluid dynamics model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Parametric Modeling and Economic Analysis of a 2MW th 3-Stream sCO 2 Heat Exchanger

Here, this paper presents the design and cost optimization of a novel 2MW th 3-stream sCO2 plate-fin heat exchanger. This heat exchanger design is unique in that it uses reduced metal oxide particle-to-sCO2 heat exchanger for cost-effective energy storage applications. The design uses low velocity, laminar air as the re-oxidizing reactant to transfer the heat of the re-oxidizing reaction to a sCO2 power loop. The design of the heat exchanger is based on a 2-D, 3-fluid plate/fin heat transfer model. The model parameterizes the size, shape, and number of passages of the heat exchanger to calculate the temperature profile, pressure drop, and fluid velocities of all three fluids. Global heat exchanger parameters such as the effectiveness and total heat transferred to the sCO2 are then calculated for overall performance. Due to the value and increased use of sCO2 heat exchangers in power cycles, a cost model of the system based on the unique high temperature/high pressure operating conditions was created using quotes from reference projects and market analysis. These quoted air-to-sCO2 heat exchangers are then processed using multiple weighting factors pertinent to heat exchanger design, including heat exchanger type, maximum temperature, differential pressures, fluids, duty, and more. These factors are then used in an exponential function in order to generate a parameterized cost curve. The design and cost of the heat exchanger are then optimized using the SMPSO genetic algorithm in Python. The optimization objectives for the system are to maximize the overall system effectiveness, including an air recuperator for preheating, and to minimize unit costs. Additional constraints are added to the system for the sCO2 and air pressure drops, air velocity to reduce particle entrainment, and the length and volume of the heat exchanger.

Cost Model↗

DOE OSTI2011

Reliable automatic alignment of tomographic projection data by passive auto-focus

Purpose: The authors present a robust algorithm that removes the blurring and double-edge artifacts in high-resolution computed tomography (CT) images that are caused by misaligned scanner components. This alleviates the time-consuming process of physically aligning hardware, which is of particular benefit if components are moved or swapped frequently. Methods: The proposed method uses the experimental data itself for calibration. A parameterized model of the scanner geometry is constructed and the parameters are varied until the sharpest 3D reconstruction is found. The concept is similar to passive auto-focus algorithms of digital optical instruments. The parameters are used to remap the projection data from the physical detector to a virtual aligned detector. This is followed by a standard reconstruction algorithm, namely the Feldkamp algorithm. Feldkamp et al.[J. Opt. Soc. Am. A 1, 612-619 (1984)]. Results: An example implementation is given for a rabbit liver specimen that was collected with a circular trajectory. The optimal parameters were determined in less computation time than that for a full reconstruction. The example serves to demonstrate that (a) sharpness is an appropriate measure for projection alignment, (b) our parameterization is sufficient to characterize misalignments for cone-beam CT, and (c) the procedure determines parameter values with sufficient precision to remove the associated artifacts. Conclusions: The algorithm is fully tested and implemented for regular use at The Australian National University micro-CT facility for both circular and helical trajectories. It can in principle be applied to more general imaging geometries and modalities. It is as robust as manual alignment but more precise since we have quantified the effect of misalignment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2012

Development of array-type prompt gamma measurement system for in vivo range verification in proton therapy

Purpose: In vivo range verification is one of the most important parts of proton therapy to fully utilize its benefits delivering high radiation dose to tumor, while sparing the normal tissue with the so-called Bragg peak. Currently, however, range verification method is not used in clinics. The purpose of the present study is to optimize and evaluate the configuration of an array-type prompt gamma measurement system on determining distal dose edge for in vivo range verification of proton therapy. Methods: To effectively measure the prompt gammas against the background gammas, the Monte Carlo simulations with the MCNPX code were employed in optimizing the configuration of the measurement system, and the Monte Carlo method was also used to understand the effect of the background gammas, mainly neutron capture gammas, in the measured gamma distribution. To reduce the effect of the background gammas, the optimized energy window of 4-10 MeV in measuring the prompt gammas was employed. A parameterized source was used to maximize computation speed in the optimization study. A simplified test measurement system, using only one detector moving from one measurement location to the next, was constructed and applied to therapeutic proton beams of 80-220 MeV. For accurate determination of the distal dose edge, the sigmoidal curve-fitting method was applied to the measured distributions of the prompt gammas, and then, the location of the half-value between the maximum and minimum value in the curve-fitting was determined as the distal dose edge and compared with the beam range assessed by the proton dose distribution. Results: The parameterized source term employed in optimization process improved the calculation speed by up to {approx}300 times. The optimization study indicates that an array-type measurement system with 3, 2, 2, and 150 mm for scintillator thickness, slit width, septal thickness, and slit length, respectively, can effectively measure the prompt gamma distributions minimizing the contribution of background gammas. The present results show that a few hundred counts of prompt gammas can be easily obtained by measuring 10 s at each measurement location for proton beams of {approx}4 nA. The distal dose edges determined by the prompt gamma distribution are 5.45, 14.73, and 27.74 cm for proton beams of 5.17 (80 MeV), 14.99 (150 MeV), and 27.38 (220 MeV) cm, respectively. Conclusions: The results show that the array-type measurement system can measure prompt gamma distributions from a therapeutic proton beam within a short measurement time, and that the distal dose edge can be determined within a few millimeters of error without using any sophisticated analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

WE-A-17A-04: Development of An Ultra-Fast Monte Carlo Dose Engine for High Dose Rate Brachytherapy

Purpose: Current clinical brachytherapy dose calculations are based on AAPM TG43 guidelines, which approximate the patient geometry as a large water phantom. This ignores heterogeneities and tends to overestimate skin dose. Although Monte Carlo (MC) dose calculations have been recognized as the most accurate method, its associated long computational time is a major bottleneck for routine clinical applications. This work aims to develop a GPU-based ultra-fast MC dose engine (gBMC) for HDR brachytherapy to provide clinical users with accurate sub-minute dose calculations. Methods: Standard photon transport with discrete events including: Compton scattering, Rayleigh scattering and Photoelectric effect, was implemented. Secondary electrons were transported under the continuous slowing down approximation. To reduce the GPU thread divergence, photons and electrons were separately transported. Transport of photons was grouped according to energy. The source model in gBMC can be either a phase-space file generated using Geant4 or a parameterized source model. This dose engine was validated against TG43 in a water phantom and against Geant4 calculations in heterogeneous patient geometries. Results: A phase space file was generated for the Varian VS2000 Ir-192 source. In a water phantom, the calculated radial dose function was within 0.6% of the TG43 calculations for radial distances from 1 cm to 20 cm. The anisotropy functions were within 1% for radial distances from 1 cm to 20 cm except for polar angles larger than 173°. Local point-dose differences were within 2%. In a Mammosite breast cancer case with 22 dwell locations, gBMC and Geant4 isodose lines compared well. The computation time was about 28 seconds using the phase-space file source and 20 seconds using the parameterized source to simulate 1 billion particles, yielding less than 1% statistical uncertainty. Conclusion: The gBMC dose engine makes it possible to use fast and accurate MC dose calculations for clinical work.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

SU-E-P-59: A Graphical Interface for XCAT Phantom Configuration, Generation and Processing

Purpose: To design a comprehensive open-source, publicly available, graphical user interface (GUI) to facilitate the configuration, generation, processing and use of the 4D Extended Cardiac-Torso (XCAT) phantom. Methods: The XCAT phantom includes over 9000 anatomical objects as well as respiratory, cardiac and tumor motion. It is widely used for research studies in medical imaging and radiotherapy. The phantom generation process involves the configuration of a text script to parameterize the geometry, motion, and composition of the whole body and objects within it, and to generate simulated PET or CT images. To avoid the need for manual editing or script writing, our MATLAB-based GUI uses slider controls, drop-down lists, buttons and graphical text input to parameterize and process the phantom. Results: Our GUI can be used to: a) generate parameter files; b) generate the voxelized phantom; c) combine the phantom with a lesion; d) display the phantom; e) produce average and maximum intensity images from the phantom output files; f) incorporate irregular patient breathing patterns; and f) generate DICOM files containing phantom images. The GUI provides local help information using tool-tip strings on the currently selected phantom, minimizing the need for external documentation. The DICOM generation feature is intended to simplify the process of importing the phantom images into radiotherapy treatment planning systems or other clinical software. Conclusion: The GUI simplifies and automates the use of the XCAT phantom for imaging-based research projects in medical imaging or radiotherapy. This has the potential to accelerate research conducted with the XCAT phantom, or to ease the learning curve for new users. This tool does not include the XCAT phantom software itself. We would like to acknowledge funding from MRA, Varian Medical Systems Inc.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-113: Inherent Functional Dependence of Spinal Cord Doses of Variable Irradiated Volumes in Spine SBRT

Purpose: Spinal cord tolerance for SBRT has been recommended for the maximum point dose level or at irradiated volumes such as 0.35 mL or 10% of contoured volumes. In this study, we investigated an inherent functional relationship that associates these dose surrogates for irradiated spinal cord volumes of up to 3.0 mL. Methods: A hidden variable termed as Effective Dose Radius (EDR) was formulated based on a dose fall-off model to correlate dose at irradiated spinal cord volumes ranging from 0 mL (point maximum) to 3.0 mL. A cohort of 15 spine SBRT cases was randomly selected to derive an EDR-parameterized formula. The mean prescription dose for the studied cases was 21.0±8.0 Gy (range, 10–40Gy) delivered in 3±1 fractions with target volumes of 39.1 ± 70.6 mL. Linear regression and variance analysis were performed for the fitting parameters of variable EDR values. Results: No direct correlation was found between the dose at maximum point and doses at variable spinal cord volumes. For example, Pearson R{sup 2} = 0.643 and R{sup 2}= 0.491 were obtained when correlating the point maximum dose with the spinal cord dose at 1 mL and 3 mL, respectively. However, near perfect correlation (R{sup 2} ≥0.99) was obtained when corresponding parameterized EDRs. Specifically, Pearson R{sup 2}= 0.996 and R{sup 2} = 0.990 were obtained when correlating EDR (maximum point dose) with EDR (dose at 1 mL) and EDR(dose at 3 mL), respectively. As a result, high confidence level look-up tables were established to correlate spinal cord doses at the maximum point to any finite irradiated volumes. Conclusion: An inherent functional relationship was demonstrated for spine SBRT. Such a relationship unifies dose surrogates at variable cord volumes and proves that a single dose surrogate (e.g. point maximum dose) is mathematically sufficient in constraining the overall spinal cord dose tolerance for SBRT.

60 APPLIED LIFE SCIENCES↗