Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Based Procedures”

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 775 records · Page 43

A Survey of Colormaps in Visualization

Colormaps are a vital method for users to gain insights into data in a visualization. With a good choice of colormaps, users are able to acquire information in the data more effectively and efficiently. In this survey, we attempt to provide readers with a comprehensive review of colormap generation techniques and provide readers a taxonomy which is helpful for finding appropriate techniques to use for their data and applications. Specifically, we first briefly introduce the basics of color spaces including color appearance models. In the core of our paper, we survey colormap generation techniques, including the latest advances in the field by grouping these techniques into four classes: procedural methods, user-study based methods, rule-based methods, and data-driven methods; we also include a section on methods that are beyond pure data comprehension purposes. We then classify colormapping techniques into a taxonomy for readers to quickly identify the appropriate techniques they might use. Furthermore, a representative set of visualization techniques that explicitly discuss the use of colormaps is reviewed and classified based on the nature of the data in these applications. Lastly, our paper is also intended to be a reference of colormap choices for readers when they are faced with similar data and/or tasks.

97 MATHEMATICS AND COMPUTING↗

Resonance scattering and radiation force calculations for an elastic cylinder using the translational addition theorem for cylindrical wave functions

The standard Resonance Scattering Theory (RST) of plane waves is extended for the case of any two-dimensional (2D) arbitrarily-shaped monochromatic beam incident upon an elastic cylinder with arbitrary location using an exact methodology based on Graf’s translational addition theorem for the cylindrical wave functions. The analysis is exact as it does not require numerical integration procedures. The formulation is valid for any cylinder of finite size and material that is immersed in a nonviscous fluid. Partial-wave series expansions (PWSEs) for the incident, internal and scattered linear pressure fields are derived, and the analysis is further extended to obtain generalized expressions for the on-axis and off-axis acoustic radiation force components. The wave-fields are expressed using generalized PWSEs involving the beam-shape coefficients (BSCs) and the scattering coefficients of the cylinder. The off-axial BSCs are expressed analytically in terms of an infinite PWSE with emphasis on the translational offset distance d. Numerical computations are considered for a zeroth-order quasi-Gaussian beam chosen as an example to illustrate the analysis. Acoustic resonance scattering directivity diagrams are calculated by subtracting an appropriate background from the expression of the scattered pressure field. In addition, computations for the radiation force exerted on an elastic cylinder centered on the axis of wave propagation of the beam, and shifted off-axially are analyzed and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sparse thermal data for cellular automata modeling of grain structure in additive manufacturing

Grain growth in the wake of the melt pool formed during alloy-based additive manufacturing (AM) is complex and multifaceted, depending on parameters governing heat transport, fluid flow, and solidification itself. Cellular automata (CA) models have proven effective in providing computationally efficient and physically sound predictions of grain structure for several AM problems, but their efficiency is tied to the performance of heat transport models. CA models use only a small portion of the problem's temperature data (near the moving melt pool boundary), and much of the CA calculations do not affect the final result due to re-melting of material. Coupling of and communication between heat transport and solidification models, and eliminating operations irrelevant towards final grain structure prediction, will be necessary for using these methods for efficient simulation of large parts. In this work, we introduce a procedure of decoupling the CA from temperature field simulation, using files of relevant temperature data written by the heat transport model. This approach is validated against the standard coupling approach using data obtained through the computational fluid dynamics software OpenFOAM. Negligible differences are seen in grain size, volume, and texture distributions for multilayer simulation of test problems, while the quantity of temperature data for these test problems was reduced by four orders of magnitude (from 100s of GB to 10s of MB) and the code performance sped up by a factor of around 50. Furthermore, variability in microstructure as a function of cell size, substrate, time step, and nucleation parameters is studied, and it is found that cell sizes less than or equal to 1.67 μm and sufficiently small time steps yield Αstatistically equivalent microstructures. Finally, a potential use case for this CA approach—the layer-wise convergence in grain structure starting from extremes in initial grain size—is examined. This approach's ability to simulate expected trends in nucleation and epitaxial grain growth for large regions of microstructure, simulated independently of heat transport models themselves, should prove useful for investigation of various microstructure uncertainties and prediction of part-scale experimental results.

36 MATERIALS SCIENCE↗

DEVELOPMENT OF A VERIFICATION BENCHMARK SUITE FOR THE ATR REACTOR PHYSICS UPGRADE

The result of a numerical algorithm is only as good as the effort made toward its verification. Comprehensive computer codes, especially for application to nuclear reactors, require in-depth assessment of the numerical methods used and their computational implementation. This presentation describes a suite of semi-analytical benchmarks for the ATR Reactor Physics Upgrade as part of the ATR Life Extension Project (LEP). Several neutron transport codes including deterministic codes HELIOS, NEWT and probabilistic codes MCNP6, SERPANT and MC21 are upgrade candidates. To properly evaluate code performance capabilities, code verification is an essential element. A credible verification program requires education to guide proper application, training in proper input preparation and output expectation and verification to instill confidence in the coded numerical procedures. The subset of the benchmark suite to be described here is primarily of the 1D transport variety (however, a 2D case does indeed exist) with the intent to provide unquestionably high order accurate semi-analytical numerical solutions to the transport equation. We base the analytical solutions on analytical Fourier series and discrete ordinates representations. Here we present two 1D benchmarks based on Fourier series.

99 GENERAL AND MISCELLANEOUS↗

Universal Workflow Language and Software Enable Geometric Learning and FAIR Scientific Protocol Reporting

Written language and conventional data structures for representing scientific procedures suffer from low process detail, often fail to accurately represent protocols, and lack universality. New strategies for the handling of experimental data are needed to provide viable process information for both humans and machines. In this work, we present the universal workflow language (UWL) and interface (UWLi). UWL is a findable, accessible, interoperable, and reusable (FAIR)-compatible, graph-based data architecture that can capture arbitrary scientific procedures through workflow representation, and UWLi is an accompanying software package for building, manipulating, and interpreting UWL entries. The UWL format was found to be highly effective in identifying deficiencies in the reported process details of high-impact, peer-reviewed scientific journals, and in simulated scenarios, the graph format was shown to be more effective than conventional methods in predictively modeling the outcome of diverse scientific protocols. Implementation of UWL could enable more accurate scientific communication and more impactful process datasets.

14 SOLAR ENERGY↗

Coupled-cluster based basis sets for valence correlation calculations

Novel basis sets are generated that target the description of valence correlation in atoms H through Ar. The new contraction coefficients are obtained according to the Atomic Natural Orbital (ANO) procedure from CCSD(T) (coupled-cluster singles and doubles with perturbative triples correction) density matrices starting from the primitive functions of Dunning et al. [J. Chem. Phys. 90, 1007 (1989); ibid. 98, 1358 (1993); ibid. 100, 2975 (1993)] (correlation consistent polarized valence X-tuple zeta, cc-pVXZ). The exponents of the primitive Gaussian functions are subject to uniform scaling in order to ensure satisfaction of the virial theorem for the corresponding atoms. These new sets, named ANO-VT-XZ (Atomic Natural Orbital Virial Theorem X-tuple Zeta), have the same number of contracted functions as their cc-pVXZ counterparts in each subshell. The performance of these basis sets is assessed by the evaluation of the contraction errors in four distinct computations: correlation energies in atoms, probing the density in different regions of space via 〈r{sup n}〉 (−3 ≤ n ≤ 3) in atoms, correlation energies in diatomic molecules, and the quality of fitting potential energy curves as measured by spectroscopic constants. All energy calculations with ANO-VT-QZ have contraction errors within “chemical accuracy” of 1 kcal/mol, which is not true for cc-pVQZ, suggesting some improvement compared to the correlation consistent series of Dunning and co-workers.

74 ATOMIC AND MOLECULAR PHYSICS↗

Distribution System Secondary Circuit Parameter Estimation for Model Calibration

To analyze and coordinate the operation of distribution systems with rapidly increasing amounts of PV, more accurate distribution system models are required, especially for the distribution system secondary (low-voltage) circuits down to the point of common coupling for distributed PV. There is a growing need for automated procedures to calibrate the distribution system secondary circuit models that are typically either not modeled at all or are modeled with a lower level of detail than the better modeled medium-voltage systems. This report presents an accurate, flexible, and computationally efficient method to use measurement data to estimate secondary circuit series impedance parameters in existing utility feeder models. The parameter estimation method assumes well-modeled primary circuit models, known secondary circuit topologies, and AMI active power, and reactive power measurements at all the loads in the secondary circuit. The method also requires AMI voltage measurement at most of the loads in the secondary circuit but can handle loads that do not have voltage measurements. No existing secondary circuit model information is needed, except for topology. The method is based on the well-known linearized voltage drop approximation and linear regression. The performance of the method is demonstrated on a three-phase test circuit with ten different secondary circuit topologies and on the Georgia Tech campus distribution system with AMI data. The developed method can be utilized to improve existing utility feeder models for more accurate analysis and operation with ubiquitous distributed PV interconnected on the low-voltage circuits.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Operating Hardware Impact on the Heat Transfer Properties of Windows

Despite significant advancements in fenestration technology in the last two decades, the thermal transmittance of fenestration products is still significantly higher than that of walls. This corresponds to 60% of the total energy loss of a modern building envelope through the windows. Hence, further development and improvements of fenestration products are necessary. Increasingly stringent codes and standards for fenestration stimulate industry to work on improved solutions. Thus, it is crucial that assessment techniques are able to account for innovations accurately. The thermal effects of non-continuous hardware in window frames are currently ignored by international rating procedures. A preliminary investigation conducted by our team showed significant performance degradation in two of the three out-opening casement profiles caused by the presence of operating hardware. Frames with the structure made of vinyl and fiberglass consist of many air cavities that are penetrated by operating hardware made of highly conductive materials. In these frames, in order to have an accurate assessment, it may be required to employ three-dimensional modeling due to the convective nature of heat transfer within the cavities. However, in this study, we demonstrate that the three-dimensional (3D) effects of non-continuous hardware can be approximated accurately with simpler two-dimensional (2D) simulations. We then develop a simplified model based on weighted average capable of replacing the time- and computation-intensive 3D simulations with 2D simulations and validate it against market available frames and their corresponding hardware. Validation results show that our approximation technique results in discrepancies lower than 0.05 W/(m2K), or 3% of the total thermal transmittance. Thus, we conclude that simplified 2D simulation models may be used for predicting hardware impact in window frames with reasonable accuracy. As windows and glazing structures are becoming ever better thermally insulated, it is becoming even more important to be able to model the impact of the operating hardware on the total thermal performance in order to design the best windows possible and not let the operating hardware ruin an otherwise well-proven design, which is hence addressed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Empirical Validation of Building Energy Modeling using Flexible Research Platform

Recently, efforts have been made to provide reliable empirical data for ANSI/ASHRAE Standard 140, Standard Method of Test for the Evaluation of Building Energy Analysis Computer Programs, to enable improved accuracy of building energy model (BEM) engines and improved characterization of their accuracy. The motivation for this effort is that the use of reliable empirical validation data sets in the evaluation of building energy modeling tools will lead to more consistent and validated simulation engines across all software vendors. This would expedite the use of building energy modeling in designing new buildings and retrofitting existing buildings, which delivers more energy efficient buildings.As part of a three-year multi-lab empirical validation project sponsored by U.S. DOE, this research project generated cooling season test plans by reviewing ASHRAE Standard 140, and the tests were performed based on the test plan. Finally, the experimental data sets were compared with the EnergyPlus model to demonstrate the validation procedure.

Im, Piljae↗

Progressive collapse analysis of stiffened panel structure with consideration of actual welding distortion and residual stresses

To date, inclusion of welding distortion and residual stresses in determining the strength of marine structures is based on empirical formulations due to welding complexity. In this study, a more feasible numerical model was developed from the perspective of inherent strain theory. The accuracy of this simplified numerical model was confirmed by transient thermo-elastic–plastic finite element analysis using a novel dynamic mesh refining and iterative substructure method which is 50 times faster than conventional finite element method. Welding distortion by the two numerical approaches was validated with experimental measurement. The collapse analyses indicated that the inherent strain model predicted a compressive strength and load-displacement curve with high accuracy, as validated by the transient thermo-mechanical model. The global stiffness and the ultimate strength of the welded structure would be overestimated by about 22% if no welding effects were considered. The proposed computation procedure for determining the ultimate strength of welded structures is accurate, simple, and practical and enables improvement of ultimate strength through optimization of welding joint design to reduce the welding residual stress and distortion.

Huang, Hui↗

Infections in patients using ventricular-assist devices: Comparison of the diagnostic performance of {sup 18}F-FDG PET/CT scan and leucocyte-labeled scintigraphy

Background: The usage of left-ventricular-assist device (LVAD) is increasing in patients presenting with advanced heart failure. However, device-related infections are a challenge to recognize and to treat, with an important morbidity and mortality rate. The role of nuclear medicine imaging remains not well established for LVAD infections. The present study compared the accuracy of positron emission tomography/computed tomography with {sup 18}F-fludeoxyglucose ({sup 18}F-FDG PET/CT) and radiolabeled leucocyte scintigraphy for the diagnosis of infections in patients supported with a continuous-flow LVAD. Methods: From a prospectively maintained database, we retrospectively analyzed the diagnostic performance of radiolabeled leucocyte scintigraphy and {sup 18}F-FDG PET/CT in 24 patients who had a LVAD with a suspected device-related infection. Both examinations were routinely performed in all patients. Infection was assessed by the International Society for Heart and Lung Transplantation criteria. Results: Twenty-four patients were included: 15 had a specific VAD infection (5 cardiac-LVAD and 10 driveline), 6 had a VAD-related infection, while 3 patients had a non-VAD-related infection. Sensitivity, specificity, positive predictive value, negative predictive value, and accuracy were 95.2%, 66.7%, 95.2%, 66.7%, and 91.6%, respectively, for {sup 18}F-FDG-PET; and 71.4%, 100%, 100%, 33.3%, and 75%, respectively, for leucocyte scintigraphy. {sup 18}F-FDG PET/CT showed significantly higher sensitivity (P = 0.01) than leucocyte scintigraphy. Conclusion: {sup 18}F-FDG PET/CT and radiolabeled leucocyte scintigraphy single-photon emission computed tomography carry high performance in the diagnostic of LVAD infections. {sup 18}F-FDG PET/CT shows significantly higher sensitivity and could be proposed as first-line nuclear medicine procedure. [French] French Abstract: Contexte: L’utilisation du dispositif d’assistance ventriculaire gauche (DAVG) est en utilisation croissante chez les patients avec insuffisance cardiaque avancée. Cependant, les infections liées au dispositif sont un défi à reconnaître et à traiter en raison d’ un taux de morbidité et de mortalité important. Le rôle de l’imagerie de médecine nucléaire n’est pas encore bien établi pour les infections secondaires au DAVG. La présente étude compare la précision de la tomographie par émission de positrons / tomodensitométrie au 18F-fluorodéoxyglucose (18F-FDG TEP/TDM) et la scintigraphie leucocytaire radio-marquée pour le diagnostic d’infections chez des patients porteurs par un DAVG à flux continu. Méthodes: A partir d’une base de données complétée de manière prospective nous avons analysé rétrospectivement les performances diagnostiques de la scintigraphie leucocytaire radio-marquée et de la TEP / TDM au 18F-FDG chez 24 patients ayant un DAVG avec une infection suspectée liée à l’appareillage. Les deux types d’examens scintigraphiques ont été systématiquement réalisé chez tous les patients. Le diagnostic d’infection fut établi selon les critères ISHLT. Résultats: Parmi les 24 patients inclus, 15 avaient une infection spécifique de l’appareillage AVG (5 au niveau de la pompe et/ou canule d’insertion myocardique) et 10 au niveau du driveline, 6 avaient une infection indirectement liée à l’appareillage AVG, tandis que 3 patients avaient une infection non liée au DAVG. La sensibilité, la spécificité, la valeur prédictive positive, la valeur prédictive négative et la précision étaient de 95,2%, 66,7%, 95,2%, 66,7% et 91,6%, respectivement, pour le 18F-FDG-PET et 71,4%, 100%, 100%, 33,3 %, et 75%, respectivement, pour la scintigraphie leucocytaire. La TEP / TDM au 18F-FDG a montré une sensitivité significativement plus élevée (p = 0,01) que la scintigraphie leucocytaire. Conclusion: 18F-FDG TEP/TDM et scintigraphie leucocytaire radio-marquée SPECT s’avèrent très efficaces pour le diagnostic des infections à DVAG. La TEP/TDM au 18F-FDG montre une sensibilité significativement plus élevée et pourrait être proposée en tant que la procédure de médecine nucléaire de premier choix. [Spanish] Spanish Abstract: Antecedentes: El uso de dispositivos de asistencia ventricular izquierda (DAVI) está incrementando en pacientes con insuficiencia cardiaca avanzada. Sin embargo, las infecciones relacionadas al dispositivo representan un reto diagnóstico y de tratamiento, con una importante tasa de morbilidad y mortalidad. El papel de los estudios de imagen de medicina nuclear en las infecciones de DAVI no se ha establecido claramente. El presente estudio comparó la precisión de la tomografía por emisión de positrones/tomografía computada con {sup 18}F-fluorodesoxiglucosa ({sup 18}F-FDG PET/CT) y la gammagrafía con leucocitos marcados para el diagnóstico de infecciones en pacientes utilizando un DAVI de flujo continuo. Métodos: Desde una base de datos obtenida de manera prospectiva, analizamos de manera retrospectiva el rendimiento diagnóstico de la gammagrafía de leucocitos marcados y del {sup 18}F-FDG PET/CT en 24 pacientes quienes tenían un DAVI con sospecha de infección relacionada al dispositivo. Ambos estudios fueron realizados rutinariamente en todos los pacientes. Se determinó la presencia de infección mediante los criterios de la ISHTL. Resultados: Se incluyeron veinticuatro pacientes, de los cuales 15 tuvieron una infección específica de DAV (5 con infecciones DAVI, y 10 con infección del driveline), 6 tuvieron una infección relacionada al DAV, y 3 pacientes tuvieron una infección no relacionada al DAV. La sensibilidad, especificidad, valor predictivo positivo, valor predictivo negativo y precisión diagnóstica fueron 95.2%, 66.7%, 95.2%, 66.7% y 91.6%, respectivamente para {sup 18}F-FDG PET/CT, y 71.4%, 100%, 100%, 33.3% y 75%, respectivamente para la gammagrafía de leucocitos marcados. El {sup 18}F-FDG PET/CT mostró una significativa mayor sensibilidad (p=0.01) que la gammagrafía con leucocitos marcados. Conclusión: El {sup 18}F-FDG PET/CT y la gammagrafía de leucocitos marcados tienen un alto desempeño en el diagnóstico de infecciones de DAVI. El {sup 18}F-FDG PET/CT muestra una significativa mayor sensibilidad y podría ser propuesto como el procedimiento de medicina nuclear de primera línea. [Chinese] Chinese Abstract: 背景: 左心室辅助装置 (LVAD) 越来越多的应用于重度心力衰竭患者。然而, 设备相关的感染有着重要的发病率和死亡率, 并且对其诊断和治疗均是一项挑战。对于 LVAD 相关感染, 核医学影像可扮演的角色一直没有完善。本研究针对装备有连续流式 LVAD 的感染病人,对比了 {sup 18}F-FDG PET/CT 造影术和放射性标记白细胞闪烁扫描法对感染诊断的准确性。 方法: 从一个前瞻性的数据库, 入选了24名疑似设备感染的 LVAD 病人, 并回顾性的分析了放射性标记白细胞闪烁扫描法和 {sup 18}F-FDG PET/CT 的诊断表现。所有病人都常规进行了上述两种检查。利用 ISHLT 标准对感染情况进行评估。 结果: 在入选的 24 名病人中, 15 名有明确的 VAD 感染 (5 名心脏的 LVAD 和 10 名动力传动系统装置), 6 名有 VAD 相关的感染,另外3名有 VAD 无关的感染。对于 {sup 18}F-FDG PET/CT 诊断方法,其敏感性,特异性,阳性预测值,阴性预测值和准确率分别为: 95.2%, 66.7%, 95.2%, 66.7%, 和 91.6%; 而对于白细胞闪烁扫描诊断方法, 其敏感性, 特异性, 阳性预测值,阴性预测值和准确率分别为: 71.4%, 100%, 100%, 33.3% 和 75%。相比于白细胞闪烁扫描, {sup 18}F-FDG PET/CT 明显展现出较高的敏感性 (p=0.01)。 结论: {sup 18}F-FDG PET/CT 和白细胞闪烁扫描 SPECT 对于 LVAD 感染的诊断具有很好的表现。{sup 18}F-FDG PET/CT 展现出较高的敏感性,可以作为核医学的首要检查。.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Convergence Analysis for an Online Data-Driven Feedback Control Algorithm

This paper presents convergence analysis of a novel data-driven feedback control algorithm designed for generating online controls based on partial noisy observational data. The algorithm comprises a particle filter-enabled state estimation component, estimating the controlled system’s state via indirect observations, alongside an efficient stochastic maximum principle-type optimal control solver. By integrating weak convergence techniques for the particle filter with convergence analysis for the stochastic maximum principle control solver, we derive a weak convergence result for the optimization procedure in search of optimal data-driven feedback control. Numerical experiments are performed to validate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗

GradDFT. A software library for machine learning enhanced density functional theory

Density functional theory (DFT) stands as a cornerstone method in computational quantum chemistry and materials science due to its remarkable versatility and scalability. Yet, it suffers from limitations in accuracy, particularly when dealing with strongly correlated systems. To address these shortcomings, recent work has begun to explore how machine learning can expand the capabilities of DFT: an endeavor with many open questions and technical challenges. In this work, we present GradDFT a fully differentiable JAX-based DFT library, enabling quick prototyping and experimentation with machine learning-enhanced exchange–correlation energy functionals. GradDFT employs a pioneering parametrization of exchange–correlation functionals constructed using a weighted sum of energy densities, where the weights are determined using neural networks. Moreover, GradDFT encompasses a comprehensive suite of auxiliary functions, notably featuring a just-in-time compilable and fully differentiable self-consistent iterative procedure. To support training and benchmarking efforts, we additionally compile a curated dataset of experimental dissociation energies of dimers, half of which contain transition metal atoms characterized by strong electronic correlations. The software library is tested against experimental results to study the generalization capabilities of a neural functional across potential energy surfaces and atomic species, as well as the effect of training data noise on the resulting model accuracy.

Chemistry↗

A Methodology for the Rapid Development of Scalable HPC Data Services

There is growing evidence in the scientifc computing community that parallel file systems are not sufficient for all HPC storage workloads. This realization has motivated extensive research in new storage system designs. The question of which alternative we should turn towards implies that there could be a single answer satisfying a wide range of very diverse applications. We argue that such a generic solution does not exist. Instead, custom data services should be designed and tailored to the needs of specific applications on specific hardware. Furthermore, custom data services should be designed in close collaboration with users. In this paper, we present a methodology for the rapid development of such data services. This methodology promotes the design of reusable building blocks that can be composed together efficiently through a runtime based on high-performance threading, tasking, and remote procedure calls. We illustrate the success of our methodology by showcasing three examples of data services designed from the same building blocks, yet targeting entirely different applications.

Dorier, Mattieu↗

Experimental and computing strategies in advanced material characterization problems

The mechanical characterization of materials relies more and more often on sophisticated experimental methods that permit to acquire a large amount of data and, contemporarily, to reduce the invasiveness of the tests. This evolution accompanies the growing demand of non-destructive diagnostic tools that assess the safety level of components in use in structures and infrastructures, for instance in the strategic energy sector. Advanced material systems and properties that are not amenable to traditional techniques, for instance thin layered structures and their adhesion on the relevant substrates, can be also characterized by means of combined experimental-numerical tools elaborating data acquired by full-field measurement techniques. In this context, parameter identification procedures involve the repeated simulation of the laboratory or in situ tests by sophisticated and usually expensive non-linear analyses while, in some situation, reliable and accurate results would be required in real time. The effectiveness and the filtering capabilities of reduced models based on decomposition and interpolation techniques can be profitably used to meet these conflicting requirements. This communication intends to summarize some results recently achieved in this field by the author and her co-workers. The aim is to foster further interaction between engineering and mathematical communities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine learning based job status prediction in scientific clusters

Large high-performance computing systems are built with increasing number of components with more CPU cores, more memory, and more storage space. At the same time, scientific applications have been growing in complexity. Together, they are leading to more frequent unsuccessful job statuses on HPC systems. From measured job statuses, 23.4% of CPU time was spent to the unsuccessful jobs. Here, we set out to study whether these unsuccessful job statuses could be anticipated from known job characteristics. To explore this possibility, we have developed a job status prediction method for the execution of jobs on scientific clusters. The Random Forests algorithm was applied to extract and characterize the patterns of unsuccessful job statuses. Experimental results show that our method can predict the unsuccessful job statuses from the monitored ongoing job executions in 99.8% the cases with 83.6% recall and 94.8% precision. Lastly, this prediction accuracy can be sufficiently high that it can be used to mitigation procedures of predicted failures.

97 MATHEMATICS AND COMPUTING↗

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

Multilayered triangular and quadrilateral flat shell elements based on the Refined Zigzag Theory

The paper presents a class of C 0 -continuous, flat shell elements based on the Refined Zigzag Theory (RZT) for the analysis of multilayered and curved composite and sandwich structures. In this study, the use of the interdependent interpolation strategy allows eliminating the shear-locking phenomenon and introducing the drilling rotation necessary to complete the set of classical nodal degrees of freedom (three displacements and three rotations). Additional kinematic variables are present in RZT, the zigzag rotations around the in-plane axes that measure the normal distortion typical of multilayered structures. An additional “drilling” zigzag rotation is therefore included among the nodal degrees of freedom in order to properly model curved and built-up structures. A stabilization procedure is adopted to suppress spurious zero-energy modes. A three-node triangular and a four-node quadrilateral flat shell element are formulated with 9 degrees of freedom per node. Example problems involving flat and curved multilayered structures are presented and discussed in order to assess the accuracy and convergence properties of the presented elements. Both static response predictions and free vibrations analyses are considered and the comparison is made with analytic RZT solutions, high-fidelity 3D finite element models and FSDT-based flat shell elements.

97 MATHEMATICS AND COMPUTING↗