Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model benchmark”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,747 records · Page 97

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Preparation of excited states for nuclear dynamics on a quantum computer

We study two different methods to prepare excited states on a quantum computer, a key initial step to study dynamics within linear response theory. The first method uses unitary evolution for a short time T = O(√1 - F) to approximate the action of an excitation operator Ô with fidelity F and success probability P ≈ 1 – F. The second method probabilistically applies the excitation operator using the Linear Combination of Unitaries (LCU) algorithm. We benchmark these techniques on emulated and real quantum devices, using a toy model for thermal neutron-proton capture. Despite its larger memory footprint, the LCU-based method is efficient even on current generation noisy devices and can be implemented at a lower gate cost than a naive analysis would suggest. Here, these findings show that quantum techniques designed to achieve good asymptotic scaling on fault tolerant quantum devices might also provide practical benefits on devices with limited connectivity and gate fidelity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

β-Delayed One and Two Neutron Emission Probabilities Southeast of 132 Sn and the Odd-Even Systematics in r-Process Nuclide Abundances

Here, the β-delayed one- and two-neutron emission probabilities (P 1n and P 2n ) of 20 neutron-rich nuclei with N≥82 have been measured at the RIBF facility of the RIKEN Nishina Center. P 1n of 130,131 Ag, 133,134 Cd, 135,136 In, and 138,139 Sn were determined for the first time, and stringent upper limits were placed on P2n for nearly all cases. β-delayed two-neutron emission (β2n) was unambiguously identified in 133 Cd and 135,136 In, and their P 2n were measured. Weak β2n was also detected from 137,138 Sn. Our results highlight the effect of the N=82 and Z=50 shell closures on β-delayed neutron emission probability and provide stringent benchmarks for newly developed macroscopic-microscopic and self-consistent global models with the inclusion of a statistical treatment of neutron and γ emission. The impact of our measurements on r-process nucleosynthesis was studied in a neutron star merger scenario. Our P 1n and P 2n have a direct impact on the odd-even staggering of the final abundance, improving the agreement between calculated and observed Solar System abundances. The odd isotope fraction of Ba in r-process-enhanced (r-II) stars is also better reproduced using our new data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MedNAS: Multiscale Training-Free Neural Architecture Search for Medical Image Analysis

Deep neural networks have demonstrated impressive results in medical image analysis, but designing suitable architectures for each specific task is expertise dependent and time consuming. Neural architecture search (NAS) offers an effective means of discovering architectures. It has been highly successful in numerous applications, particularly in natural image classification. Yet, medical images possess unique characteristics, such as small regions and a wide variety of lesion sizes, that differentiate them from natural images. Furthermore, most current NAS methods struggle with high computational costs, especially when dealing with high-resolution image datasets. In this article, we present a novel evolutionary NAS method called multiscale training-free neural architecture search (MSTF-NAS) to address these challenges. Specifically, to accommodate the broad range of lesion region sizes in disease diagnosis, we develop a new reduction cell search space that enables the search algorithm to explicitly identify the optimal scale combination for multiscale feature extraction. Further, to overcome the issue of high computational costs, we utilize training-free indicators as performance measures for candidate architectures, which allows us to search for the optimal architecture more efficiently. More specifically, by considering the capability and simplicity of various networks, we formulate a multiobjective optimization problem that involves two training-free indicators and model complexity for candidate architectures. Extensive experiments on a large medical image benchmark and a publicly available breast cancer detection dataset are conducted. The empirical results demonstrate that our MSTF-NAS outperforms both human-designed architectures and current state-of-the-art NAS algorithms on both datasets, indicating the effectiveness of our proposed method.

60 APPLIED LIFE SCIENCES↗

Feature Extraction for Pipeline Defects Inspection Based Upon Distributed Acoustic Fiber Optic Sensing Data

Fiber-optic distributed acoustic sensing (DAS) is becoming an increasingly important tool for real-time monitoring of energy and civil infrastructure structural health such as pipelines. We present a systematic theoretical study of the potential for DAS to be directly coupled with guided ultrasonic waves typically used in conventional acoustic non-destructive evaluation (NDE) methods for real-time pipeline health monitoring. We are referring to this innovative new NDE technique as ultrasonic guided wave and optical fiber sensor fusion. In the practical application of DAS coupled with guided ultrasonic waves, the structural design of (1) the specific guided waves excited, (2) the physical installation of the acoustic transducers and the fiber optic sensors, and (3) the functional performance specifications (gauge length, sensitivity, Etc.) of fiber optic DAS have an important influence on overall capabilities of the monitoring system. Meanwhile, physics-based analysis of acoustic waves is still a challenge due to the complex nature of the Lamb wave when it propagates, scatters, and disperses in the presence of structural defects. In this work, we simulate carbon steel pipes relevant for oil and gas pipeline applications with diameters of approximately 6-12” and wall thickness of 0.5” as the objects to be monitored. By establishing and implementing these capabilities, we seek to pursue an in-depth study on structural parameter optimization of DAS network, measurement range, and signal processing with an ultimate goal of increasing the sensitivity and efficacy of DAS to defect identification for various modes of corrosion expected in practice. To study the characteristics of scattered acoustic waves and performance of DAS for defect identification, we simulated the response of DAS for multiple pipe structures, defect types, and DAS sensor network configuration using finite element software Ansys, then the properties of signal response are extracted to construct defect-sensitive features. The raw data simulated, and the associated features extracted can ultimately be utilized as annotated training data to benchmark various designs for DAS applications, guided acoustic excitation sources, and learning model parameters to enhance early detection of potentially problematic defects.

Pipeline Defects Inspection, Fiber-optic sensors, ↗

TwoFold: Highly accurate structure and affinity prediction for protein-ligand complexes from sequences

We describe our development of ab initio protein-ligand binding pose prediction models based on transformers and binding affinity prediction models based on the neural tangent kernel (NTK). Folding both protein and ligand, the TwoFold models achieve efficient and quality predictions matching state-of-the-art implementations while additionally reconstructing protein structures. In conclusion, solving NTK models points to a new use case for highly optimized linear solver benchmarking codes on HPC.

60 APPLIED LIFE SCIENCES↗

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗

The Polar System Analysis Package (Ver. 1 Specifications)

The Polar System Analysis Package (POLeSTAr) is a modern data fusion tool being designed to benchmark,evaluate and compare the polar physics and biogeochemistry of Earth System Models (ESMs). It is being written in Python to make use of libraries including Xarray, Matplotlib, SciPy, NumPy and PyTorch within algorithms developed specifically for polar system modeling. Initial POLeSTAr algorithms will be adapted from the MATLAB Ridgepack package that was designed to analyze the Model for Prediction Across Scales (MPAS) sea ice component of the Energy Exascale Earth System Model (E3SM) and sea ice model CICE within the Community Earth System Model (CESM).

97 MATHEMATICS AND COMPUTING↗

Modeling Transverse Space Charge effects in IOTA with pyORBIT

The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.

43 PARTICLE ACCELERATORS↗

Stability Analysis of Power Systems with High Penetration of State-of-the-Art Inverter Technologies

With the increasing level of inverter-based resources (IBRs) in modern power systems, this paper presents a small-signal stability analysis for power systems comprising synchronous generators (SGs) and IBRs. Four types of inverter controls are considered: two grid-following (GFL) controls, with or without grid support functions; droop-based grid-forming (GFM) controls; and virtual oscillator control-based GFM. We also analyze the impact of STATCOM and synchronous condensers on system stability to assess their role in the energy mix transition. With the small-signal dynamic behavior of the major technologies modeled, this paper provides stringent stability assessments using the IEEE 39-bus benchmark system modified to simulate future power systems. The exhaustive test cases allow for (a) assessing the impacts of different types and controls of generation and supplementary grid assets, as well as system inertia and line impedance on grid stability, and (b) elucidating pathways for the stabilization of IBR-dominated power systems. The analysis also indicates that future power systems can be stabilized with only a fraction of the total generation as voltage sources without SGs or significant system inertia if they are well distributed. This study provides insights into future power system operations with a high level of IBRs that can also be used for planning and operation studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The boundary layer on compressor cascade blades

The flow field about an airfoil in cascade at a Reynolds number of 5 x 10 to the 5th power is described. Hot wire and laser anemometry are combined with flow visualization techniques in order to obtain detailed flow data (e.g., boundary layer profiles, points of separation, and the transition zone) on a cascade of relatively highly loaded blades. Benchmark data is provided for the evaluation of current and future predictive models, in this way aiding in the compressor design process.

Deutsch, S.↗

The boundary layer on compressor cascade blades

The flow field about an airfoil in a cascade at chord Reynolds number (R sub C) near 50,000. The program is experimental and combines laser Doppler anemometry (LDA) with flow visualization techniques in order to obtain detailed flow data (e.g., boundary layer profiles, points of separation and the transition zone) on a cascade of highly-loaded compressor blades. The information provided is to serve as benchmark data for the evaluation of current and future compressor cascade predictive models, in this way aiding in the compressor design process. The completed pressure surface mean velocity profiles, as well as two detailed near wake velocity profiles, all at a single incidence angle are provided.

Deutsch, S.↗

The boundary layer on compressor cascade blades

The purpose of NASA Research Grant NSG-3264 is to characterize the flowfield about an airfoil in a cascade at chord Reynolds number(R sub C)near 5 x 10 to the 5th power. The program is experimental and combines laser Doppler velocimeter (LDV) measurements with flow visualization techniques in order to obtain detailed flow data, e.g., boundary layer profiles, points of separation and the transition zone, on a cascade of highly-loaded compressor blades. The information provided by this study is to serve as benchmark data for the evaluation of current and future compressor cascade predictive models, in this way aiding in the compressor design process. Summarized is the research activity for the period 1 December 1985 through 1 June 1986. Progress made from 1 June 1979 through 1 December 1985 is presented. Detailed measurements have been completed at the initial cascade angle of 53 deg. (incidence angle 5 degrees). A three part study, based on that data, has been accepted as part of the 1986 Gas Turbine Conference and will be submitted for subsequent journal publication. Also presented are data for a second cascade angle of 45 deg (an incidence angle of 3 degrees).

Deutsch, S.↗

Flexible missile autopilot design studies with PC-MATLAB/386

Development of a responsive, high-bandwidth missile autopilot for airframes which have structural modes of unusually low frequency presents a challenging design task. Such systems are viable candidates for modern, state-space control design methods. The PC-MATLAB interactive software package provides an environment well-suited to the development of candidate linear control laws for flexible missile autopilots. The strengths of MATLAB include: (1) exceptionally high speed (MATLAB's version for 80386-based PC's offers benchmarks approaching minicomputer and mainframe performance); (2) ability to handle large design models of several hundred degrees of freedom, if necessary; and (3) broad extensibility through user-defined functions. To characterize MATLAB capabilities, a simplified design example is presented. This involves interactive definition of an observer-based state-space compensator for a flexible missile autopilot design task. MATLAB capabilities and limitations, in the context of this design task, are then summarized.

Ruth, Michael J.↗

A flux-limited treatment for the conductive evaporation of spherical interstellar gas clouds

In this work, we present and analyze a new analytic solution for the saturated (flux-limited) thermal evaporation of a spherical cloud. This work is distinguished from earlier analytic studies by allowing the thermal conductivity to change continuously from a diffusive to a saturated form, in a manner usually employed only in numerical calculations. This closed form solution will be of interest as a computational benchmark. Using our calculated temperature profiles and mass-loss rates, we model the thermal evaporation of such a cloud under typical interstellar medium (ISM) conditions, with some restrictions. We examine the ionization structure of the cloud-ISM interface and evaluate column densities of carbon, nitrogen, oxygen, neon, and silicon ions toward the cloud. In accord with other investigations, we find that ionization equilibrium is far from satisfied under the assumed conditions. Since the inclusion of saturation effects in the heat flux narrows the thermal interface relative to its classical structure, we also find that saturation effects tend to lower predicted column densities.

Dalton, William W.↗

Research of Stochastic Robustness: Results and conclusions

With stochastic robustness we are creating tools to design robust compensators for practical systems. During this year, the stochastic robustness research achieved the following results: refined the search tools needed for synthesis; successfully designed robust compensators for the American Controls Conference benchmark problem; and successfully designed robust compensators for a nonlinear hypersonic aircraft model with uncertainties in 28 parameters.

Marrison, Chris↗

Some Comments on the Status of Shell Theory at the End of the 20th Century: Complaints and Correctives

This review is divided into complaints and correctives. Complaints are directed at: sloppy refereeing and editing; authors who fail to read or acknowledge what others have done; the mis-naming or mis-crediting of results; the misunderstanding and misuse of the Kirchhoff hypothesis; inflated claims of accuracy based on overly-simplified benchmark problems; the failure to appreciate the inherent errors in various shell models; the failure to appreciate that the physical response and mathematical structure of shell theory are fundamentally different from 3-dimensional elasticity; and the irrelevance of Cosserat-type theories. Correctives include a simple, straight-forward derivation of a general nonlinear dynamic shell theory with the following features: (1) the equations of motion and kinematics (and those of thermodynamics, if desired) are exact consequences of their 3-dimensional counterparts; (2) there are no asymptotic or series expansions through the thickness; (3) all approximations (including the Kirchhoff Hypothesis) occur in the constitutive relations; (4) in static problems, there is a mixed form of the governing equations involving a mixed-energy density and exhibiting remnants of the well-known static-geometric duality of linear theory which is numerically robust because the limiting cases of nonlinear membrane theory and inextensional bending theory fall out naturally. (These latter two special cases are known to produce numerical nightmares unless treated with great care); and (5) all equations may be expressed in coordinate-free form (although, sometimes, a hybrid form is shown to be superior).

Simmonds, James G.↗

WIND Validation Cases: Computational Study of Thermally-perfect Gases

The ability of the WIND Navier-Stokes code to predict the physics of multi-species gases is investigated in support of future high-speed, high-temperature propulsion applications relevant to NASA's Space Transportation efforts. Three benchmark cases are investigated to evaluate the capability of the WIND chemistry model to accurately predict the aerodynamics of multi-species chemically non-reacting (frozen) gases. Case 1 represents turbulent mixing of sonic hydrogen and supersonic vitiated air. Case 2 consists of heated and unheated round supersonic jet exiting to ambient. Case 3 represents 2-D flow through a converging-diverging Mach 2 nozzle. For Case 1, the WIND results agree fairly well with experimental results and that significant mixing occurs downstream of the hydrogen injection point. For Case 2, the results show that the Wilke and Sutherland viscosity laws gave similar results, and the available SST turbulence model does not predict round supersonic nozzle flows accurately. For Case 3, results show that experimental, frozen, and 1-D gas results agree fairly well, and that frozen, homogeneous, multi-species gas calculations can be approximated by running in perfect gas mode while specifying the mixture gas constant and Ratio of Specific Heats.

DalBello, Teryn↗