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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Grid generation for general 2-D regions using hyperbolic equations

A method for applying a hyperbolic grid generation scheme to the construction of meshes in general 2-D regions has been developed. This approach, which follows the theory developed by Steger and Chaussee (1980) and the algorithm outlined by Kinsey and Barth (1984), is based on improving local grid control. This is accomplished by adding an angle control source term to the equations and using a new algorithm for computing the volume source term. These modifications lead to superior methods for fixing the 'local' problems of hyperbolic grid generation, namely, propagation of initial discontinuities and formation of grid shocks (crossing grid lines). More importantly, a method for solving the global problem of constraining the grid with more than one boundary (internal grid generation) has been developed. These algorithms have been implemented in an interactive grid generation program and the results for several geometries are presented and discussed.

Cordova, Jeffrey Q.↗
DOE OSTI2017

$L^1$ penalization of volumetric dose objectives in optimal control of PDEs

This work is concerned with a class of PDE-constrained optimization problems that are motivated by an application in radiotherapy treatment planning. Here the primary design objective is to minimize the volume where a functional of the state violates a prescribed level, but prescribing these levels in the form of pointwise state constraints leads to infeasible problems. We therefore propose an alternative approach based on L 1 penalization of the violation that is also applicable when state constraints are infeasible. We establish well-posedness of the corresponding optimal control problem, derive first-order optimality conditions, discuss convergence of minimizers as the penalty parameter tends to infinity, and present a semismooth Newton method for their efficient numerical solution. Finally, the performance of this method for a model problem is illustrated and contrasted with an alternative approach based on (regularized) state constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SU-F-T-335: Piecewise Uniform Dose Prescription and Optimization Based On PET/CT Images

Purpose: In intensity modulated radiation therapy (IMRT), the tumor target volume is given a uniform dose prescription, which does not consider the heterogeneous characteristics of tumor such as hypoxia, clonogen density, radiosensitivity, tumor proliferation rate and so on. Our goal is to develop a nonuniform target dose prescription method which can spare organs at risk (OARs) better and does not decrease the tumor control probability (TCP). Methods: We propose a piecewise uniform dose prescription (PUDP) based on PET/CT images of tumor. First, we propose to delineate biological target volumes (BTV) and sub-biological target volumes (sub-BTVs) by our Hierarchical Mumford-Shah Vector Model based on PET/CT images of tumor. Then, in order to spare OARs better, we make the BTV mean dose minimized while restrict the TCP to a constant. So, we can get a general formula for determining an optimal dose prescription based on a linearquadratic model (LQ). However, this dose prescription is high heterogeneous, it is very difficult to deliver by IMRT. Therefore we propose to use the equivalent uniform dose (EUD) in each sub-BTV as its final dose prescription, which makes a PUDP for the BTV. Results: We have evaluated the IMRT planning of a patient with nasopharyngeal carcinoma respectively using PUDP and UDP. The results show that the highest and mean doses inside brain stem are 48.425Gy and 19.151Gy respectively when the PUDP is used for IMRT planning, while they are 52.975Gy and 20.0776Gy respectively when the UDP is used. Both of the resulting TCPs(0.9245, 0.9674) are higher than the theoretical TCP(0.8739), when 70Gy is delivered to the BTV. Conclusion: Comparing with the UDP, the PUDP can spare the OARs better while the resulting TCP by PUDP is not significantly lower than by UDP. This work was supported in part by National Natural Science Foundation of China undergrant no.61271382 and by the foundation for construction of scientific project platform forthe cancer hospital of Hunan province.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1990

Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction

High order accurate finite-volume schemes for solving the Euler equations of gasdynamics are developed. Central to the development of these methods are the construction of a k-exact reconstruction operator given cell-averaged quantities and the use of high order flux quadrature formulas. General polygonal control volumes (with curved boundary edges) are considered. The formulations presented make no explicit assumption as to complexity or convexity of control volumes. Numerical examples are presented for Ringleb flow to validate the methodology.

Barth, Timothy J.↗
NASA NTRS2012

Method and apparatus for electrokinetic transport

Controlled electrokinetic transport of constituents of liquid media can be achieved by connecting at least two volumes containing liquid media with at least one dielectric medium with opposing dielectric surfaces in direct contact with said liquid media, and establishing at least one conduit across said dielectric medium, with a conduit inner surface surrounding a conduit volume and at least a first opening and a second opening opposite to the first opening. The conduit is arranged to connect two volumes containing liquid media and includes a set of at least three electrodes positioned in proximity of the inner conduit surface. A power supply is arranged to deliver energy to the electrodes such that time-varying potentials inside the conduit volume are established, where the superposition of said potentials represents at least one controllable traveling potential well that can travel between the opposing conduit openings.

James, Patrick Ismail↗
DOE OSTI2017

All-Materials-Inclusive Flash Spark Plasma Sintering

A new flash (ultra-rapid) spark plasma sintering method applicable to various materials systems, regardless of their electrical resistivity, is developed. A number of powders ranging from metals to electrically insulative ceramics have been successfully densified resulting in homogeneous microstructures within sintering times of 8–35 s. A finite element simulation reveals that the developed method, providing an extraordinary fast and homogeneous heating concentrated in the sample’s volume and punches, is applicable to all the different samples tested. The utilized uniquely controllable flash phenomenon is enabled by the combination of the electric current concentration around the sample and the confinement of the heat generated in this area by the lateral thermal contact resistance. The presented new method allows: extending flash sintering to nearly all materials, controlling sample shape by an added graphite die, and an energy efficient mass production of small and intermediate size objects. This approach represents also a potential venue for future investigations of flash sintering of complex shapes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Kimberlina 1.2 Velocity Models and Seismic Data

Kimberlina 1.2 Velocity model and synthetic seismic data, produced in collaboration of teams at the National Energy Technology Laboratory, Los Alamos National Laboratory, and Lawrence Livermore National Laboratory through the National Risk Assessment Partnership. Data is associated with the following publication: Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019. The Kimberlina 1.2 Velocity models were produced by Zan Wang, Robert Dilmore, William Harbert, and Lianjie Huang at NETL. The following citations are directly related to the creation of the velocity models: Wang, Z. Harbert, W., Dilmore, R., Huang, L. Modeling of time-lapse seismic monitoring using CO2 leakage simulations for a model CO2 storage site with realistic geology: Application in assessment of early leak-detection capabilities. International Journal of Greenhouse Gas Control. V. 76, September 2018, Pages 39-52. https://doi.org/10.1016/j.ijggc.2018.06.011 Wang, Z., Dilmore, R., Harbert, W. Inferring CO2 saturation from synthetic surface seismic and downhole monitoring data using machine learning for leakage detection at CO2 sequestration sites. International Journal of Greenhouse Gas Control, V. 100, September 2020. https://doi.org/10.1016/j.ijggc.2020.103115 The velocity models were built based on the Kimberlina 1.2 aquifer impact data which is associated with the following publications: Buscheck, T.A., Mansoor, K., Yang, X., Wainwright, H., and Carroll, S. (2019). Downhole pressure and chemical monitoring for CO2 and brine leak detection in aquifers above a CO2 storage reservoir. International Journal of Greenhouse Gas Control. 91. 102812. 10.1016/j.ijggc.2019.102812. Xianjin Yang, Thomas A. Buscheck, Kayyum Mansoor, Zan Wang, Kai Gao, Lianjie Huang, Delphine Appriou, Susan A. Carroll, Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, Volume 90, 2019, 102803, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2019.102803 The synthetic seismic data was produced by Youzuo Lin and team at LANL, and are associated with the following citations: Jordan, P. D., and J. L. Wagoner. Characterizing Construction of Existing Wells to a CO2 Storage Target: The Kimberlina Site, California. Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019.

CO2 Leakage↗
NASA NTRS2002

Spectral (Finite) Volume Method for Conservation Laws on Unstructured Grids II: Extension to Two Dimensional Scalar Equation

The framework for constructing a high-order, conservative Spectral (Finite) Volume (SV) method is presented for two-dimensional scalar hyperbolic conservation laws on unstructured triangular grids. Each triangular grid cell forms a spectral volume (SV), and the SV is further subdivided into polygonal control volumes (CVs) to supported high-order data reconstructions. Cell-averaged solutions from these CVs are used to reconstruct a high order polynomial approximation in the SV. Each CV is then updated independently with a Godunov-type finite volume method and a high-order Runge-Kutta time integration scheme. A universal reconstruction is obtained by partitioning all SVs in a geometrically similar manner. The convergence of the SV method is shown to depend on how a SV is partitioned. A criterion based on the Lebesgue constant has been developed and used successfully to determine the quality of various partitions. Symmetric, stable, and convergent linear, quadratic, and cubic SVs have been obtained, and many different types of partitions have been evaluated. The SV method is tested for both linear and non-linear model problems with and without discontinuities.

Wang, Z. J.↗
NASA NTRS2011

Acoustic Testing of Flight Hardware Using Loudspeakers: How Much do We Know About This Method of Testing?

Loudspeakers have been used for acoustic qualification of spacecrafts, reflectors, solar panels, and other acoustically responsive structures for more than a decade. Even though a lot of hardware has been acoustic tested using this method, the nature of the acoustic field generated by controlling an ensemble of speakers with and without the hardware in the test volume has not been thoroughly investigated. Limited measurements from some of the recent speaker tests used to qualify flight hardware have indicated significant spatial variation of the acoustic field within the test volume. Also structural responses have been reported to differ when similar tests were performed using reverberant chambers. Unlike the reverberant chamber acoustic test, for which the acoustic field in most chambers is known to be diffuse except below several tens of Hz where acoustic standing waves and large spatial variations exist, the characteristics of the acoustic field within the speaker test volume has not been quantified. It has only been recently that a detailed acoustic field characterization of speaker testing has been made at Jet Propulsion Laboratory (JPL) with involvement of various organizations. To address the impact of non-uniform acoustic field on structures, a series of acoustic tests were performed using a flat panel and a 3-ft cylinder exposed to the field controlled by speakers and repeated in a reverberant chamber. The analysis of the data from this exercise reveals that there are significant differences both in the acoustic field and in the structural responses. In this paper the differences between the two methods are reviewed in some detail and the over- or under-testing of articles that could pose un-anticipated structural and flight qualification issues are discussed. A framework for discussing the validity of the speaker acoustic testing method with the current control system and a path forward for improving it will be provided.

Acoustics↗
NASA NTRS1987

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗
NASA NTRS1988

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗
NASA NTRS2011

Linear Approximation to Optimal Control Allocation for Rocket Nozzles with Elliptical Constraints

In this paper we present a straightforward technique for assessing and realizing the maximum control moment effectiveness for a launch vehicle with multiple constrained rocket nozzles, where elliptical deflection limits in gimbal axes are expressed as an ensemble of independent quadratic constraints. A direct method of determining an approximating ellipsoid that inscribes the set of attainable angular accelerations is derived. In the case of a parameterized linear generalized inverse, the geometry of the attainable set is computationally expensive to obtain but can be approximated to a high degree of accuracy with the proposed method. A linear inverse can then be optimized to maximize the volume of the true attainable set by maximizing the volume of the approximating ellipsoid. The use of a linear inverse does not preclude the use of linear methods for stability analysis and control design, preferred in practice for assessing the stability characteristics of the inertial and servoelastic coupling appearing in large boosters. The present techniques are demonstrated via application to the control allocation scheme for a concept heavy-lift launch vehicle.

Orr, Jeb S.↗
DOE OSTI2016

Technical Note: Cortical thickness and density estimation from clinical CT using a prior thickness-density relationship

Purpose: Cortical thickness and density are critical components in determining the strength of bony structures. Computed tomography (CT) is one possible modality for analyzing the cortex in 3D. In this paper, a model-based approach for measuring the cortical bone thickness and density from clinical CT images is proposed. Methods: Density variations across the cortex were modeled as a function of the cortical thickness and density, location of the cortex, density of surrounding tissues, and imaging blur. High resolution micro-CT data of cadaver proximal femurs were analyzed to determine a relationship between cortical thickness and density. This thickness-density relationship was used as prior information to be incorporated in the model to obtain accurate measurements of cortical thickness and density from clinical CT volumes. The method was validated using micro-CT scans of 23 cadaver proximal femurs. Simulated clinical CT images with different voxel sizes were generated from the micro-CT data. Cortical thickness and density were estimated from the simulated images using the proposed method and compared with measurements obtained using the micro-CT images to evaluate the effect of voxel size on the accuracy of the method. Then, 19 of the 23 specimens were imaged using a clinical CT scanner. Cortical thickness and density were estimated from the clinical CT images using the proposed method and compared with the micro-CT measurements. Finally, a case-control study including 20 patients with osteoporosis and 20 age-matched controls with normal bone density was performed to evaluate the proposed method in a clinical context. Results: Cortical thickness (density) estimation errors were 0.07 ± 0.19 mm (−18 ± 92 mg/cm{sup 3}) using the simulated clinical CT volumes with the smallest voxel size (0.33 × 0.33 × 0.5 mm{sup 3}), and 0.10 ± 0.24 mm (−10 ± 115 mg/cm{sup 3}) using the volumes with the largest voxel size (1.0 × 1.0 × 3.0 mm{sup 3}). A trend for the cortical thickness and density estimation errors to increase with voxel size was observed and was more pronounced for thin cortices. Using clinical CT data for 19 of the 23 samples, mean errors of 0.18 ± 0.24 mm for the cortical thickness and 15 ± 106 mg/cm{sup 3} for the density were found. The case-control study showed that osteoporotic patients had a thinner cortex and a lower cortical density, with average differences of −0.8 mm and −58.6 mg/cm{sup 3} at the proximal femur in comparison with age-matched controls (p-value < 0.001). Conclusions: This method might be a promising approach for the quantification of cortical bone thickness and density using clinical routine imaging techniques. Future work will concentrate on investigating how this approach can improve the estimation of mechanical strength of bony structures, the prevention of fracture, and the management of osteoporosis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Low temperature synthesis of hierarchical TiO 2 nanostructures for high performance perovskite solar cells by pulsed laser deposition

A promising way to advance perovskite solar cells is to improve the quality of the electron transport material e.g., titanium dioxide (TiO 2 ) in a direction that increases electron transport and extraction. Although dense TiO 2 films are easily grown in solution, efficient electron extraction suffers due to a lack of interfacial contact area with the perovskite. Conversely, mesoporous films do offer high surface-area-to-volume ratios, thereby promoting efficient electron extraction, but their morphology is relatively difficult to control via conventional solution synthesis methods. Here, a pulsed laser deposition method was used to assemble TiO 2 nanoparticles into TiO 2 hierarchical nanoarchitectures having the anatase crystal structure, and prototype solar cells employing these structures yielded power conversion efficiencies of ~ 14%. Our approach demonstrates a way to grow high aspect-ratio TiO 2 nanostructures for improved interfacial contact between TiO 2 and perovskite materials, leading to high electron-hole pair separation and electron extraction efficiencies for superior photovoltaic performance. In addition, compared to conventional solution-processed TiO 2 films that require 500 °C to obtain a good crystallinity, our relatively low temperature (300 °C) TiO 2 processing method may promote reduced energy-consumption during device fabrication as well as enable compatibility with various flexible polymer substrates.

14 SOLAR ENERGY↗
NASA NTRS1983

Derivation of the fundamental equation of sound generated by moving aerodynamic surfaces

Ffowcs Williams and Hawkings (1969) based their derivation of the fundamental equation of the sound generated by arbitrarily moving aerodynamic surfaces on the study of mass and momentum balance of a control volume imbedding a mathematical surface(s) exactly corresponding to real surface(s). These investigators also sketched an alternative method, employing generalized functions, for its derivation. This latter method, which was later developed by Farassat (1975), is purely mathematical and formal. Goldstein (1976) used the free-space Green function to produce an implicit derivation of the Ffowcs Williams and Hawkings equation. In the study presented here, Lowson's (1965) concept of moving point singularities is generalized to moving surface singularities, and a new derivation is given of the fundamental equation. The derivation is based on topological considerations of the underlying space, the fluid medium, and the integral properties of the Dirac delta function.

Aggarwal, H. R.↗
NASA NTRS1992

On-orbit cryogenic fluid transfer research at NASA Lewis Research Center

A summary of research into on-orbit cryogenic fluid transfer at the NASA Lewis Research Center (LRC) is presented. Variable test parameters and liquid injection configurations elucidated the conditions necessary for a successful transfer of liquid hydrogen by the no-vent fill method. The model is based on conservation of mass and a first-law energy balance for a control volume. The ullage, the bulk liquid, and the tank wall are each represented by a single node. The magnitude of the maximum receiver tank pressure was found to be dependent on the liquid inlet temperature, the inlet mass flow rate, and the initial temperature. A comparison of the test data and the analytical results for a no-vent fill test with a small receiver tank is presented.

Taylor, W. J.↗
NASA NTRS2017

Influence of Containment on the Growth of Germanium-Silicon in Microgravity

A series of Ge(sub 1-x)Si(sub x) crystal growth experiments are planned to be conducted in the Low Gradient Furnace (LGF) onboard the International Space Station. The primary objective of the research is to determine the influence of containment on the processing-induced defects and impurity incorporation in germanium-silicon alloy crystals. A comparison will be made between crystals grown by the normal and 'detached' Bridgman methods and the ground-based float zone technique. 'Detached' or 'dewetted' Bridgman growth is similar to regular Bridgman growth in that most of the melt is in contact with the crucible wall, but the crystal is separated from the wall by a small gap, typically of the order of 10-100 microns. A meniscus bridges this gap between the top of the crystal and the crucible wall. Theoretical models indicate that an important parameter governing detachment is the pressure differential across this meniscus. An experimental method has been developed to control this pressure differential in microgravity that does not require connection of the ampoule volume to external gases or changes in the temperature profile during growth. Experiments will be conducted with positive, negative or zero pressure differential across the meniscus. Characterization results of ground-based experiments, including etch pit density, synchrotron white beam X-ray topography and double axis X-ray diffraction will also be described.

Volz, M. P.↗