Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function Estimation”

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,495 records · Page 83

Limits of applicability of a two-temperature model under nonuniform heating of metal by an ultrashort laser pulse

The heating of metals (silver and aluminium) by ultrashort laser pulses is analysed proceeding from a spatially nonuniform kinetic equation for the electron distribution function. The electron subsystem thermalisation is estimated in a wide range of absorbed pulse energy density. The limits of applicability are determined for the two-temperature model. (interaction of laser radiation with matter)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Toward in vivo lung's tissue incompressibility characterization for tumor motion modeling in radiation therapy

Purpose: A novel technique is proposed to characterize lung tissue incompressibility variation during respiration. Estimating lung tissue incompressibility parameter variations resulting from air content variation throughout respiration is critical for computer assisted tumor motion tracking. Continuous tumor motion is a major challenge in lung cancer radiotherapy, especially with external beam radiotherapy. If not accounted for, this motion may lead to areas of radiation overdosage for normal tissue. Given the unavailability of imaging modality that can be used effectively for real-time lung tumor tracking, computer assisted approach based on tissue deformation estimation can be a good alternative. This approach involves lung biomechanical model where its fidelity depends on input tissue properties. This investigation shows that considering variable tissue incompressibility parameter is very important for predicting tumor motion accurately, hence improving the lung radiotherapy outcome. Methods: First, an in silico lung phantom study was conducted to demonstrate the importance of employing variable Poisson's ratio for tumor motion predication. After it was established that modeling this variability is critical for accurate tumor motion prediction, an optimization based technique was developed to estimate lung tissue Poisson's ratio as a function of respiration cycle time. In this technique, the Poisson's ratio and lung pressure value were varied systematically until optimal values were obtained, leading to maximum similarity between acquired and simulated 4D CT lung images. This technique was applied in an ex vivo porcine lung study where simulated images were constructed using the end exhale CT image and deformation fields obtained from the lung's FE modeling of each respiration time increment. To model the tissue, linear elastic and Marlow hyperelastic material models in conjunction with variable Poisson's ratio were used. Results: The phantom study showed that the tumor motion trajectory and its final locations obtained from simulations with and without considering tissue incompressibility variation were very different. For example, tumor displacements in the z direction were -11.23 and -38.10 mm obtained with the Marlow hyperelastic material model in conjunction with constant and variable Poisson's ratio, respectively. By comparing the acquired 4D-CT image sequence of the porcine lung with their image sequence counterparts obtained from the hyperelastic model with constant and variable Poisson's ratio, it was shown that using variable tissue incompressibility reduced errors significantly in tumor motion prediction. Conclusions: This investigation demonstrates the importance of incompressibility variation estimation and utilization for accurate tumor tracking in computer assisted lung external beam radiation therapy. An optimization framework was developed to estimate a Poisson's ratio function in terms of respiration cycle time using experimental image data of the lung. Utilizing this function along with respiratory system FE modeling may lead to more effective tumor targeting, hence potentially improving the outcome of lung external beam radiation therapy techniques. This is particularly true for stereotactic body radiation therapy where only one or a few fraction treatments are applied, precluding the possibility of averaging out dosimetric deviations introduced by the respiratory motion.

60 APPLIED LIFE SCIENCES↗

Estimates of cellular mutagenesis from cosmic rays

A parametric track structure model is used to estimate the cross section as a function of particle velocity and charge for mutations at the hypoxanthine guanine phosphoribosyl transferase (HGPRT) locus in human fibroblast cell cultures. Experiments that report the fraction of mutations per surviving cell for human lung and skin fibroblast cells indicate small differences in the mutation cross section for these two cell lines when differences in inactivation rates between these cell lines are considered. Using models of cosmic ray transport, the mutation rate at the HGPRT locus is estimated for cell cultures in space flight and rates of about 2 to 10 x 10(exp -6) per year are found for typical spacecraft shielding. A discussion of how model assumptions may alter the predictions is also presented.

Cucinotta, Francis A.↗

Vertical Photon Transport in Cloud Remote Sensing Problems

Photon transport in plane-parallel, vertically inhomogeneous clouds is investigated and applied to cloud remote sensing techniques that use solar reflectance or transmittance measurements for retrieving droplet effective radius. Transport is couched in terms of weighting functions which approximate the relative contribution of individual layers to the overall retrieval. Two vertical weightings are investigated, including one based on the average number of scatterings encountered by reflected and transmitted photons in any given layer. A simpler vertical weighting based on the maximum penetration of reflected photons proves useful for solar reflectance measurements. These weighting functions are highly dependent on droplet absorption and solar/viewing geometry. A superposition technique, using adding/doubling radiative transfer procedures, is derived to accurately determine both weightings, avoiding time consuming Monte Carlo methods. Superposition calculations are made for a variety of geometries and cloud models, and selected results are compared with Monte Carlo calculations. Effective radius retrievals from modeled vertically inhomogeneous liquid water clouds are then made using the standard near-infrared bands, and compared with size estimates based on the proposed weighting functions. Agreement between the two methods is generally within several tenths of a micrometer, much better than expected retrieval accuracy. Though the emphasis is on photon transport in clouds, the derived weightings can be applied to any multiple scattering plane-parallel radiative transfer problem, including arbitrary combinations of cloud, aerosol, and gas layers.

Platnick, S.↗

Magnetoacoustic tomography with magnetic induction for high-resolution bioimepedance imaging through vector source reconstruction under the static field of MRI magnet

Purpose: Magnetoacoustic tomography with magnetic induction (MAT-MI) is an imaging modality to reconstruct the electrical conductivity of biological tissue based on the acoustic measurements of Lorentz force induced tissue vibration. This study presents the feasibility of the authors' new MAT-MI system and vector source imaging algorithm to perform a complete reconstruction of the conductivity distribution of real biological tissues with ultrasound spatial resolution. Methods: In the present study, using ultrasound beamformation, imaging point spread functions are designed to reconstruct the induced vector source in the object which is used to estimate the object conductivity distribution. Both numerical studies and phantom experiments are performed to demonstrate the merits of the proposed method. Also, through the numerical simulations, the full width half maximum of the imaging point spread function is calculated to estimate of the spatial resolution. The tissue phantom experiments are performed with a MAT-MI imaging system in the static field of a 9.4 T magnetic resonance imaging magnet. Results: The image reconstruction through vector beamformation in the numerical and experimental studies gives a reliable estimate of the conductivity distribution in the object with a ∼1.5 mm spatial resolution corresponding to the imaging system frequency of 500 kHz ultrasound. In addition, the experiment results suggest that MAT-MI under high static magnetic field environment is able to reconstruct images of tissue-mimicking gel phantoms and real tissue samples with reliable conductivity contrast. Conclusions: The results demonstrate that MAT-MI is able to image the electrical conductivity properties of biological tissues with better than 2 mm spatial resolution at 500 kHz, and the imaging with MAT-MI under a high static magnetic field environment is able to provide improved imaging contrast for biological tissue conductivity reconstruction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SIMULTANEOUS ESTIMATION OF PHOTOMETRIC REDSHIFTS AND SED PARAMETERS: IMPROVED TECHNIQUES AND A REALISTIC ERROR BUDGET

We seek to improve the accuracy of joint galaxy photometric redshift estimation and spectral energy distribution (SED) fitting. By simulating different sources of uncorrected systematic errors, we demonstrate that if the uncertainties in the photometric redshifts are estimated correctly, so are those on the other SED fitting parameters, such as stellar mass, stellar age, and dust reddening. Furthermore, we find that if the redshift uncertainties are over(under)-estimated, the uncertainties in SED parameters tend to be over(under)-estimated by similar amounts. These results hold even in the presence of severe systematics and provide, for the first time, a mechanism to validate the uncertainties on these parameters via comparison with spectroscopic redshifts. We propose a new technique (annealing) to re-calibrate the joint uncertainties in the photo-z and SED fitting parameters without compromising the performance of the SED fitting + photo-z estimation. This procedure provides a consistent estimation of the multi-dimensional probability distribution function in SED fitting + z parameter space, including all correlations. While the performance of joint SED fitting and photo-z estimation might be hindered by template incompleteness, we demonstrate that the latter is “flagged” by a large fraction of outliers in redshift, and that significant improvements can be achieved by using flexible stellar populations synthesis models and more realistic star formation histories. In all cases, we find that the median stellar age is better recovered than the time elapsed from the onset of star formation. Finally, we show that using a photometric redshift code such as EAZY to obtain redshift probability distributions that are then used as priors for SED fitting codes leads to only a modest bias in the SED fitting parameters and is thus a viable alternative to the simultaneous estimation of SED parameters and photometric redshifts.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical Properties for Flows with Unbounded Roof Function, Including the Lorenz Attractor

For geometric Lorenz attractors (including the classical Lorenz attractor) we obtain a greatly simplified proof of the central limit theorem which applies also to the more general class of codimension two singular hyperbolic attractors. We also obtain the functional central limit theorem and moment estimates, as well as iterated versions of these results. A consequence is deterministic homogenisation (convergence to a stochastic differential equation) for fast-slow dynamical systems whenever the fast dynamics is singularly hyperbolic of codimension two.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffractometric method for determining the degree of crystallinity of materials

A new method for determining the degree of crystallinity of a material from X-ray diffraction data has been developed. The method is based on estimating the rate of change in function I = f(2θ) in the entire range of scattering angles. A calculation is performed using the ratio of the integral modulus of the first derivative of intensity with respect to angle 2θ to the integral area under the diffraction pattern curve. The method was tested on two substances with known amorphous and crystalline components. A linear relationship is revealed between the specified ratio of crystalline and amorphous parts and the calculated crystallinity index. The proposed method allows one to estimate impartially and compare the degree of crystallinity for samples of different nature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Diffusion-Model-Assisted Supervised Learning of Generative Models for Density Estimation

Here, we present a supervised learning framework of training generative models for density estimation. Generative models, including generative adversarial networks (GANs), normalizing flows, and variational auto-encoders (VAEs), are usually considered as unsupervised learning models, because labeled data are usually unavailable for training. Despite the success of the generative models, there are several issues with the unsupervised training, e.g., requirement of reversible architectures, vanishing gradients, and training instability. To enable supervised learning in generative models, we utilize the score-based diffusion model to generate labeled data. Unlike existing diffusion models that train neural networks to learn the score function, we develop a training-free score estimation method. This approach uses mini-batch-based Monte Carlo estimators to directly approximate the score function at any spatial-temporal location in solving an ordinary differential equation (ODE), corresponding to the reverse-time stochastic differential equation (SDE). This approach can offer both high accuracy and substantial time savings in neural network training. Once the labeled data are generated, we can train a simple, fully connected neural network to learn the generative model in the supervised manner. Compared with existing normalizing flow models, our method does not require the use of reversible neural networks and avoids the computation of the Jacobian matrix. Compared with existing diffusion models, our method does not need to solve the reverse-time SDE to generate new samples. As a result, the sampling efficiency is significantly improved. We demonstrate the performance of our method by applying it to a set of 2D datasets as well as real data from the University of California Irvine (UCI) repository.

97 MATHEMATICS AND COMPUTING↗

Pre-Commercial Demonstration of Direct Non-Iterative State Estimator (DNSE+)

The objective of this DOE co-funded project is to demonstrate functionality and performance of a production-grade Direct Non-iterative State Estimator (DNSE) through its integration with NYPA’s Energy Management System (EMS) and an enhanced Real Time Dynamics Monitoring System (RTDMS) synchrophasor platform from Electric Power Group (EPG). A production-grade DNSE aims to overcome a major obstacle to the operational use of synchrophasor data in Synchro-Phasor Management Systems (SPMS) by providing SPMS applications with a consistent and a complete synchrophasor data foundation in the same way that a traditional EMS State Estimator (SE) provides to EMS applications. Specifically, the DNSE uses synchrophasor measurements and Supervisory Control and Data Acquisition (SCADA) measurements, together with the complete power system model, to simultaneously obtain the complete state of the entire utility’s operating model at rates that are close to the synchrophasor data rates. The full system coverage property is what distinguishes the DNSE from the linear synchrophasor SEs (which cover only parts of the operating model that are visible through synchrophasor measurements), and DNSE differs from traditional EMS SE in that the solution is obtained without iterations, without initial values (hence the “direct” in the name), and at a much higher execution rate. This presentation provides an overview of the project objectives, a brief description of the main tasks including the development work to productize and optimize DNSE from its current MATLAB code base, RTDMS enhancement for interfacing with the DNSE, factory acceptance testing, field installation, testing and demonstration, as well as the current status of the project execution.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Personalized Assessment of kV Cone Beam Computed Tomography Doses in Image-guided Radiotherapy of Pediatric Cancer Patients

Purpose: To develop a quantitative method for the estimation of kV cone beam computed tomography (kVCBCT) doses in pediatric patients undergoing image-guided radiotherapy. Methods and Materials: Forty-two children were retrospectively analyzed in subgroups of different scanned regions: one group in the head-and-neck and the other group in the pelvis. Critical structures in planning CT images were delineated on an Eclipse treatment planning system before being converted into CT phantoms for Monte Carlo simulations. A benchmarked EGS4 Monte Carlo code was used to calculate three-dimensional dose distributions of kVCBCT scans with full-fan high-quality head or half-fan pelvis protocols predefined by the manufacturer. Based on planning CT images and structures exported in DICOM RT format, occipital-frontal circumferences (OFC) were calculated for head-and-neck patients using DICOMan software. Similarly, hip circumferences (HIP) were acquired for the pelvic group. Correlations between mean organ doses and age, weight, OFC, and HIP values were analyzed with SigmaPlot software suite, where regression performances were analyzed with relative dose differences (RDD) and coefficients of determination (R{sup 2}). Results: kVCBCT-contributed mean doses to all critical structures decreased monotonically with studied parameters, with a steeper decrease in the pelvis than in the head. Empirical functions have been developed for a dose estimation of the major organs at risk in the head and pelvis, respectively. If evaluated with physical parameters other than age, a mean RDD of up to 7.9% was observed for all the structures in our population of 42 patients. Conclusions: kVCBCT doses are highly correlated with patient size. According to this study, weight can be used as a primary index for dose assessment in both head and pelvis scans, while OFC and HIP may serve as secondary indices for dose estimation in corresponding regions. With the proposed empirical functions, it is possible to perform an individualized quantitative dose assessment of kVCBCT scans.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Single transverse spin asymmetry of prompt photon production

We study the single transverse spin asymmetry of prompt photon production in high energy proton–proton scattering. We include the contributions from both the direct and fragmentation photons. While the asymmetry for direct photon production receives only the Sivers type of contribution, the asymmetry for fragmentation photons receives both the Sivers and Collins types of contributions. We make a model calculation for quark-to-photon Collins function, which is then used to estimate the Collins asymmetry for fragmentation photons. We find that the Collins asymmetry for fragmentation photons is very small, thus the single transverse spin asymmetry of prompt photon production is mainly coming from the Sivers asymmetry in direct and fragmentation photons. We make predictions for the prompt photon spin asymmetry at RHIC energy, and emphasize the importance of such a measurement. Furthermore, the asymmetry of prompt photon production can provide a good measurement for the important twist-three quark–gluon correlation function, which is urgently needed in order to resolve the “sign mismatch” puzzle.

Gamberg, Leonard↗

Nuclear kinetic density from ab initio theory

The nuclear kinetic density is one of many fundamental, nonobservable quantities in density functional theory (DFT) dependent on the nonlocal nuclear density. Often, approximations may be made when computing the density that may result in spurious contributions in other DFT quantities. With the ability to compute the nonlocal nuclear density from ab initio wave functions, it is now possible to estimate effects of such spurious contributions. Herein, we derive the kinetic density using ab initio nonlocal scalar one-body nuclear densities computed within the no-core shell model (NCSM) approach, utilizing two- and three-nucleon chiral interactions as the sole input. The ability to compute translationally invariant nonlocal densities allows us to gauge the impact of the spurious center-of-mass (c.m.) contributions in DFT quantities, such as the kinetic density, and provide ab initio insight into refining energy density functionals. The nonlocal nuclear densities are derived from the NCSM one-body densities calculated in second quantization. We present a review of c.m. contaminated and translationally invariant nuclear densities. We then derive an analytic expression for the kinetic density using these nonlocal densities, producing an ab initio kinetic density. The ground-state nonlocal densities of He 4 , 6 , 8 , C 12 , and O 16 are used to compute the kinetic densities of the aforementioned nuclei. The impact of c.m. removal techniques in the density are discussed and compared to a procedure applied in DFT. The results of this work can be extended to other fundamental quantities in DFT. The use of a general nonlocal density allows for the calculation of fundamental quantities taken as input in theories such as DFT. This allows benchmarking c.m. removal procedures and provides a bridge for comparison between ab initio and DFT many-body techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗

Estimates of explosive yield

Explosive yield as a function of propellant type and accident conditions was discussed along with the characteristics of propellant explosions. Three types of fuel-hypergolic, liquid oxygen-hydrocarbon, and liquid oxygen-liquid-hydrogen were considered. Data are also given on three modes (failure modes) of mixing: confinement by missile, confinement by ground surface, and high velocity impact.

Source record↗

Estimation of an Optimal Stimulus Amplitude for Using Vestibular Stochastic Stimulation to Improve Balance Function

Sensorimotor changes such as postural and gait instabilities can affect the functional performance of astronauts when they transition across different gravity environments. We are developing a method, based on stochastic resonance (SR), to enhance information transfer by applying non-zero levels of external noise on the vestibular system (vestibular stochastic resonance, VSR). The goal of this project was to determine optimal levels of stimulation for SR applications by using a defined vestibular threshold of motion detection.

Goel, R.↗

Isopiestic Measurements of Water Adsorption on Uranium Bearing Solids - 1900w

The research focus was to investigate the stability of the multicomponent uranium-bearing precipitates mimicking those created in sediments as a result of alkaline ammonia gas treatment in the vadose zone. Gravimetric measurements of moisture uptake on uranium-bearing precipitates as a function of increasing relative humidity can provide reasonable estimates of the deliquescence point of solid phases or its components. Water moisture uptake by uranium-bearing precipitates as a function of increasing relative humidity was measured at 25 deg. C via the isopiestic method. The conditions for each isotherm were calculated by a division of the amount of adsorbed moisture by the amount of dry precipitate samples. This method can provide estimates of water quantity retained within precipitates conducted in the range of water activities. The desorption portion of the isotherm was obtained by progressively decreasing the water activity or relative pressure in the isopiestic chamber while monitoring the equilibrium amount of moisture sorbed at each relative pressure. The amount of water adsorbed by multicomponent precipitates at a{sub w} = 0.67 ranged between 0.63 and 0.77 g water/g precipitate. Higher water activities or relative humidity values eventually lead to the deliquescence of the solid or its components and the formation of a liquid phase. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Features of ZnO Nanocrystal Size Distribution in a Polymer Matrix

X-ray diffraction patterns of a photocurable epoxy acryl-based polymer composite with addition of ZnO nanocrystals are analyzed. We recover the size distribution function of the ZnO nanocrystals and estimate the median size of the nanocrystals (∼40–50 nm) from a detailed analysis of the most intensive peak in the X-ray diffraction pattern.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗