Search NASA⌕ Search

SEARCH · Search NASA

Results for “CASE METHOD”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2011

Linear entropy in quantum phase space

We calculate the quantum Renyi entropy in a phase-space representation for either fermions or bosons. This can also be used to calculate purity and fidelity, or the entanglement between two systems. We show that it is possible to calculate the entropy from sampled phase-space distributions in normally ordered representations, although this is not possible for all quantum states. We give an example of the use of this method in an exactly soluble thermal case. The quantum entropy cannot be calculated at all using sampling methods in classical symmetric (Wigner) or antinormally ordered (Husimi) phase spaces, due to inner-product divergences. The preferred method is to use generalized Gaussian phase-space methods, which utilize a distribution over stochastic Green's functions. We illustrate this approach by calculating the reduced entropy and entanglement of bosonic or fermionic modes coupled to a time-evolving, non-Markovian reservoir.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Validation of an in-vivo proton beam range check method in an anthropomorphic pelvic phantom using dose measurements

Purpose: In-vivo dosimetry and beam range verification in proton therapy could play significant role in proton treatment validation and improvements. In-vivo beam range verification, in particular, could enable new treatment techniques one of which could be the use of anterior fields for prostate treatment instead of opposed lateral fields as in current practice. This paper reports validation study of an in-vivo range verification method which can reduce the range uncertainty to submillimeter levels and potentially allow for in-vivo dosimetry. Methods: An anthropomorphic pelvic phantom is used to validate the clinical potential of the time-resolved dose method for range verification in the case of prostrate treatment using range modulated anterior proton beams. The method uses a 3 × 4 matrix of 1 mm diodes mounted in water balloon which are read by an ADC system at 100 kHz. The method is first validated against beam range measurements by dose extinction measurements. The validation is first completed in water phantom and then in pelvic phantom for both open field and treatment field configurations. Later, the beam range results are compared with the water equivalent path length (WEPL) values computed from the treatment planning system XIO. Results: Beam range measurements from both time-resolved dose method and the dose extinction method agree with submillimeter precision in water phantom. For the pelvic phantom, when discarding two of the diodes that show sign of significant range mixing, the two methods agree with ±1 mm. Only a dose of 7 mGy is sufficient to achieve this result. The comparison to the computed WEPL by the treatment planning system (XIO) shows that XIO underestimates the protons beam range. Quantifying the exact XIO range underestimation depends on the strategy used to evaluate the WEPL results. To our best evaluation, XIO underestimates the treatment beam range between a minimum of 1.7% and maximum of 4.1%. Conclusions: Time-resolved dose measurement method satisfies the two basic requirements, WEPL accuracy and minimum dose, necessary for clinical use, thus, its potential for in-vivo protons range verification. Further development is needed, namely, devising a workflow that takes into account the limits imposed by proton range mixing and the susceptibility of the comparison of measured and expected WEPLs to errors on the detector positions. The methods may also be used for in-vivo dosimetry and could benefit various proton therapy treatments.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1950

Application of radial-equilibrium condition to axial-flow compressor and turbine design

Basic general equations governing the three-dimensional compressible flow of gas through a compressor or turbine are given in terms of total enthalpy, entropy, and velocity components of the gas. Two methods of solution are obtained for the simplified, steady axially symmetric flow; one involves the use of a number of successive planes normal to the axis of the machine and short distances apart, and the other involves only three stations for a stage in which an appropriate radial-flow path is used. Methods of calculation for the limiting cases of zero and infinite blade aspect ratios and an approximate method of calculation for finite blade aspect ratio are also given. In these methods, the blade loading and the shape of the annular passage wall may be arbitrarily specified.

Wu, Chung-Hua↗
NASA NTRS2003

Launch Vehicle Design and Optimization Methods and Priority for the Advanced Engineering Environment

NASA's Advanced Engineering Environment (AEE) is a research and development program that will improve collaboration among design engineers for launch vehicle conceptual design and provide the infrastructure (methods and framework) necessary to enable that environment. In this paper, three major technical challenges facing the AEE program are identified, and three specific design problems are selected to demonstrate how advanced methods can improve current design activities. References are made to studies that demonstrate these design problems and methods, and these studies will provide the detailed information and check cases to support incorporation of these methods into the AEE. This paper provides background and terminology for discussing the launch vehicle conceptual design problem so that the diverse AEE user community can participate in prioritizing the AEE development effort.

Rowell, Lawrence F.↗
NASA NTRS1988

General image method in a plane-layered elastostatic medium

The general-image method presently used to obtain the elastostatic fields in plane-layered media relies on the use of potentials in order to represent elastic fields. For the case of a single interface, this method yields the displacement field in closed form, and is applicable to antiplane, plane, and three-dimensional problems. In the case of multiplane interfaces, the image method generates the displacement fields in terms of infinite series whose convergences can be accelerated to improve method efficiency.

Fares, N.↗
DOE OSTI2014

One-parameter families of supersymmetric isospectral potentials from Riccati solutions in function composition form

In the context of supersymmetric quantum mechanics, we define a potential through a particular Riccati solution of the composition form (F∘f)(x)=F(f(x)) and obtain a generalized Mielnik construction of one-parameter isospectral potentials when we use the general Riccati solution. Some examples for special cases of F and f are given to illustrate the method. An interesting result is obtained in the case of a parametric double well potential generated by this method, for which it is shown that the parameter of the potential controls the heights of the localization probability in the two wells, and for certain values of the parameter the height of the localization probability can be higher in the smaller well. -- Highlights: •Function-composition generalization of parametric isospectral potentials is presented. •Mielnik one-parameter family of harmonic potentials is obtained as a particular case. •Graphical discussion of regular and singular regions in the parameter space is given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1959

The Stabilizing Effectiveness of Conical Flares on Bodies with Conical Noses

An analysis is presented of published results of force tests on 80 cone-cylinder-flare configurations at Mach numbers of 2.18, 2.81, and 4.04. The contributions, excluding interference effects, of the cone-cylinder bodies to the over-all normal force derivatives have been removed by means of the second-order shock-expansion method, and the normal force derivatives at zero angle of attack due to the flares alone are shown. The results from a wide variety of configurations are correlated by plotting ratios of the normal force derivatives of the flares to the normal force derivatives of cones having the same included angle. Comparisons are made of the experimental normal force results with the normal force derivatives obtained by assuming conical flow over the flares and with those obtained by use of the second-order shock-expansion method. The comparisons show that use of the second-order shock-expansion method is generally the superior of the two, and in most cases gives values of the normal force derivatives of the flares which agree very well with the experimental results. Centers of pressure of the flares are presented and comparisons are made with results obtained from the theories mentioned. In general, the comparisons show that the assumption of conical flow over the flares is comparable to use of the second-order shock-expansion method in determining the centers of pressure, and in many cases both methods give values which agree closely with the experimental results.

Kirk, Donn B.↗
DOE OSTI2013

A BAYESIAN METHOD FOR THE ANALYSIS OF THE DUST EMISSION IN THE FAR-INFRARED AND SUBMILLIMETER

We present a method, based on Bayesian statistics, to fit the dust emission parameters in the far-infrared and submillimeter wavelengths. The method estimates the dust temperature and spectral emissivity index, plus their relationship, properly taking into account the statistical and systematic uncertainties. We test it on three sets of simulated sources detectable by the Herschel Space Observatory in the PACS and SPIRE spectral bands (70-500 {mu}m), spanning over a wide range of dust temperatures. The simulated observations are a one-component interstellar medium and two two-component sources, both warm (H II regions) and cold (cold clumps (CCs)). We first define a procedure to identify the better model, then we recover the parameters of the model and measure their physical correlations by means of a Markov Chain Monte Carlo algorithm adopting multi-variate Gaussian priors. In this process, we assess the reliability of the model recovery and of parameter estimation. We conclude that the model and parameters are properly recovered only under certain circumstances and that false models may be derived in some cases. We applied the method to a set of 91 starless CCs in an interarm region of the Galactic plane with low star formation activity, observed by Herschel in the Hi-GAL survey. Our results are consistent with a temperature-independent spectral index.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1963

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗
DOE OSTI2019

Dynamic Mode Decomposition for Compressive System Identification

Dynamic mode decomposition has emerged as a leading technique to identify spatiotemporal coherent structures from high-dimensional data, benefiting from a strong connection to nonlinear dynamical systems via the Koopman operator. In this work, we integrate and unify two recent innovations that extend DMD to systems with actuation [Proctor et al., 2016] and systems with heavily subsampled measurements [Brunton et al., 2015]. When combined, these methods yield a novel framework for compressive system identification [code is publicly available at: https://github.com/zhbai/cDMDc]. It is possible to identify a low-order model from limited input-output data and reconstruct the associated full-state dynamic modes with compressed sensing, adding interpretability to the state of the reduced-order model. Moreover, when full-state data is available, it is possible to dramatically accelerate downstream computations by first compressing the data. We demonstrate this unified framework on two model systems, investigating the effects of sensor noise, different types of measurements (e.g., point sensors, Gaussian random projections, etc.), compression ratios, and different choices of actuation (e.g., localized, broadband, etc.). In the first example, we explore this architecture on a test system with known low-rank dynamics and an artificially inflated state dimension. The second example consists of a real-world engineering application given by the fluid flow past a pitching airfoil at low Reynolds number. This example provides a challenging and realistic test-case for the proposed method, and results demonstrate that the dominant coherent structures are well characterized despite actuation and heavily subsampled data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Extended FFT-based micromechanical formulation to consider general non-periodic boundary conditions

Here, this paper presents a new approach for applying non-periodic boundary conditions in the context of FFT-based methods to solve micromechanical problems in heterogeneous solids. The domain of the original problem is extended to satisfy the periodicity requirements at the boundary of the extended domain. The velocity constraint on the boundary of the original domain is replaced by a corresponding constraint on the velocity gradient in the extended volume, and a two-level augmented Lagrangian method is used to enforce the constraint. The proposed method is implemented as an extension of the large-strain elasto-viscoplastic FFT-based (LS-EVPFFT) model of Zecevic et al. (2022). The proposed method is verified in the cases of fully imposed velocity boundary conditions and mixed velocity/traction-free boundary conditions. The accuracy and convergence of the method are studied next, followed by applications to bending and indentation of polycrystals that illustrate the extended capabilities of the proposed formulation.

36 MATERIALS SCIENCE↗
DOE OSTI2015

SU-E-T-25: A Dosimetric Comparison of Three-Dimension Conformal and Intensity-Modulated Radiation Therapy in Esophageal Cancer

Purpose: dosimetric comparison of 3DCRT and IMRT in 9 esophageal cancer. The aim of this paper is to know which of these two techniques is dosimetrically more favorable dosimetrically at both the CTV coverage and dose obtained in the relevant organs at risk, in this case, lungs and heart, as the spinal cord received in all cases below 45 Gy. Methods: we chose 9 patients from our center (CHN) with the same type of esophageal cancer and in which the prescribed dose was the same, 54 Gy. For these treatments we have used the same fields and the same angles (AP (0 °), OPD (225°–240°) and OPI (125°–135°)).All plans have been implemented using Eclipse (version 11.0) with AAA( Analytical Anisotropic Algorithm )(Version 11.0.31). Results: To analyze the coverage of the CTV, we have evaluated the D99% and found that the average dose received by 99% of CTV with IMRT is 53.8 ± 0.4 Gy (99.6% of the prescribed dose) and the mean value obtained with 3DCRT is 52.3 ± 0.6 Gy (96.8% of the prescribed dose).The last data analyzed was the D2% of PTV, a fact that gives us information on the maximum dose received by our PTV. D2% of the PTV for IMRT planning is 55.4 ± 0.4 Gy (102.6% of the prescribed dose) and with 3DCRT is 56.8 ± 0.7 Gy (105.2% of the prescribed dose).All parameters analyzed at risk organs (V30, V40, V45 and V50 for the case of heart and V5, V10, V15 and V20 for the case of the lungs) provide us irradiated volume percentages lower in IMRT than 3DCRT. Conclusion: IMRT provides a considerable improvement in the coverage of the CTV and the doses to organs at risk.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-642: PTV Is the Voxel-Wise Worst-Case of CTV in Prostate Photon Therapy

Purpose: To examine the adequacy of the planning target volume (PTV) dose distribution as the worst-case representation of clinical target volume (CTV) dose distribution in prostate volumetric-modulated arc therapy (VMAT) plans. Methods: Ten intact prostate cancer cases treated by VMAT at our institution were randomly selected. Isocenter was shifted in the three cardinal directions by a displacement equal to the PTV expansion on the CTV (±3 mm) for a total of six shifted plans per original plan. Rotationally-perturbed plans were generated with a couch rotation of ±1° to simulate patient yaw. The eight perturbed dose distributions were recalculated in the treatment planning system using the same, fixed fluence map as the original plan. The voxel-wise worst-case CTV dose distribution was constructed from the minimum value per voxel from the eight perturbed doses. The resulting dose volume histograms (DVH) were evaluated for statistical correlation between the worst-case CTV and nominal PTV dose distributions based on D95% by Wilcoxon signed-rank test with significance level p ≤ 0.05. Results: Inspection demonstrates the PTV DVH in the nominal dose distribution is bounded by the CTV DVH in the worst-case dose distribution. Comparison of D95% for the two dose distributions by Wilcoxon signed-rank test gives p = 0.131. Therefore the null hypothesis cannot be rejected since the difference in median values is not statistically significant. Conclusion: The assumption that the nominal dose distribution for PTV represents the worst-case dose distribution for CTV appears valid for the ten plans under examination. Although the worst-case dose distribution is unphysical since the dose per voxel is chosen independently, it serves as a lower bound for the possible CTV coverage. Furthermore, this is consistent with the unphysical nature of the PTV. Minor discrepancies between the two dose distributions are expected since the dose cloud is not strictly static. Funding Support: 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↗
DOE OSTI2021

A Frequency-Domain-Based Algorithm for Detecting Microseismicity Using Dense Surface Seismic Arrays

We propose a new frequency-domain-based algorithm for detecting small-magnitude seismic events using dense surface seismic arrays. Our proposed method takes advantage of the high energy carried by S waves, and approximate known source locations, which are used to rotate the horizontal components to obtain the maximum amplitude. By surrounding the known source area with surface geophones, we achieve a favorable geometry for locating the detected seismic events with the backprojection method. To test our new detection method, we used a dense circular array, consisting of 151 5 Hz three-component geophones, over a 5 km aperture that was in operation at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) in southcentral Utah. We apply the new detection method during a small-scale test injection phase at FORGE, and during an aftershock sequence of an Mw 4.1 earthquake located ~30 km north of the geophone array, within the Black Rock volcanic field. We are able to detect and locate microseismic events (Mw<0) during injections, despite the high level of anthropogenic activity, and several aftershocks that are missing from the regional catalog. By comparing our method with known algorithms that operate both in the time and frequency domain, we show that our proposed method performs better in the case of the FORGE injection monitoring, and equally well for the off-array aftershock sequence. Our new method has the potential to improve microseismic event detections even in extremely noisy environments, and the proposed location scheme serves as a direct discriminant between true and false detections.

Geochemistry & Geophysics↗
NASA NTRS1977

A general method for calculating three-dimensional compressible laminar and turbulent boundary layers on arbitrary wings

The method described utilizes a nonorthogonal coordinate system for boundary-layer calculations. It includes a geometry program that represents the wing analytically, and a velocity program that computes the external velocity components from a given experimental pressure distribution when the external velocity distribution is not computed theoretically. The boundary layer method is general, however, and can also be used for an external velocity distribution computed theoretically. Several test cases were computed by this method and the results were checked with other numerical calculations and with experiments when available. A typical computation time (CPU) on an IBM 370/165 computer for one surface of a wing which roughly consist of 30 spanwise stations and 25 streamwise stations, with 30 points across the boundary layer is less than 30 seconds for an incompressible flow and a little more for a compressible flow.

Cebeci, T.↗
DOE OSTI2017

FY17 Status Report on the Initial EPP Finite Element Analysis of Grade 91 Steel

This report describes a modification to the elastic-perfectly plastic (EPP) strain limits design method to account for cyclic softening in Gr. 91 steel. The report demonstrates that the unmodified EPP strain limits method described in current ASME code case is not conservative for materials with substantial cyclic softening behavior like Gr. 91 steel. However, the EPP strain limits method can be modified to be conservative for softening materials by using softened isochronous stress-strain curves in place of the standard curves developed from unsoftened creep experiments. The report provides softened curves derived from inelastic material simulations and factors describing the transformation of unsoftened curves to a softened state. Furthermore, the report outlines a method for deriving these factors directly from creep/fatigue tests. If the material softening saturates the proposed EPP strain limits method can be further simplified, providing a methodology based on temperature-dependent softening factors that could be implemented in an ASME code case allowing the use of the EPP strain limits method with Gr. 91. Finally, the report demonstrates the conservatism of the modified method when applied to inelastic simulation results and two bar experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A nonparametric method for producing isolines of bivariate exceedance probabilities

We present a method for drawing isolines indicating regions of equal joint exceedance probability for bivariate data. The method relies on bivariate regular variation, a dependence framework widely used for extremes. The method we utilize for characterizing dependence in the tail is largely nonparametric. The extremes framework enables drawing isolines corresponding to very low exceedance probabilities and may even lie beyond the range of the data; such cases would be problematic for standard nonparametric methods. Furthermore, we extend this method to the case of asymptotic independence and propose a procedure which smooths the transition from hidden regular variation in the interior to the first-order behavior on the axes. We propose a diagnostic plot for assessing the isoline estimate and choice of smoothing, and a bootstrap procedure to visually assess uncertainty.

97 MATHEMATICS AND COMPUTING↗