Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set methods”

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,279 records · Page 71

What Makes You Hold on to That Old Car? Joint Insights From Machine Learning and Multinomial Logit on Vehicle-Level Transaction Decisions

What makes you hold on to that old car? While the vast majority of household vehicles are still powered by conventional internal combustion engines, the progress of adopting emerging vehicle technologies will critically depend on how soon the existing vehicles are transacted out of the household fleet. Leveraging a nationally representative longitudinal data set, the Panel Study of Income Dynamics, this study examines how household decisions to dispose of or replace a given vehicle are: 1) influenced by the vehicle’s attributes, 2) mediated by households’ concurrent socio-demographic and economic attributes, and 3) triggered by key life cycle events. Coupled with a newly developed machine learning interpretation tool, TreeExplainer, we demonstrate an innovative use of machine learning models to augment traditional logit modeling to both generate behavioral insights and improve model performance. We find the two gradient-boosting-based methods, CatBoost and LightGBM, are the best performing machine learning models for this problem. The multinomial logistic model can achieve similar performance levels after its model specification is informed by TreeExplainer. Both machine learning and multinomial logit models suggest that while older vehicles are more likely to be disposed of or replaced than newer ones, such probability decreases as the vehicles serve the family longer. Pickup trucks and sport utility vehicles are less likely to be disposed of or replaced than cars, and leased vehicles are more likely to be transacted than owned vehicles. We find that married families, families with higher education levels, homeowners, and older families tend to keep their vehicles longer. Life events such as childbirth, residential relocation, and change of household composition and income are found to increase vehicle disposal and/or replacement. We provide additional insights on the timing of vehicle replacement or disposal, in particular, the presence of children and childbirth events are more strongly associated with vehicle replacement among younger parents.

33 ADVANCED PROPULSION SYSTEMS↗

Bayesian data assimilation for stochastic multiscale models of transport in porous media.

We investigate Bayesian techniques that can be used to reconstruct field variables from partial observations. In particular, we target fields that exhibit spatial structures with a large spectrum of lengthscales. Contemporary methods typically describe the field on a grid and estimate structures which can be resolved by it. In contrast, we address the reconstruction of grid-resolved structures as well as estimation of statistical summaries of subgrid structures, which are smaller than the grid resolution. We perform this in two different ways (a) via a physical (phenomenological), parameterized subgrid model that summarizes the impact of the unresolved scales at the coarse level and (b) via multiscale finite elements, where specially designed prolongation and restriction operators establish the interscale link between the same problem defined on a coarse and fine mesh. The estimation problem is posed as a Bayesian inverse problem. Dimensionality reduction is performed by projecting the field to be inferred on a suitable orthogonal basis set, viz. the Karhunen-Loeve expansion of a multiGaussian. We first demonstrate our techniques on the reconstruction of a binary medium consisting of a matrix with embedded inclusions, which are too small to be grid-resolved. The reconstruction is performed using an adaptive Markov chain Monte Carlo method. We find that the posterior distributions of the inferred parameters are approximately Gaussian. We exploit this finding to reconstruct a permeability field with long, but narrow embedded fractures (which are too fine to be grid-resolved) using scalable ensemble Kalman filters; this also allows us to address larger grids. Ensemble Kalman filtering is then used to estimate the values of hydraulic conductivity and specific yield in a model of the High Plains Aquifer in Kansas. Strong conditioning of the spatial structure of the parameters and the non-linear aspects of the water table aquifer create difficulty for the ensemble Kalman filter. We conclude with a demonstration of the use of multiscale stochastic finite elements to reconstruct permeability fields. This method, though computationally intensive, is general and can be used for multiscale inference in cases where a subgrid model cannot be constructed.

54 ENVIRONMENTAL SCIENCES↗

Respiratory Amplitude Guided 4-Dimensional Magnetic Resonance Imaging

Purpose: To evaluate the feasibility of prospectively guiding 4-dimensional (4D) magnetic resonance imaging (MRI) image acquisition using triggers at preselected respiratory amplitudes to achieve T{sub 2} weighting for abdominal motion tracking. Methods and Materials: A respiratory amplitude-based triggering system was developed and integrated into a commercial turbo spin echo MRI sequence. Initial feasibility tests were performed on healthy human study participants. Four respiratory states, the middle and the end of inhalation and exhalation, were used to trigger 4D MRI image acquisition of the liver. To achieve T{sub 2} weighting, the echo time and repetition time were set to 75 milliseconds and 4108 milliseconds, respectively. Single-shot acquisition, together with parallel imaging and partial k-space imaging techniques, was used to improve image acquisition efficiency. 4D MRI image sets composed of axial or sagittal slices were acquired. Results: Respiratory data measured and logged by the MRI scanner showed that the triggers occurred at the appropriate respiratory levels. Liver motion could be easily observed on both 4D MRI image datasets by sensing either the change of liver in size and shape (axial) or diaphragm motion (sagittal). Both 4D MRI image datasets were T{sub 2}-weighted as expected. Conclusions: This study demonstrated the feasibility of achieving T{sub 2}-weighted 4D MRI images using amplitude-based respiratory triggers. With the aid of the respiratory amplitude-based triggering system, the proposed method is compatible with most MRI sequences and therefore has the potential to improve tumor-tissue contrast in abdominal tumor motion imaging.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Removal of Historic Low-Level Radioactive Sediment from the Port Hope Harbour - 13314

At the Port Hope Harbour, located on the north shore of Lake Ontario, the presence of low-level radioactive sediment, resulting from a former radium and uranium refinery that operated alongside the Harbour, currently limits redevelopment and revitalization opportunities. These waste materials contain radium-226, uranium, arsenic and other contaminants. Several other on-land locations within the community of Port Hope are also affected by the low-level radioactive waste management practices of the past. The Port Hope Project is a community initiated undertaking that will result in the consolidation of an estimated 1.2 million cubic metres of the low-level radioactive waste from the various sites in Port Hope into a new engineered above ground long-term waste management facility. The remediation of the estimated 120,000 m{sup 3} of contaminated sediments from the Port Hope Harbour is one of the more challenging components of the Port Hope Project. Following a thorough review of various options, the proposed method of contaminated sediment removal is by dredging. The sediment from the dredge will then be pumped as a sediment-water slurry mixture into geo-synthetic containment tubes for dewatering. Due to the hard substrate below the contaminated sediment, the challenge has been to set performance standards in terms of low residual surface concentrations that are attainable in an operationally efficient manner. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Radio-Frequency Tank Eigenmode Sensor for Propellant Quantity Gauging

Although there are several methods for determining liquid level in a tank, there are no proven methods to quickly gauge the amount of propellant in a tank while it is in low gravity or under low-settling thrust conditions where propellant sloshing is an issue. Having the ability to quickly and accurately gauge propellant tanks in low-gravity is an enabling technology that would allow a spacecraft crew or mission control to always know the amount of propellant onboard, thus increasing the chances for a successful mission. The Radio Frequency Mass Gauge (RFMG) technique measures the electromagnetic eigenmodes, or natural resonant frequencies, of a tank containing a dielectric fluid. The essential hardware components consist of an RF network analyzer that measures the reflected power from an antenna probe mounted internal to the tank. At a resonant frequency, there is a drop in the reflected power, and these inverted peaks in the reflected power spectrum are identified as the tank eigenmode frequencies using a peak-detection software algorithm. This information is passed to a pattern-matching algorithm, which compares the measured eigenmode frequencies with a database of simulated eigenmode frequencies at various fill levels. A best match between the simulated and measured frequency values occurs at some fill level, which is then reported as the gauged fill level. The database of simulated eigenmode frequencies is created by using RF simulation software to calculate the tank eigenmodes at various fill levels. The input to the simulations consists of a fairly high-fidelity tank model with proper dimensions and including internal tank hardware, the dielectric properties of the fluid, and a defined liquid/vapor interface. Because of small discrepancies between the model and actual hardware, the measured empty tank spectra and simulations are used to create a set of correction factors for each mode (typically in the range of 0.999 1.001), which effectively accounts for the small discrepancies. These correction factors are multiplied to the modes at all fill levels. By comparing several measured modes with the simulations, it is possible to accurately gauge the amount of propellant in the tank. An advantage of the RFMG approach of applying computer simulations and a pattern-matching algorithm is that the Although there are several methods for determining liquid level in a tank, there are no proven methods to quickly gauge the amount of propellant in a tank while it is in low gravity or under low-settling thrust conditions where propellant sloshing is an issue. Having the ability to quickly and accurately gauge propellant tanks in low-gravity is an enabling technology that would allow a spacecraft crew or mission control to always know the amount of propellant onboard, thus increasing the chances for a successful mission. The Radio Frequency Mass Gauge (RFMG) technique measures the electromagnetic eigenmodes, or natural resonant frequencies, of a tank containing a dielectric fluid. The essential hardware components consist of an RF network analyzer that measures the reflected power from an antenna probe mounted internal to the tank. At a resonant frequency, there is a drop in the reflected power, and these inverted peaks in the reflected power spectrum are identified as the tank eigenmode frequencies using a peak-detection software algorithm. This information is passed to a pattern-matching algorithm, which compares the measured eigenmode frequencies with a database of simulated eigenmode frequencies at various fill levels. A best match between the simulated and measured frequency values occurs at some fill level, which is then reported as the gauged fill level. The database of simulated eigenmode frequencies is created by using RF simulation software to calculate the tank eigenmodes at various fill levels. The input to the simulations consists of a fairly high-fidelity tank model with proper dimensions and including internal tank harare, the dielectric properties of the fluid, and a defined liquid/vapor interface. Because of small discrepancies between the model and actual hardware, the measured empty tank spectra and simulations are used to create a set of correction factors for each mode (typically in the range of 0.999 1.001), which effectively accounts for the small discrepancies. These correction factors are multiplied to the modes at all fill levels. By comparing several measured modes with the simulations, it is possible to accurately gauge the amount of propellant in the tank. An advantage of the RFMG approach of applying computer simulations and a pattern-matching algorithm is that the

Zimmerli, Gregory A.↗

International Spine Radiosurgery Consortium Consensus Guidelines for Target Volume Definition in Spinal Stereotactic Radiosurgery

Purpose: Spinal stereotactic radiosurgery (SRS) is increasingly used to manage spinal metastases. However, target volume definition varies considerably and no consensus target volume guidelines exist. This study proposes consensus target volume definitions using common scenarios in metastatic spine radiosurgery. Methods and Materials: Seven radiation oncologists and 3 neurological surgeons with spinal radiosurgery expertise independently contoured target and critical normal structures for 10 cases representing common scenarios in metastatic spine radiosurgery. Each set of volumes was imported into the Computational Environment for Radiotherapy Research. Quantitative analysis was performed using an expectation maximization algorithm for Simultaneous Truth and Performance Level Estimation (STAPLE) with kappa statistics calculating agreement between physicians. Optimized confidence level consensus contours were identified using histogram agreement analysis and characterized to create target volume definition guidelines. Results: Mean STAPLE agreement sensitivity and specificity was 0.76 (range, 0.67-0.84) and 0.97 (range, 0.94-0.99), respectively, for gross tumor volume (GTV) and 0.79 (range, 0.66-0.91) and 0.96 (range, 0.92-0.98), respectively, for clinical target volume (CTV). Mean kappa agreement was 0.65 (range, 0.54-0.79) for GTV and 0.64 (range, 0.54-0.82) for CTV (P<.01 for GTV and CTV in all cases). STAPLE histogram agreement analysis identified optimal consensus contours (80% confidence limit). Consensus recommendations include that the CTV should include abnormal marrow signal suspicious for microscopic invasion and an adjacent normal bony expansion to account for subclinical tumor spread in the marrow space. No epidural CTV expansion is recommended without epidural disease, and circumferential CTVs encircling the cord should be used only when the vertebral body, bilateral pedicles/lamina, and spinous process are all involved or there is extensive metastatic disease along the circumference of the epidural space. Conclusions: This report provides consensus guidelines for target volume definition for spinal metastases receiving upfront SRS in common clinical situations.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

High-Grade Glioma Radiation Therapy Target Volumes and Patterns of Failure Obtained From Magnetic Resonance Imaging and {sup 18}F-FDOPA Positron Emission Tomography Delineations From Multiple Observers

Purpose: The objective of this study was to compare recurrent tumor locations after radiation therapy with pretreatment delineations of high-grade gliomas from magnetic resonance imaging (MRI) and 3,4-dihydroxy-6-[{sup 18}F]fluoro-L-phenylalanine ({sup 18}F-FDOPA) positron emission tomography (PET) using contours delineated by multiple observers. Methods and Materials: Nineteen patients with newly diagnosed high-grade gliomas underwent computed tomography (CT), gadolinium contrast-enhanced MRI, and {sup 18}F-FDOPA PET/CT. The image sets (CT, MRI, and PET/CT) were registered, and 5 observers contoured gross tumor volumes (GTVs) using MRI and PET. Consensus contours were obtained by simultaneous truth and performance level estimation (STAPLE). Interobserver variability was quantified by the percentage of volume overlap. Recurrent tumor locations after radiation therapy were contoured by each observer using CT or MRI. Consensus recurrence contours were obtained with STAPLE. Results: The mean interobserver volume overlap for PET GTVs (42% ± 22%) and MRI GTVs (41% ± 22%) was not significantly different (P=.67). The mean consensus volume was significantly larger for PET GTVs (58.6 ± 52.4 cm{sup 3}) than for MRI GTVs (30.8 ± 26.0 cm{sup 3}, P=.003). More than 95% of the consensus recurrence volume was within the 95% isodose surface for 11 of 12 (92%) cases with recurrent tumor imaging. Ten (91%) of these cases extended beyond the PET GTV, and 9 (82%) were contained within a 2-cm margin on the MRI GTV. One recurrence (8%) was located outside the 95% isodose surface. Conclusions: High-grade glioma contours obtained with {sup 18}F-FDOPA PET had similar interobserver agreement to volumes obtained with MRI. Although PET-based consensus target volumes were larger than MRI-based volumes, treatment planning using PET-based volumes may not have yielded better treatment outcomes, given that all but 1 recurrence extended beyond the PET GTV and most were contained by a 2-cm margin on the MRI GTV.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

WE-FG-207B-07: Feasibility of Low Dose Lung Cancer Screening with a Whole-Body Photon Counting CT: First Human Results

Purpose: To evaluate the feasibility of using a whole-body photon counting detector (PCD) CT scanner for low dose lung cancer screening compared to a conventional energy integrating detector (EID) system. Methods: Radiation dose-matched EID and PCD scans of the COPDGene 2 phantom and 2 human volunteers were acquired. Phantom images were acquired at different radiation dose levels (CTDIvol: 3.0, 1.5, and 0.75 mGy) and different tube voltages (120, 100, and 80 kVp), while human images were acquired at vendor recommended low-dose lung cancer screening settings. EID and PCD images were compared for quantitative Hounsfield unit accuracy, noise levels, and contrast-to-noise ratios (CNR) for detection of ground-glass nodules (GGNs) and emphysema. Results: The PCD Hounsfield unit accuracy was better for water at all scan parameters, and for lung, GGN and emphysema equivalent regions of interest (ROIs) at 1.5 and 0.75 mGy. PCD attenuation accuracy was more consistent for all scan parameters (all P<0.01), while Hounsfield units for lung, GGN and emphysema ROIs changed significantly for EID with decreasing dose (all P<0.001). PCD showed lower noise levels at the lowest dose setting at 120, 100 and 80 kVp (15.2±0.3 vs 15.8±0.2, P=0.03; 16.1±0.3 vs 18.0±0.4, P=0.003; and 16.1±0.3 vs 17.9±0.3, P=0.001, respectively), resulting in superior CNR for the detection of GGNs and emphysema at 100 and 80 kVp. Significantly lower PCD noise levels were confirmed in volunteer images. Conclusion: PCD provided better Hounsfield unit accuracy for lung, ground-glass, and emphysema-equivalent foams at 1.5 and 0.75 mGy with less variability than EID. Additionally, PCD showed less noise, and higher CNR at 0.75 mGy for both 100 and 80 kVp. PCD technology may help reduce radiation exposure in lung cancer screening while maintaining diagnostic quality.

60 APPLIED LIFE SCIENCES↗

Measurements and predictions of flyover and static noise of a TF30 afterburning turbofan engine

The noise of the TF30 afterburning turbofan engine in an F-111 airplane was determined from static (ground) and flyover tests. A survey was made to measure the exhaust temperature and velocity profiles for a range of power settings. Comparisons were made between predicted and measured jet mixing, internal, and shock noise. It was found that the noise produced at static conditions was dominated by jet mixing noise, and was adequately predicted by current methods. The noise produced during flyovers exhibited large contributions from internally generated noise in the forward arc. For flyovers with the engine at nonafterburning power, the internal noise, shock noise, and jet mixing noise were accurately predicted. During flyovers with afterburning power settings, however, additional internal noise believed to be due to the afterburning process was evident; its level was as much as 8 decibels above the nonafterburning internal noise. Power settings that produced exhausts with inverted velocity profiles appeared to be slightly less noisy than power settings of equal thrust that produced uniform exhaust velocity profiles both in flight and in static testing.

Burcham, F. W., Jr.↗

Measurements and predictions of flyover and static noise of an afterburning turbofan engine in an F-111 airplane

The noise of the TF30 afterburning turbofan engine in an F-111 airplane was determined from static (ground) and flyover tests. Exhaust temperatures and velocity profiles were measured for a range of power settings. Comparisons were made between predicted and measured jet mixing, internal, and shock noise. It was found that the noise produced at static conditions was dominated by jet mixing noise, and was adequately predicted by current methods. The noise produced during flyovers exhibited large contributions from internally generated noise in the forward arc. For flyovers with the engine at nonafterburning power, the internal noise, shock noise, and jet mixing noise were accurately predicted. During flyovers with afterburning power settings, however, additional internal noise believed to be due to the afterburning process was evident; its level was as much as 8 decibels above the nonafterburning internal noise.

Burcham, F. W., Jr.↗

Objective measurement of human tolerance to +G sub z acceleration stress

The efficacy of a new objective technique using a transcutaneous Doppler flowmeter to monitor superficial temporal artery blood flow velocity during acceleration was investigated. The results were correlated with current objective and subjective G tolerance end points. In over 1300 centrifuge runs, retrograde eye level blood flow leading to total flow cessation was consistently recorded and preceded visual field deterioration leading to blackout by 3 to 23 seconds. The new method was successfully applied as an objective indication of tolerance in a variety of test situations including evaluation of g-suits, straining maneuvers, and 13 deg, 45 deg and 65 deg set back angles.

Rositano, S. A.↗

A baroclinic quasigeostrophic open ocean model

A baroclinic quasigeostrophic open ocean model is presented, calibrated by a series of test problems, and demonstrated to be feasible and efficient for application to realistic mid-oceanic mesoscale eddy flow regimes. Two methods of treating the depth dependence of the flow, a finite difference method and a collocation method, are tested and intercompared. Sample Rossby wave calculations with and without advection are performed with constant stratification and two levels of nonlinearity, one weaker than and one typical of real ocean flows. Using exact analytical solutions for comparison, the accuracy and efficiency of the model is tabulated as a function of the computational parameters and stability limits set; typically, errors were controlled between 1 percent and 10 percent RMS after two wave periods. Further Rossby wave tests with realistic stratification and wave parameters chosen to mimic real ocean conditions were performed to determine computational parameters for use with real and simulated data. Finally, a prototype calculation with quasiturbulent simulated data was performed successfully, which demonstrates the practicality of the model for scientific use.

Miller, R. N.↗

Changes on the ice plain of Ice Stream B and Ross Ice Shelf

During the 1970's and 1980's, nearly 200 stations from which accurate, three dimensional position fixes have been obtained from TRANSIT satellites were occupied throughout the Ross Ice Shelf. We have transformed the elevations obtained by satellite altimetry to the same geodetic datum, and then applied a second transformation to reduce the geodetic heights to elevations above mean sea level using the GEM-10C geoidal height. On the IGY Ross Ice Shelf traverse between Oct. 1957 and Feb. 1958, an accurate method of barometric altimetry was used on a loop around the ice shelf that was directly tied to the sea at both ends of the travel route, thus providing absolute elevations. Comparisons of the two sets of data at 32 station pairs on floating ice show a mean difference of 0 +/- 1 m. The elevation data were also compared with theoretical values of elevations for a hydrostatically floating ice shelf. The mean difference between theoretical and measured values of elevations is -2 +/- 1 m.

Shabtaie, Sion↗

Radio frequency tank eigenmode sensor for propellant quantity gauging

A method for measuring the quantity of fluid in a tank may include the steps of selecting a match between a measured set of electromagnetic eigenfrequencies and a simulated plurality of sets of electromagnetic eigenfrequencies using a matching algorithm, wherein the match is one simulated set of electromagnetic eigenfrequencies from the simulated plurality of sets of electromagnetic eigenfrequencies, and determining the fill level of the tank based upon the match.

Zimmerli, Gregory A.↗

Recommendations on Evidence and Process for Certification of Learning-enabled Components in Aerospace Systems

This report primarily identifies a collection of relevant and necessary evidence for assurance of machine learnt components (MLCs)—also known as learning-enabled components—integrated into aircraft systems, and gives preliminary suggestions on the elements of a certification process that invoke the identified evidence. The main focus is on feedforward neural networks that are static and trained offline through supervised learning. A brief background on the generic elements of the lifecycle of an MLC is given to contextualize the assurance considerations and, consequently, the evidence that is relevant and necessary to support certification. At the level of an MLC, those considerations relate to: (i) the consistency and correctness of MLC contributions to system functions in the context of a validated functional intent; and (ii) the absence of MLC contributions to aircraft-level failure conditions. At an ML model level, confidence in model and data properties contribute to assurance of the containing MLC, in particular: (a) generalizability and robustness of models, in the presence of inputs not previously seen during training, disturbances to inputs, and unexpected inputs; and (b) valid data, i.e., data that are at least representative, relevant, complete, and accurate. Evidence for the above span the elements of the ML lifecycle, and includes, at a minimum, lifecycle artifacts that pertain to: (1) properties of requirements capturing functional intent, safety constraints, and aspects of the intended use and operating environment; (2) model performance, model complexity and design, and algorithm choice; (3) achievement of required performance at the levels of a trained model during model development, a trained model after model development is complete, and a trained model that is transformed into an executable equivalent; (4) model implementation aspects necessary for transforming a trained model into the executable equivalent; (5) integration of the executable trained model into the containing MLC, and eventually the larger system; and, (6) lastly, the verification and validation (V&V) of each of the above. Such V&V lifecycle artifacts themselves include: aspects of coverage, e.g., of various levels of requirements by the input space of the model and the data; traceability (where applicable); application of formal methods for property specification, analysis, and checking. Examples of evidence generation methods and tools further ground the discussion on what constitutes evidence, and the contribution to assurance during certification. The identified assurance considerations and supporting evidence is not a comprehensive set. Additionally, neither what should be considered as sufficient evidence relative to the assigned criticality of an MLC, nor how criticality ought to be determined and adjusted, have been considered in this report. However, suggestions are made for potential activities of the ML lifecycle that are aimed at providing confidence that an MLC can be relied upon when integrated into its containing (aircraft) system. Those activities are proposed as candidate elements of a certification process for MLCs. The main purpose of this report to inform regulatory guidance and consensus standards that may be used to meet the safety intent of the applicable regulations.

Aviation safety↗

Operator Workload and Task Allocation in m:N Operational Architectures of Uncrewed Aerial Systems

Uncrewed aerial systems (UAS) show promise in urban air transport, package delivery, and emergency services. UAS efficiency can be significantly improved by having fewer operators (m) manage a greater number of vehicles (N), or the m:N architecture of operation. The current study investigates how workload affects operators’ task-allocation decision-making and potential effects of two crucial human factors: trust and self-confidence. In the context of a simulated UAS package-delivery task, 10 participants with expertise in UAS operation were recruited. Each participant reported their preferred task-allocation strategy for a set of five subtasks while watching two sets of videos with different workload levels. Perceived workload, trust, and self-confidence were also measured after each video session. Overall, participants indicated a preference for automation for most of the subtasks under the delivery mission. Trust, rather than workload and self-confidence, played a significant role in experts’ decisions of task-allocation and assignment methods. Higher trust led to higher preference for automation.

workload↗

Optimizing a Multiple Right-Hand Side Dslash Kernel for Intel Knights Corner

Lattice Quantum Chromo-dynamics (LQCD) stands unique as the only computationally tractable, non-perturbative, and model-independent quantum field theory of the strong nuclear force. LQCD achieves tractability by discretizing our familiar 4-dimensional space-time onto a hypercubic lattice. The core component of the linear solver used in LQCD is the Wilson-Dirac operator whose finite difference component is called the Dslash operator. Dslash is a nearest neighbor stencil operator summing matrix-vector multiplications over lattice points. The arithmetic intensity of a regular Dslash operator acting on a single vector in single precision is between 0.88-2.29 (depending on the tricks employed) making it memory bandwidth bound on most current architectures, including Xeon Phi. Maximizing the performance of Dslash on Xeon Phi requires deliberate exploitation of its architectural features. One key feature of Xeon Phi is its powerful vector processing unit with lengthy vector registers. There is a class of LQCD computations which aims to solve the Dirac equation for the same lattice gauge fields multiple times. In this case one can gain additional reuse of the gauge field, and potentially simplified vectorization, through applying Dslash to several vectors at the same time; often referred to as Multiple Right Hand Side (MRHS) Dslash. Previous work in this area focused on applying this approach to the Highly Improved Staggered operator, and in cases to the multi-Grid algorithm for Wilson clover fermions. Here, we present two approaches to MRHS Dslash: a vector register blocking approach and one using the software package QPhiX with a custom Code Generator for low-level intrinsics. We explore the Knights Corner instruction set, memory hierarchy, and threading capabilities. For standard lattice sizes used in testing, we observed significant speedups, with sustained performance of over 700 GFLOPS (single precision) in one instance. We achieved up to 27% of theoretical peak performance compared to a maximum of 12% obtained by the previous single RHS method in QPhiX.

Joó, Bálint↗

Analysis of CDC Bioreactor Internal Thermal Measurements and Sample Coupon Temperatures

The Center for Disease Control and Prevention (CDC) bioreactor is an integral laboratory tool for the Environmental Control and Life Support Systems (ECLSS) biofilm formation and growth research program. Critical to this research is the need to adjust and maintain various surface temperatures of the coupons housed within the CDC Bioreactor. The purpose of this study was to provide quantitative temperature gradient information when the CDC Bioreactor was operating according to several process scenarios. Two primary process parameters were evaluated. For the first set of test parameters, the liquid level was maintained at 350 mL, with an inlet flowrate of 0.1 mL/min, 1 mL/min, 10 mL/min. The liquid was allowed to gravity drain out of the outlet spout. For the second set of test parameters, the liquid level within the reactor was maintained at 550 mL, with an inlet flow of 0.1 mL/min, 1 mL/min, 10 mL/min and draining intermittently controlled to 0.8 mL/min to maintain the 550 mL level. Due to the placement of the thermocouple in the reactor, a difference in temperature occurred between the coupon surfaces and target Bioreactor temperature when operated according to the first set of test parameters. When the reactor was operated according to the second set of parameters, which resulted in the thermocouple being submerged, the temperature gradient was eliminated. The results demonstrated minimal temperature gradient between the top and bottom coupon surfaces for coupons placed in a single rod within the CDC Bioreactor for both sets of test parameters evaluated. The collection of this information helped to explain previous ECLSS biofilm formation test runs, along with providing guidance on best operating practices for future ECLSS experiments. The placement of the thermocouple also helps to explain the challenge of achieving and maintaining bulk liquid temperatures when biofilm is grown according to the standardized test methods.

CDC Bioreactor↗