Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Soft computing”

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 685 records · Page 38

SU-E-I-74: Image-Matching Technique of Computed Tomography Images for Personal Identification: A Preliminary Study Using Anthropomorphic Chest Phantoms

Purpose: Fingerprints, dental impressions, and DNA are used to identify unidentified bodies in forensic medicine. Cranial Computed tomography (CT) images and/or dental radiographs are also used for identification. Radiological identification is important, particularly in the absence of comparative fingerprints, dental impressions, and DNA samples. The development of an automated radiological identification system for unidentified bodies is desirable. We investigated the potential usefulness of bone structure for matching chest CT images. Methods: CT images of three anthropomorphic chest phantoms were obtained on different days in various settings. One of the phantoms was assumed to be an unidentified body. The bone image and the bone image with soft tissue (BST image) were extracted from the CT images. To examine the usefulness of the bone image and/or the BST image, the similarities between the two-dimensional (2D) or threedimensional (3D) images of the same and different phantoms were evaluated in terms of the normalized cross-correlation value (NCC). Results: For the 2D and 3D BST images, the NCCs obtained from the same phantom assumed to be an unidentified body (2D, 0.99; 3D, 0.93) were higher than those for the different phantoms (2D, 0.95 and 0.91; 3D, 0.89 and 0.80). The NCCs for the same phantom (2D, 0.95; 3D, 0.88) were greater compared to those of the different phantoms (2D, 0.61 and 0.25; 3D, 0.23 and 0.10) for the bone image. The difference in the NCCs between the same and different phantoms tended to be larger for the bone images than for the BST images. These findings suggest that the image-matching technique is more useful when utilizing the bone image than when utilizing the BST image to identify different people. Conclusion: This preliminary study indicated that evaluating the similarity of bone structure in 2D and 3D images is potentially useful for identifying of an unidentified body.

60 APPLIED LIFE SCIENCES↗

In situ tumor model for longitudinal in silico imaging trials

Abstract Objective.In this article, we introduce a computational model for simulating the growth of breast cancer lesions accounting for the stiffness of surrounding anatomical structures.Approach.In our model, ligaments are classified as the most rigid structures while the softer parts of the breast are occupied by fat and glandular tissues As a result of these variations in tissue elasticity, the rapidly proliferating tumor cells are met with differential resistance. It is found that these cells are likely to circumvent stiffer terrains such as ligaments, instead electing to proliferate preferentially within the more yielding confines of the breast’s soft topography. By manipulating the interstitial tumor pressure in direct proportion to the elastic constants of the tissues surrounding the tumor, this model thus creates the potential for realizing a database of unique lesion morphology sculpted by the distinctive topography of each local anatomical infrastructure. We modeled the growth of simulated lesions within volumes extracted from fatty breast models, developed by Graffet alwith a resolution of 50μm generated with the open-source and readily available Virtual Imaging Clinical Trials for Regulatory Evaluation (VICTRE) imaging pipeline. To visualize and validate the realism of the lesion models, we leveraged the imaging component of the VICTRE pipeline, which replicates the siemens mammomat inspiration mammography system in a digital format. This system was instrumental in generating digital mammogram (DM) images for each breast model containing the simulated lesions.Results.By utilizing the DM images, we were able to effectively illustrate the imaging characteristics of the lesions as they integrated with the anatomical backgrounds. Our research also involved a reader study that compared 25 simulated DM regions of interest (ROIs) with inserted lesions from our models with DM ROIs from the DDSM dataset containing real manifestations of breast cancer. In general the simulation time for the lesions was approximately 2.5 hours, but it varied depending on the lesion’s local environment.Significance.The lesion growth model will facilitate and enhance longitudinal in silico trials investigating the progression of breast cancer.

Engineering↗

Computerized analysis of coronary artery disease: Performance evaluation of segmentation and tracking of coronary arteries in CT angiograms

Purpose: The authors are developing a computer-aided detection system to assist radiologists in analysis of coronary artery disease in coronary CT angiograms (cCTA). This study evaluated the accuracy of the authors’ coronary artery segmentation and tracking method which are the essential steps to define the search space for the detection of atherosclerotic plaques. Methods: The heart region in cCTA is segmented and the vascular structures are enhanced using the authors’ multiscale coronary artery response (MSCAR) method that performed 3D multiscale filtering and analysis of the eigenvalues of Hessian matrices. Starting from seed points at the origins of the left and right coronary arteries, a 3D rolling balloon region growing (RBG) method that adapts to the local vessel size segmented and tracked each of the coronary arteries and identifies the branches along the tracked vessels. The branches are queued and subsequently tracked until the queue is exhausted. With Institutional Review Board approval, 62 cCTA were collected retrospectively from the authors’ patient files. Three experienced cardiothoracic radiologists manually tracked and marked center points of the coronary arteries as reference standard following the 17-segment model that includes clinically significant coronary arteries. Two radiologists visually examined the computer-segmented vessels and marked the mistakenly tracked veins and noisy structures as false positives (FPs). For the 62 cases, the radiologists marked a total of 10191 center points on 865 visible coronary artery segments. Results: The computer-segmented vessels overlapped with 83.6% (8520/10191) of the center points. Relative to the 865 radiologist-marked segments, the sensitivity reached 91.9% (795/865) if a true positive is defined as a computer-segmented vessel that overlapped with at least 10% of the reference center points marked on the segment. When the overlap threshold is increased to 50% and 100%, the sensitivities were 86.2% and 53.4%, respectively. For the 62 test cases, a total of 55 FPs were identified by radiologist in 23 of the cases. Conclusions: The authors’ MSCAR-RBG method achieved high sensitivity for coronary artery segmentation and tracking. Studies are underway to further improve the accuracy for the arterial segments affected by motion artifacts, severe calcified and noncalcified soft plaques, and to reduce the false tracking of the veins and other noisy structures. Methods are also being developed to detect coronary artery disease along the tracked vessels.

60 APPLIED LIFE SCIENCES↗

SU-G-206-05: A Comparison of Head Phantoms Used for Dose Determination in Imaging Procedures

Purpose: To determine similarities and differences between various head phantoms that might be used for dose measurements in diagnostic imaging procedures. Methods: We chose four frequently used anthropomorphic head phantoms (SK-150, PBU-50, RS-240T and Alderson Rando), a computational patient phantom (Zubal) and the CTDI head phantom for comparison in our study. We did a CT scan of the head phantoms using the same protocol and compared their dimensions and CT numbers. The scan data was used to calculate dose values for each of the phantoms using EGSnrc Monte Carlo software. An .egsphant file was constructed to describe these phantoms using a Visual C++ program for DOSXYZnrc/EGSnrc simulation. The lens dose was calculated for a simulated CBCT scan using DOSXYZnrc/EGSnrc and the calculated doses were validated with measurements using Gafchromic film and an ionization chamber. Similar calculations and measurements were made for PA radiography to investigate the attenuation and backscatter differences between these phantoms. We used the Zubal phantom as the standard for comparison since it was developed based on a CT scan of a patient. Results: The lens dose for the Alderson Rando phantom is around 9% different than the Zubal phantom, while the lens dose for the PBU-50 phantom was about 50% higher, possibly because its skull thickness and the density of bone and soft tissue are lower than anthropometric values. The lens dose for the CTDI phantom is about 500% higher because of its totally different structure. The entrance dose profiles are similar for the five anthropomorphic phantoms, while that for the CTDI phantom was distinctly different. Conclusion: The CTDI and PBU-50 head phantoms have substantially larger lens dose estimates in CBCT. The other four head phantoms have similar entrance dose with backscatter hence should be preferred for dose measurement in imaging procedures of the head. Partial support from NIH Grant R01-EB002873 and Toshiba Medical Systems Corp.

60 APPLIED LIFE SCIENCES↗

The seven-gluon amplitude in multi-Regge kinematics beyond leading logarithmic accuracy

We present an all-loop dispersion integral, well-defined to arbitrary logarithmic accuracy, describing the multi-Regge limit of the 2 → 5 amplitude in planar N = 4 super Yang-Mills theory. It follows from factorization, dual conformal symmetry and consistency with soft limits, and specifically holds in the region where the energies of all produced particles have been analytically continued. After promoting the known symbol of the 2-loop N -particle MHV amplitude in this region to a function, we specialize to N = 7, and extract from it the next-to-leading order (NLO) correction to the BFKL central emission vertex, namely the building block of the dispersion integral that had not yet appeared in the well-studied six-gluon case. As an application of our results, we explicitly compute the seven-gluon amplitude at next-to-leading logarithmic accuracy through 5 loops for the MHV case, and through 3 and 4 loops for the two independent NMHV helicity configurations, respectively.

1/N expansion↗

Revisiting the multi-monopole point of SU( N ) $\mathcal{N}$ = 2 gauge theory in four dimensions

Motivated by applications to soft supersymmetry breaking, we revisit the expansion of the Seiberg-Witten solution around the multi-monopole point on the Coulomb branch of pure SU(N) $\mathcal{N}$ = 2 gauge theory in four dimensions. At this point N - 1 mutually local magnetic monopoles become massless simultaneously, and in a suitable duality frame the gauge couplings logarithmically run to zero. We explicitly calculate the leading threshold corrections to this logarithmic running from the Seiberg-Witten solution by adapting a method previously introduced by D’Hoker and Phong. We compare our computation to existing results in the literature; this includes results specific to SU(2) and SU(3) gauge theories, the large-N results of Douglas and Shenker, as well as results obtained by appealing to integrable systems or topological strings. We find broad agreement, while also clarifying some lingering inconsistencies. Finally, we explicitly extend the results of Douglas and Shenker to finite N, finding exact agreement with our first calculation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the soft limit of the large scale structure power spectrum: UV dependence

We derive a non-perturbative equation for the large scale structure power spectrum of long-wavelength modes. Thereby, we use an operator product expansion together with relations between the three-point function and power spectrum in the soft limit. The resulting equation encodes the coupling to ultraviolet (UV) modes in two time-dependent coefficients, which may be obtained from response functions to (anisotropic) parameters, such as spatial curvature, in a modified cosmology. We argue that both depend weakly on fluctuations deep in the UV. As a byproduct, this implies that the renormalized leading order coefficient(s) in the effective field theory (EFT) of large scale structures receive most of their contribution from modes close to the non-linear scale. Consequently, the UV dependence found in explicit computations within standard perturbation theory stems mostly from counter-term(s). We confront a simplified version of our non-perturbative equation against existent numerical simulations, and find good agreement within the expected uncertainties. Our approach can in principle be used to precisely infer the relevance of the leading order EFT coefficient(s) using small volume simulations in an 'anisotropic separate universe' framework. Our results suggest that the importance of these coefficient(s) is a ∼ 10% effect, and plausibly smaller.

79 ASTRONOMY AND ASTROPHYSICS↗

Compressive Phase Contrast Tomography

When x-rays penetrate soft matter, their phase changes more rapidly than their amplitude. Interference effects visible with high brightness sources creates higher contrast, edge enhanced images. When the object is piecewise smooth (made of big blocks of a few components), such higher contrast datasets have a sparse solution. We apply basis pursuit solvers to improve SNR, remove ring artifacts, reduce the number of views and radiation dose from phase contrast datasets collected at the Hard X-Ray Micro Tomography Beamline at the Advanced Light Source. We report a GPU code for the most computationally intensive task, the gridding and inverse gridding algorithm (non uniform sampled Fourier transform).

47 OTHER INSTRUMENTATION↗

TU-H-207A-03: CT Hounsfield Unit Accuracy: Effect of Beam Hardening On Phantom and Clinical Whole-Body CT Images

Purpose: To assess the effect of beam hardening on measured CT HU values. Methods: An anthropomorphic knee phantom was scanned with the CT component of a GE Discovery 690 PET/CT scanner (120kVp, 300mAs, 40?0.625mm collimation, pitch=0.984, FOV=500mm, matrix=512?512) with four different scan setups, each of which induces different degrees of beam hardening by introducing additional attenuation media into the field of view. Homogeneous voxels representing “soft tissue” and “bone” were segmented by HU thresholding followed by a 3D morphological erosion operation which removes the non-homogenous voxels located on the interface of thresholded tissue mask. HU values of segmented “soft tissue” and “bone” were compared.Additionally, whole-body CT data with coverage from the skull apex to the end of toes were retrospectively retrieved from seven PET/CT exams to evaluate the effect of beam hardening in vivo. Homogeneous bone voxels were segmented with the same method previously described. Total In-Slice Attenuation (TISA) for each CT slice, defined as the summation of HU values over all voxels within a CT slice, was calculated for all slices of the seven whole-body CT datasets and evaluated against the mean HU values of homogeneous bone voxels within that slice. Results: HU values measured from the phantom showed that while “soft tissue” HU values were unaffected, added attenuation within the FOV caused noticeable decreases in the measured HU values of “bone” voxels. A linear relationship was observed between bone HU and TISA for slices of the torso and legs, but not of the skull. Conclusion: Beam hardening effect is not an issue of concern for voxels with HU in the soft tissue range, but should not be neglected for bone voxels. A linear relationship exists between bone HU and the associated TISA in non-skull CT slices, which can be exploited to develop a correction strategy.

60 APPLIED LIFE SCIENCES↗

Compiler-Directed Soft Error Detection and Recovery to Avoid DUE and SDC via Tail-DMR

This article presents Clover, a compiler-directed soft error detection and recovery scheme for lightweight soft error resilience. The compiler carefully generates soft-error-tolerant code based on idempotent processing without explicit checkpoints. During program execution, Clover relies on a small number of acoustic wave detectors deployed in the processor to identify soft errors by sensing the wave made by a particle strike. To cope with DUEs (detected unrecoverable errors) caused by the sensing latency of error detection, Clover leverages a novel selective instruction duplication technique called tail-DMR (dual modular redundancy) that provides a region-level error containment. Once a soft error is detected by either the sensors or the tail-DMR, Clover takes care of the error as in the case of exception handling. To recover from the error, Clover simply redirects program control to the beginning of the code region where the error is detected. The experimental results demonstrate that the average runtime overhead is only 26%, which is a 75% reduction compared to that of the state-of-the-art soft error resilience technique. In addition, this article evaluates an alternative technique called tail-wait, comparing it to Clover. According to the evaluation with the different processor configurations and the various error detection latencies, Clover turns out to be a superior technique, achieving 1.06 to 3.49× speedup over the tail-wait.

Computer Science↗

Exploiting Social Media Sensor Networks through Novel Data Fusion Techniques

Unprecedented amounts of data are continuously being generated by sensors (“hard” data) and by humans (“soft” data), and this data needs to be exploited to its full potential. The first step in exploiting this data is determine how the hard and soft data are related to each other. In this project we fuse hard and soft data, using the attributes of each (e.g., time and space), to gain more information about interesting events. Next, we attempt to use social networking textual data to predict the present (i.e., predict that an interesting event is occurring and details about the event) using data mining, machine learning, natural language processing, and text analysis techniques.

97 MATHEMATICS AND COMPUTING↗

Scalar-graviton amplitudes and celestial holography

We compute scattering amplitudes involving one massive scalar and two, three, or four gravitons. We show that when the conformal dimension of the massive scalar is set to zero, the resulting celestial correlators depend only on the coordinates of the gravitons. Such correlators of gravitons are well-defined and do not suffer from divergences associated with the Mellin transform of usual graviton amplitudes. Moreover, they are non-distributional and take the form of standard CFT correlators. We show that they are consistent with the usual OPEs but the statement of the soft theorem is modified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Phonon localization in binary alloys with diagonal and off-diagonal disorder: A cluster Green's function approach

We report the development and application of a method for carrying out computational investigations of the effects of mass and force-constant (FC) disorder on phonon spectra. The method is based on the recently developed typical medium dynamical cluster approach, which is a Green's function approach. Excellent quantitative agreement with previous exact diagonalization results establishes the veracity of the method. Application of the method to a model system of binary mass and an FC-disordered system leads to several findings. A narrow resonance, significantly below the van Hove singularity, that has been termed the boson peak, is seen to emerge for low soft particle concentrations. We show, using the typical phonon spectrum, that the states constituting the boson peak cross over from being completely localized to being extended as a function of increasing soft particle concentration. In general, an interplay of mass and FC disorder is found to be cooperative in nature, enhancing phonon localization over all frequencies. However, for a certain range of frequencies, and depending on material parameters, FC disorder can delocalize the states that were localized by mass disorder, and vice versa. Modeling vacancies as weakly bonded sites with vanishing mass, we find that vacancies, even at very low concentrations, are extremely effective in localizing phonons. Thus, inducing vacancies is proposed as a promising route for efficient thermoelectrics. Lastly, we use model parameters corresponding to the alloy system, Ni 1- x Pt x , and we show that mass disorder alone is insufficient to explain the pseudogap in the phonon spectrum; the concomitant presence of FC disorder is necessary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Prediction of Alpha-Particle-Immune Gate-All-Around Field-Effect Transistors (GAA-FET) Based SRAM Design

Alpha particles are known to be a major source of particles creating soft errors in semiconductor devices, such as content flipping in Static Random-Access Memory (SRAM). Recent advancements in transistor nodes have led to the introduction of Gate-All-Around Field Effect Transistors (GAA-FETs), which have better gate control, thus better electrostatics. Moreover, the introduction of bottom dielectric isolation (BDI) eliminates substrate leakage and thus is expected to enhance its radiation hardness. It is thus important to explore if one can design an SRAM that is completely radiation-hard to alpha particles. In this paper, using 3D Technology Computer-Aided-Design (TCAD) simulations, we show that it is possible to design an SRAM using GAA-FET technology so that it is immune to single alpha particle radiation error. In other words, with the design, there will be no single-event upset (SEU) due to alpha particles. We first use ab initio calculations in PHITS to show that there is a maximum linear energy transfer (LET), LET max , for the alpha particle in Si and Si x Ge 1-x . Based on that, by de signing a sub-7nm GAA-FET-based SRAM with BDI, we show that the SRAM does not flip even if the particle strike is in the worst-case scenario for LET > LET max .

42 ENGINEERING↗

Mechanical Response of Porcine Liver Tissue under High Strain Rate Compression

In automobile accidents, abdominal injuries are often life-threatening yet not apparent at the time of initial injury. The liver is the most commonly injured abdominal organ from this type of trauma. In contrast to current safety tests involving crash dummies, a more detailed, efficient approach to predict the risk of human injuries is computational modelling and simulations. Further, the development of accurate computational human models requires knowledge of the mechanical properties of tissues in various stress states, especially in high-impact scenarios. In this study, a polymeric split-Hopkinson pressure bar (PSHPB) was utilized to apply various high strain rates to porcine liver tissue to investigate its material behavior during high strain rate compression. Liver tissues were subjected to high strain rate impacts at 350, 550, 1000, and 1550 s –1 . Tissue directional dependency was also explored by PSHPB testing along three orthogonal directions of liver at a strain rate of 350 s –1 . Histology of samples from each of the three directions was performed to examine the structural properties of porcine liver. Porcine liver tissue showed an inelastic and strain rate-sensitive response at high strain rates. The liver tissue was found lacking directional dependency, which could be explained by the isotropic microstructure observed after staining and imaging. Furthermore, finite element analysis (FEA) of the PSHPB tests revealed the stress profile inside liver tissue and served as a validation of PSHPB methodology. The present findings can assist in the development of more accurate computational models of liver tissue at high-rate impact conditions allowing for understanding of subfailure and failure mechanisms.

60 APPLIED LIFE SCIENCES↗

The LPM effect in sequential bremsstrahlung 2: factorization

The splitting processes of bremsstrahlung and pair production in a medium are coherent over large distances in the very high energy limit, which leads to a suppression known as the Landau-Pomeranchuk-Migdal (LPM) effect. In this paper, we continue analysis of the case when the coherence lengths of two consecutive splitting processes overlap (which is important for understanding corrections to standard treatments of the LPM effect in QCD), avoiding soft-gluon approximations. In particular, this paper analyzes the subtle problem of how to precisely separate overlapping double splitting (e.g. overlapping double bremsstrahlung) from the case of consecutive, independent bremsstrahlung (which is the case that would be implemented in a Monte Carlo simulation based solely on single splitting rates). As an example of the method, we consider the rate of real double gluon bremsstrahlung from an initial gluon with various simplifying assumptions (thick media; q approximation; large N c ; and neglect for the moment of processes involving 4-gluon ver-tices) and explicitly compute the correction Δ dΓ/dx dy due to overlapping formation times.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

TH-A-BRF-03: Evaluation of Synthetic CTs Generated Using MR-SIM Data

Purpose: To describe and evaluate a novel algorithm for generating synthetic CT images from MR-SIM data for dose calculations in MR-only treatment planning. Methods: A voxel-based weighted summation method was implemented to generate synthetic CT (synCT) images. MR data were acquired using Philips 1.0T Panorama high-field open MR-SIM. Retrospective patient data from seven prostate patients and one brain patient (three lesions) enrolled in an IRB-approved study were used. 3D T1-weighted fast field echo and 3D T2-weighted turbo spin echo sequences were utilized for all patients. A 3D balanced turbo field echo sequence using spectral presaturation with inversion recovery was acquired for prostate patients, but 3D ultra-short echo time (UTE)-DIXON was instead acquired for the brain patient to amplify bone signal for semi-automatic bone segmentation. Weight optimization was performed using a training subset of patients. HU value differences between planning CT and synCTs were analyzed using mean absolute error (MAE). Original patient CT-based treatment plans were mapped onto synCTs, dose was recalculated using original leaf motion and MU values, and DRRs were generated. Dosevolume metrics and gamma analysis were used for dosimetric evaluation. Results: Average whole-body MAE of synCTs across all patients was 75+12 HU. In prostate cancer patients, average HU difference between planning and synCTs was 0.9±1.0% for soft tissue structures and 4.3±2.5% for bony structures. DRRs were generated from synCTs and qualitatively showed good geometric agreement with planning CT-generated DRRs. D99, mean dose, and maximum dose to CTV calculated using the synCT remained within 1.2% of planning CT-based dose calculations. All gamma analysis evaluated at 2%/2mm dose difference/distance to agreement) pass rates were greater than 95% with an average of 99.9±0.1% for prostate patients and 98.4±2.2% for three brain lesions. Conclusion: SynCTs were generated with clinically acceptable accuracy comparable to planning CTs, enabling dose computations for MR-only simulation. Research supported in part by a grant from Philips HealthCare (Best, Netherlands)

60 APPLIED LIFE SCIENCES↗

Frequency- and Time-Domain Methods in Soil-Structure Interaction Analysis

Soil-structure interaction (SSI) analysis in the nuclear industry is currently performed using linear codes that function in the frequency domain. There is a consensus that these frequency-domain codes give reasonably accurate results for low-intensity ground motions that result in almost linear response. For higher intensity ground motions, which may result in nonlinear response in the soil, structure or at the vicinity of the foundation, the adequacy of frequency-domain codes is unproven. Nonlinear analysis, which is only possible in the time domain, is theoretically more appropriate in such cases. These methods are available but are rarely used due to the large computational requirements and a lack of experience with analysts and regulators. This paper presents an assessment of the linear frequency-domain code, SASSI, which is widely used in the nuclear industry, and the time-domain commercial finite-element code, LS-DYNA, for SSI analysis. The assessment involves benchmarking the SSI analysis procedure in LS-DYNA against SASSI for linearly elastic models. After affirming that SASSI and LS-DYNA result in almost identical responses for these models, they are used to perform nonlinear SSI analyses of two structures founded on soft soil. An examination of the results shows that, in spite of using identical material properties, the predictions of frequency- and time-domain codes are significantly different in the presence of nonlinear behavior such as gapping and sliding of the foundation.

42 ENGINEERING↗