Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error vector”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

COBE nonspinning attitude propagation

The Cosmic Background Explorer (COBE) spacecraft will exhibit complex attitude motion consisting of a spin rate of approximately -0.8 revolution per minute (rpm) about the x-axis and simultaneous precession of the spin axis at a rate of one revolution per orbit (rpo) about the nearly perpendicular spacecraft-to-Sun vector. The effect of the combined spinning and precession is to make accurate attitude propagation difficult and the 1-degree (3 sigma) solution accuracy goal problematic. To improve this situation, an intermediate reference frame is introduced, and the angular velocity divided into two parts. The nonspinning part is that which would be observed if there were no rotation about the X-axis. The spinning part is simply the X-axis component of the angular velocity. The two are propagated independently and combined whenever the complete attitude is needed. This approach is better than the usual one-step method because each of the two angular velocities look nearly constant in their respective reference frames. Since the angular velocities are almost constant, the approximations made in discrete time propagation are more nearly true. To demonstrate the advantages of this nonspinning method, attitude is propagated as outlined above and is then compared with the results of the one-step method. Over the 100-minute COBE orbit, the one-step error grows to several degrees while the nonspinning error remains negligible.

Chu, D.↗
NASA NTRS2016

Polarized Radiative Transfer of a Cirrus Cloud Consisting of Randomly Oriented Hexagonal Ice Crystals: The 3 x 3 Approximation for Non-Spherical Particles

The reflection and transmission of polarized light for a cirrus cloud consisting of randomly oriented hexagonal columns were calculated by two very different vector radiative transfer models. The forward peak of the phase function for the ensemble-averaged ice crystals has a value of order 6 x 10(exp 3) so a truncation procedure was used to help produce numerically efficient yet accurate results. One of these models, the Vectorized Line-by-Line Equivalent model (VLBLE), is based on the doubling- adding principle, while the other is based on a vector discrete ordinates method (VDISORT). A comparison shows that the two models provide very close although not entirely identical results, which can be explained by differences in treatment of single scattering and the representation of the scattering phase matrix. The relative differences in the reflected I and Q Stokes parameters are within 0.5 for I and within 1.5 for Q for all viewing angles. In 1971 Hansen showed that for scattering by spherical particles the 3 x 3 approximation is sufficient to produce accurate results for the reflected radiance I and the degree of polarization (DOP), and he conjectured that these results would hold also for non-spherical particles. Simulations were conducted to test Hansen's conjecture for the cirrus cloud particles considered in this study. It was found that the 3 x 3 approximation also gives accurate results for the transmitted light, and for Q and U in addition to I and DOP. For these non-spherical ice particles the 3 x 3 approximation leads to an absolute error 2 x 10(exp -6) for the reflected and transmitted I, Q and U Stokes parameters. Hence, it appears to be an excellent approximation, which significantly reduces the computational complexity and burden required for multiple scattering calculations.

Polarization↗
NASA NTRS2004

Recent Developments and Status of the Langley Single Vector Balance Calibration System (SVS)

The Langley Research Center (LaRC) Single Vector Balance Calibration System (SVS) was first introduced in 2000 by Peter Parker. The SVS combines the Design of Experiments (DOE) methodology with a novel load application system. Since that time three systems have been designed and developed with different load range capabilities (ranging from 2 pounds to 3,000 pounds). Approximately fifteen balances have been calibrated and their data compared to conventional techniques. This paper will present results of these comparisons, based on the mathematical models and accuracies, and discuss differences that were observed. In addition, changes in the implementation of the initial load schedules developed using DOE will be highlighted. One of the principles behind DOE is randomization. The initial loading schedules used to date have been randomized in the traditional DOE sense but not for repeat calibrations or experiments. Implementation of this randomization within blocks and its impact on data quality will be reviewed. Areas of potential future development will be presented which include changes in the centers to include loads with the force position system in the pure error estimates.

Jones, Shirley M.↗
DOE OSTI2014

Noise Decoherence and Errors from Entanglement-function Theory for Quantum Computing

A significant problem in quantum computing is the development of physical realizations of algorithms that are robust against noise. One way to examine and mitigate noise would be to simulate large sets of qubits coupling to the external environment on classical computers. This is extremely challenging as quantum information processing is in some sense tied to computing resources that scale exponentially with the number of computing elements (qubits). In this LDRD, we set the foundation for a computational framework potentially allowing simulations of 1000s of qubits vs. 10s now possible. Exact wave-function-based methods demand exponentially increasing resources with system size. The method proposed, entanglement-functional theory (EFT), requires vastly fewer resources. The crucial step is to map the information contained in the wave-functions into a simpler object with associated 1.) auxiliary gate operations and 2.) entanglement functionals of this object. This is similar to the Time-dependent Density Functional Theory (TDDFT) approach that has revolutionized chemistry and materials science. Instead of dealing with the exponentially large wave-function, EFT works with a polynomially large set of projections (the density) that are easily manipulated through unitary operations. For a given set of quantum gates, an isomorphism exists that relates the sequence of events to the time-dependent density. A system of entangled qubits can be simulated at drastically reduced cost relative to existing state-of-the-art vector-state simulation codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

TH-A-BRF-08: Deformable Registration of MRI and CT Images for MRI-Guided Radiation Therapy

Purpose: To evaluate the quality of a commercially available MRI-CT image registration algorithm and then develop a method to improve the performance of this algorithm for MRI-guided prostate radiotherapy. Methods: Prostate contours were delineated on ten pairs of MRI and CT images using Eclipse. Each pair of MRI and CT images was registered with an intensity-based B-spline algorithm implemented in Velocity. A rectangular prism that contains the prostate volume was partitioned into a tetrahedral mesh which was aligned to the CT image. A finite element method (FEM) was developed on the mesh with the boundary constraints assigned from the Velocity generated displacement vector field (DVF). The resultant FEM displacements were used to adjust the Velocity DVF within the prism. Point correspondences between the CT and MR images identified within the prism could be used as additional boundary constraints to enforce the model deformation. The FEM deformation field is smooth in the interior of the prism, and equal to the Velocity displacements at the boundary of the prism. To evaluate the Velocity and FEM registration results, three criteria were used: prostate volume conservation and center consistence under contour mapping, and unbalanced energy of their deformation maps. Results: With the DVFs generated by the Velocity and FEM simulations, the prostate contours were warped from MRI to CT images. With the Velocity DVFs, the prostate volumes changed 10.2% on average, in contrast to 1.8% induced by the FEM DVFs. The average of the center deviations was 0.36 and 0.27 cm, and the unbalance energy was 2.65 and 0.38 mJ/cc3 for the Velocity and FEM registrations, respectively. Conclusion: The adaptive FEM method developed can be used to reduce the error of the MIbased registration algorithm implemented in Velocity in the prostate region, and consequently may help improve the quality of MRI-guided radiation therapy.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗
DOE OSTI2026

Leveraging machine learning to enhance aerosol classification using Single-Particle Mass Spectrometry

Advancing automated classification of atmospheric aerosols from Single-Particle Mass Spectrometry (SPMS) data remains challenging due to overlapping ion signatures, compositional diversity, and limited labeled data. This study evaluates supervised and semi-supervised learning frameworks to enhance aerosol identification by jointly leveraging labeled and unlabeled spectra. Four models were compared: a supervised Support Vector Machine (SVM), a self-training SVM, a stacked autoencoder classifier, and a stacked autoencoder trained using a temporal-ensembling Mean Teacher approach. All models achieved high and stable accuracies (90.0 %–91.1 %), surpassing previous results on the same dataset (87 %) and matching the performance of state-of-the-art deep learning methods. Despite small global metric differences (≤ 1 %), semi-supervised variants yielded up to 5 %–10 % improvements for compositionally rare particle types – such as soot (0.77 % of spectra, F1-score: 0.93–0.97) and hazelnut pollen (0.98 % of spectra, F1-score: 0.97–1.00) – equating to roughly ∼ 187 additional correctly classified spectra. These gains are scientifically significant, as such rare particles exert disproportionate influence on radiative absorption and ice nucleation processes; their improved detection reduces modeled uncertainties in aerosol absorption optical depth and mixed-phase cloud ice nucleation rates. The models' residual misclassifications (≈ 9 %) largely arise from true spectral overlap among chemically adjacent species (e.g., Na- vs. K-feldspar, coated vs. uncoated feldspars), reflecting physical compositional continuity rather than algorithmic error. Collectively, these findings demonstrate that leveraging unlabeled data to learn robust spectral representations and refine classification enhances both fidelity and interpretability, bridging data-driven analysis with aerosol–climate process understanding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Overview of new MAST physics in anticipation of first results from MAST Upgrade

We report that the mega amp spherical tokamak (MAST) was a low aspect ratio device (R/a = 0.85/0.65 ~ 1.3) with similar poloidal cross-section to other medium-size tokamaks. The physics programme concentrates on addressing key physics issues for the operation of ITER, design of DEMO and future spherical tokamaks by utilising high resolution diagnostic measurements closely coupled with theory and modelling to significantly advance our understanding. An empirical scaling of the energy confinement time that favours higher power, lower collisionality devices is consistent with gyrokinetic modelling of electron scale turbulence. Measurements of ion scale turbulence with beam emission spectroscopy and gyrokinetic modelling in up-down symmetric plasmas find that the symmetry of the turbulence is broken by flow shear. Near the non-linear stability threshold, flow shear tilts the density fluctuation correlation function and skews the fluctuation amplitude distribution. Results from fast particle physics studies include the observation that sawteeth are found to redistribute passing and trapped fast particles injected from neutral beam injectors in equal measure, suggesting that resonances between the m = 1 perturbation and the fast ion orbits may be playing a dominant role in the fast ion transport. Measured D–D fusion products from a neutron camera and a charged fusion product detector are 40% lower than predictions from TRANSP/NUBEAM, highlighting possible deficiencies in the guiding centre approximation. Modelling of fast ion losses in the presence of resonant magnetic perturbations (RMPs) can reproduce trends observed in experiments when the plasma response and charge-exchange losses are accounted for. Measurements with a neutral particle analyser during merging-compression start-up indicate the acceleration of ions and electrons. Transport at the plasma edge has been improved through reciprocating probe measurements that have characterised a geodesic acoustic mode at the edge of an ohmic L-mode plasma and particle-in-cell modelling has improved the interpretation of plasma potential estimates from ball-pen probes. The application of RMPs leads to a reduction in particle confinement in L-mode and H-mode and an increase in the core ionization source. The ejection of secondary filaments following type-I ELMs correlates with interactions with surfaces near the X-point. Simulations of the interaction between pairs of filaments in the scrape-off layer suggest this results in modest changes to their velocity, and in most cases can be treated as moving independently. A stochastic model of scrape-off layer profile formation based on the superposition of non-interacting filaments is in good agreement with measured time-average profiles. Transport in the divertor has been improved through fast camera imaging, indicating the presence of a quiescent region devoid of filament near the X-point, extending from the separatrix to $ψ_n$ ~ 1.02. Simulations of turbulent transport in the divertor show that the angle between the divertor leg on the curvature vector strongly influences transport into the private flux region via the interchange mechanism. Coherence imaging measurements show counter-streaming flows of impurities due to gas puffing increasing the pressure on field lines where the gas is ionised. MAST Upgrade is based on the original MAST device, with substantially improved capabilities to operate with a Super-X divertor to test extended divertor leg concepts. SOLPS-ITER modelling predicts the detachment threshold will be reduced by more than a factor of 2, in terms of upstream density, in the Super-X compared with a conventional configuration and that the radiation front movement is passively stabilised before it reaches the X-point. 1D fluid modelling reveals the key role of momentum and power loss mechanisms in governing detachment onset and evolution. Analytic modelling indicates that long legs placed at large major radius, or equivalently low at the target compared with the X-point are more amenable to external control. With MAST Upgrade experiments expected in 2019, a thorough characterisation of the sources of the intrinsic error field has been carried out and a mitigation strategy developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A hybrid reconstruction algorithm for fast and accurate 4D cone-beam CT imaging

Purpose: 4D cone beam CT (4D-CBCT) has been utilized in radiation therapy to provide 4D image guidance in lung and upper abdomen area. However, clinical application of 4D-CBCT is currently limited due to the long scan time and low image quality. The purpose of this paper is to develop a new 4D-CBCT reconstruction method that restores volumetric images based on the 1-min scan data acquired with a standard 3D-CBCT protocol. Methods: The model optimizes a deformation vector field that deforms a patient-specific planning CT (p-CT), so that the calculated 4D-CBCT projections match measurements. A forward-backward splitting (FBS) method is invented to solve the optimization problem. It splits the original problem into two well-studied subproblems, i.e., image reconstruction and deformable image registration. By iteratively solving the two subproblems, FBS gradually yields correct deformation information, while maintaining high image quality. The whole workflow is implemented on a graphic-processing-unit to improve efficiency. Comprehensive evaluations have been conducted on a moving phantom and three real patient cases regarding the accuracy and quality of the reconstructed images, as well as the algorithm robustness and efficiency. Results: The proposed algorithm reconstructs 4D-CBCT images from highly under-sampled projection data acquired with 1-min scans. Regarding the anatomical structure location accuracy, 0.204 mm average differences and 0.484 mm maximum difference are found for the phantom case, and the maximum differences of 0.3–0.5 mm for patients 1–3 are observed. As for the image quality, intensity errors below 5 and 20 HU compared to the planning CT are achieved for the phantom and the patient cases, respectively. Signal-noise-ratio values are improved by 12.74 and 5.12 times compared to results from FDK algorithm using the 1-min data and 4-min data, respectively. The computation time of the algorithm on a NVIDIA GTX590 card is 1–1.5 min per phase. Conclusions: High-quality 4D-CBCT imaging based on the clinically standard 1-min 3D CBCT scanning protocol is feasible via the proposed hybrid reconstruction algorithm.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

First-Principles Quantum Transport Modeling of Spin-Transfer and Spin-Orbit Torques in Magnetic Multilayers

A unified approach for computing (i) spin-transfer torque in magnetic trilayers like spin valve and magnetic tunnel junction, where injected charge current flows perpendicularly to interfaces, and (ii) spin-orbit torque in magnetic bilayers of the type ferromagnet/spin-orbit-coupled material, where injected charge current flows parallel to the interface, was reviewed. The experimentally explored and technologically relevant spin-orbit-coupled materials include 5d heavy metals, topological insulators, Weyl semimetals, and transition metal dichalcogenides. This approach requires to construct the torque operator for a given Hamiltonian of the device and the steady-state nonequilibrium density matrix, where the latter is expressed in terms of the nonequilibrium Green’s functions and split into three contributions. Tracing these contributions with the torque operator automatically yields field-like and damping-like components of spin-transfer torque or spin-orbit torque vector, which is particularly advantageous for spin-orbit torque where the direction of these components depends on the unknown-in-advance orientation of the current-driven nonequilibrium spin density in the presence of spin-orbit coupling. Illustrative examples are provided by computing spin-transfer torque in a one-dimensional toy model of a magnetic tunnel junction and realistic Co/Cu/Co spin valve, both of which are described by first-principles Hamiltonians obtained from noncollinear density functional theory calculations, as well as by computing spin-orbit torque in a ferromagnetic layer described by a tight-binding Hamiltonian which includes spin-orbit proximity effect within ferromagnetic monolayers assumed to be generated by the adjacent monolayer transition metal dichalcogenide. In addition, it is shown here how spin-orbit proximity effect, quantified by computing (via first-principles retarded Green’s function) spectral functions and spin textures on monolayers of realistic ferromagnetic material like Co in contact with heavy metal or monolayer transition metal dichalcogenide, can be tailored to enhance the magnitude of spin-orbit torque. Furthermore, errors made in the calculation of spin-transfer torque are quantified when using Hamiltonian from collinear density functional theory, with rigidly rotated magnetic moments to create noncollinear magnetization configurations, instead of proper (but computationally more expensive) self-consistent Hamiltonian obtained from noncollinear density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Online adaptive basis refinement and compression for reduced-order models via vector-space sieving

In many applications, projection-based reduced-order models (ROMs) have demonstrated the ability to provide rapid approximate solutions to high-fidelity full-order models (FOMs). However, there is no a priori assurance that these approximate solutions are accurate; their accuracy depends on the ability of the low-dimensional trial basis to represent the FOM solution. As a result, ROMs can generate inaccurate approximate solutions, e.g., when the FOM solution at the online prediction point is not well represented by training data used to construct the trial basis. To address this fundamental deficiency of standard model-reduction approaches, this work proposes a novel online-adaptive mechanism for efficiently enriching the trial basis in a manner that ensures convergence of the ROM to the FOM, yet does not incur any FOM solves. The mechanism is based on the previously proposed adaptive $h$-refinement method for ROMs, but improves upon this work in two crucial ways. First, the proposed method enables basis refinement with respect to any orthogonal basis (not just the Kronecker basis), thereby generalizing the refinement mechanism and enabling it to be tailored to the physics characterizing the problem at hand. Second, the proposed method provides a fast online algorithm for periodically compressing the enriched basis via an efficient proper orthogonal decomposition (POD) method, which does not incur any operations that scale with the FOM dimension. We report these two features allow the proposed method to serve as (1) a failsafe mechanism for ROMs, as the method enables the ROM to satisfy any prescribed error tolerance online (even in the case of inadequate training), and (2) an efficient online basis-adaptation mechanism, as the combination of basis enrichment and compression enables the basis to adapt online while controlling its dimension.

42 ENGINEERING↗
DOE OSTI2025

Bounding entanglement entropy with Clifford double cosets

Following on our previous work studying the orbits of quantum states under Clifford circuits via reachability graphs, we introduce contracted graphs whose vertices represent classes of quantum states with the same entropy vector. These contracted graphs represent the double cosets of the Clifford group, where the left cosets are built from the stabilizer subgroup of the starting state and the right cosets are built from the entropy-preserving operators. We study contracted graphs for stabilizer states, as well as 𝑊 states and Dicke states, discussing how the diameter of a state's contracted graph constrains the entropic diversity of its two-qubit Clifford orbit. We derive an upper bound on the number of entropy vectors that can be generated using any 𝑛-qubit Clifford circuit, for any quantum state. Here, we speculate on the holographic implications for the relative proximity of gravitational duals of states within the same Clifford orbit. Although we concentrate on how entropy evolves under the Clifford group, our double-coset formalism, and thus the contracted graph picture, is extendable to generic gate sets and generic state properties.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Numerically approximating centrality for graph ranking guarantees

Many real-world datasets can be represented as graphs. Using iterative solvers to approximate graph centrality measures allows us to obtain a ranking vector on the nodes of the graph, consisting of a number for each vertex in the graph identifying its relative importance. In this study the centrality measures we use are Katz Centrality and PageRank. Given an approximate solution, we use the residual to accurately estimate how much of the ranking matches the ranking given by the exact solution. Using probabilistic matrix norms, we obtain bounds on the accuracy of the approximation compared to the exact solution with respect to the highly ranked nodes and apply numerical analysis to the computation of centrality with iterative methods. This relates the numerical accuracy of the linear solver to the data analysis accuracy of finding the correct ranking. In particular, we answer the question of which pairwise rankings are reliable given an approximate solution to the linear system. Experiments on many real-world undirected and directed networks up to several million vertices and several hundred million edges validate our theory and show that we are able to accurately estimate large portions of the approximation. We also analyze the difference between global centrality scores and personalized scores (w.r.t. specific seed vertices). By analyzing convergence error, we develop confidence in the ranking schemes of data mining. We show we are able to accurately guarantee ranking of vertices with an approximation to centrality metrics faster than current methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2010

Dalitz Plot Analysis of B{sup -{yields}}D{sup +{pi}-{pi}-}

We present a Dalitz plot analysis of B{sup -{yields}}D{sup +{pi}-{pi}-} decays, based on a sample of about 383 million {Upsilon}(4S){yields}BB-bar decays collected by the BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC. The analysis has been published previously in [1]. We measure the inclusive branching fraction of the three-body decay to be B(B{sup -{yields}}D{sup +{pi}+{pi}-}) = (1.08{+-}0.03(stat){+-}0.05(syst))x10{sup -3}. We observe the established D{sub 2}*{sup 0} and confirm the existence of D{sub 0}*{sup 0} in their decays to D{sup +{pi}-}, where the D{sub 2}*{sup 0} and D{sub 0}*{sup 0} are the 2{sup +} and 0{sup +} cu-bar P-wave states, respectively. We measure the masses and widths of D{sub 2}*{sup 0} and {Gamma}D{sub 0}*{sup 0} to be: mD{sub 2}* = (2460.4{+-}1.2{+-}1.2{+-}1.9) MeV/c{sup 2}, {Gamma}D{sub 2}* = (41.8{+-}2.5{+-}2.1{+-}2.0) MeV, mD{sub 0}* = (2297{+-}8{+-}5{+-}19) MeV/c{sup 2}, {Gamma}D{sub 0}* = (273{+-}12{+-}17{+-}45) MeV. The stated errors reflect the statistical and systematic uncertainties, and the uncertainty related to the assumed composition of signal events and the theoretical model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

WE-G-BRF-04: Robust Real-Time Volumetric Imaging Based On One Single Projection

Purpose: Real-time volumetric imaging is highly desirable to provide instantaneous image guidance for lung radiation therapy. This study proposes a scheme to achieve this goal using one single projection by utilizing sparse learning and a principal component analysis (PCA) based lung motion model. Methods: A patient-specific PCA-based lung motion model is first constructed by analyzing deformable vector fields (DVFs) between a reference image and 4DCT images at each phase. At the training stage, we “learn” the relationship between the DVFs and the projection using sparse learning. Specifically, we first partition the projections into patches, and then apply sparse learning to automatically identify patches that best correlate with the principal components of the DVFs. Once the relationship is established, at the application stage, we first employ a patchbased intensity correction method to overcome the problem of different intensity scale between the calculated projection in the training stage and the measured projection in the application stage. The corrected projection image is then fed to the trained model to derive a DVF, which is applied to the reference image, yielding a volumetric image corresponding to the projection. We have validated our method through a NCAT phantom simulation case and one experiment case. Results: Sparse learning can automatically select those patches containing motion information, such as those around diaphragm. For the simulation case, over 98% of the lung region pass the generalized gamma test (10HU/1mm), indicating combined accuracy in both intensity and spatial domain. For the experimental case, the average tumor localization errors projected to the imager are 0.68 mm and 0.4 mm on the axial and tangential direction, respectively. Conclusion: The proposed method is capable of accurately generating a volumetric image using one single projection. It will potentially offer real-time volumetric image guidance to facilitate lung radiotherapy. This study is supported in part by NIH (1R01CA154747-01), The Core Technology Research in Strategic Emerging Industry, Guangdong, China(2011A081402003)

60 APPLIED LIFE SCIENCES↗
NASA NTRS2005

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗
NASA NTRS1998

Experiment 1: Orbital Processing of High-Quality Zn-Alloyed CdTe Compound Semiconductors

Zinc-alloyed cadmium telluride (Zn:CdTe) crystals were grown in unit gravity and in microgravity for comparative analysis. Two crystals were grown on USML-1/STS-50 in 1992, and two additional crystals were grown on USML-2/STS-73 in 1995. The Crystal Growth Furnace (CGF) in the seeded Bridgman-Stockbarger crystal growth geometry was utilized on both missions. Crystals grown on USML-1/STS-50 were found to have solidified with partial wall contact due to the near-absence of the hydrostatic pressure in microgravity, a residual g-vector that was not axial, and the nonwetting sample/ampoule wetting conditions. Crystals grown on USML-2/STS-73 included: a sample/ampoule identical to the USML-1/STS-73 sample/ampoule, with the addition of a restraining spring to simulate hydrostatic pressure internally, and a tapered ampoule which accomplished 2.2 cm of crystal growth without wall contact. Infrared transmission of all ground and flight samples was found to be 63-66 percent, very close to the theoretical 66 percent, suggesting good stoichiometric control. Infrared microscopy confirmed that the primary precipitates were Te and their size (1-10 microns) and density suggested that the flight- and ground-base samples experienced similar thermal histories. Longitudinal macrosegregation, calculated using scaling analysis, was predicted to be low. Nearly diffusion controlled growth was achieved even in unit gravity and macrosegregation data could be fit with a diffusion controlled model. Radial segregation was monitored and was found to vary with fraction solidified, particularly through the shoulder region, where the sample cross section was varying significantly. It was also disturbed in the flight samples in regions where asymmetric wall contact was noted. In regions where a steady state was established, the radial segregation was invariant within our experimental measurement error. Flight samples were found to be much higher in structural perfection than samples processed in unit gravity under identical growth conditions. In regions where solidification had occurred without wall contact, the free surfaces evidenced virtually no twinning, although twins appeared in the flight samples in regions of wall contact and were pervasive in the ground samples. These results were confirmed using optical microscopy and synchrotron x-ray white beam topography. Full-width half-maximum rocking curve widths (FWHM), recorded in arc-seconds, were significantly reduced from 20 a-s (1 g) to 9 a-s (micro-g) for the best regions of the crystals. The 9 a-s (FWHM) rocking curve value in the unconfined flight samples equals the best value reported terrestrially for this material. The ground samples exhibited a fully developed (111)[110] dislocation mosaic structure, whereas dislocations within the flight samples were discrete and no mosaic structure was evident. The defect density was quantitatively reduced from 75,000 (1 g) to 800 (micro-g) +/-50 percent. Dislocation etch pit density results were confirmed using transmission synchrotron white beam and monochromated beam topography. The low defect density is thought to have resulted from the near absence of hydrostatic pressure, which allowed the molten boule to solidify with little or no wall contact. This minimized the transfer of hoop stresses during solidification and post-solidification processing.

Larson, D. J., Jr.↗
NASA NTRS2002

The 1999 Hector Mine Earthquake, Southern California: Vector Near-Field Displacements from ERS InSAR

Two components of fault slip are uniquely determined from two line-of-sight (LOS) radar interferograms by assuming that the fault-normal component of displacement is zero. We use this approach with ascending and descending interferograms from the ERS satellites to estimate surface slip along the Hector Mine earthquake rupture. The LOS displacement is determined by visually counting fringes to within 1 km of the outboard ruptures. These LOS estimates and uncertainties are then transformed into strike- and dip-slip estimates and uncertainties; the transformation is singular for a N-S oriented fault and optimal for an E-W oriented fault. In contrast to our previous strike-slip estimates, which were based only on a descending interferogram, we now find good agreement with the geological measurements, except at the ends of the rupture. The ascending interferogram reveals significant west-sidedown dip-slip (approximately 1.0 m) which reduces the strike-slip estimates by 1 to 2 m, especially along the northern half of the rupture. A spike in the strike-slip displacement of 6 m is observed in central part of the rupture. This large offset is confirmed by subpixel cross correlation of features in the before and after amplitude images. In addition to strike slip and dip slip, we identify uplift and subsidence along the fault, related to the restraining and releasing bends in the fault trace, respectively. Our main conclusion is that at least two look directions are required for accurate estimates of surface slip even along a pure strike-slip fault. Models and results based only on a single look direction could have major errors. Our new estimates of strike slip and dip slip along the rupture provide a boundary condition for dislocation modeling. A simple model, which has uniform slip to a depth of 12 km, shows good agreement with the observed ascending and descending interferograms.

Sandwell, David T.↗