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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Effects of size reduction on microstructural, optical, vibrational, magnetic and photocatalytic properties of ZnO nanocrystals

Highlights: • ZnO nanocrystals of different sizes 4–19 nm, thus different S/V ratios synthesized. • Lattice strain and lattice volume increase with size reduction. • Raman spectra show enhanced surface effects with phonon mode disappearance. • ESR shows exchange interaction; VSM shows surface originated ferromagnetism. • Photocatalysis is less affected with size reduction and surface enhancement. - Abstract: ZnO nanocrystals were synthesized in the size range 4 nm–19 nm by chemical co-precipitation method and extensively characterized by various techniques. Size control was achieved by changing the concentration of ME that acted as the capping agent. Lattice strain and lattice volume were found to increase with the decrease in crystallite size that could be attributed to an enhanced influence of surface atoms. Raman spectroscopy provided evidence of enhancement of surface influence and reduction of crystallinity with the observation of decrease in E{sub 2}{sup high} phonon mode intensity with a decrease in crystallite size. Electron spin resonance spectra revealed the presence of exchange interaction as well as presence of lattice disorder. Vibrating sample magnetometry results showed that M-H curves exhibit hysteresis loops which confirmed the ferromagnetic nature of the nanocrystals and it could be attributed to the interaction between localized electron spins resulting from oxygen vacancies at the surface. Blue shift of the UV–Visible absorption spectra was observed with decrease in crystallite size, which could be attributed to quantum confinement effect. Room temperature photoluminescence spectra were recorded under UV light excitation, which showed emission of different colours (light violet, violet, blue and green) in the visible region, attributable to the presence of various types of defects such as zinc interstitial, zinc vacancy and oxygen vacancy in the ZnO nanocrystals. Photodegradation of aqueous methylene blue by the ZnO nanocrystals was performed under simulated sunlight irradiation. There was a very small increase in efficiency of photodegradation with size reduction, which showed that there was no significant change in surface defects with size reduction in the nanocrystals despite the surface enhancement. A mechanism for photocatalysis involving defect states has been proposed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2013

High-Dose and Extended-Field Intensity Modulated Radiation Therapy for Early-Stage NK/T-Cell Lymphoma of Waldeyer's Ring: Dosimetric Analysis and Clinical Outcome

Purpose: To assess the dosimetric benefit, treatment outcome, and toxicity of high-dose and extended-field intensity modulated radiation therapy (IMRT) in patients with early-stage NK/T-cell lymphoma of Waldeyer's ring (WR-NKTCL). Methods and Materials: Thirty patients with early-stage WR-NKTCL who received extended-field IMRT were retrospectively reviewed. The prescribed dose was 50 Gy to the primary involved regions and positive cervical lymph nodes (planning target volume requiring radical irradiation [PTV{sub 50}]) and 40 Gy to the negative cervical nodes (PTV{sub 40}). Dosimetric parameters for the target volume and critical normal structures were evaluated. Locoregional control (LRC), overall survival (OS), and progression-free survival (PFS) were calculated using the Kaplan-Meier method. Results: The median mean doses to the PTV{sub 50} and PTV{sub 40} were 53.2 Gy and 43.0 Gy, respectively. Only 1.4% of the PTV{sub 50} and 0.9% of the PTV{sub 40} received less than 95% of the prescribed dose, indicating excellent target coverage. The average mean doses to the left and right parotid glands were 27.7 and 28.4 Gy, respectively. The 2-year OS, PFS, and LRC rates were 71.2%, 57.4%, and 87.8%. Most acute toxicities were grade 1 to 2, except for grade ≥3 dysphagia and mucositis. The most common late toxicity was grade 1-2 xerostomia, and no patient developed any ≥grade 3 late toxicities. A correlation between the mean dose to the parotid glands and the degree of late xerostomia was observed. Conclusions: IMRT achieves excellent target coverage and dose conformity, as well as favorable survival and locoregional control rates with acceptable toxicities in patients with WR-NKTCL.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2013

Ordered nanoporous carbon for increasing CO{sub 2} capture

Ordered nanoporous carbons (ONCs) were prepared using a soft-templating method. The prepared ONCs materials were subjected to a controlled carbonization temperature over the temperature range, 700-1000 Degree-Sign C, to increase the specific surface area and total pore volume of ordered nanoporous carbon followed by carbonization of the phenolic resin. ONCs materials synthesized at various carbonization temperatures were used as adsorbents to improve the CO{sub 2} adsorption efficiency. The surface properties of the ONCs materials were examined by X-ray photoelectron spectroscopy. The structural properties of the ONCs materials were analyzed by X-ray diffraction. The textural properties of the ONCs materials were examined using the N{sub 2}/77 K adsorption isotherms according to the Brunauer-Emmett-Teller equation. The CO{sub 2} adsorption capacity was measured by CO{sub 2} isothermal adsorption at 298 K/30 bar and 298 K/1 bar. The carbonization temperature was found to have a major effect on the CO{sub 2} adsorption capacity, resulting from the specific surface area and total pore volumes of the ONCs materials. - Graphical abstract: This schematic diagram described synthesis of ONCs. Highlights: Black-Right-Pointing-Pointer ONCs materials can be prepared readily using the direct-triblock-copolymer-templating method. Black-Right-Pointing-Pointer The distributions show that prominent development can be observed around the micro-pore region. Black-Right-Pointing-Pointer The soft-templating method provides opportunities for controlling the pore structure of ONCs. Black-Right-Pointing-Pointer From thermal power plants for CO2 capture by adsorption technology, is a new direction.

36 MATERIALS SCIENCE↗

NASA NTRS2001

A Study of Green's Function Methods Applied to Space Radiation Protection

The purpose of this research was to study the propagation of galactic ions through various materials. Galactic light ions result from the break up of heavy ion particles and their propagation through materials is modeled using the one-dimensional Boltzmann equation. When ions enter materials there can occur (i) the interaction of ions with orbital electrons which causes ionization within the material and (ii) ions collide with atoms causing production of secondary particles which penetrate deeper within the material. These processes are modeled by a continuum model. The basic idea is to place a control volume within the material and examine the change in ion flux across this control volume. In this way on can derive the basic equations for the transport of light and heavy ions in matter. Green's function perturbation methods can then be employed to solve the resulting equations using energy dependent nuclear cross sections.

Heinbockel, John H.↗

DOE OSTI2012

Tumor Volume-Adapted Dosing in Stereotactic Ablative Radiotherapy of Lung Tumors

Purpose: Current stereotactic ablative radiotherapy (SABR) protocols for lung tumors prescribe a uniform dose regimen irrespective of tumor size. We report the outcomes of a lung tumor volume-adapted SABR dosing strategy. Methods and Materials: We retrospectively reviewed the outcomes in 111 patients with a total of 138 primary or metastatic lung tumors treated by SABR, including local control, regional control, distant metastasis, overall survival, and treatment toxicity. We also performed subset analysis on 83 patients with 97 tumors treated with a volume-adapted dosing strategy in which small tumors (gross tumor volume <12 mL) received single-fraction regimens with biologically effective doses (BED) <100 Gy (total dose, 18-25 Gy) (Group 1), and larger tumors (gross tumor volume {>=}12 mL) received multifraction regimens with BED {>=}100 Gy (total dose, 50-60 Gy in three to four fractions) (Group 2). Results: The median follow-up time was 13.5 months. Local control for Groups 1 and 2 was 91.4% and 92.5%, respectively (p = 0.24) at 12 months. For primary lung tumors only (excluding metastases), local control was 92.6% and 91.7%, respectively (p = 0.58). Regional control, freedom from distant metastasis, and overall survival did not differ significantly between Groups 1 and 2. Rates of radiation pneumonitis, chest wall toxicity, and esophagitis were low in both groups, but all Grade 3 toxicities developed in Group 2 (p = 0.02). Conclusion: A volume-adapted dosing approach for SABR of lung tumors seems to provide excellent local control for both small- and large-volume tumors and may reduce toxicity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2021

Demonstration of Dimensional and Microstructural Postirradiation Examination Capabilities on MiniFuel Disk Specimens

This work demonstrates two new MiniFuel postirradiation examination (PIE) capabilities, laser profilometry and serial sectioning, for future use on alpha-U specimens for the upcoming Mini-99 irradiation campaign. Representative U–Mo disk specimens—3 mm in diameter and approximately 0.75 mm in thickness—were subjected to dimensional measurements using standard calipers, as well as a new Keyence CL-3000 series laser profilometry system designed to measure the surface of both sides of the disk simultaneously. The laser profilometry system captured thickness gradients across the specimens as large as 0.2 mm. Additionally, because of the limitations inherent to the use of calipers, caliper measurements consistently overestimated the initial volume of MiniFuel disk specimens compared to the laser profilometry method by up to 5%. However, the results show that optimal data collection and analysis procedures are crucial to avoiding significant artifacts during data collection that can lead to large errors (approximately 10%). Serial sectioning using a combination of ion milling and electron imaging was also demonstrated for the analysis of micron-sized control volumes. Nanometer- to micron-scale porosity was resolved on a representative U–Mo specimen before irradiation, even when the surface of the specimen appeared free of porosity. Future efforts to relate microstructural features (porosity, fission product distributions, and so on) to macroscale effects (dimensional changes, thermal conductivity degradation) require accurate characterization methods such as those demonstrated in this work. These advanced PIE capabilities will be applied to optimized alpha-U specimens for the upcoming Mini-99 neutron irradiation campaign.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

DOE OSTI2024

Comparing the MELCOR Aerosol Deposition Model with Exact Analytical Solutions

Exact analytical solutions are presented for the evolution of the aerosol particle mass density function in a control volume for particle deposition due to gravitational settling, thermophoresis, and diffusion. The solutions are for arbitrary initial mass density functions and are applied for an initial lognormal density function. Integration of these solutions provides the suspended mass in the control volume as a function of time. These solutions serve as an exact benchmark to assess the accuracy of numerical methods. For the numerical algorithm used in MELCOR, excellent agreement is obtained for gravitational settling, diffusive deposition, and thermophoretic deposition for the suspended aerosol mass. In all cases, the default number of discrete particle size bins of 10 is shown to converge, with hardly any advantage to using 20 size bins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗

NASA NTRS1988

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗

NASA NTRS1989

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

Woods, C. M.↗

DOE OSTI2019

MELT POOL MONITORING FOR CONTROL AND DATA ANALYTICS IN LARGE-SCALE METAL ADDITIVE MANUFACTURING

Laser-wire based Directed Energy Deposition (DED) processes are being developed at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility, in collaboration with GKN Aerospace, for large volume metal additive manufacturing of metallic structures. The technology has the potential to deliver reduced costs and lead times when compared to conventional methods of manufacturing in several industries. While complex structures are being produced with relatively high deposition rates and at near-net shape, issues persist with achieving consistent geometric accuracy and thermal stability. A significant research effort is focused on developing coordinated, multi-modal sensing and control for addressing these issues. An introduction to this large-volume metal process will be provided, followed by investigations into real-time thermal monitoring capabilities that support control and data analytics frameworks. Studies focused on melt pool monitoring via infrared thermography, thermal and geometric effects on melt pool size measurement, and interactions between process variables and thermal monitoring metrics are presented.

Gibson, Brian↗

DOE OSTI2016

TH-CD-209-04: Fuzzy Robust Optimization in Intensity-Modulated Proton Therapy Planning to Account for Range and Patient Setup Uncertainties

Purpose: We propose to apply a robust optimization model based on fuzzy-logic constraints in the intensity-modulated proton therapy (IMPT) planning subject to range and patient setup uncertainties. The purpose is to ensure the plan robustness under uncertainty and obtain the best trade-off between tumor dose coverage and organ-at-risk(OAR) sparing. Methods: Two IMPT plans were generated for 3 head-and-neck cancer patients: one used the planning target volume(PTV) method; the other used the fuzzy robust optimization method. In the latter method, nine dose distributions were computed - the nominal one and one each for ±3mm setup uncertainties along three cardinal axes and for ±3.5% range uncertainty. For tumors, these nine dose distributions were explicitly controlled by adding hard constraints with adjustable parameters. For OARs, fuzzy constraints that allow the dose to vary within a certain range were used so that the tumor dose distribution was guaranteed by minimum compromise of that of OARs. We rendered this model tractable by converting the fuzzy constraints to linear constraints. The plan quality was evaluated using dose-volume histogram(DVH) indices such as tumor dose coverage(D95%), homogeneity(D5%-D95%), plan robustness(DVH band at D95%), and OAR sparing like D1% of brain and D1% of brainstem. Results: Our model could yield clinically acceptable plans. The fuzzy-logic robust optimization method produced IMPT plans with comparable target dose coverage and homogeneity compared to the PTV method(unit: Gy[RBE]; average[min, max])(CTV D95%: 59 [52.7, 63.5] vs 53.5[46.4, 60.1], CTV D5% - D95%: 11.1[5.3, 18.6] vs 14.4[9.2, 21.5]). It also generated more robust plans(CTV DVH band at D95%: 3.8[1.2, 5.6] vs 11.5[6.2, 16.7]). The parameters of tumor constraints could be adjusted to control the tradeoff between tumor coverage and OAR sparing. Conclusion: The fuzzy-logic robust optimization generates superior IMPT with minimum compromise of OAR sparing. This research was supported by the National Cancer Institute Career Developmental Award K25CA168984, by the Fraternal Order of Eagles Cancer Research Fund Career Development Award, by The Lawrence W. and Marilyn W. Matteson Fund for Cancer Research, by Mayo Arizona State University Seed Grant, and by The Kemper Marley Foundation. eRA Person ID(s) for the Principal Investigator: 11017970 (Research Supported by National Institutes of Health)

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Removing of Core Constructions of the MR and RFT Reactors in NRC 'Kurchatov Institute' - 16054

In frame of the FTP 'Providing nuclear and radiation safety of Russia in 2008 and for the period till 2015' works of a decommissioning of the MR and RFT research reactors are carried out in National Research Center 'Kurchatov Institute'. Within these works, the equipment of 9 loop-back facilities of the MR reactor was dismantled. In the conditions, which are most approached to natural, one researched a heat, hydrodynamic and strength characteristics of elements of cores and the capital equipment of the majority of power reactors at these facilities. In 2014-2015 dismantle of core constructions of the MR reactor was carried out and dismantle of a graphite laying of the RFT reactor is begun. The beryllium blocks, which were constructional elements of the core, and the graphite blocks playing a reactor reflector role were also taken, packed into cases or containers and removed: part on an accumulative platform, and part are placed in storage of high level radioactive waste of the Center for decay and further sending for long storage. Works were performed by remote controlled mechanisms, control of their performance was exercised by means of remote diagnostics: Gammavizor, the scanning spectrometer systems. The following methods and the organizational principles were applied to providing of radiation safety of these works: - the constant radiation control including individual control of a dose of external radiation, control of volume activity of radionuclides in air, annual control of internal radiation doses of the personnel; - dust suppression and fixing of radioactive contamination by dispersion of polymeric structures; - carrying out decontamination on a shift basis rooms with use of detergents; - use of the remote controlled mechanisms equipped with hydraulic scissors, nippers, application of the diamond are sharp that allowed to perform the majority of technological operations remotely, using 'cold' methods of fragmentation; - development on each type of works of the operation project and schedules. As a result, of works of 2011-2014 of the decommissioning of the MR and RFT research reactors more than 630 t of the equipment were dismantled. From them more than 455 t of radioactive waste was. The total activity of radwaste exceeded 300.0 GBq. From the pool of the reactor, about 160 channels of loop-back installations and near 230 beryllium and graphite are taken, fragmented and packed. The total activity on radionuclides {sup 60}Co and {sup 137}Cs of them was near 55 TBq. From the territory of the center about 2550 m{sup 3} of firm and 4600 m{sup 3} of liquid radwaste are removed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2012

A method for deriving a 4D-interpolated balanced planning target for mobile tumor radiotherapy

Purpose: Tumor control and normal tissue toxicity are strongly correlated to the tumor and normal tissue volumes receiving high prescribed dose levels in the course of radiotherapy. Planning target definition is, therefore, crucial to ensure favorable clinical outcomes. This is especially important for stereotactic body radiation therapy of lung cancers, characterized by high fractional doses and steep dose gradients. The shift in recent years from population-based to patient-specific treatment margins, as facilitated by the emergence of 4D medical imaging capabilities, is a major improvement. The commonly used motion-encompassing, or internal-target volume (ITV), target definition approach provides a high likelihood of coverage for the mobile tumor but inevitably exposes healthy tissue to high prescribed dose levels. The goal of this work was to generate an interpolated balanced planning target that takes into account both tumor coverage and normal tissue sparing from high prescribed dose levels, thereby improving on the ITV approach. Methods: For each 4DCT dataset, 4D deformable image registration was used to derive two bounding targets, namely, a 4D-intersection and a 4D-composite target which minimized normal tissue exposure to high prescribed dose levels and maximized tumor coverage, respectively. Through definition of an ''effective overlap volume histogram'' the authors derived an ''interpolated balanced planning target'' intended to balance normal tissue sparing from prescribed doses with tumor coverage. To demonstrate the dosimetric efficacy of the interpolated balanced planning target, the authors performed 4D treatment planning based on deformable image registration of 4D-CT data for five previously treated lung cancer patients. Two 4D plans were generated per patient, one based on the interpolated balanced planning target and the other based on the conventional ITV target. Plans were compared for tumor coverage and the degree of normal tissue sparing resulting from the new approach was quantified. Results: Analysis of the 4D dose distributions from all five patients showed that while achieving tumor coverage comparable to the ITV approach, the new planning target definition resulted in reductions of lung V{sub 10}, V{sub 20}, and V{sub 30} of 6.3% {+-} 1.7%, 10.6% {+-} 3.9%, and 12.9% {+-} 5.5%, respectively, as well as reductions in mean lung dose, mean dose to the GTV-ring and mean heart dose of 8.8% {+-} 2.5%, 7.2% {+-} 2.5%, and 10.6% {+-} 3.6%, respectively. Conclusions: The authors have developed a simple and systematic approach to generate a 4D-interpolated balanced planning target volume that implicitly incorporates the dynamics of respiratory-organ motion without requiring 4D-dose computation or optimization. Preliminary results based on 4D-CT data of five previously treated lung patients showed that this new planning target approach may improve normal tissue sparing without sacrificing tumor coverage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

NASA NTRS1991

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗

NASA NTRS1993

Aerodynamic control of NASP-type vehicles through Vortex manipulation. Volume 1: Static water tunnel tests

Water tunnel tests were conducted on a NASP-type configuration to evaluate different pneumatic Forebody Vortex Control (FVC) methods. Flow visualization and yawing moment measurements were performed at angles of attack from 0 deg to 30 deg. The pneumatic techniques tested included jet and slot blowing. In general, blowing can be used efficiently to manipulate the forebody vortices at angles of attack greater than 20 deg. These vortices are naturally symmetric up to alpha = 25 deg and asymmetric between 25 deg and 30 deg angle of attack. Results indicate that tangential aft jet blowing is the most promising method for this configuration. Aft jet blowing produces a yawing moment towards the blowing side and the trends with blowing rate are well behaved. The size of the nozzle is not the dominant factor in the blowing process; the change in the blowing 'momentum,' i.e., the product of the mass flow rate and the velocity of the jet, appears to be the important parameter in the water tunnel (incompressible and unchoked flow at the nozzle exit). Forward jet blowing is very unpredictable and sensitive to mass flow rate changes. Slot blowing (with the exception of very low blowing rates) acts as a flow 'separator'; it promotes early separation on the blow side, producing a yawing moment toward the non-blowing side for the C(sub mu) range investigated.

Suarez, Carlos J.↗

NASA NTRS1996

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗