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 361 records · Page 20

Control rod treatment for FEM based radiation transport methods

This paper presents a novel control rod treatment for finite element method (FEM) based reactor physics simulations. The new method falls into the class of homogenization methods which treat material interfaces within mesh elements by defining a representative, homogenized nuclear cross section set. In contrast to existing homogenization based cusping treatments that attempt to approximate flux-volume weighting of the cross sections, we propose to exactly integrate terms appearing in the variational form of the transport or diffusion equation over mixed material elements. This goal is achieved by defining smoothly varying cross sections in elements that straddle material boundaries. It is shown that the smoothly varying cross sections can be obtained by projecting the piecewise set of cross sections onto the set of Legendre polynomials up to a finite order. Implementation of smoothly varying cross sections is supported within the Rattlesnake radiation multiphysics code that is currently under development at the Idaho National Laboratory. Homogenization methods traditionally suffer from cusping effects that originate from overestimation of absorption cross sections in strongly neutron-absorbing materials as found in control rods. We demonstrate that the presented method is (1) highly effective in removing cusping effects and (2) provides accurate results based on the C5G7-TD benchmark, a time-dependent version of the C5G7 benchmark.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-throughput continuous flow synthesis of nickel nanoparticles for the catalytic hydrodeoxygenation of guaiacol

The translation of batch chemistries to high-throughput continuous flow methods dresses scaling, automation, and reproducibility concerns associated with the implementation of colloidally prepared nanoparticle (NP) catalysts for industrial catalytic processes. Nickel NPs were synthesized by the high-temperature amine reduction of a Ni2+ precursor using a continuous millifluidic (mF) flow method, achieving yields greater than 60%. The resulting Ni NP catalysts were compared against catalysts prepared in a batch reaction under conditions analogous to the continuous flow conditions with respect to total reaction volume, time, and temperature and by traditional incipient wetness (IW) impregnation for the hydrodeoxygenation (HDO) of guaiacol under ex situ catalytic fast pyrolysis conditions. Compared to the IW method, the colloidally prepared NPs displayed increased morphological control and narrowed size distributions, and the NPs prepared by both methods showed similar size, shape, and crystallinity. The Ni NP catalyst synthesized by the continuous flow method exhibited similar H-adsorption site densities, site-time yields, and selectivities towards deoxygenated products as compared to the analogous batch reaction, and outperformed the IW catalyst with respect to higher selectivity to lower oxygen content products and a 6.9-fold slower deactivation rate. These results demonstrate the utility of synthesizing colloidal Ni NP catalysts using continuous flow methods while maintaining the catalytic properties displayed by the batch equivalent. Finally, this methodology can be extended to other catalytically relevant base metals for the high-throughput synthesis of metal NPs for the catalytic production of biofuels.

09 BIOMASS FUELS↗

Laser 3D printing of highly compacted protonic ceramic electrolyzer stack

Solid oxide electrolysis cells (SOECs) for H 2 production is a core technology for H2@scale. However, its conventional manufacturing methods adapted from the manufacture of solid oxide fuel cells (SOFC) need high cost, especially for small-volume production. The emerging laser 3D printing (L3DP) technology with computer-aided 3D printing and computer-controlled laser processing can fulfill layer-by-layer digital shaping and rapid in-situ consolidating feedstock into complicated geometries. L3DP, a promising additive manufacturing (AM) technology, has achieved significant success in manufacturing plastic and metal parts, which is currently attracting considerable attention for the cost-effective, rapid, and flexible manufacturing of heterogeneous multilayered ceramic devices (e.g., SOECs, SOFCs, and solid-state batteries). It is believed that the L3DP can integrate the advantages of selective consolidation of heterogeneous layers, accurate control of layer microstructures, high processing heat efficiency, high processing speed, high stack compactness, high stack design flexibility, and low stack sealing area for cost-effective, rapid, and flexible manufacturing of SOECs to meet DOE’s electrolyzer target.

08 HYDROGEN↗

Piezoelectric pushers for active vibration control of rotating machinery

The active control of rotordynamic vibrations and stability by magnetic bearings and electromagnetic shakers was discussed extensively in the literature. These devices, though effective, are usually large in volume and add significant weight to the stator. The use of piezoelectric pushers may provide similar degrees of effectiveness in light, compact packages. Analyses are contained which extend quadratic regulator, pole placement and derivative feedback control methods to the prescribed displacement character of piezoelectric pushers. The structural stiffness of the pusher is also included in the theory. Tests are currently being conducted at NASA Lewis Research Center with piezoelectric pusher-based active vibration control. Results performed on the NASA test rig as preliminary verification of the related theory are presented.

Palazzolo, A. B.↗

Computation of solid/liquid phase change including free convection - Comparison with data

A computational model is presented for solid/liquid phase-change energy transport including free convection fluid flow in the liquid phase. The computational model considers the velocity components of all nonliquid control volumes to be zero but fully solves the coupled mass-momentum problem within the liquid. The thermal energy model includes the entire domain and employs an enthalpy-like model and a recently developed method for handling the phase-change interface nonlinearity. Convergence studies are performed and comparisons made with experimental data for two different problems. Grid independence is achieved, and the comparison with experimental data indicates excellent quantitative prediction of the melt fraction evolution. Qualitative data are also provided as velocity vector and isotherm plots. The computational costs incurred are quite low by comparison with other models.

Schneider, G. E.↗

Prediction of cascade performance using an incompressible Navier-Stokes technique

A fully elliptic, control volume solution of the two-dimensional incompressible Navier-Stokes equations for the prediction of cascade performance over a wide incidence range is presented. The numerical technique is based on a new pressure substitution method. A Poisson equation is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. The analysis of a double circular arc compressor cascade is presented, and the results are compared with the available experimental data at various incidence angles. Good agreement is obtained for the blade pressure distribution, boundary layer and wake profiles, skin friction coefficient, losses and outlet angles. Turbulence effects are simulated by the Low-Reynolds-Number version of the k-epsilon turbulence model.

Hobson, G. V.↗

Numerical Formulation for the Prediction of Solid/Liquid Change of a Binary Alloy

A computational model is presented for the prediction of solid/liquid phase change energy transport including the influence of free convection fluid flow in the liquid phase region. The computational model considers the velocity components of all non-liquid phase change material control volumes to be zero but fully solves the coupled mass-momentum problem within the liquid region. The thermal energy model includes the entire domain and uses an enthalpy like model and a recently developed method for handling the phase change interface nonlinearity. Convergence studies are performed and comparisons made with experimental data for two different problem specifications. The convergence studies indicate that grid independence was achieved and the comparison with experimental data indicates excellent quantitative prediction of the melt fraction evolution. Qualitative data is also provided in the form of velocity vector diagrams and isotherm plots for selected times in the evolution of both problems. The computational costs incurred are quite low by comparison with previous efforts on solving these problems.

Schneider, G. E.↗

Recent refinements and increased capabilities in balloon vertical performance analysis

The NASA Thermal Trajectory analysis model for predicting the vertical performance of balloons is described in terms of two critical refinements and recent results. The model employs heat-balance equations and a vertical equation of motion, and revised methods are employed to predict maximum balloon volume, balloon shape, and the continuous computation of gas loss. Also incorporated are modeling techniques for assessing the solar radiation input that can affect the load tapes, and autoballast control algorithms are employed. The model permits the evaluation of the presence of cap layers as well as the effects of small leaks and programmed venting. The revised model is shown to give good predictions of balloon performance during flight in terms of altitude changes, descent rates, gas-mass flow, and leak potential. The model is found to be more versatile due to autoballasting and leak analysis, and good prediction are possible for the float phase and descent of current balloon shapes.

Conrad, G. R.↗

{sup 18}F-Fluorodeoxyglucose/Positron Emission Tomography Predicts Patterns of Failure After Definitive Chemoradiation Therapy for Locally Advanced Non-Small Cell Lung Cancer

Background: We previously reported that pretreatment positron emission tomography (PET) identifies lesions at high risk for progression after concurrent chemoradiation therapy (CRT) for locally advanced non-small cell lung cancer (NSCLC). Here we validate those findings and generate tumor control probability (TCP) models. Methods: We identified patients treated with definitive, concurrent CRT for locally advanced NSCLC who underwent staging {sup 18}F-fluorodeoxyglucose/PET/computed tomography. Visible hypermetabolic lesions (primary tumors and lymph nodes) were delineated on each patient's pretreatment PET scan. Posttreatment imaging was reviewed to identify locations of disease progression. Competing risks analyses were performed to examine metabolic tumor volume (MTV) and radiation therapy dose as predictors of local disease progression. TCP modeling was performed to describe the likelihood of local disease control as a function of lesion size. Results: Eighty-nine patients with 259 hypermetabolic lesions (83 primary tumors and 176 regional lymph nodes) met the inclusion criteria. Twenty-eight patients were included in our previous report, and the remaining 61 constituted our validation cohort. The median follow-up time was 22.7 months for living patients. In 20 patients, the first site of progression was a primary tumor or lymph node treated with radiation therapy. The median time to progression for those patients was 11.5 months. Data from our validation cohort confirmed that lesion MTV predicts local progression, with a 30-month cumulative incidence rate of 23% for lesions above 25 cc compared with 4% for lesions below 25 cc (P=.008). We found no evidence that radiation therapy dose was associated with local progression risk. TCP modeling yielded predicted 30-month local control rates of 98% for a 1-cc lesion, 94% for a 10-cc lesion, and 74% for a 50-cc lesion. Conclusion: Pretreatment FDG-PET identifies lesions at risk for progression after CRT for locally advanced NSCLC. Strategies to improve local control should be tested on high-risk lesions, and treatment deintensification for low-risk lesions should be explored.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

WE-AB-BRB-01: Development of a Probe-Format Graphite Calorimeter for Practical Clinical Dosimetry: Numerical Design Optimization, Prototyping, and Experimental Proof-Of-Concept

Purpose: In this work, the feasibility of performing absolute dose to water measurements using a constant temperature graphite probe calorimeter (GPC) in a clinical environment is established. Methods: A numerical design optimization study was conducted by simulating the heat transfer in the GPC resulting from irradiation using a finite element method software package. The choice of device shape, dimensions, and materials was made to minimize the heat loss in the sensitive volume of the GPC. The resulting design, which incorporates a novel aerogel-based thermal insulator, and 15 temperature sensitive resistors capable of both Joule heating and measuring temperature, was constructed in house. A software based process controller was developed to stabilize the temperatures of the GPC’s constituent graphite components to within a few 10’s of µK. This control system enables the GPC to operate in either the quasi-adiabatic or isothermal mode, two well-known, and independent calorimetry techniques. Absorbed dose to water measurements were made using these two methods under standard conditions in a 6 MV 1000 MU/min photon beam and subsequently compared against TG-51 derived values. Results: Compared to an expected dose to water of 76.9 cGy/100 MU, the average GPC-measured doses were 76.5 ± 0.5 and 76.9 ± 0.5 cGy/100 MU for the adiabatic and isothermal modes, respectively. The Monte Carlo calculated graphite to water dose conversion was 1.013, and the adiabatic heat loss correction was 1.003. With an overall uncertainty of about 1%, the most significant contributions were the specific heat capacity (type B, 0.8%) and the repeatability (type A, 0.6%). Conclusion: While the quasi-adiabatic mode of operation had been validated in previous work, this is the first time that the GPC has been successfully used isothermally. This proof-of-concept will serve as the basis for further study into the GPC’s application to small fields and MRI-linac dosimetry. This work has been supported in part by the CREATE Medical Physics Research Training Network of the Natural Sciences and Engineering Research Council (NSERC) grant 432290, NSERC grants RGPIN 298191 & 435608-13, Canadian Institutes of Health Research doctoral scholarship GSD-121793. This work has also been supported by Sun Nuclear Corporation.

60 APPLIED LIFE SCIENCES↗

A program for calculating load coefficient matrices utilizing the force summation method, L218 (LOADS). Volume 2: Supplemental system design and maintenance document

The LOADS computer program L218 which calculates dynamic load coefficient matrices utilizing the force summation method is described. The load equations are derived for a flight vehicle in straight and level flight and excited by gusts and/or control motions. In addition, sensor equations are calculated for use with an active control system. The load coefficient matrices are calculated for the following types of loads: (1) translational and rotational accelerations, velocities, and displacements; (2) panel aerodynamic forces; (3) net panel forces; and (4) shears, bending moments, and torsions.

Anderson, L. R.↗

A program for calculating load coefficient matrices utilizing the force summation method, L218 (LOADS). Volume 1: Engineering and usage

The LOADS program L218, a digital computer program that calculates dynamic load coefficient matrices utilizing the force summation method, is described. The load equations are derived for a flight vehicle in straight and level flight and excited by gusts and/or control motions. In addition, sensor equations are calculated for use with an active control system. The load coefficient matrices are calculated for the following types of loads: translational and rotational accelerations, velocities, and displacements; panel aerodynamic forces; net panel forces; shears and moments. Program usage and a brief description of the analysis used are presented. A description of the design and structure of the program to aid those who will maintain and/or modify the program in the future is included.

Miller, R. D.↗

SU-E-T-574: Novel Chance-Constrained Optimization in Intensity-Modulated Proton Therapy Planning to Account for Range and Patient Setup Uncertainties

Purpose: We propose to apply a probabilistic framework, namely chanceconstrained optimization, in the intensity-modulated proton therapy (IMPT) planning subject to range and patient setup uncertainties. The purpose is to hedge against the influence of uncertainties and improve robustness of treatment plans. Methods: IMPT plans were generated for a typical prostate patient. Nine dose distributions are computed — the nominal one and one each for ±5mm setup uncertainties along three cardinal axes and for ±3.5% range uncertainty. These nine dose distributions are supplied to the solver CPLEX as chance constraints to explicitly control plan robustness under these representative uncertainty scenarios with certain probability. This probability is determined by the tolerance level. We make the chance-constrained model tractable by converting it to a mixed integer optimization problem. The quality of plans derived from this method is evaluated using dose-volume histogram (DVH) indices such as tumor dose homogeneity (D5% – D95%) and coverage (D95%) and normal tissue sparing like V70 of rectum, V65, and V40 of bladder. We also compare the results from this novel method with the conventional PTV-based method to further demonstrate its effectiveness Results: Our model can yield clinically acceptable plans within 50 seconds. The chance-constrained optimization produces IMPT plans with comparable target coverage, better target dose homogeneity, and better normal tissue sparing compared to the PTV-based optimization [D95% CTV: 67.9 vs 68.7 (Gy), D5% – D95% CTV: 11.9 vs 18 (Gy), V70 rectum: 0.0 % vs 0.33%, V65 bladder: 2.17% vs 9.33%, V40 bladder: 8.83% vs 21.83%]. It also simultaneously makes the plan more robust [Width of DVH band at D50%: 2.0 vs 10.0 (Gy)]. The tolerance level may be varied to control the tradeoff between plan robustness and quality. Conclusion: The chance-constrained optimization generates superior IMPT plan compared to the PTV-based optimization with explicit control of plan robustness. NIH/NCI K25CA168984, Eagles Cancer Research Career Development, The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, Mayo ASU Seed Grant, and The Kemper Marley Foundation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A Versatile Microfluidic Device for Automating Synthetic Biology

New microbes are being engineered that contain the genetic circuitry, metabolic pathways, and other cellular functions required for a wide range of applications such as producing biofuels, biobased chemicals, and pharmaceuticals. Although currently available tools are useful in improving the synthetic biology process, further improvements in physical automation would help to lower the barrier of entry into this field. We present an innovative microfluidic platform for assembling DNA fragments with 10× lower volumes (compared to that of current microfluidic platforms) and with integrated region-specific temperature control and on-chip transformation. Integration of these steps minimizes the loss of reagents and products compared to that with conventional methods, which require multiple pipetting steps. For assembling DNA fragments, we implemented three commonly used DNA assembly protocols on our microfluidic device: Golden Gate assembly, Gibson assembly, and yeast assembly (i.e., TAR cloning, DNA Assembler). Here, we demonstrate the utility of these methods by assembling two combinatorial libraries of 16 plasmids each. Each DNA plasmid is transformed into Escherichia coli or Saccharomyces cerevisiae using on-chip electroporation and further sequenced to verify the assembly. Lastly, we anticipate that this platform will enable new research that can integrate this automated microfluidic platform to generate large combinatorial libraries of plasmids and will help to expedite the overall synthetic biology process.

59 BASIC BIOLOGICAL SCIENCES↗

Non-Aqueous Electrochemical Fluorination of Used Nuclear Fuel as an Advanced Separation Process

Development of efficient and environmentally benign methods to reprocess used nuclear fuel (UNF) will enable technologies for a nuclear renaissance. Electrochemical and fluorination methods for reprocessing UNF have been proposed, but combinations of electrochemical and fluorination methods have not been investigated. Electrochemical fluorination may reduce waste volumes compared to the main large-scale aqueous methods for processing used nuclear fuel. A non-aqueous electrochemical fluorination reprocessing method has been developed and demonstrated that enables gas phase uranium recovery while allowing for potential control of the reaction using a single process. Thermodynamic modeling showed a galvanic reaction between UNF and a fluorinating agent, such as NF 3 , in a molten fluoride electrolyte was possible and could selectively fluorinate U to UF 6 . To verify the model results, a reactor system for electrochemical fluorination was constructed along with vessels for product collection and analysis. Initial trials of the electrochemical reaction were performed and characterized by several electrochemical methods. The electrochemical oxidation of U using NF 3 as an oxidizing agent was demonstrated and UF 6 as a reaction product was detected.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Subplane Decusping for BWRs in MPACT

Control blade cusping can introduce significant error in boiling water reactor (BWR) calculations with MPACT when blade tips fall partway within an axial method of characteristics (MOC) plane, requiring homogenization of controlled and uncontrolled regions. This work implements subplane decusping for BWRs in MPACT by enabling BWR-compatible subplane coarse mesh finite difference (CMFD) and extending the existing decusping framework to represent between-assembly control blades that insert from the bottom of the core. The method resolves axial heterogeneity on a refined subplane mesh in the low-order solve and uses the resulting subplane flux shape to form flux-volume homogenized transport cross sections for the partially rodded MOC plane. The capability is evaluated using a single physics General Electric (GE)-14 assembly and a multiphysics Peach Bottom Unit 2 Type 1 assembly (PB2T1A) with thermal hydraulic feedback. In both cases, coarse axial meshes with and without subplane decusping are compared against fine mesh reference solutions over the full range of blade withdrawal positions. Subplane decusping reduces maximum/average eigenvalue errors from 3,369/220 pcm to 303/36 pcm for GE-14 and from 6,689/585 pcm to 1,582/109 pcm for PB2T1A. Additionally, it reduces pin power root mean square errors from 5.2%/0.8% to 1.8%/0.2% for GE-14 and from 6.9%/0.6% to 2.3%/0.1% for PB2T1A. These results demonstrate an effective, practical correction for BWR blade cusping in VERA-MPACT.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗

Communication: Control of chemical reactions using electric field gradients

We examine theoretically a new idea for spatial and temporal control of chemical reactions. When chemical reactions take place in a mixture of solvents, an external electric field can alter the local mixture composition, thereby accelerating or decelerating the rate of reaction. The spatial distribution of electric field strength can be non-trivial and depends on the arrangement of the electrodes producing it. In the absence of electric field, the mixture is homogeneous and the reaction takes place uniformly in the reactor volume. When an electric field is applied, the solvents separate and the reactants are concentrated in the same phase or separate to different phases, depending on their relative miscibility in the solvents, and this can have a large effect on the kinetics of the reaction. This method could provide an alternative way to control runaway reactions and to increase the reaction rate without using catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗