Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “control volume methods”

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 685 records · Page 38

Controlling sign problems in spin models using tensor renormalization

We consider the sign problem for classical spin models at complex $\beta =1/g_0^2$ on $L\times L$ lattices. We show that the tensor renormalization group method allows reliable calculations for larger Im$\beta$ than the reweighting Monte Carlo method. For the Ising model with complex $\beta$ we compare our results with the exact Onsager-Kaufman solution at finite volume. The Fisher zeros can be determined precisely with the TRG method. We check the convergence of the TRG method for the O(2) model on $L\times L$ lattices when the number of states $D_s$ increases. We show that the finite size scaling of the calculated Fisher zeros agrees very well with the Kosterlitz-Thouless transition assumption and predict the locations for larger volume. The location of these zeros agree with Monte Carlo reweighting calculation for small volume. The application of the method for the O(2) model with a chemical potential is briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

System and method for computed axial lithography (CAL) for 3D additive manufacturing

In one aspect the present disclosure relates to a system for forming a three dimensional (3D) object from a volume of photo-responsive material contained within a material container. A plurality of optical projection subsystems are adapted to be arranged at least partially circumferentially around a material container containing the volume of photo-responsive material, and are controlled by a controller. The optical projection subsystems direct optical projections at a plurality of angles θ through the volume of photo-responsive material. The optical projections are further directed about a z axis extending through the volume of photo-responsive material. The plurality of optical subsystems are further controlled by the controller to provide each of the optical projections with a calculated three-dimensional intensity distribution acting over a fixed temporal exposure period, which is sufficient to at least one of cure or remove selected portions of the volume of photo-responsive material, and leave other portions unmodified, to form a desired 3D part.

Kelly, Brett↗

Accumulated Delivered Dose Response of Stereotactic Body Radiation Therapy for Liver Metastases

Purpose: To determine whether the accumulated dose using image guided radiation therapy is a stronger predictor of clinical outcomes than the planned dose in stereotactic body radiation therapy (SBRT) for liver metastases. Methods and Materials: From 2003 to 2009, 81 patients with 142 metastases were treated in institutional review board–approved SBRT studies (5-10 fractions). Patients were treated during free breathing (with or without abdominal compression) or with controlled exhale breath-holding. SBRT was planned on a static exhale computed tomography (CT) scan, and the minimum planning target volume dose to 0.5 cm{sup 3} (minPTV) was recorded. The accumulated minimum dose to the 0.5 cm{sup 3} gross tumor volume (accGTV) was calculated after performing dose accumulation from exported image guided radiation therapy data sets registered to the planning CT using rigid (2-dimensional MV/kV orthogonal) or deformable (3-dimensional/4-dimensional cone beam CT) image registration. Univariate and multivariate Cox regression models assessed the factors influencing the time to local progression (TTLP). Hazard ratios for accGTV and minPTV were compared using model goodness-of-fit and bootstrapping. Results: Overall, the accGTV dose exceeded the minPTV dose in 98% of the lesions. For 5 to 6 fractions, accGTV doses of >45 Gy were associated with 1-year local control of 86%. On univariate analysis, the cancer subtype (breast), smaller tumor volume, and increased dose were significant predictors for improved TTLP. The dose and volume were uncorrelated; the accGTV dose and minPTV dose were correlated and were tested separately on multivariate models. Breast cancer subtype, accGTV dose (P<.001), and minPTV dose (P=.02) retained significance in the multivariate models. The univariate hazard ratio for TTLP for 5-Gy increases in accGTV versus minPTV was 0.67 versus 0.74 (all patients; 95% confidence interval of difference 0.03-0.14). Goodness-of-fit testing confirmed the accGTV dose as a stronger dose–response predictor than the minPTV dose. Conclusions: The accGTV dose is a better predictor of TTLP than the minPTV dose for liver metastasis SBRT. The use of modern image guided radiation therapy in future analyses of dose–response outcomes should increase the concordance between the planned and delivered doses.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

One-pot synthesis of co-substituted manganese oxide nanosheets and physical properties of lamellar aggregates

Graphical abstract: Display Omitted Highlights: ► Solid solution nanosheets, (Mn{sub 1−x}Co{sub x})O{sub 2}, synthesized via facile one-pot process. ► The structural characterization of nanosheets revealing a single (Mn,Co)O{sub 2} layer and the solubility limit as x ∼ 0.20. ► The invariant charge density of the layer upon Co substitution. ► Systematic dependence of magnetic and optical properties of the lamellar aggregates. -- Abstract: Co-substituted manganese oxide nanosheets, (Mn{sub 1−x}Co{sub x})O{sub 2} have been synthesized in the form of a colloidal suspension via a simple one-pot method. Substitution effects on the structural, optical absorption, and magnetic properties are investigated for the nanosheets and their lamellar aggregates. The composition of the (Mn{sub 1−x}Co{sub x})O{sub 2} nanosheets can be controlled continuously by adjusting the molar ratio of the starting materials. The solubility limit is x ∼ 0.20 based on the cell volume. In the 0.00 ≤ x ≤ 0.20 range, the band gap energy, magnetic moment, and Weiss temperature change systematically with x. The charge density of the (Mn,Co)O{sub 2} layer is independent of x (i.e., [(Mn,Co)O{sub 2}]{sup 0.2−}) and the cobalt ions are trivalent in low-spin state.

36 MATERIALS SCIENCE↗

SU-E-T-454: Impact of Calculation Grid Size On Dosimetry and Radiobiological Parameters for Head and Neck IMRT

Purpose: IMRT has become standard of care for complex treatments to optimize dose to target and spare normal tissues. However, the impact of calculation grid size is not widely known especially dose distribution, tumor control probability (TCP) and normal tissue complication probability (NTCP) which is investigated in this study. Methods: Ten head and neck IMRT patients treated with 6 MV photons were chosen for this study. Using Eclipse TPS, treatment plans were generated for different grid sizes in the range 1–5 mm for the same optimization criterion with specific dose-volume constraints. The dose volume histogram (DVH) was calculated for all IMRT plans and dosimetric data were compared. ICRU-83 dose points such as D2%, D50%, D98%, as well as the homogeneity and conformity indices (HI, CI) were calculated. In addition, TCP and NTCP were calculated from DVH data. Results: The PTV mean dose and TCP decreases with increasing grid size with an average decrease in mean dose by 2% and TCP by 3% respectively. Increasing grid size from 1–5 mm grid size, the average mean dose and NTCP for left parotid was increased by 6.0% and 8.0% respectively. Similar patterns were observed for other OARs such as cochlea, parotids and spinal cord. The HI increases up to 60% and CI decreases on average by 3.5% between 1 and 5 mm grid that resulted in decreased TCP and increased NTCP values. The number of points meeting the gamma criteria of ±3% dose difference and ±3mm DTA was higher with a 1 mm on average (97.2%) than with a 5 mm grid (91.3%). Conclusion: A smaller calculation grid provides superior dosimetry with improved TCP and reduced NTCP values. The effect is more pronounced for smaller OARs. Thus, the smallest possible grid size should be used for accurate dose calculation especially in H and N planning.

60 APPLIED LIFE SCIENCES↗

TU-F-CAMPUS-J-01: Dosimetric Effects of HU Changes During the Course of Proton Therapy for Lung Cancer

Purpose: To characterize the changes in Hounsfield unit (HU) in lung radiotherapy with proton beams during the course of treatment and to study the effect on the proton plan dose distribution. Methods: Twenty consecutive patients with non-small cell lung cancer treated with proton radiotherapy who underwent multiple CT scans including the planning CT and weekly verification CTs were studied. HU histograms were computed for irradiated lung volumes in beam paths for all scans using the same treatment plan. Histograms for un-irradiated lung volume were used as control to characterize inter-scan variations. HU statistics were calculated for both irradiated and un-irradiated lung volumes for each patient scan. Further, multiple CT scans based on the same planning CT were generated by replacing the HU of the lung based on the verification CT scans HU values. Using the same beam arrangement, we created plans for each of the altered CT scans to study the dosimetric effect using the dose volume histogram. Results: Lung HU decreased for irradiated lung volume during the course of radiotherapy. The magnitude of this change increased with total irradiation dose. On average, HU changed by −53.8 in the irradiated volume. This change resulted in less than 0.5mm of beam overshoot in tissue for every 1cm beam traversed in the irradiated lung. The dose modification is about +3% for the lung, and less than +1% for the primary tumor. Conclusion: HU of the lung decrease throughout the course of radiation therapy. This change results in a beam overshoot (e.g. 3mm for 6cm of lung traversed) and causes a small dose modification in the overall plan. However, this overshoot does not affect the quality of plans since the margins used in planning, based on proton range uncertainty, are greater. HU needs to change by 150 units before re-planning is warranted.

60 APPLIED LIFE SCIENCES↗

Reverse-Mode Algorithmic Differentiation of an OpenMP-Parallel Compressible Flow Solver

y Reverse-mode algorithmic differentiation (AD) is an established method for obtaining adjoint derivatives of computer simulation applications. In computational fluid dynamics (CFD), adjoint derivatives of a cost function output such as drag or lift with respect to design parameters such as surface coordinates or geometry control points are a key ingredient for shape optimization, uncertainty quantification and flow control. The computational cost of CFD applications and their derivatives makes it essential to use high-performance computing hardware efficiently, including multi- and many-core architectures. Nevertheless, OpenMP is not supported in most AD tools, and previously shown methods achieve poor scalability of the derivative code. We present the AD of an OpenMP-parallelized finite volume compressible flow solver for unstructured meshes. Our approach enables us to reuse the parallelization of the original code in the computation of adjoint derivatives. The method works by identifying code segments that can be differentiated in reverse-mode without changing their memory access pattern. The OpenMP parallelization is integrated into the derivative code during the build process in a way that is robust to modifications of the original code and independent of the OpenMP support of the differentiation tool. We show the scalability of our adjoint CFD solver on test cases ranging from thousands to millions of finite volume mesh cells on CPUs with up to 16 threads as well as on an Intel XeonPhi card with 236 threads. We demonstrate that our approach is more practical to implement for production-sized CFD codes and produces more efficient adjoint derivative code than previously shown AD methods.

Huckelheim, Jan↗

Near-IR absorption in high-purity photothermorefractive glass and holographic optical elements: measurement and application for high-energy lasers

Volume Bragg gratings (VBGs) in photothermorefractive (PTR) glass are widely used for laser beam control including high-power laser systems. Among them, spectral beam combining based on VBGs is one of the most promising. Achieving 100+ kW of combined laser beams requires the development of PTR glass and VBGs with an extremely low absorption coefficient and therefore methods of its measurement. This paper describes the calorimetric method that was developed for measuring a low absorption coefficient in PTR glass and VBGs. It is based on transmission monitoring of the intrinsic Fabry-Perot interferometer produced by the plane-parallel surfaces of the measured optical elements when heated by high-power laser radiation. An absorption coefficient at 1085 nm as low as 5x10{sup -5} cm{sup -1} is demonstrated in pristine PTR glass while an absorption coefficient as low as 1x10{sup -4} cm{sup -1} is measured in high-efficiency reflecting Bragg gratings with highest purity. The actual level of absorption in PTR glass allows laser beam control at the 10 kW level, while the 100 kW level would require active cooling and/or decreasing the absorption in PTR Bragg gratings to a value similar to that in virgin PTR glass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Fifth International Workshop on Ice Nucleation phase 2 (FIN-02): laboratory intercomparison of ice nucleation measurements

Abstract. The second phase of the Fifth International Ice Nucleation Workshop (FIN-02)involved the gathering of a large number of researchers at the KarlsruheInstitute of Technology's Aerosol Interactions and Dynamics of the Atmosphere(AIDA) facility to promote characterization and understanding of icenucleation measurements made by a variety of methods used worldwide.Compared to the previous workshop in 2007, participation was doubled,reflecting a vibrant research area. Experimental methods involved sampling ofaerosol particles by direct processing ice nucleation measuring systems fromthe same volume of air in separate experiments using different ice nucleatingparticle (INP) types, and collections of aerosol particle samples ontofilters or into liquid for sharing amongst measurement techniques thatpost-process these samples. In this manner, any errors introduced bydifferences in generation methods when samples are shared across laboratorieswere mitigated. Furthermore, as much as possible, aerosol particle sizedistribution was controlled so that the size limitations of different methodswere minimized. The results presented here use data from the workshop toassess the comparability of immersion freezing measurement methods activatingINPs in bulk suspensions, methods that activate INPs in condensation and/orimmersion freezing modes as single particles on a substrate, continuous flowdiffusion chambers (CFDCs) directly sampling and processing particles wellabove water saturation to maximize immersion and subsequent freezing ofaerosol particles, and expansion cloud chamber simulations in which liquidcloud droplets were first activated on aerosol particles prior to freezing.The AIDA expansion chamber measurements are expected to be the closestrepresentation to INP activation in atmospheric cloud parcels in thesecomparisons, due to exposing particles freely to adiabatic cooling. The different particle types used as INPs included the minerals illite NX andpotassium feldspar (K-feldspar), two natural soil dusts representative of arable sandy loam(Argentina) and highly erodible sandy dryland (Tunisia) soils, respectively,and a bacterial INP (Snomax ® ). Consideredtogether, the agreement among post-processed immersion freezing measurementsof the numbers and fractions of particles active at different temperaturesfollowing bulk collection of particles into liquid was excellent, withpossible temperature uncertainties inferred to be a key factor in determiningINP uncertainties. Collection onto filters for rinsing versus directly intoliquid in impingers made little difference. For methods that activatedcollected single particles on a substrate at a controlled humidity at orabove water saturation, agreement with immersion freezing methods was good inmost cases, but was biased low in a few others for reasons that have not beenresolved, but could relate to water vapor competition effects. AmongstCFDC-style instruments, various factors requiring (variable) highersupersaturations to achieve equivalent immersion freezing activation dominatethe uncertainty between these measurements, and for comparison with bulkimmersion freezing methods. When operated above water saturation to includeassessment of immersion freezing, CFDC measurements often measured at orabove the upper bound of immersion freezing device measurements, but oftenunderestimated INP concentration in comparison to an immersion freezingmethod that first activates all particles into liquid droplets prior tocooling (the PIMCA-PINC device, or Portable Immersion Mode Cooling chAmber–Portable Ice Nucleation Chamber), and typically slightly underestimated INPnumber concentrations in comparison to cloud parcel expansions in the AIDAchamber; this can be largely mitigated when it is possible to raise therelative humidity to sufficiently high values in the CFDCs, although this isnot always possible operationally. Correspondence of measurements of INPs among direct sampling andpost-processing systems varied depending on the INP type. Agreement was bestfor Snomax ® particles in the temperature regimecolder than −10 ∘ C, where their ice nucleation activity is nearlymaximized and changes very little with temperature. At temperatures warmer than −10 ∘ C, Snomax ® INP measurements (allvia freezing of suspensions) demonstrated discrepancies consistent withprevious reports of the instability of its protein aggregates that appear tomake it less suitable as a calibration INP at these temperatures. ForArgentinian soil dust particles, there was excellent agreement across allmeasurement methods; measures ranged within 1 order of magnitude for INPnumber concentrations, active fractions and calculated active site densitiesover a 25 to 30 ∘ C range and 5 to 8 orders of correspondingmagnitude change in number concentrations. This was also the case for alltemperatures warmer than −25 ∘ C in Tunisian dust experiments. Incontrast, discrepancies in measurements of INP concentrations or active sitedensities that exceeded 2 orders of magnitude across a broad range of temperaturemeasurements found at temperatures warmer than −25 ∘ C in a previous study werereplicated for illite NX. Discrepancies also exceeded 2 orders of magnitude attemperatures of −20 to −25 ∘ C for potassium feldspar (K-feldspar), but these coincidedwith the range of temperatures at which INP concentrations increase rapidly atapproximately an order of magnitude per 2 ∘ C cooling forK-feldspar. These few discrepancies did not outweigh the overall positive outcomes of theworkshop activity, nor the future utility of this data set or future similarefforts for resolving remaining measurement issues. Measurements of the samematerials were repeatable over the time of the workshop and demonstratedstrong consistency with prior studies, as reflected by agreement of databroadly with parameterizations of different specific or general (e.g., soildust) aerosol types. The divergent measurements of the INP activity of illiteNX by direct versus post-processing methods were not repeated for otherparticle types, and the Snomax ® datademonstrated that, at least for a biological INP type, there is no expectedmeasurement bias between bulk collection and direct immediately processedfreezing methods to as warm as −10 ∘ C. Since particle size rangeswere limited for this workshop, it can be expected that for atmosphericpopulations of INPs, measurement discrepancies will appear due to thedifferent capabilities of methods for sampling the full aerosol sizedistribution, or due to limitations on achieving sufficient watersupersaturations to fully capture immersion freezing in direct processinginstruments. Overall, this workshop presents an improved picture of presentcapabilities for measuring INPs than in past workshops, and providesdirection toward addressing remaining measurement issues.

54 ENVIRONMENTAL SCIENCES↗

Suitability of Gray Water for Hydroponic Crop Production Following Biological and Physical Chemical and Biological Subsystems

The water present in waste streams from a human habitat must be recycled in Controlled Ecological Life Support Systems (CELSS) to limit resupply needs and attain self-sufficiency. Plants play an important role in providing food, regenerating air, and producing purified water via transpiration. However, we have shown that the surfactants present in hygiene waste water have acute toxic effects on plant growth (Bubenheim et al. 1994; Greene et al., 1994). These phytotoxic affects can be mitigated by allowing the microbial population on the root surface to degrade the surfactant, however, a significant suppression (several days) in crop performance is experienced prior to reaching sub-toxic surfactant levels and plant recovery. An effective alternative is to stabilize the microbial population responsible for degradation of the surfactant on an aerobic bioreactor and process the waste water prior to utilization in the hydroponic solution (Wisniewski and Bubenheim, 1993). A sensitive bioassay indicates that the surfactant phytotoxicity is suppressed by more than 90% within 5 hours of introduction of the gray water to the bioreactor; processing for more than 12 hours degrades more than 99% of the phytotoxin. Vapor Compression Distillation (VCD) is a physical / chemical method for water purification which employees sequential distillation steps to separate water from solids and to volatilize contaminants. The solids from the waste water are concentrated in a brine and the pure product water (70 - 90% of the total waste water volume depending on operating conditions) retains non of the phytotoxic effects. Results of the bioassay were used to guide evaluations of the suitability of recovered gray water following biological and VCD processing for hydroponic lettuce production in controlled environments. Lettuce crops were grown for 28 days with 100% of the input water supplied with recovered water from the biological processor or VCD. When compared with the growth of plants in control hydroponic solution containing pure deionized water, no growth difference could be measured resulting from any of the recovered water treatments. Both biological treatment and VCD offer alternative technology approaches to recovering water from waste streams appropriate for input into a crop production system. A high level of crop performance (food, air, and water production) can be maintained with either processor; selection decisions can be based on other factors regarding system integration.

Bubenheim, David L.↗

Stereotactic Ablative Radiation Therapy for Centrally Located Early Stage or Isolated Parenchymal Recurrences of Non-Small Cell Lung Cancer: How to Fly in a “No Fly Zone”

Purpose: We extended our previous experience with stereotactic ablative radiation therapy (SABR; 50 Gy in 4 fractions) for centrally located non-small cell lung cancer (NSCLC); explored the use of 70 Gy in 10 fractions for cases in which dose-volume constraints could not be met with the previous regimen; and suggested modified dose-volume constraints. Methods and Materials: Four-dimensional computed tomography (4DCT)-based volumetric image-guided SABR was used for 100 patients with biopsy-proven, central T1-T2N0M0 (n=81) or isolated parenchymal recurrence of NSCLC (n=19). All disease was staged with positron emission tomography/CT; all tumors were within 2 cm of the bronchial tree, trachea, major vessels, esophagus, heart, pericardium, brachial plexus, or vertebral body. Endpoints were toxicity, overall survival (OS), local and regional control, and distant metastasis. Results: At a median follow-up time of 30.6 months, median OS time was 55.6 months, and the 3-year OS rate was 70.5%. Three-year cumulative actuarial local, regional, and distant control rates were 96.5%, 87.9%, and 77.2%, respectively. The most common toxicities were chest-wall pain (18% grade 1, 13% grade 2) and radiation pneumonitis (11% grade 2 and 1% grade 3). No patient experienced grade 4 or 5 toxicity. Among the 82 patients receiving 50 Gy in 4 fractions, multivariate analyses showed mean total lung dose >6 Gy, V{sub 20} >12%, or ipsilateral lung V{sub 30} >15% to independently predict radiation pneumonitis; and 3 of 9 patients with brachial plexus D{sub max} >35 Gy experienced brachial neuropathy versus none of 73 patients with brachial D{sub max} <35 Gy (P=.001). Other toxicities were analyzed and new dose-volume constraints are proposed. Conclusions: SABR for centrally located lesions produces clinical outcomes similar to those for peripheral lesions when normal tissue constraints are respected.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Intracranial Effects of Artificial Gravity: A 3T MRI Study

INTRODUCTION Spaceflight associated neuro-ocular syndrome (SANS) is characterized by the development of optic disc edema, posterior globe flattening, choroidal/retinal folds and hyperopic refractive errors1. SANS is hypothesized to be a result of headward fluid shifts that invariably occurs in the microgravity environment. As a countermeasure, artificial gravity (AG) through centrifugation has been proposed to reduce this headward fluid shift, however there is no current proof of benefit. The goal of this study was to determine if the application of AG can prevent or reduce known changes in brain volumetry, internal carotid artery (ICA) stroke volume and cerebral spinal fluid (CSF) flow velocity that occurs during simulated chronic headward fluid shift using head down tilt bed rest (HDTBR) methodology2 as an indicator of countermeasure efficacy. METHODS Healthy volunteers were recruited for an IRB approved HDTBR study performed at the German Aerospace Center in Cologne, Germany. Strict six-degree HDTBR was used as a spaceflight analog to induce a continuous headward fluid shift. HDTBR was carried out for 60 days for all subjects. Short-arm centrifugation was utilized to generate AG equating to ~0.3g of acceleration at the level of the eye. The subjects were divided equally into three groups: NoAG (control; n=8), daily intermittent AG (6 x 5 min iAG; n=8), and daily continuous 30 min (cAG; n=8). All studies were performed on a single dedicated 3T MRI Scanner. Pulse-gated MRI phase-contrast flow imaging was used to quantify ICA stroke volume and peak-to-peak CSF flow velocity in the mid cerebral aqueduct. 3D-SPGR was acquired for volumetric segmentation of the brain and CSF spaces. MRI acquisitions were obtained at baseline (BDC), 14 days into HDTBR (HDTBR14), 52 days into HDTBR (HDTBR52) and 3-5 days after HDTBR (recovery, R+3/5).The data were analyzed by the mixed model, which included intervention and time (BDC, HDTBR 14, HDTBR 52, R+3/5) as the fixed effects and included subject as the random effect.RESULTS24 healthy subject volunteers (16 men, 8 women, mean age = 33 years ± 9 [standard deviation] and mean BMI = 24.3 kg/m2 ± 2.0) successfully completed all phases of the study. Strict six-degree HDTBR was characterized by progressive and statistically significant (p<.01) increases in mean combined brain and CSF volumes and mean aqueductal CSF peak-to-peak flow velocity, as well as statistically significant (p<.01) progressive decrease in mean ICA stroke volume from baseline to 52 days post intervention (Figs. 1-3). Compared to baseline, only combined brain and CSF volumes did not return to baseline values in the recovery period (p=NS). Neither iAG nor cAG exerted any significant effects on the measured MRI brain parameters as compared to HDTBR alone (p=NS). CONCLUSION Our results indicate that HDTBR at 6-degrees was effective in producing alterations in ICA stroke volume, aqueductal CSF flow velocity, and combined brain and CSF volumetric change that is associated with chronic headward fluid shift. Short duration, 30-min daily exposure to either iAG or cAG appears to be insufficient in preventing or reducing the effects of chronic HDTBR and thus may not be a suitable countermeasure as currently deployed. AG protocol modifications, including increased duration and magnitude of exposure, should be considered for future research.

L A Kramer↗

Plasmon-Mediated Intramolecular Methyl Migration with Nanoscale Spatial Control

Plasmonic materials interact strongly with light to focus and enhance electromagnetic radiation down to nanoscale volumes. Due to this localized confinement, materials that support localized surface plasmon resonances are capable of driving energetically unfavorable chemical reactions. In certain cases, the plasmonic nanostructures are able to preferentially catalyze the formation of specific photoproducts, which offers an opportunity for the development of solar-driven chemical synthesis. Here, using plasmonic environments, we report inducing an intramolecular methyl migration reaction, forming 4-methylpyridine from N-methylpyridinium. Using both experimental and computational methods, we were able to confirm the identity of the N-methylpyridinium by making spectral comparisons against possible photoproducts. This reaction involves breaking a C–N bond and forming a new C–C bond, highlighting the ability of plasmonic materials to drive complex and selective reactions. Additionally, we observe that the product yield depends strongly on optical illumination conditions. This is likely due to steric hindrance in specific regions on the nanostructured plasmonic substrate, providing an optical handle for driving plasmonic catalysis with spatial specificity. Furthermore, this work adds yet another class of reactions accessible by surface plasmon excitation to the ever-growing library of plasmon-mediated chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Subtask 1.1 – Advanced Characterization of Unconventional Oil and Gas Reservoirs to Enhance CO2 Storage Resource Estimates

The Energy & Environmental Research Center (EERC) developed methods and approaches for advanced characterization of unconventional oil and gas reservoirs to enhance CO2 storage resource estimates using samples collected from organic-rich shales of the Bakken Formation. The project was funded under U.S. Department of Energy (DOE) Cooperative Agreement No. DE-FE0024233, with nonfederal cost share provided by Hitachi High Technologies America. The project participants included the EERC, Hitachi High Technologies America, the team at the National Energy Technology Laboratory’s (NETL’s) computed tomography (CT) scanning laboratory in Morgantown, West Virginia, and several NETL staff who have worked on CO2 storage resource estimation. Investigations on the use of subsurface geologic formations for CO2 storage have been ongoing for many years. In the last decade, improvements in directional drilling technology and hydraulic fracturing of organic-rich shale and other tight rock formations have opened new opportunities for CO2 storage, with the potential cobenefit of enhancing gas and oil recovery. Gas-bearing, organic-rich shales have been of particular interest as a carbon storage target because their ability to store large volumes of CO2 is implied by their ability to store methane; however, estimating the CO2 storage resource of organic-rich shales and other unconventional reservoirs is challenging because of the lack of detailed geologic and petrophysical data that are needed to develop improved volumetric equations. According to NETL’s method for estimating regional-scale CO2 storage resources, the key factors that control CO2 storage in geologic formations include void space within induced fractures, natural fractures and matrix pores, as well as the occurrence and distribution of organic matter (OM) and clay minerals. A key goal of this effort was to use field emission scanning electron microscopy (FESEM) analysis of fine-grain unconventional reservoir rocks coupled with advanced image analysis and machine learning techniques to quantify the occurrence and distribution of porosity, OM, and clay minerals. By better understanding the potential migration pathways of CO2 within tight rocks and quantifying the occurrence of sorptive components (OM and clays) along potential flow pathways, methods to estimate the CO2 storage resource of such rocks can be refined. The Bakken petroleum system (Bakken) was used as the source of rock samples for this effort because it includes organic-rich shales with contrasting degrees of thermal maturity as well as tight, nonshale clastic rocks with varying clay mineral content. To better understand the key components of interest in tight rock samples, the EERC worked closely with Hitachi to optimize the FESEM settings necessary to resolve the sample mineralogy, porosity and OM content at micro- to nanoscale resolution. The EERC also worked closely with Hitachi to apply and improve its advanced mineral identification and characterization system (AMICS) software so that it is better suited for mineral identification at the high resolution necessary to resolve features of interest in fine-grained rocks such as organic-rich shales. The EERC also applied advanced image analysis software packages, including Ilastik, to classify and quantify the key components of interest within the samples using FESEM imagery. These advanced techniques allowed for quantification and differentiation of the void space within fractures, pores and OM, as well as a technique to determine and quantify the key mineral and OM components surrounding pore spaces. Being able to better quantify the void space and the potential accessibility of CO2 to sorptive minerals or OM within tight rocks allows for improved estimation of the CO2 storage resource potential of organic-rich shales and other unconventional formations. The information derived from the image analysis techniques was used to refine NETL’s method for estimating regional-scale CO2 storage resources in collaboration with staff from NETL. An additional goal of this work was to better understand rates of CO2 permeation, flow pathways, and sorption in organic-rich shales and nonshale reservoir rocks and to evaluate potential changes in CO2 flow over time. The EERC conducted long-duration (several weeks) CO2 flow-through testing at reservoir conditions using Bakken shale and nonshale samples to assess CO2 permeation rates and retention/sorption within core plugs to represent bulk matrix behavior. Conventional isothermal adsorption testing over a range of pressure conditions was also performed on pulverized Bakken samples to better assess potential maximum sorptive capacity of the samples. The results of this testing illustrated the higher sorptive capacity of the Bakken shales when compared to the sorptive capacity of the nonshale reservoir rock, which contains very little OM when compared to the Bakken shales. The results also illustrate the differences in permeation behavior as a function of fracture flow versus matrix flow. In an attempt to visualize CO2 permeation and flow behavior in Bakken shale samples over time, the EERC collaborated with NETL’s CT Lab in Morgantown, West Virginia. A variety of tests were performed at the NETL lab on unfractured and fractured samples of the Bakken shale. The unfractured Bakken shale samples did not result in CO2 migration within the samples over the testing duration; thus imagery of CO2 migration pathways within the matrix of Bakken samples was not obtained. The results of the flow-through testing within fractured samples illustrate changes in CO2 permeation rates over time, with a trend of decreasing CO2 flow over time. The results of this effort suggest that FESEM analysis coupled with advanced image analysis and machine learning techniques are effective in characterizing and quantifying the mineralogy, porosity, fractures, and OM content of 2-D imagery at the high resolution needed to identify features of interest in tight rocks. In addition, the image analysis techniques allow for identification and quantification of different porosity types as well as the key mineral components and/or OM that CO2 would contact in the reservoir as it migrates through pores and fracture networks. This is key information needed to better estimate the CO2 storage resource potential of organic-rich shales and other tight reservoirs. Additional work is needed to determine the sorptive capacity of the individual sample components (OM and clay minerals) at reservoir conditions so that this approach and the refined CO2 storage resource equation can be applied to other shale and unconventional reservoirs. In addition, subsequent work should be performed to better understand the number of images needed to use stereological approaches to estimate the volume of sample components using 2-D imagery, to explore machine learning techniques to differentiate between induced vs. natural fractures, and to estimate the connectivity of pores and fracture networks using 2-D imagery. Being able to derive this information at the microscale would facilitate scaling up of the sample characterization data used for reservoir modeling and simulation to estimate the CO2 storage resource potential of tight formations.

03 NATURAL GAS↗

Shape Optimization for Control and Isolation of Structural Vibrations in Aerospace and Defense Applications

Among the main challenges in shape optimization is the coupling of Finite Element Method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multi-physics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. There are two mathematically equivalent approaches to computing the shape derivative: the volume method, and the boundary method. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. Prior implementations of shape derivatives at Sandia have been based on the volume method. We introduce the strip method, which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. The development of the strip method also offers us the opportunity to share some lessons learned about implementing the volume method and boundary method, to show shape optimization results on problems of interest, and to begin addressing the other main challenges at hand: constraints on optimized shapes, and their interplay with optimization algorithms.

97 MATHEMATICS AND COMPUTING↗

Winter Simulation Conference, Miami Beach, Fla., December 4-6, 1978, Proceedings. Volumes 1 & 2

The papers report on the various aspects of simulation such as random variate generation, simulation optimization, ranking and selection of alternatives, model management, documentation, data bases, and instructional methods. Simulation studies in a wide variety of fields are described, including system design and scheduling, government and social systems, agriculture, computer systems, the military, transportation, corporate planning, ecosystems, health care, manufacturing and industrial systems, computer networks, education, energy, production planning and control, financial models, behavioral models, information systems, and inventory control.

Highland, H. J.↗

A new method for analyzing and visualizing plasma simulations using a phase-space tessellation

We apply a novel phase-space interpolation technique referred to as the simplex-in-cell (SIC) method to analyze two- and three-dimensional particle-in-cell (PIC) simulations of electromagnetic plasmas. SIC relies on a discretization of the initial phase-space distribution function into simplices, which allows an approximation to the full, continuously defined distribution function to be constructed at any later time in the simulation. This allows densities, currents, and even full momentum distribution functions to be measured at any point in the simulation domain without averaging over control volumes. The SIC approach applies to any PIC simulation for which a tessellation of the initial particle distribution can be constructed. In this study, we use outputs from standard PIC simulations of the Weibel instability and compare physical quantities such as charge and current densities calculated in postprocessing using SIC and standard particle deposits. Using 2D simulations with 1–65 536 particles-per-cell, we find that SIC eliminates discrete particle noise and in some cases can reach a given noise level using ~1000 times fewer simulation particles than with standard particle deposition schemes. In regions of low density, such as between current filaments, SIC is able to capture small amplitude features even with fewer particles than gridpoints due to the deformable nature of the SIC volume elements. Here, by calculating momentum distributions, we show how SIC can capture low density tails in the spectrum using far fewer particles than are necessary for standard particle deposits. We calculate the charge density on spatial grids of increasing resolution to demonstrate the ability of SIC to reveal fine-scale details that are not accessible with standard particle deposits. Finally, we show how SIC can be extended to 3D and give an example of its use to calculate the charge density from 3D PIC simulations of the Weibel instability. These results motivate the future implementation of SIC directly in the simulation force calculation for a novel low-noise electromagnetic plasma simulation method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Prognostic Factors for Survival in Patients Treated With Stereotactic Radiosurgery for Recurrent Brain Metastases After Prior Whole Brain Radiotherapy

Purpose: To evaluate prognostic factors for survival after stereotactic radiosurgery (SRS) for new, progressive, or recurrent brain metastases (BM) after prior whole brain radiotherapy (WBRT). Methods and Materials: Patients treated between 1991 and 2007 with Gamma Knife SRS for BM after prior WBRT were retrospectively reviewed. Potential prognostic factors were analyzed overall and by primary site using univariate and stepwise multivariate analyses and recursive partitioning analysis, including age, Karnofsky performance status (KPS), primary tumor control, extracranial metastases, number of BM treated, total SRS target volume, and interval from WBRT to SRS. Results: A total of 310 patients were analyzed, including 90 breast, 113 non-small-cell lung, 31 small-cell lung, 42 melanoma, and 34 miscellaneous patients. The median age was 56, KPS 80, number of BM treated 3, and interval from WBRT to SRS 8.1 months; 76% had controlled primary tumor and 60% had extracranial metastases. The median survival was 8.4 months overall and 12.0 vs. 7.9 months for single vs. multiple BM treated (p = 0.001). There was no relationship between number of BM and survival after excluding single-BM patients. On multivariate analysis, favorable prognostic factors included age <50, smaller total target volume, and longer interval from WBRT to SRS in breast cancer patients; smaller number of BM, KPS >60, and controlled primary in non-small-cell lung cancer patients; and smaller total target volume in melanoma patients. Conclusions: Among patients treated with salvage SRS for BM after prior WBRT, prognostic factors appeared to vary by primary site. Although survival time was significantly longer for patients with a single BM, the median survival time of 7.9 months for patients with multiple BM seems sufficiently long for salvage SRS to appear to be worthwhile, and no evidence was found to support the use of a cutoff for number of BM appropriate for salvage SRS.

62 RADIOLOGY AND NUCLEAR MEDICINE↗