Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Robustness”

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 973 records · Page 54

Hybrid super homogenization and discontinuity factor method for continuous finite element diffusion

This paper presents a novel homogenization equivalence technique aiming to simultaneously leveragethe simplicity of the Super Homogenization (SPH) method to reproduce reference reaction rates andthe ability to preserve reference leakage rates at desired surfaces through Discontinuity Factors (DF).The need for this new class of methods arises from the inability of the current state-of-the-art SPH technology to properly reproduce the reactor multiplication factor for problems with significant leakage. This work shows that this defect lies in the use of normalization factors in the SPH algorithm: while they are necessary with purely reflecting problems to ensure uniqueness of the solution, they introduce homogenization inconsistencies if at least one vacuum boundary condition is present. Two solutions to this problem are presented in this work: (i) to simply remove the normalization factor or (ii) to introduce additional degrees of freedom in the form of DFs in such a way that the normalization factors can still constrain the problem. While the former clearly offers unrivaled simplicity and can lead to very satisfactory results for leakage-dominated cores or if spatial restriction of the SPH regions is applied – as demonstrated for the High Temperature Test Reactor, the latter is more robust and does not require to scale the fission terms with the reference eigenvalue.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Peer-to-Peer Communication Trade-Offs for Smart Grid Applications: Preprint

Peer-to-peer energy management systems for smart grids require developers to consider the trade-offs between the amount of communication traffic generated and the quality and speed of convergence of the control algorithms that are deployed. Employing a fully connected communication causes messages to scale exponentially with the number of nodes, while using a sparse connectivity causes less information dissemination leading to degradation of the algorithm performance. The best communication topology for a particular application lies somewhere in between and often requires empirical evaluation by application designers. Existing methods do not put focus on the needs for smart grid applications, which is information dissemination throughout the network and they do not provide a flexible solution for application developers to prototype and deploy different topologies without modifying the application code. This paper introduces a configurable virtual communication topology framework TopLinkMgr, allowing users to specify any chosen communication topology and deploy peer-to-peer applications using it. It also introduces a self-adaptive, fault-tolerant topology management algorithm, Bounded Path Dissemination that can ensure the dissemination of information to all peers within a specified threshold for a sparsely connected topology. Experiments show that the algorithm improves on convergence speed and accuracy over state-of-the-art methods and is also robust against node failures. The results indicate the possibility of achieving a close-to optimal convergence without overloading the network allowing the realization of peer-to-peer control platforms covering larger and more complex power systems.

Bounded Path Dissemination↗

MPD thruster technology

MPD (MagnetoPlasmaDynamic) thrusters demonstrated between 2000 and 7000 seconds specific impulse at efficiencies approaching 40 percent, and were operated continuously at power levels over 500 kW. These demonstrated capabilities, combined with the simplicity and robustness of the thruster, make them attractive candidates for application to both unmanned and manned orbit raising, lunar, and planetary missions. To date, however, only a limited number of thruster configurations, propellants, and operating conditions were studied. The present status of MPD research is reviewed, including developments in the measured performance levels and electrode erosion rates. Theoretical studies of the thruster dynamics are also described. Significant progress was made in establishing empirical scaling laws, performance and lifetime limitations and in the development of numerical codes to simulate the flow field and electrode processes.

Myers, Roger M.↗

Adaptive digital beamforming for a CDMA mobile communications payload

In recent years, Spread-Spectrum Code Division Multiple Access (CDMA) has become a very popular access scheme for mobile communications due to a variety of reasons: excellent performance in multipath environments, high scope for frequency reuse, graceful degradation near saturation, etc. In this way, a CDMA system can support simultaneous digital communication among a large community of relatively uncoordinated users sharing a given frequency band. Nevertheless, there are also important problems associated with the use of CDMA. First, in a conventional CDMA scheme, the signature sequences of asynchronous users are not orthogonal and, as the number of active users increases, the self-noise generated by the mutual interference between users considerably degrades the performance, particularly in the return link. Furthermore, when there is a large disparity in received powers - due to differences in slant range or atmospheric attenuation - the non-zero cross-correlation between the signals gives rise to the so-called near-far problem. This leads to an inefficient utilization of the satellite resources and, consequently, to a drastic reduction in capacity. Several techniques were proposed to overcome this problem, such as Synchronized CDMA - in which the signature sequences of the different users are quasi-orthogonal - and power control. At the expense of increased network complexity and user coordination, these techniques enable the system capacity to be restored by equitably sharing the satellite resources among the users. An alternative solution is presented based upon the use of time-reference adaptive digital beamforming on board the satellite. This technique enables a high number of independently steered beams to be generated from a single phased array antenna, which automatically track the desired user signal and null the unwanted interference source. In order to use a time-reference adaptive antenna in a communications system, the main challenge is to obtain a reference signal highly correlated with the desired user signal and uncorrelated with the interferences. CDMA lends itself very easily to the generation of such a reference signal, thanks to the a priori knowledge of the user's signature sequence. First, the integration of an adaptive antenna in an asynchronous CDMA system is analyzed. The adaptive antenna system can provide increased interference rejection - much higher than that afforded by the code alone - and, since CDMA is mainly interference limited, any reduction in interference converts directly and linearly into an increase in capacity. Analyses and computer simulations are presented that show how an asynchronous CDMA system incorporating adaptive beamforming can provide at least as much capacity as a synchronous system. More importantly, the proposed concept allows the near-far effect to be mitigated without requiring a tight coordination of the users in terms of transmitted power control or network synchronization. The system is extremely robust to the near-far effect because the signals reaching the satellite from directions other than that of the desired user - which are likely to have different power levels - are adaptively canceled by the antenna. Finally, a payload architecture is presented that illustrates the practical implementation of this concept. This digital payload architecture demonstrates that with the advent of high performance CMOS digital processing, the on-board implementation of complex DSP techniques - in particular Digital Beamforming - has become possible, being most attractive for Mobile Satellite Communications.

Munoz-Garcia, Samuel G.↗

Robustness of Modeling of Out-of-Service Gas Mechanical Face Seal

Gas lubricated mechanical face seal are ubiquitous in many high performance applications such as compressors and gas turbines. The literature contains various analyses of seals having orderly face patterns (radial taper, waves, spiral grooves, etc.). These are useful for design purposes and for performance predictions. However, seals returning from service (or from testing) inevitably contain wear tracks and warped faces that depart from the aforementioned orderly patterns. Questions then arise as to the heat generated at the interface, leakage rates, axial displacement and tilts, minimum film thickness, contact forces, etc. This work describes an analysis of seals that may inherit any (i.e., random) face pattern. A comprehensive computer code is developed, based upon the Newton- Raphson method, which solves for the equilibrium of the axial force and tilting moments that are generated by asperity contact and fluid film effects. A contact mechanics model is incorporated along with a finite volume method that solves the compressible Reynolds equation. Results are presented for a production seal that has sustained a testing cycle.

Green, Itzhak↗

Characterizing Wheel-Soil Interaction Loads Using Meshfree Finite Element Methods: A Sensitivity Analysis for Design Trade Studies

A wheel experiencing sinkage and slippage events poses a high risk to planetary rover missions as evidenced by the mobility challenges endured by the Mars Exploration Rover (MER) project. Current wheel design practice utilizes loads derived from a series of events in the life cycle of the rover which do not include (1) failure metrics related to wheel sinkage and slippage and (2) performance trade-offs based on grouser placement/orientation. Wheel designs are rigorously tested experimentally through a variety of drive scenarios and simulated soil environments; however, a robust simulation capability is still in development due to myriad of complex interaction phenomena that contribute to wheel sinkage and slippage conditions such as soil composition, large deformation soil behavior, wheel geometry, nonlinear contact forces, terrain irregularity, etc. For the purposes of modeling wheel sinkage and slippage at an engineering scale, meshfree nite element approaches enable simulations that capture su cient detail of wheel-soil interaction while remaining computationally feasible. This study implements the JPL wheel-soil benchmark problem in the commercial code environment utilizing the large deformation modeling capability of Smooth Particle Hydrodynamics (SPH) meshfree methods. The nominal, benchmark wheel-soil interaction model that produces numerically stable and physically realistic results is presented and simulations are shown for both wheel traverse and wheel sinkage cases. A sensitivity analysis developing the capability and framework for future ight applications is conducted to illustrate the importance of perturbations to critical material properties and parameters. Implementation of the proposed soil-wheel interaction simulation capability and associated sensitivity framework has the potential to reduce experimentation cost and improve the early stage wheel design proce

terramechanics↗

Development of Neutron Energy Spectral Signatures for Passive Monitoring of Spent Nuclear Fuels in Dry Cask Storage

Demand for spent nuclear fuel dry casks as an interim storage solution has increased globally and the IAEA has expressed a need for robust safeguards and verification technologies for ensuring the continuity of knowledge and the integrity of radioactive materials inside spent fuel casks. Existing research has been focusing on “fingerprinting” casks based on count rate statistics to represent radiation emission signatures. The current research aims to expand to include neutron energy spectral information as part of the fuel characteristics. First, spent fuel composition data are taken from the Next Generation Safeguards Initiative Spent Fuel Libraries, representative for Westinghouse 17 x 17 PWR assemblies. The ORIGEN-S code then calculates the spontaneous fission and (α,n) emissions for individual fuel rods, followed by detailed MCNP simulations of neutrons transported through the fuel assemblies. A comprehensive database of neutron energy spectral profiles is to be constructed, with different enrichment, burn-up, and cooling time conditions. The end goal is to utilize the computational spent fuel library, predictive algorithm, and a pressurized 4 He scintillator to verify the spent fuel assemblies inside a cask. This work identifies neutron spectral signatures that correlate with the cooling time of spent fuel. Both the total and relative contributions from spontaneous fission and (α,n) change noticeably with respect to cooling time, due to the relatively short half-life (18 years) of the major neutron source 244 Cm. Identification of this and other neutron spectral signatures allows the characterization of spent nuclear fuels in dry cask storage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

QuantumNAT: quantum noise-aware training with noise injection, quantization and normalization

Parameterized Quantum Circuits (PQC) are promising towards quantum advantage on near-term quantum hardware. However, due to the large quantum noises (errors), the performance of PQC models has a severe degradation on real quantum devices. Take Quantum Neural Network (QNN) as an example, the accuracy gap between noise-free simulation and noisy results on IBMQ-Yorktown for MNIST-4 classification is over 60%. Existing noise mitigation methods are general ones without leveraging unique characteristics of PQC; on the other hand, existing PQC work does not consider noise effect. To this end, we present QuantumNAT, a PQC-specific framework to perform noise-aware optimizations in both training and inference stages to improve robustness. We experimentally observe that the effect of quantum noise to PQC measurement outcome is a linear map from noise-free outcome with a scaling and a shift factor. Motivated by that, we propose post-measurement normalization to mitigate the feature distribution differences between noise-free and noisy scenarios. Furthermore, to improve the robustness against noise, we propose noise injection to the training process by inserting quantum error gates to PQC according to realistic noise models of quantum hardware. Lastly, post-measurement quantization is introduced to quantize the measurement outcomes to discrete values, achieving the denoising effect. Extensive experiments on 8 classification tasks using 6 quantum devices demonstrate that QuantumNAT improves accuracy by up to 43%, and achieves over 94% 2-class, 80% 4-class, and 34% 10-class classification accuracy measured on real quantum computers. The code for construction and noise-aware training of PQC is available in the TorchQuantum library.

Wang, Hanrui↗

BISON Contact Algorithm Improvements in Support of Pellet Cladding Mechanical Interaction Modeling

The thermal and mechanical behavior of the gap between the fuel and cladding plays an extremely important role in determining the response of light water reactor (LWR) fuel under irradiation. This is especially true during situations when pellet-cladding mechanical interactions (PCMI) have a significant effect on the behavior of the fuel/cladding system. The BISON code, which is being used as the fuel performance simulation tool for the CASL program, employs contact detection and enforcement algorithms to realistically simulate the thermal and mechanical interactions across the fuel/cladding interaction in LWR fuel. These contact algorithms permit accurate representation of these interactions even when there is significant relative movement between the two surfaces, which is typically the case in many regions of an LWR fuel rod. This report summarizes work done during the 2017 fiscal year to further improve the accuracy and robustness of the contact enforcement algorithms in BISON. This work will benefit all BISON simulations, but will particularly benefit simulations scenarios where PCMI is of interest.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Deoptimization of FMDV P1 Region Results in Robust Serotype-Independent Viral Attenuation

Foot-and-mouth disease (FMD), caused by the FMD virus (FMDV), is a highly contagious disease of cloven-hoofed livestock that can have severe economic impacts. Control and prevention strategies, including the development of improved vaccines, are urgently needed to effectively control FMD outbreaks in endemic settings. Previously, we employed two distinct strategies (codon pair bias deoptimization (CPD) and codon bias deoptimization (CD)) to deoptimize various regions of the FMDV serotype A subtype A12 genome, which resulted in the development of an attenuated virus in vitro and in vivo, inducing varying levels of humoral responses. In the current study, we examined the versatility of the system by using CPD applied to the P1 capsid coding region of FMDV serotype A subtype, A24, and another serotype, Asia1. Viruses carrying recoded P1 (A24-P1Deopt or Asia1-P1Deopt) exhibited different degrees of attenuation (i.e., delayed viral growth kinetics and replication) in cultured cells. Studies in vivo using a mouse model of FMD demonstrated that inoculation with the A24-P1Deopt and Asia1-P1Deopt strains elicited a strong humoral immune response capable of offering protection against challenge with homologous wildtype (WT) viruses. However, different results were obtained in pigs. While clear attenuation was detected for both the A24-P1Deopt and Asia1-P1Deopt strains, only a limited induction of adaptive immunity and protection against challenge was detected, depending on the inoculated dose and serotype deoptimized. Our work demonstrates that while CPD of the P1 coding region attenuates viral strains of multiple FMDV serotypes/subtypes, a thorough assessment of virulence and induction of adaptive immunity in the natural host is required in each case in order to finely adjust the degree of deoptimization required for attenuation without affecting the induction of protective adaptive immune responses.

60 APPLIED LIFE SCIENCES↗

Petascale solvers for anisotropic PDEs in atmospheric modelling on GPU clusters

Memory bound applications such as solvers for large sparse systems of equations remain a challenge for GPUs. Fast solvers should be based on numerically efficient algorithms and implemented such that global memory access is minimised. To solve systems with trillions ( O ( 10 12 ) ) unknowns the code has to make efficient use of several million individual processor cores on large GPU clusters. We describe the multi-GPU implementation of two algorithmically optimal iterative solvers for anisotropic PDEs which are encountered in (semi-) implicit time stepping procedures in atmospheric modelling. In this application the condition number is large but independent of the grid resolution and both methods are asymptotically optimal, albeit with different absolute performance. In particular, an important constant in the discretisation is the CFL number; only the multigrid solver is robust to changes in this constant. We parallelise the solvers and adapt them to the specific features of GPU architectures, paying particular attention to efficient global memory access. We achieve a performance of up to 0.78 PFLOPs when solving an equation with 0.55 · 10 12 unknowns on 16384 GPUs; this corresponds to about 3% of the theoretical peak performance of the machine and we use more than 40% of the peak memory bandwidth with a Conjugate Gradient (CG) solver. Although the other solver, a geometric multigrid algorithm, has a slightly worse performance in terms of FLOPs per second, overall it is faster as it needs less iterations to converge; the multigrid algorithm can solve a linear PDE with half a trillion unknowns in about one second.

Computer Science↗

Opening the reheating box in multifield inflation

The robustness of multi-field inflation to the physics of reheating is investigated. In order to carry out this study, reheating is described in detail by means of a formalism which tracks the evolution of scalar fields and perfect fluids in interaction (the inflatons and their decay products). This framework is then used to establish the general equations of motion of the background and perturbative quantities controlling the evolution of the system during reheating. Next, these equations are solved exactly by means of a new numerical code. Moreover, new analytical techniques, allowing us to interpret and approximate these solutions, are developed. As an illustration of a physical prediction that could be affected by the micro-physics of reheating, the amplitude of non-adiabatic perturbations in double inflation is considered. It is found that ignoring the fine-structure of reheating, as usually done in the standard approach, can lead to differences as big as ~ 50%, while our semi-analytic estimates can reduce this error to ~ 10%. Finally, we conclude that, in multi-field inflation, tracking the perturbations through the details of the reheating process is important and, to achieve good precision, requires the use of numerical calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

Improving Incremental Balance in the GSI 3DVAR Analysis System

The Gridpoint Statistical Interpolation (GSI) analysis system is a unified global/regional 3DVAR analysis code that has been under development for several years at the National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center. It has recently been implemented into operations at NCEP in both the global and North American data assimilation systems (GDAS and NDAS). An important aspect of this development has been improving the balance of the analysis produced by GSI. The improved balance between variables has been achieved through the inclusion of a Tangent Linear Normal Mode Constraint (TLNMC). The TLNMC method has proven to be very robust and effective. The TLNMC as part of the global GSI system has resulted in substantial improvement in data assimilation both at NCEP and at the NASA Global Modeling and Assimilation Office (GMAO).

Errico, Ronald M.↗

Simultaneous Development and Robust Optimization of a Microstructure Dependent Material

Recent microstructure characterization techniques combined with Symbolic Regression(SR)analysis has been proven to generate white box plasticity models well suited for incorporation into FEA software.The current work builds upon those efforts and demonstrates the applicability of Sequential Monte-Carlo (SMC) methods within SR analysis to condense model development and robust optimization into a single, co-dependent process. In this project, SMC methods provide a mechanism through which the observed microstructure features and associated variability can be incorporated into the discovery phase of model development and simultaneously recover approximate parameter distributions through SR analysis. The demonstration utilized a data set consisting of tensile test results from a limited number of sample specimens with corresponding EBSD data from which microstructure features were characterized.The maximum threshold stress model in the Visco-Plastic Self-Consistent (VPSC) code developed by Los Alamos National Laboratories was calibrated using mechanical test data.Synthetic volume elements with statistically equivalent microstructure were generated with DREAM3Dbased on the observed EBSD data. VPSC was used to simulate the corresponding tensile test response for each of the synthetic volume elements. The simulated microstructure and tensile test data was used astraining datafor SMC-SR algorithm and the resulting model was validated with data from the original empirical data set.

Karl Garbrecht↗

Tuning the spectral emittance of α-SiC open-cell foams up to 1300 K with their macro porosity

A simple and robust analytical model is used to finely predict the spectral emittance under air up to 1300 K of α-SiC open-cell foams constituted of optically thick struts. The model integrates both the chemical composition and the macro-porosity and is valid only if foams have volumes higher than their Representative Elementary Volumes required for determining their emittance. Infrared emission spectroscopy carried out on a doped silicon carbide single crystal associated to homemade numerical tools based on 3D meshed images (Monte Carlo Ray Tracing code, foam generator) make possible to understand the exact role of the cell network in emittance. Finally, one can tune the spectral emittance of α-SiC foams up to 1300 K by simply changing their porosity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-E-P-05: Electronic Brachytherapy: A Physics Perspective On Field Implementation

Purpose: We want to summarize our experience implementing a successful program of electronic brachytherapy at several dermatology clinics with the help of a cloud based software to help us define the key program parameters and capture physics QA aspects. Optimally developed software helps the physicist in peer review and qualify the physical parameters. Methods: Using the XOFT™ Axxent™ electronic brachytherapy system in conjunction with a cloud-based software, a process was setup to capture and record treatments. It was implemented initially at about 10 sites in California. For dosimetric purposes, the software facilitated storage of the physics parameters of surface applicators used in treatment and other source calibration parameters. In addition, the patient prescription, pathology and other setup considerations were input by radiation oncologist and the therapist. This facilitated physics planning of the treatment parameters and also independent check of the dwell time. From 2013–2014, nearly1500 such calculation were completed by a group of physicists. A total of 800 patients with multiple lesions have been treated successfully during this period. The treatment log files have been uploaded and documented in the software which facilitated physics peer review of treatments per the standards in place by AAPM and ACR. Results: The program model was implemented successfully at multiple sites. The cloud based software allowed for proper peer review and compliance of the program at 10 clinical sites. Dosimtery was done on 800 patients and executed in a timely fashion to suit the clinical needs. Accumulated physics data in the software from the clinics allows for robust analysis and future development. Conclusion: Electronic brachytherapy implementation experience from a quality assurance perspective was greatly enhanced by using a cloud based software. The comprehensive database will pave the way for future developments to yield superior physics outcomes.

07 ISOTOPE AND RADIATION SOURCES↗

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗

A fast, matrix-based method to perform omnidirectional pressure integration

Abstract Experimentally-measured pressure fields play an important role in understanding many fluid dynamics problems. Unfortunately, pressure fields are difficult to measure directly with non-invasive, spatially resolved diagnostics, and calculations of pressure from velocity have proven sensitive to error in the data. Omnidirectional line integration methods are usually more accurate and robust to these effects as compared to implicit Poisson equations, but have seen slower uptake due to the higher computational and memory costs, particularly in 3D domains. This paper demonstrates how omnidirectional line integration approaches can be converted to a matrix inversion problem. This novel formulation uses an iterative approach so that the boundary conditions are updated each step, preserving the convergence behavior of omnidirectional schemes while also keeping the computational efficiency of Poisson solvers. This method is implemented in Matlab and also as a GPU-accelerated code in CUDA-C++. The behavior of the new method is demonstrated on 2D and 3D synthetic and experimental data. Three-dimensional grid sizes of up to 125 million grid points are tractable with this method, opening exciting opportunities to perform volumetric pressure field estimation from 3D PIV measurements.

42 ENGINEERING↗