Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data Intensive 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 883 records · Page 49

PINN surrogate of Li-ion battery models for parameter inference, Part II: Regularization and application of the pseudo-2D model

Bayesian parameter inference is useful to improve Li-ion battery diagnostics and can help formulate battery aging models. However, it is computationally intensive and cannot be easily repeated for multiple cycles, multiple operating conditions, or multiple replicate cells. To reduce the computational cost of Bayesian calibration, numerical solvers for physics-based models can be replaced with faster surrogates. A physics-informed neural network (PINN) is developed as a surrogate for the pseudo-2D (P2D) battery model calibration. For the P2D surrogate, additional training regularization was needed as compared to the PINN single-particle model (SPM) developed in Part I. Both the PINN SPM and P2D surrogate models are exercised for parameter inference and compared to data obtained from a direct numerical solution of the governing equations. A parameter inference study highlights the ability to use these PINNs to calibrate scaling parameters for the cathode Li diffusion and the anode exchange current density. By realizing computational speed-ups of ~2250x for the P2D model, as compared to using standard integrating methods, the PINN surrogates enable rapid state-of-health diagnostics. Finally, in the low-data availability scenario, the testing error was estimated to ~2 mV for the SPM surrogate and ~10 mV for the P2D surrogate which could be mitigated with additional data.

25 ENERGY STORAGE↗

Nuclear Resonance Vibrational Spectroscopy Definition of O 2 Intermediates in an Extradiol Dioxygenase: Correlation to Crystallography and Reactivity

The extradiol dioxygenases are a large subclass of mononuclear nonheme Fe enzymes that catalyze the oxidative cleavage of catechols distal to their OH groups. These enzymes are important in bioremediation, and there has been significant interest in understanding how they activate O 2 . The extradiol dioxygenase homoprotocatechuate 2,3-dioxygenase (HPCD) provides an opportunity to study this process, as two O 2 intermediates have been trapped and crystallographically defined using the slow substrate 4-nitrocatechol (4NC): a side-on Fe–O 2 –4NC species and a Fe–O 2 –4NC peroxy bridged species. Also with 4NC, two solution intermediates have been trapped in the H200N variant, where H200 provides a second-sphere hydrogen bond in the wild-type enzyme. While the electronic structure of these solution intermediates has been defined previously as Fe III -superoxo-catecholate and Fe III -peroxy-semiquinone, their geometric structures are unknown. Nuclear resonance vibrational spectroscopy (NRVS) is an important tool for structural definition of nonheme Fe–O 2 intermediates, as all normal modes with Fe displacement have intensity in the NRVS spectrum. In this study, NRVS is used to define the geometric structure of the H200N–4NC solution intermediates in HPCD as an end-on Fe III -superoxo-catecholate and an end-on Fe III -hydroperoxo-semiquinone. Parallel calculations are performed to define the electronic structures and protonation states of the crystallographically defined wild-type HPCD–4NC intermediates, where the side-on intermediate is found to be a Fe III -hydroperoxo-semiquinone. The assignment of this crystallographic intermediate is validated by correlation to the NRVS data through computational removal of H200. While the side-on hydroperoxo semiquinone intermediate is computationally found to be nonreactive in peroxide bridge formation, it is isoenergetic with a superoxo catecholate species that is competent in performing this reaction. In conclusion, this study provides insight into the relative reactivities of Fe III -superoxo and Fe III -hydroperoxo intermediates in nonheme Fe enzymes and into the role H200 plays in facilitating extradiol catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic measurements and simulations of airborne picolitre-droplet coalescence in holographic optical tweezers

We report studies of the coalescence of pairs of picolitre aerosol droplets manipulated with holographic optical tweezers, probing the shape relaxation dynamics following coalescence by simultaneously monitoring the intensity of elastic backscattered light (EBL) from the trapping laser beam (time resolution on the order of 100 ns) while recording high frame rate camera images (time resolution <10 μs). The goals of this work are to: resolve the dynamics of droplet coalescence in holographic optical traps; assign the origin of key features in the time-dependent EBL intensity; and validate the use of the EBL alone to precisely determine droplet surface tension and viscosity. For low viscosity droplets, two sequential processes are evident: binary coalescence first results from the overlap of the optical traps on the time scale of microseconds followed by the recapture of the composite droplet in an optical trap on the time scale of milliseconds. As droplet viscosity increases, the relaxation in droplet shape eventually occurs on the same time scale as recapture, resulting in a convoluted evolution of the EBL intensity that inhibits quantitative determination of the relaxation time scale. Droplet coalescence was simulated using a computational framework to validate both experimental approaches. The results indicate that time-dependent monitoring of droplet shape from the EBL intensity allows for robust determination of properties such as surface tension and viscosity. Finally, the potential of high frame rate imaging to examine the coalescence of dissimilar viscosity droplets is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffraction pattern simulation of cellulose fibrils using distributed and quantized pair distances

Intensity simulation of X-ray scattering from large twisted cellulose molecular fibrils is important in understanding the impact of chemical or physical treatments on structural properties such as twisting or coiling. This paper describes a highly efficient method for the simulation of X-ray diffraction patterns from complex fibrils using atom-type-specific pair-distance quantization. Pair distances are sorted into arrays which are labelled by atom type. Histograms of pair distances in each array are computed and binned and the resulting population distributions are used to represent the whole pair-distance data set. These quantized pair-distance arrays are used with a modified and vectorized Debye formula to simulate diffraction patterns. This approach utilizes fewer pair distances in each iteration, and atomic scattering factors are moved outside the iteration since the arrays are labelled by atom type. As a result, this algorithm significantly reduces the computation time while maintaining the accuracy of diffraction pattern simulation, making possible the simulation of diffraction patterns from large twisted fibrils in a relatively short period of time, as is required for model testing and refinement.

09 BIOMASS FUELS↗

SU-F-J-217: Accurate Dose Volume Parameters Calculation for Revealing Rectum Dose-Toxicity Effect Using Deformable Registration in Cervical Cancer Brachytherapy: A Pilot Study

Purpose: To study the feasibility of employing deformable registration methods for accurate rectum dose volume parameters calculation and their potentials in revealing rectum dose-toxicity between complication and non-complication cervical cancer patients with brachytherapy treatment. Method and Materials: Data from 60 patients treated with BT including planning images, treatment plans, and follow-up clinical exam were retrospectively collected. Among them, 12 patients complained about hematochezia were further examined with colonoscopy and scored as Grade 1–3 complication (CP). Meanwhile, another 12 non-complication (NCP) patients were selected as a reference group. To seek for potential gains in rectum toxicity prediction when fractional anatomical deformations are account for, the rectum dose volume parameters D0.1/1/2cc of the selected patients were retrospectively computed by three different approaches: the simple “worstcase scenario” (WS) addition method, an intensity-based deformable image registration (DIR) algorithm-Demons, and a more accurate, recent developed local topology preserved non-rigid point matching algorithm (TOP). Statistical significance of the differences between rectum doses of the CP group and the NCP group were tested by a two-tailed t-test and results were considered to be statistically significant if p < 0.05. Results: For the D0.1cc, no statistical differences are found between the CP and NCP group in all three methods. For the D1cc, dose difference is not detected by the WS method, however, statistical differences between the two groups are observed by both Demons and TOP, and more evident in TOP. For the D2cc, the CP and NCP cases are statistically significance of the difference for all three methods but more pronounced with TOP. Conclusion: In this study, we calculated the rectum D0.1/1/2cc by simple WS addition and two DIR methods and seek for gains in rectum toxicity prediction. The results favor the claim that accurate dose deformation and summation tend to be more sensitive in unveiling the dose-toxicity relationship. This work is supported in part by grant from VARIAN MEDICAL SYSTEMS INC, the National Natural Science Foundation of China (no 81428019 and no 81301940), the Guangdong Natural Science Foundation (2015A030313302)and the 2015 Pearl River S&T Nova Program of Guangzhou (201506010096).

60 APPLIED LIFE SCIENCES↗

Expectation-maximization algorithm for amplitude estimation of saturated optical transient signals

This paper addresses parameter estimation for an optical transient signal when the received data has been right-censored. We develop an expectation-maximization (EM) algorithm to estimate the amplitude of a Poisson intensity with a known shape in the presence of additive background counts, where the measurements are subject to saturation effects. We then compare the results of our algorithm with those of an EM algorithm that is unaware of the censoring.

47 OTHER INSTRUMENTATION↗

National Smart Manufacturing Strategic Plan: To Facilitate More Rapid Development, Deployment and Adoption of Smart Manufacturing Technologies

Smart manufacturing technologies provide real-time data and insight to improve the productivity, efficiency, and competitiveness of U.S. manufacturing, creating the potential for new jobs in the manufacturing sector. These technologies can support U.S. manufacturers’ ability to increase throughput and energy efficiency, and decrease waste, defects, and costs. A wide range of manufacturing and industrial subsectors, particularly energy-intensive and energy-dependent industries, have the potential to benefit from smart manufacturing technologies. There are technical improvements still needed to reduce the costs and barriers (trained workforce and upskilling, software-hardware integration, cost, and technical barriers to deployment of advanced sensors, computing, and communication technologies for existing manufacturing assets) to the adoption of smart manufacturing technologies, especially to small and medium-sized manufacturers, and subsequently, increase the overall adoption rate of smart manufacturing technologies. This report outlines the Department of Energy’s (DOE) strategic plan to accelerate the development and implementation of smart manufacturing technologies in the United States and the actions DOE has taken to use these technologies in smart manufacturing for the United States. DOE’s plan to facilitate more rapid development, deployment and widespread adoption of smart manufacturing technologies derive from the 2018 Strategy for American Leadership in Advanced Manufacturing. The strategy outlines a vision for America’s leadership in advanced manufacturing including the development of intelligent manufacturing systems to optimize manufacturing facilities and support the manufacturing of transformative materials. DOE’s Clean Energy Smart Manufacturing Innovation Institute (CESMII), a Manufacturing USA Institute, focuses on accelerating the development and adoption of advanced sensors, controls, platforms, and models needed for smart manufacturing. The objective is to enhance U.S. manufacturing productivity and global competitiveness through the research and development of these technologies. Smart manufacturing has the potential to improve the overall performance, energy productivity, and efficiency of manufacturing, fostering the economic competitiveness of the manufacturing sector.

99 GENERAL AND MISCELLANEOUS↗

A spherical harmonic transform approach to the indexing of electron back-scattered diffraction patterns

A new approach is proposed for the indexing of electron back-scattered diffraction (EBSD) patterns. The algorithm employs a spherical master EBSD pattern and computes its cross-correlation with a back-projected experimental pattern using the spherical harmonic transform (SHT). This approach is significantly faster than the recent dictionary indexing algorithm, but shares the latter’s robustness against noise. The underlying theory is presented, followed by example applications, one on a series of Ni EBSD data sets recorded with decreasing signal-to-noise ratio, the other on a large shot-peened Al data set. The dependence of indexing speed and memory usage on the SHT bandwidth is explored. The speed gains of the new algorithm are achieved by executing real-valued Fast Fourier Transforms, explicitly incorporating crystallographic symmetry in the cross-correlation computation, and using efficient loop ordering to improve the caching behavior. The algorithm produces a cross-correlation array in the zyz Euler space; an orientation refinement procedure is proposed based on analytical derivatives of the Wigner d functions. As a result, the new approach can be applied to any diffraction modality for which the scattered intensity can be represented on a spherical surface.

36 MATERIALS SCIENCE↗

Water at the graphene–substrate interface: interaction with short laser pulses

We have investigated the role of adsorption water in the local transformation of multilayer graphene deposited on an oxidised silicon substrate, which was exposed to nanosecond low-intensity focused laser radiation with a wavelength of 532 nm in the air. Experimental data obtained for a laser energy density E = 0.04 J cm{sup -2} suggest that the formation of micropits (craters) is a consequence of the multipulse removal of the layer of a water adsorbate, which is intercalated between graphene and the substrate, from the zone of laser irradiation of the graphene sheet. The energy threshold of graphene damage in the regions devoid of water was found to be higher in comparison with the initial one (0.058 against 0.048 J cm{sup -2}). According to computer simulations of the heating dynamics of the sample and the heat distribution in the substrate – adsorbate – graphene multilayer system, at energy densities corresponding to the experimental ones the water adsorbate layer heats to a temperature sufficiently high to form an increased-pressure vapour cavity under the graphene film. (interaction of laser radiation with matter)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-F-J-94: Development of a Plug-in Based Image Analysis Tool for Integration Into Treatment Planning

Purpose: To extend the functionality of a commercial treatment planning system (TPS) to support (i) direct use of quantitative image-based metrics within treatment plan optimization and (ii) evaluation of dose-functional volume relationships to assist in functional image adaptive radiotherapy. Methods: A script was written that interfaces with a commercial TPS via an Application Programming Interface (API). The script executes a program that performs dose-functional volume analyses. Written in C#, the script reads the dose grid and correlates it with image data on a voxel-by-voxel basis through API extensions that can access registration transforms. A user interface was designed through WinForms to input parameters and display results. To test the performance of this program, image- and dose-based metrics computed from perfusion SPECT images aligned to the treatment planning CT were generated, validated, and compared. Results: The integration of image analysis information was successfully implemented as a plug-in to a commercial TPS. Perfusion SPECT images were used to validate the calculation and display of image-based metrics as well as dose-intensity metrics and histograms for defined structures on the treatment planning CT. Various biological dose correction models, custom image-based metrics, dose-intensity computations, and dose-intensity histograms were applied to analyze the image-dose profile. Conclusion: It is possible to add image analysis features to commercial TPSs through custom scripting applications. A tool was developed to enable the evaluation of image-intensity-based metrics in the context of functional targeting and avoidance. In addition to providing dose-intensity metrics and histograms that can be easily extracted from a plan database and correlated with outcomes, the system can also be extended to a plug-in optimization system, which can directly use the computed metrics for optimization of post-treatment tumor or normal tissue response models. Supported by NIH - P01 - CA059827.

60 APPLIED LIFE SCIENCES↗

High-resolution broadband spectroscopy using externally dispersed interferometry at the Hale telescope: Part 1, data analysis and results

High-resolution broadband spectroscopy at near-infrared wavelengths (950 to 2450 nm) has been performed using externally dispersed interferometry (EDI) at the Hale telescope at Mt. Palomar. Observations of stars were performed with the “TEDI” interferometer mounted within the central hole of the 200-in. primary mirror in series with the comounted TripleSpec near-infrared echelle spectrograph. These are the first multidelay EDI demonstrations on starlight, as earlier measurements used a single delay or laboratory sources. We demonstrate very high (10×) resolution boost, from original 2700 to 27,000 with current set of delays (up to 3 cm), well beyond the classical limits enforced by the slit width and detector pixel Nyquist limit. Significantly, the EDI used with multiple delays rather than a single delay as used previously yields an order of magnitude or more improvement in the stability against native spectrograph point spread function (PSF) drifts along the dispersion direction. We observe a dramatic (20×) reduction in sensitivity to PSF shift using our standard processing. A recently realized method of further reducing the PSF shift sensitivity to zero is described theoretically and demonstrated in a simple simulation which produces a 350× times reduction. We demonstrate superb rejection of fixed pattern noise due to bad detector pixels—EDI only responds to changes in pixel intensity synchronous to applied dithering. This part 1 describes data analysis, results, and instrument noise. Lastly, a section on theoretical photon limited sensitivity is in a companion paper, part 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structure and behavior of water-laden CH 4 /air counterflow diffusion flames

A counterflow configuration was used to measure thermal and species structure in water-vapor diluted nonpremixed methane–air flames. The motivation is to understand the chemical and thermal effects that water has when it is introduced as a diluent into the fuel side. Here, this work is relevant to combustion processes where water is incorporated naturally in the fuel; e.g., methane hydrates, and when water is added intentionally for emission reduction such as in flares and H 2 O/fuel emulsions combustion. Experimental data are compared to 1-D computations. The agreement is generally very good, but the one dimensional counterflow diffusion model overpredicts flame temperature and major radical, OH, concentration very near extinction in highly diluted H 2 O–methane/air diffusion flames. Changes in flame position, flame width, and peak temperature with the addition of water were measured. Flame temperatures were measured with thin filament pyrometry. OH-PLIF is used to characterize the flame reaction zone with water dilution; the OH distribution, flame position and thickness from the OH-PLIF images were measured. The results show that the OH intensity and reaction zone thickness decreases with the increase in water. Predictions and experiments demonstrate that water mainly acts thermally to lower the flame temperature until extinction. The OH maximum intensity shifts towards the air side of the counterflow burner with water addition. OH is also measured with CO 2 dilution of the fuel stream, and the results are compared with H 2 O addition, including comparisons with the OH molar peak predictions obtained using the GRI 3.0 mechanism and the CHEMKIN Pro one dimensional counterflow model. The study indicates that water’s chemical effects are to change the production and depletion of OH, H and O radicals, especially near extinction. Chemical kinetics simulation of the flame demonstrates good agreement in OH and flame temperatures over a wide range of dilution away from extinction, particularly for CO 2 . Lastly, an over prediction of the water carrying capacity near extinction is found for highly water-diluted flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SmartArc-Based Volumetric Modulated Arc Therapy for Oropharyngeal Cancer: A Dosimetric Comparison With Both Intensity-Modulated Radiation Therapy and Helical Tomotherapy

Purpose: To investigate the roles of volumetric modulated arc therapy with SmartArc (VMAT-S), intensity-modulated radiation therapy (IMRT), and helical tomotherapy (HT) for oropharyngeal cancer using a simultaneous integrated boost (SIB) approach. Methods and Materials: Eight patients treated with IMRT were selected at random. Plans were computed for both IMRT and VMAT-S (using Pinnacle TPS for an Elekta Infinity linac) along with HT. A three-dose level prescription was used to deliver 70 Gy, 63 Gy, and 58.1 Gy to regions of macroscopic, microscopic high-risk, and microscopic low-risk disease, respectively. All doses were given in 35 fractions. Comparisons were performed on dose-volume histogram data, monitor units per fraction (MU/fx), and delivery time. Results: VMAT-S target coverage was close to that achieved by IMRT, but inferior to HT. The conformity and homogeneity within the PTV were improved for HT over all strategies. Sparing of the organs at risk (OAR) was achieved with all modalities. VMAT-S (along with HT) shortened delivery time (mean, -38%) and reduced MU/fx (mean, -28%) compared with IMRT. Conclusion: VMAT-S represents an attractive solution because of the shorter delivery time and the lower number of MU/fx compared with IMRT. However, in this complex clinical setting, current VMAT-S does not appear to provide any distinct advantage compared with helical tomotherapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An accelerated framework for predicting creep rupture lifetimes in engineering alloys

Confidently predicting high-temperature deformation, including creep and creep rupture, is paramount for the design and commercialization of candidate materials for advanced nuclear energy systems. To accelerate creep quantification, we introduce a framework that enables rapid, cost-effective, and reliable prediction of creep rupture lifetimes, minimizing reliance on time-intensive bulk creep testing. Unlike conventional creep analysis, which requires extensive time and resources, our method leverages a maximum of four short-term bulk creep tests as training data for prediction. This framework combines high-throughput nanoindentation up to 700 °C with these targeted bulk tests to inform our creep rupture model in order to predict rupture lifetimes. The strong agreement between our predictions and conventional experimental data demonstrates the effectiveness of our approach for accelerated creep analysis and lifetime prediction of structural components in high-temperature applications. Our multi-pronged approach motivates further integration of computational tools and advanced instrumentation to establish a universal framework for understanding high-temperature material responses.

36 MATERIALS SCIENCE↗

Consistency in the Sum Frequency Generation Intensity and Phase Vibrational Spectra of the Air/Neat Water Interface

Tremendous progresses have been made in quantitative understanding and interpretation of the hydrogen bonding and ordering structure at the air/water interface since the first sum-frequency generation vibrational spectroscopy (SFG-VS) measurement on the neat air/water interface by Q. Du et al. in 1993 (PRL, 70, 2312-2316, 1993.). However, there are still disagreements and controversies on the consistency between the different experiment measurements and the theoretical computational results. One critical problem lies in the inconsistency between the SFG-VS intensity measurements and the recently developed SFG-VS phase spectra measurements of the neat air/water interface, which has inspired various theoretical efforts trying to understand them. In this report, the reliability of the SFG-VS intensity spectra of the neat air/water interface is to be quantitatively examined, and the sources of possible inaccuracies in the SFG-VS phase spectral measurement is to be discussed based on the non-resonant SHG phase measurement results. The conclusion is that the SFG-VS intensity spectra data from different laboratories are now quantitatively converging and in agreement with each other, and the possible inaccuracies and inconsistencies in the SFG-VS phase spectra measurements need to be carefully examined against the properly corrected phase standard.

08 HYDROGEN↗

Reconstructing three-dimensional protein crystal intensities from sparse unoriented two-axis X-ray diffraction patterns

Recently, there has been a growing interest in adapting serial microcrystallography (SMX) experiments to existing storage ring (SR) sources. For very small crystals, however, radiation damage occurs before sufficient numbers of photons are diffracted to determine the orientation of the crystal. The challenge is to merge data from a large number of such `sparse' frames in order to measure the full reciprocal space intensity. To simulate sparse frames, a dataset was collected from a large lysozyme crystal illuminated by a dim X-ray source. The crystal was continuously rotated about two orthogonal axes to sample a subset of the rotation space. With the EMC algorithm [expand–maximize–compress; Loh & Elser (2009). Phys. Rev. E , 80 , 026705], it is shown that the diffracted intensity of the crystal can still be reconstructed even without knowledge of the orientation of the crystal in any sparse frame. Moreover, parallel computation implementations were designed to considerably improve the time and memory scaling of the algorithm. The results show that EMC-based SMX experiments should be feasible at SR sources.

Chemistry↗

Clinical evaluations of an amplitude-based binning algorithm for 4DCT reconstruction in radiation therapy

Purpose: Phase-binning algorithms are commonly utilized in 4DCT image reconstruction for characterization of tumor or organ shape and respiration motion, but breathing irregularities occurring during 4DCT acquisition can cause considerable image distortions. Recently, amplitude-binning algorithms have been evaluated as a potential improvement to phase-binning algorithms for 4DCT image reconstruction. The purpose of this study was to evaluate the performance of the first commercially available on-line retrospective amplitude-binning algorithm for comparison to the traditional phase-binning algorithm. Methods: Both phantom and clinical data were used for evaluation. A phantom of known geometry was mounted on a 4D motion platform programmed with seven respiratory waves (two computer generated and five patient trajectories) and scanned with a Philips Brilliance Big bore 16-slice CT simulator. 4DCT images were reconstructed using commercial amplitude- and phase-binning algorithms. Image quality of the amplitude- and phase-binned image sets was compared by evaluation of shape and volume distortions in reconstructed images. Clinical evaluations were performed on 64 4DCT patient image sets in a blinded review process. The amplitude- and phase-binned 4DCT maximum intensity projection (MIP) images were further evaluated for 28 stereotactic body radiation therapy (SBRT) cases of total 64 cases. A preliminary investigation of the effects of respiratory amplitude and pattern irregularities on motion artifact severity was conducted. Results: The phantom experiments illustrated that, as expected, maximum inhalation occurred at the 0% amplitude and maximum exhalation occurred at the 50% amplitude of the amplitude-binned 4DCT image sets. The phantom shape distortions were more severe in the images reconstructed from the phase-binning algorithm. In the clinical study, compared to the phase-binning algorithm, the amplitude-binning algorithm yielded fewer or less severe motion artifacts in 37.5% of the cases (24/64), comparable artifacts in 54.7% of the cases (35/64), and slightly greater artifacts in 7.8% of the cases (5/64). Evaluation of SBRT cases demonstrated that the reconstructed tumor sizes and locations were comparable in 96% (1/28) of the MIP image pairs generated from both amplitude- and phase-binning algorithms. In this case the amplitude-binned image set rendered a smaller tumor size, which was likely due to very shallow respiratory amplitudes occurring over several breathing cycles. Conclusions: Overall, the amplitude-binning algorithm for 4DCT reconstruction reduced the severity of tumor distortion and image artifacts compared to the phase-binning algorithm. However, the full range of motion may not be characterized using amplitude-binning algorithms. Despite superior performance, amplitude binning can still be susceptible to motion artifacts caused by large variations in amplitude of respiratory waves.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Development of a Virtual CFR Engine Model for Knocking Combustion Analysis

Knock is a major bottleneck to achieving higher thermal efficiency in spark ignition (SI) engines. The overall tendency to knock is highly dependent on fuel anti-knock quality as well as engine operating conditions. It is, therefore, critical to gain a better understanding of fuel-engine interactions in order to develop robust knock mitigation strategies. In the present work, a numerical model based on three-dimensional (3-D) computational fluid dynamics (CFD) was developed to capture knock in a Cooperative Fuel Research (CFR) engine. For combustion modeling, a hybrid approach incorporating the G-equation model to track turbulent flame propagation, and a homogeneous reactor multi-zone model to predict end-gas auto-ignition ahead of the flame front and post-flame oxidation in the burned zone, was employed. In addition, a hybrid methodology was implemented wherein a laminar flame speed lookup table generated a priori from a chemical kinetic mechanism could be used to provide flame speed as an input to the G-equation model, instead of using conventional empirical correlations. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations were performed for two different spark timings (STs) corresponding to non-knocking and knocking conditions, with other operating conditions kept the same as those of a standard research octane number (RON) test. Isooctane was considered as the fuel for the numerical study. Two different reduced kinetic mechanisms were employed to describe end-gas auto-ignition chemistry and to generate the flame speed lookup table. Experimental data, including intake/exhaust boundary conditions, was provided by an isooctane ST sweep study conducted in an in-house CFR engine. Moreover, cylinder wall/valve/port surface temperatures and residual gas fraction (RGF) were estimated using a well-calibrated one-dimensional (1-D) model. On the other hand, a novel methodology was also developed to analyze experimental data for the knocking case and identify the most representative cycle. For the non-knocking case, a good agreement was found between experiment and CFD simulation, with respect to cycle-averaged values of 10% burn point (CA10), 50% burn point (CA50), and peak pressure magnitude/location. The virtual CFR engine model was also demonstrated to be capable of predicting average knock characteristics for the knocking case, such as knock point, knock intensity, and energy of resonance, with good accuracy.

42 ENGINEERING↗