Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level set method”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The case considered was the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge.

Davoudzadeh, Farhad↗
NASA NTRS1990

Two- and three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction (BVI) problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The two-dimensional study considered the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The 2-D BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge. The three-dimensional interaction study considered the case of a vortex at 20 deg incidence to the blade leading edge. Although the qualitative results were similar to that of the two-dimensional interaction, details clearly showed the three-dimensional nature of the interaction process.

Davoudzadeh, F.↗
NASA NTRS2017

Cryogenic Fluid Management Technology Development Roadmaps

Advancement in Cryogenic Fluid Management (CFM) Technologies is essential for achieving NASA's future long duration missions. Propulsion systems utilizing cryogens are necessary to achieve mission success. Current State Of the Art (SOA) CFM technologies enable cryogenic propellants to be stored for several hours. However, some envisioned mission architectures require cryogens to be stored for two years or longer. The fundamental roles of CFM technologies are long term storage of cryogens, propellant tank pressure control and propellant delivery. In the presence of heat, the cryogens will "boil-off" over time resulting in excessive pressure buildup, off-nominal propellant conditions, and propellant loss. To achieve long term storage and tank pressure control, the CFM elements will intercept and/or remove any heat from the propulsion system. All functions are required to perform both with and without the presence of a gravitational field. Which CFM technologies are required is a function of the cryogens used, mission architecture, vehicle design and propellant tank size. To enable NASA's crewed mission to the Martian surface, a total of seventeen CFM technologies have been identified to support an In-Space Stage and a Lander/Ascent Vehicle. Recognizing that FY2020 includes a Decision Point regarding the In-Space Stage Architecture, a set of CFM Technology Development Roadmaps have been created identifying the current Technology Readiness Level (TRL) of each element, current technology "gaps", and existing technology development efforts. The roadmaps include a methodical approach and schedule to achieve a flight demonstration in FY2023, hence maturing CFM technologies to TRL 7 for infusion into the In-Space Stage Preliminary Design.

Stephens, J. R.↗
NASA NTRS2018

NASA's Cryogenic Fluid Management Technology Development Roadmaps

The maturation of Cryogenic Fluid Management (CFM) Technologies is essential for achieving NASA's future long duration missions. Propulsion systems utilizing cryogens are necessary to achieve NASA's exploration missions to the moon, Mars, and beyond. Current State Of the Art (SOA) CFM technologies enable cryogenic propellants to be stored for several hours prior to their use. However, some envisioned mission architectures require that cryogens to be stored for two years or longer. The fundamental roles of CFM technologies are long term storage of cryogens, propellant tank pressure control and propellant delivery. In the presence of heat, the cryogens will "boil-off" over time resulting in excessive pressure buildup, off-nominal propellant conditions for engine consumption, and propellant loss. To achieve long term storage and tank pressure control, the CFM elements will intercept and/or remove any heat from the propulsion system. All functions are required to be performed both with and without the presence of a gravitational field. Which CFM technologies are required is a function of the cryogens used, mission architecture, vehicle design and propellant tank size. To enable NASA's crewed missions beyond Low Earth Orbit, a total of twenty-seven CFM technologies have been identified to support various In-Space Stages and Lander/Ascent Vehicles. A set of CFM Technology Development Roadmaps have been created identifying the current Technology Readiness Level (TRL) of each element, current technology "gaps", and existing technology development efforts. The roadmaps include a methodical approach and schedule to achieve a flight demonstration, hence maturing CFM technologies to TRL 6/7 for infusion into the NASA's exploration elements. Additionally, a survey of the aerospace industry was completed to understand their views on the various technologies and how they would be infused. This does not cover all possible technologies, but rather those that are of interest to NASA specifically.

Cryogenic fluid management↗
NASA NTRSDate not supplied

Towards Integrated Computational Materials Engineering for Quantifying Performance Impacts of Microstructure and Defect Interactions in Powder Bed Fusion Parts

Powder bed fusion (PBF) additive manufacturing (AM) enables the creation of parts with complexity and functionality levels that were previously impossible with traditional manufacturing methods. By modifying the laser power, hatch spacing, or the numerous other processing parameters, the PBF process supports the production of a wide set of materials and geometries. However, that same process parameter design flexibility causes the process-design space of PBF to be massive and expensive to explore experimentally. Another challenge is quality variation across a build. As a part is being built, geometric variance between locations, such as at a thin-wall section vs. the bulk material, may cause the specified processing parameters to no longer be acceptable for producing defect-free printing. Furthermore, if the processing parameters deviate during the print process, it is difficult to assess if the part will still perform satisfactorily. Integrated Computational Materials Engineering (ICME) provides a way to understand and address these various challenges. This talk will present advancements in process-structure simulations of PBF at NASA Langley Research Center. The ability to simulate grain-scale PBF microstructures using the Physically Based Monte Carlo method will be demonstrated and compared to experimental measurements. Techniques for simulating three-dimensional lack-of-fusion and keyhole porosity defects based on the specific processing conditions and approaches for integrating the two porosity prediction techniques alongside the computational microstructure evolution models will be shown. Finally, the integration of simulated PBF microstructures, embedded process defects, and crystal plasticity finite element models to elucidate the interaction of porosity and microstructure on micromechanical fields will be demonstrated. These integrated techniques demonstrate an example of using ICME to relate processing to performance for PBF AM materials. With continued maturity, it is hoped that such ICME approaches will lead to next-generation computational-materials supported qualification and certification of AM parts.

Additive manufacturing↗
DOE OSTI2020

Assessing United States County-Level Exposure for Research on Tropical Cyclones and Human Health

Tropical cyclone epidemiology can be advanced through exposure assessment methods that are comprehensive and consistent across space and time, as these facilitate multiyear, multistorm studies. Further, an understanding of patterns in and between exposure metrics that are based on specific hazards of the storm can help in designing tropical cyclone epidemiological research. a) Provide an open-source data set for tropical cyclone exposure assessment for epidemiological research; and b) investigate patterns and agreement between county-level assessments of tropical cyclone exposure based on different storm hazards. We created an open-source data set with data at the county level on exposure to four tropical cyclone hazards: peak sustained wind, rainfall, flooding, and tornadoes. The data cover all eastern U.S. counties for all land-falling or near-land Atlantic basin storms, covering 1996–2011 for all metrics and up to 1988–2018 for specific metrics. We validated measurements against other data sources and investigated patterns and agreement among binary exposure classifications based on these metrics, as well as compared them to use of distance from the storm’s track, which has been used as a proxy for exposure in some epidemiological studies. Our open-source data set was typically consistent with data from other sources, and we present and discuss areas of disagreement and other caveats. Over the study period and area, tropical cyclones typically brought different hazards to different counties. Therefore, when comparing exposure assessment between different hazard-specific metrics, agreement was usually low, as it also was when comparing exposure assessment based on a distance-based proxy measurement and any of the hazard-specific metrics. Our results provide a multihazard data set that can be leveraged for epidemiological research on tropical cyclones, as well as insights that can inform the design and analysis for tropical cyclone epidemiological research.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Complete collision data set for electrons scattering on molecular hydrogen and its isotopologues: I. Fully vibrationally-resolved electronic excitation of H 2 $X^1Σ^+_g)$

Here, we present a comprehensive set of vibrationally-resolved cross sections for electron-impact electronic excitation of molecular hydrogen suitable for implementation in collisional-radiative models. The adiabatic-nuclei molecular convergent close-coupling method is used to calculate cross sections for excitation of all bound vibrational levels and dissociative excitation of the B 1 Σ u + , C 1 Π u , E F 1 Σ g + , B ′ 1 Σ u + , G K 1 Σ g + , I 1 Π g , J 1 Δ g , D 1 Π u , H 1 Σ g + , b 3 Σ u + , c 3 Π u , a 3 Σ g + , e 3 Σ u + , d 3 Π u , h 3 Σ g + , g 3 Σ g + , i 3 Π g , and j 3 Δ g electronic states from all –14 bound vibrational levels of the ground electronic ($\text{X}$ $^1&#x3A3^+_g$) state. The data set consists of cross sections from threshold to 500 eV for over 5000 transitions, representing all possible electronic and vibrational transitions between the state and the –3 singlet and triplet states (where refers to the united-atoms-limit principle quantum number). The cross sections are presented in graphical form and provided as both numerical values and analytic fit functions in supplementary data files.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Custom Accessors: Enabling Scalable Data Ingestion, (Re-)Organization, and Analysis on Distributed Systems

The emerging class of high velocity and high volume data analytic workflows comprise interwoven data ingestion, organization, and processing stages, with ingestion and organization steps often contributing comparable or even higher computational costs than actual processing steps. Since complex workflows consist of a variety of phases that view and use data differently, being able to construct efficient, scalable, distributed data structures (arrays, vectors, sets, maps, and multi-maps) is essential and requires custom methods to extend and shrink containers, analyze and position data, and, maintain globallyconsistent meta-data. In this paper, we propose a novel datastructure access paradigm based on the concept of Accessors. At a high level, accessors are customizable callable objects that can modify the behavior of insert, read, update, and delete operations for distributed containers while preserving atomicity guarantees. Accessors provide a very clean and natural way to implement a variety of programming patterns, e.g., conditional insertion/deletion and cascading computations, which would be otherwise hard (or even impossible) to express in parallel and distributed settings without using locks. We demonstrate the practicality and usefulness of our approach with two representative use cases and study the performance of these applications on a distributed High-Performance Computing system. Our analysis highlights that our proposed abstraction allows for an effective overlapping and concurrent execution of different workflow steps (e.g., data ingestion and analysis), which in a conventional analytics pipeline would execute sequentially, contributing cumulatively to the overall latency.

Castellana, Vito G. [BATTELLE (PACIFIC NW LAB)] (O↗
DOE OSTI2019

Improving efficiency of semi‐direct møller–plesset second‐order perturbation methods through oversubscription on multiple nodes

The purpose of this work is to evaluate the efficacy of oversubscription, at the 1 n , 2 n , and 3 n levels for n physical cores, on semi‐direct MP2 methods within NWChem when using two and three Intel nodes. Semi‐direct MP2 energy and gradient calculations were performed on chemical systems ranging from 824 to 1626 basis functions using the cc‐pVDZ basis set. Wall times for semi‐direct MP2 energies were reduced by as much as 36% using two nodes and 44% using three nodes compared to no oversubscription. Total energy consumed by the CPU and DRAM was also reduced by as much as 12% using two nodes and as much as 20% using three nodes when oversubscribing. MP2 gradient wall times improved by as much as 16% using two nodes and 18% using three nodes compared to execution at the 1 n level; however, energy savings were insignificant. Intel performance‐counter data show a strong correlation between total wall time saved and less time spent in the idle state, indicating a more efficient use of the processors when oversubscribing. © 2019 Wiley Periodicals, Inc.

Fought, Ellie L.↗
DOE OSTI2016

From plane waves to local Gaussians for the simulation of correlated periodic systems

Here, we present a simple, robust, and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cutoff radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artifacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid, and water adsorption on a LiH surface, at the level of second-order Møller–Plesset perturbation theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Predicting motion sickness during parabolic flight

BACKGROUND: There are large individual differences in susceptibility to motion sickness. Attempts to predict who will become motion sick have had limited success. In the present study, we examined gender differences in resting levels of salivary amylase and total protein, cardiac interbeat intervals (R-R intervals), and a sympathovagal index and evaluated their potential to correctly classify individuals into two motion sickness severity groups. METHODS: Sixteen subjects (10 men and 6 women) flew four sets of 10 parabolas aboard NASA's KC-135 aircraft. Saliva samples for amylase and total protein were collected preflight on the day of the flight and motion sickness symptoms were recorded during each parabola. Cardiovascular parameters were collected in the supine position 1-5 days before the flight. RESULTS: There were no significant gender differences in sickness severity or any of the other variables mentioned above. Discriminant analysis using salivary amylase, R-R intervals and the sympathovagal index produced a significant Wilks' lambda coefficient of 0.36, p=0.006. The analysis correctly classified 87% of the subjects into the none-mild sickness or the moderate-severe sickness group. CONCLUSIONS: The linear combination of resting levels of salivary amylase, high-frequency R-R interval levels, and a sympathovagal index may be useful in predicting motion sickness severity.

manned↗
NASA NTRS2002

Predicting Motion Sickness During Parabolic Flight

Background: There are large individual differences in susceptibility to motion sickness. Attempts to predict who will become motion sick have had limited success. In the present study we examined gender differences in resting levels of salivary amylase and total protein, cardiac interbeat intervals (R-R intervals), and a sympathovagal index and evaluated their potential to correctly classify individuals into two motion sickness severity groups. Methods: Sixteen subjects (10 men and 6 women) flew 4 sets of 10 parabolas aboard NASA's KC-135 aircraft. Saliva samples for amylase and total protein were collected preflight on the day of the flight and motion sickness symptoms were recorded during each parabola. Cardiovascular parameters were collected in the supine position 1-5 days prior to the flight. Results: There were no significant gender differences in sickness severity or any of the other variables mentioned above. Discriminant analysis using salivary amylase, R-R intervals and the sympathovagal index produced a significant Wilks' lambda coefficient of 0.36, p= 0.006. The analysis correctly classified 87% of the subjects into the none-mild sickness or the moderate-severe sickness group. Conclusions: The linear combination of resting levels of salivary amylase, high frequency R-R interval levels, and a sympathovagal index may be useful in predicting motion sickness severity.

Harm, Deborah L.↗
DOE OSTI2018

A hybrid wavelet de-noising and Rank-Set Pair Analysis approach for forecasting hydro-meteorological time series

Highlights: • A new hybrid approach is proposed to improve forecasts of hydro-meteorological time series. • Rank-Set Pair Analysis combined with wavelet de-noising markedly improves forecasting accuracy. • The performance of the proposed approach proves best among its present competitors even when the extreme value occurs. Accurate, fast forecasting of hydro-meteorological time series is presently a major challenge in drought and flood mitigation. This paper proposes a hybrid approach, wavelet de-noising (WD) and Rank-Set Pair Analysis (RSPA), that takes full advantage of a combination of the two approaches to improve forecasts of hydro-meteorological time series. WD allows decomposition and reconstruction of a time series by the wavelet transform, and hence separation of the noise from the original series. RSPA, a more reliable and efficient version of Set Pair Analysis, is integrated with WD to form the hybrid WD-RSPA approach. Two types of hydro-meteorological data sets with different characteristics and different levels of human influences at some representative stations are used to illustrate the WD-RSPA approach. The approach is also compared to three other generic methods: the conventional Auto Regressive Integrated Moving Average (ARIMA) method, Artificial Neural Networks (ANNs) (BP-error Back Propagation, MLP-Multilayer Perceptron and RBF-Radial Basis Function), and RSPA alone. Nine error metrics are used to evaluate the model performance. Compared to three other generic methods, the results generated by WD-REPA model presented invariably smaller error measures which means the forecasting capability of the WD-REPA model is better than other models. The results show that WD-RSPA is accurate, feasible, and effective. In particular, WD-RSPA is found to be the best among the various generic methods compared in this paper, even when the extreme events are included within a time series.

58 GEOSCIENCES↗
DOE OSTI2025

Engineering Polyketide Stereocenters with Ketoreductase Domain Exchanges

Polyketide synthases (PKSs) are versatile biosynthetic megasynthases capable of producing a diverse range of natural products with many applications, including in pharmaceuticals. The stereochemical precision of PKSs makes them a powerful tool for engineering tailored, unnatural polyketides; however, modifying the stereocenters of a PKS product while maintaining production levels remains a significant challenge. In this study, we systematically tested and evaluated strategies for ketoreductase (KR) domain exchanges, the domain responsible for setting stereocenters of polyketide products. After first optimizing the method for KR exchanges, we then performed 44 KR domain exchanges on three different PKSs to obtain high production of all four stereoisomers in vivo. By testing both one- and two-module PKS systems, we investigated how downstream modules process intermediates with altered stereochemistry and found that the configuration of the α-substituents was critical for gatekeeping by the ketosynthase (KS). To overcome this constraint, we investigated two different strategies for altering the KS domain, including introducing targeted mutations in the downstream KS, and exploring boundaries in exchanging the entire functional unit from the donor PKS. Both strategies successfully modified the KS stereocontrol with distinct trade-offs; the functional unit exchange resulted in higher titer improvements, though it was more likely to break the entire PKS. This study demonstrates a comprehensive approach to successfully engineering all four stereochemical configurations in multiple PKS systems, advancing our understanding of and ability to rationally modify polyketide stereochemistry through multiple engineering strategies.

Keiser, Leah S. [Joint BioEnergy Institute (JBEI),↗
DOE OSTI2023

Fringing Analysis and Simulation for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time

Abstract The presence of fringing in astronomical Charge-Coupled Device (CCD) images will have an impact on photometric quality and measurements. Yet its impact on the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) has not been fully studied. We present a detailed study on fringing for CCDs already implemented on the Rubin Observatory LSST Camera’s focal plane. After making physical measurements and knowing the compositions, we developed a model for the e2v CCDs. We present a method to fit for the internal height variation of the epoxy layer within the sensors based on fringing measurements in a laboratory setting. This method is generic enough that it can be easily modified to work for other CCDs. Using the derived fringing model, we successfully reproduce comparable fringing amplitudes that match the observed levels in images taken by existing telescopes with different optical designs. This model is then used to forecast the expected level of fringing in a single LSST y -band sky background exposure with Rubin telescope optics in the presence of a realistic time-varying sky spectrum. The predicted fringing amplitude in LSST images ranges from 0.04% to 0.2% depending on the location of a CCD on the focal plane. We find that the predicted variation in surface brightness caused by fringing in LSST y -band sky background images is about 0.6 μ Jy arcsec − 2 , which is 40 times larger than the current measurement error. We conclude that it is necessary to include fringing correction in the Rubin’s LSST image processing pipeline.

Astronomy & Astrophysics↗
DOE OSTI2017

Recent advances in stellarator optimization

Computational optimization has revolutionized the field of stellarator design. To date, optimizations have focused primarily on optimization of neoclassical confinement and ideal MHD stability, although limited optimization of other parameters has also been performed. Here, we outline a select set of new concepts for stellarator optimization that, when taken as a group, present a significant step forward in the stellarator concept. One of the criticisms that has been leveled at existing methods of design is the complexity of the resultant field coils. Recently, a new coil optimization code—COILOPT++, which uses a spline instead of a Fourier representation of the coils,—was written and included in the STELLOPT suite of codes. The advantage of this method is that it allows the addition of real space constraints on the locations of the coils. The code has been tested by generating coil designs for optimized quasi-axisymmetric stellarator plasma configurations of different aspect ratios. As an initial exercise, a constraint that the windings be vertical was placed on large major radius half of the non-planar coils. Further constraints were also imposed that guaranteed that sector blanket modules could be removed from between the coils, enabling a sector maintenance scheme. Results of this exercise will be presented. New ideas on methods for the optimization of turbulent transport have garnered much attention since these methods have led to design concepts that are calculated to have reduced turbulent heat loss. We have explored possibilities for generating an experimental database to test whether the reduction in transport that is predicted is consistent with experimental observations. Thus, a series of equilibria that can be made in the now latent QUASAR experiment have been identified that will test the predicted transport scalings. Fast particle confinement studies aimed at developing a generalized optimization algorithm are also discussed. A new algorithm developed for the design of the scraper element on W7-X is presented along with ideas for automating the optimization approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Toward an optimal online checkpoint solution under a two-level HPC checkpoint model

The traditional single-level checkpointing method suffers from significant overhead on large-scale platforms. Hence, multilevel checkpointing protocols have been studied extensively in recent years. The multilevel checkpoint approach allows different levels of checkpoints to be set (each with different checkpoint overheads and recovery abilities), in order to further improve the fault tolerance performance of extreme-scale HPC applications. How to optimize the checkpoint intervals for each level, however, is an extremely difficult problem. In this paper, we construct an easy-to-use two-level checkpoint model. Checkpoint level 1 deals with errors with low checkpoint/recovery overheads such as transient memory errors, while checkpoint level 2 deals with hardware crashes such as node failures. Compared with previous optimization work, our new optimal checkpoint solution offers two improvements: (1) it is an online solution without requiring knowledge of the job length in advance, and (2) it shows that periodic patterns are optimal and determines the best pattern. We evaluate the proposed solution and compare it with the most up-to-date related approaches on an extreme-scale simulation testbed constructed based on a real HPC application execution. Simulation results show that our proposed solution outperforms other optimized solutions and can improve the performance significantly in some cases. Specifically, with the new solution the wall-clock time can be reduced by up to 25.3% over that of other state-of-the-art approaches. Lastly, a brute-force comparison with all possible patterns shows that our solution is always within 1% of the best pattern in the experiments.

97 MATHEMATICS AND COMPUTING↗