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 433 records · Page 24

Issues in measure-preserving three dimensional flow integrators: Self-adjointness, reversibility, and non-uniform time stepping

Properties of integration schemes for solenoidal fields in three dimensions are studied, with a focus on integrating magnetic field lines in a plasma using adaptive time stepping. It is shown that implicit midpoint (IM) and a scheme we call three-dimensional leapfrog (LF) can do a good job (in the sense of preserving KAM tori) of integrating fields that are reversible, or (for LF) have a “special divergence-free” (SDF) property. We review the notion of a self-adjoint scheme, showing that such schemes are at least second order accurate and can always be formed by composing an arbitrary scheme with its adjoint. We also review the concept of reversibility, showing that a reversible but not exactly volume-preserving scheme can lead to a fractal invariant measure in a chaotic region, although this property may not often be observable. We also show numerical results indicating that the IM and LF schemes can fail to preserve KAM tori when the reversibility property (and the SDF property for LF) of the field is broken. We discuss extensions to measure preserving flows, the integration of magnetic field lines in a plasma and the integration of rays for several plasma waves. The main new result of this paper relates to non-uniform time stepping for volume-preserving flows. We investigate two potential schemes, both based on the general method of Feng and Shang [Numer. Math. 71, 451 (1995)], in which the flow is integrated in split time steps, each Hamiltonian in two dimensions. The first scheme is an extension of the method of extended phase space, a well-proven method of symplectic integration with non-uniform time steps. This method is found not to work, and an explanation is given. The second method investigated is a method based on transformation to canonical variables for the two split-step Hamiltonian systems. This method, which is related to the method of non-canonical generating functions of Richardson and Finn [Plasma Phys. Controlled Fusion 54, 014004 (2012)], appears to work very well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Comparison of Thruster Implementation Strategies for a Deep Space Nanosatellite

Attitude control for very small spacecraft, commonly referred to as nanosatellites or CubeSats, has traditionally been carried out using reaction wheels and magnetic torquers as the primary actuators. However, as these spacecraft begin to be considered for a broader range of scientific applications, including some beyond low Earth orbit, it has become necessary to also consider thruster systems for actuation. In recent years a number of thruster designs that conform to the mass, volume, and power constraints of nanosatellites have become commercially available, including cold gas systems, pulsed plasma thrusters (PPTs), and micro-electrospray propulsion (MEP) systems. The challenge now facing the nanosatellite community is to determine which thruster solutions are appropriate for a particular application, and what the best method of control might be. This paper will compare the implementation of a cold gas system with that of an MEP or PPT system for an upcoming nanosatellite mission using a previously reported saturation-restricted control law. Results are presented for this controller both with and without a fuel-optimal thruster allocation scheme, and an assessment on incorporating these technologies in an upcoming NASA mission is offered.

Nehrenz, Matt↗

High-Dose Hypofractionated Radiation Therapy for Noncompressive Vertebral Metastases in Combination With Zoledronate: A Phase 1 Study

Introduction: Hypofractionated stereotactic radiation therapy (HSRT) for vertebral metastases gives good results in terms of local control but increases the risk of fracture in the treated volume. Preclinical and clinical studies have shown that zoledronate not only reduces the risk of fracture and stimulates osteoclastic remodeling but also increases the immune response and radiosensitivity. This study aimed to evaluate the tolerability and effectiveness of zoledronate in association with radiation therapy. Patients and Methods: We conducted a multicenter phase 1 study that combined HSRT (3 × 9 Gy) and zoledronate in patients with vertebral metastasis ( (NCT01219790)). The principal objective was the absence of spinal cord adverse reactions at 1 year. The secondary objectives were acute tolerability, the presentation of a bone event, local tumor control, pain control, progression-free survival, and overall survival. Results: Thirty patients (25 male, 5 female), median age 66 years, who were followed up for a median period of 19.2 months, received treatment for 49 vertebral metastases. A grade 3 acute mucosal adverse event occurred in 1 patient during the treatment and in 2 more at 1 month. No late neurologic adverse events were reported at 1 year. The mean pain scores diminished significantly at 1 month (1.35; P=.0125) and 3 months (0.77; P<.0001) compared with pain scores at study entry (2.49). Vertebral collapse in the irradiated zone occurred in 1 (2%) treated vertebra. Control of local disease was achieved in 94% of irradiated patients (3 local recurrences). Conclusion: The combination of zoledronate and HSRT in the treatment of vertebral metastasis is well tolerated and seems to reduce the rate of vertebral collapse, effectively relieve pain, and achieve good local tumor control with no late neurologic adverse effects.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

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↗

Dose-Escalated Intensity-Modulated Radiotherapy Is Feasible and May Improve Locoregional Control and Laryngeal Preservation in Laryngo-Hypopharyngeal Cancers

Purpose: To determine the safety and outcomes of induction chemotherapy followed by dose-escalated intensity-modulated radiotherapy (IMRT) with concomitant chemotherapy in locally advanced squamous cell cancer of the larynx and hypopharynx (LA-SCCL/H). Methods and Materials: A sequential cohort Phase I/II trial design was used to evaluate moderate acceleration and dose escalation. Patients with LA-SCCL/H received IMRT at two dose levels (DL): DL1, 63 Gy/28 fractions (Fx) to planning target volume 1 (PTV1) and 51.8 Gy/28 Fx to PTV2; DL2, 67.2 Gy/28 Fx and 56 Gy/28 Fx to PTV1 and PTV2, respectively. Patients received induction cisplatin/5-fluorouracil and concomitant cisplatin. Acute and late toxicities and tumor control rates were recorded. Results: Between September 2002 and January 2008, 60 patients (29 DL1, 31 DL2) with Stage III (41% DL1, 52% DL2) and Stage IV (52% DL1, 48% DL2) disease were recruited. Median (range) follow-up for DL1 was 51.2 (12.1-77.3) months and for DL2 was 36.2 (4.2-63.3) months. Acute Grade 3 (G3) dysphagia was higher in DL2 (87% DL2 vs. 59% DL1), but other toxicities were equivalent. One patient in DL1 required dilatation of a pharyngeal stricture (G3 dysphagia). In DL2, 2 patients developed benign pharyngeal strictures at 1 year. One underwent a laryngo-pharyngectomy and the other a dilatation. No other G3/G4 toxicities were reported. Overall complete response was 79% (DL1) and 84% (DL2). Two-year locoregional progression-free survival rates were 64.2% (95% confidence interval, 43.5-78.9%) in DL1 and 78.4% (58.1-89.7%) in DL2. Two-year laryngeal preservation rates were 88.7% (68.5-96.3%) in DL1 and 96.4% (77.7-99.5%) in DL2. Conclusions: At a mean follow-up of 36 months, dose-escalated chemotherapy-IMRT at DL2 has so far been safe to deliver. In this study, DL2 delivered high rates of locoregional control, progression-free survival, and organ preservation and has been selected as the experimental arm in a Cancer Research UK Phase III study.

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↗