Search NASA⌕ Search

SEARCH · Search NASA

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

Coupling the System Analysis Module with SAS4A/SASSYS-1

SAS4A/SASSYS-1 is a simulation tool used to perform deterministic analysis of anticipated events as well as design basis and beyond design basis accidents for advanced reactors, with an emphasis on sodium fast reactors. SAS4A/SASSYS-1 has been under development and in active use for nearly forty-five years, and is currently maintained by the U.S. Department of Energy under the Office of Advanced Reactor Technology. Although SAS4A/SASSYS-1 contains a very capable primary and intermediate system modeling component, PRIMAR-4, it also has some shortcomings: outdated data management and code structure makes extension of the PRIMAR-4 module somewhat difficult. The user input format for PRIMAR-4 also limits the number of volumes and segments that can be used to describe a given system. The System Analysis Module (SAM) is a fairly new code development effort being carried out under the U.S. DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM is being developed with advanced physical models, numerical methods, and software engineering practices; however, it is currently somewhat limited in the system components and phenomena that can be represented. For example, component models for electromagnetic pumps and multi-layer stratified volumes have not yet been developed. Nor is there support for a balance of plant model. Similarly, system-level phenomena such as control-rod driveline expansion and vessel elongation are not represented. This report documents fiscal year 2016 work that was carried out to couple the transient safety analysis capabilities of SAS4A/SASSYS-1 with the system modeling capabilities of SAM under the joint support of the ART and NEAMS programs. The coupling effort was successful and is demonstrated by evaluating an unprotected loss of flow transient for the Advanced Burner Test Reactor (ABTR) design. There are differences between the stand-alone SAS4A/SASSYS-1 simulations and the coupled SAS/SAM simulations, but these are mainly attributed to the limited maturity of the SAM development effort. The severe accident modeling capabilities in SAS4A/SASSYS-1 (sodium boiling, fuel melting and relocation) will continue to play a vital role for a long time. Therefore, the SAS4A/SASSYS-1 modernization effort should remain a high priority task under the ART program to ensure continued participation in domestic and international SFR safety collaborations and design optimizations. On the other hand, SAM provides an advanced system analysis tool, with improved numerical solution schemes, data management, code flexibility, and accuracy. SAM is still in early stages of development and will require continued support from NEAMS to fulfill its potential and to mature into a production tool for advanced reactor safety analysis. The effort to couple SAS4A/SASSYS-1 and SAM is the first step on the integration of these modeling capabilities.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Ultra-compact and efficient heat exchanger for solar desalination with unprecedented scaling resistance

The objective is to design, prototype, and test a heat exchanger technology with an efficiency improvement of more than 150% compared to state-of-the-art heat exchanger technology (UA/volume), with significantly improved fouling/scaling resistance at elevated temperatures (>120°C) when compared to existing technology where scaling is a major performance-limiting issue, as well as at lower operating-temperatures (~63-110°C). We now report our overall progress on the project. An ultra-compact, efficient heat exchanger geometry has been developed that can maximize surface-area-to-volume ratio so that proposed target (UA/volume). Advanced manufacturing techniques, such additive manufacturing, have made it possible to create complex structures cost effectively. In addition to traditional manufacturing methods, additive manufacturing has been successfully used for the development of high-performance geometries. These methods allow the ability to create custom surface textures, which are deployed to dramatically control wettability and create energy-gradient surfaces. These techniques allow fouling/scaling nucleation to be spatially controlled, such that as the fouling layer grows it will be broken and removed by flow forces.

42 ENGINEERING↗
NASA NTRS1993

An evaluation of three spatial differencing schemes for the discrete ordinates method in participating media

Three popular spatial differencing schemes for the discrete ordinates method are examined for two-dimensional Cartesian coordinates system. These are a positive, the step, and the diamond schemes. Contrary to the common belief that negative intensities will not occur when fine spatial discretization is used with the diamond scheme, under certain conditions, the diamond scheme will produce negative intensities irrespective of the number of control volums employed. The positive scheme can produce physically unrealistic trends. The diamond and positive schemes are also capable of producing physically unrealistic overshoots. In absorbing-emitting or absorbing-emitting-scattering media, grid refinement can result in negative intensities when the diamond or positive scheme is used.

Chai, John C.↗
NASA NTRS1993

Temperature dependence of protein solubility-determination, application to crystallization, and growth kinetics studies

A scintillation method was developed for determinations of the temperature dependence of the solubility, and of nucleation induction times of proteins, in 50-100 mu(l) volumes of solution. Solubility data for lysozyme and horse serum albumin were obtained for various combinations of pH and precipitant concentrations. These data and the nucleation induction information were used for dynamic crystallization control, that is, for the controlled separation of nucleation and growth stages. Individual lysozyme and horse serum albumin crystals were grown in 15-20 mu(l) solution volumes contained in x-ray capillaries. The morphology and kinetics of the growth and dissolution of lysozyme in aqueous solutions with 2.5 percent NaCl and at pH = 4.5 was studied in situ with a depth resolution of 300 A (4 unit cells) by high resolution optical microscopy and digital image processing. The bulk super- or under saturation, sigma, of the solution inside a closed growth cell was controlled by temperature. The growth habit was bound by (110) and (101) faces that grew through layer spreading, although with different growth rate dependencies on supersaturation/temperature. At sigma less than 10 (obtained at higher temperatures) growth was purely kinetic ally controlled, with impurity effects (macrostep formation and kinetic hindrance) becoming significant for sigma less than 2. At sigma greater than 10 (lower temperatures), anisotropies in the interfacial kinetics were more pronounced, with interfacial kinetics and bulk transport becoming equally important to the growth morphology. Growth rates were growth history dependent. The formation of striations (layers of irregularly incorporated solution) was unambiguously correlated with growth temperature variations. Etching exposed dislocations and various high-index faces whose growth morphologies were studied during return to the steady state growth form. Growth steps were observed to originate from two-dimensional nuclei or from outcrops of growth striations, and from dislocations that preferentially formed in growth sector boundaries.

Rosenberger, Franz↗
NASA NTRS2003

Tetrahedron Formation Control

Spacecraft flying in tetrahedron formations are excellent instrument platforms for electromagnetic and plasma studies. A minimum of four spacecraft - to establish a volume - is required to study some of the key regions of a planetary magnetic field. The usefulness of the measurements recorded is strongly affected by the tetrahedron orbital evolution. This paper considers the preliminary development of a general optimization procedure for tetrahedron formation control. The maneuvers are assumed to be impulsive and a multi-stage optimization method is employed. The stages include targeting to a fixed tetrahedron orientation, rotating and translating the tetrahedron and/or varying the initial and final times. The number of impulsive maneuvers citn also be varied. As the impulse locations and times change, new arcs are computed using a differential corrections scheme that varies the impulse magnitudes and directions. The result is a continuous trajectory with velocity discontinuities. The velocity discontinuities are then used to formulate the cost function. Direct optimization techniques are employed. The procedure is applied to the Magnetospheric Multiscale Mission (MMS) to compute preliminary formation control fuel requirements.

Guzman, Jose J.↗
NASA NTRS2014

Stability Analysis of ISS Medications

It is known that medications degrade over time, and that extreme storage conditions will hasten their degradation. The temperature and humidity conditions of the ISS have been shown to be within the ideal ranges for medication storage, but the effects of other environmental factors, like elevated exposure to radiation, have not yet been evaluated. Current operational procedures ensure that ISS medications are re-stocked before expiration, but this may not be possible on long duration exploration missions. For this reason, medications that have experienced long duration storage on the ISS were returned to JSC for analysis to determine any unusual effects of aging in the low- Earth orbit environment. METHODS Medications were obtained by the JSC Pharmacy from commercial distributors and were re-packaged by JSC pharmacists to conserve up mass and volume. All medication doses were part of the ISS crew medical kit and were transported to the International Space Station (ISS) via NASA's Shuttle Transportation System (Space Shuttle). After 568 days of storage, the medications were removed from the supply chain and returned to Earth on a Dragon (SpaceX) capsule. Upon return to Earth, medications were transferred to temperature and humidity controlled environmental chambers until analysis. Nine medications were chosen on the basis of their availability for study. The medications included several of the most heavily used by US crewmembers: 2 sleep aids, 2 antihistamines/decongestants, 3 pain relievers, an antidiarrheal and an alertness medication. Each medication was available at a single time point; analysis of the same medication at multiple time points was not possible. Because the samples examined in this study were obtained opportunistically from medical supplies, there were no control samples available (i.e. samples aged for a similar period of time on the ground); a significant limitation of this study. Medications were analyzed using the HPLC/MS methods described in the United States Pharmacopeia (USP) to measure the amount of intact active ingredient, identify degradation products and measure their amounts. Some analyses were conducted by an independent analytical laboratory, but certain (Schedule) medications could not be shipped to their facility and were analyzed at JSC. RESULTS Nine medications were analyzed with respect to active pharmaceutical ingredient (API) and degradant amounts. Results were compared to the USP requirements for API and degradants/impurities content for every FDA-approved medication. One medication met USP requirements at 5 months after its expiration date. Four of the nine (44% of those tested) medications tested met USP requirements up to 8 months post-expiration. Another 3 medications (33% of those tested) met USP guidelines 2-3 months before expiration. One medication, a compound classed by the FDA as a dietary supplement and sometimes used as a sleep aid, failed to meet USP requirements at 11 months post-expiration. CONCLUSION Analysis of each medication at a single time point provides limited information on the stability of a medication stored in particular conditions; it is not possible to predict how long a medication may be safe and effective from these data. Notwithstanding, five of the nine medications tested (56%) met USP requirements for API and degradants/impurities at least 5 months past expiration dates. The single compound that failed to meet USP requirements is not regulated as strictly as prescription medications are during manufacture; it is unknown if this medication would have met the requirements prior to flight. Notably, it was the furthest beyond its expiration date. Only more comprehensive analysis of flight-aged samples compared to appropriate ground controls will permit determination of spaceflight effects on medication stability.

Wotring, V. E.↗
DOE OSTI2024

Generation of Large-Volume High-Pressure Plasma by Spatio-Temporal Control of Space Charge

Due to the fundamental limitations of scaling, very little progress has been made towards achieving a large-volume dense non-equilibrium plasma near atmospheric pressures. Commercially available state-of-the art glow-like atmospheric plasma devices for industrial applications are either narrow tubes or wide slits with narrow openings. Often multiple sources are put together in various configurations to process larger surfaces. The traditional methods of exciting electrodes create spatially fixed electric fields. As a result, the space charge at atmospheric pressure tends to be spatially confined resulting in non-uniformity which eventually leads to instability as the discharge is scaled. Theoretical work done under this project showed that it is possible to generate a spatially rotating electric field by exciting a set of electrodes with phase staggered sinusoidal waveforms. The modeling and simulations were done using plasma fluid models. It was shown that such a field can produce a uniform plasma. At the conclusion of the project, experimental proof of concept with an eight-electrode system in various gases (Air, Helium and Argon) was demonstrated. Power measurements and spectral investigation show that the concept can be used to generate a uniform stable plasma. This plasma source has the potential of opening new applications of nonthermal plasma including combustion of carbon-free fuel. The current limitation of the proposed method in scaling to higher volume and pressure is the need for multiple high voltage amplifiers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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↗
DOE OSTI2015

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↗
DOE OSTI2015

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↗
DOE OSTI2019

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↗
DOE OSTI2018

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↗
NASA NTRS1994

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.↗
NASA NTRSDate not supplied

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↗
DOE OSTI2020

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↗
DOE OSTI2019

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↗
DOE OSTI2020

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↗