Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Quality”

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,117 records · Page 62

Model Predictive Control of Discrete-Continuous Energy Systems via Generalized Disjunctive Programming

Generalized Disjunctive Programming (GDP) provides an alternative framework to model optimization problems with both discrete and continuous variables. The key idea behind GDP involves the use of logical disjunctions to represent discrete decisions in the continuous space, and logical propositions to denote algebraic constraints in the discrete space. Compared to traditional mixed-integer programming (MIP), the inherent logic structure in GDP yields tighter relaxations that are exploited by global branch and bound algorithms to improve solution quality. In this paper, we present a general GDP model for optimal control of hybrid systems that exhibit both discrete and continuous dynamics. Specifically, we use GDP to formulate a model predictive control (MPC) model for piecewise-affine systems with implicit switching logic. As an example, the GDP-based MPC approach is used as a supervisory control to improve energy efficiency in residential buildings with binary on/off, relay-based thermostats. A simulation study is used to demonstrate the validity of the proposed approach, and the improved solution quality compared to existing MIPbased control approaches.

Bhattacharya, Arnab↗

Dual redundant arm system operational quality measures and their applications - Dynamic measures

Dual-arm dynamic operation quality measures are presented which quantify the efficiency and capability of generating Cartesian accelerations by two cooperative arms based on the analysis of dual-arm dynamic interactions. Dual-arm dynamic manipulability is defined as the efficiency of generating Cartesian accelerations under the dynamic and kinematic interactions between individual arms and an object under manipulation. The analysis of dual-arm dynamic interactions is based on the so-called Cartesian space agent model of an arm, which represents an individual arm as a force source acting upon a point mass with the effective Cartesian space arm dynamics and an environment or an object under manipulation. The Cartesian space agent model of an arm makes it possible to derive the dynamic and kinematic constraints involved in the transport, assembly and grasping modes of dual-arm cooperation. A task-oriented operational quality measure, (TOQd) is defined by evaluating dual-arm dynamic manipulability in terms of given task requirements. TOQd is used in dual-arm joint configuration optimization. Simulation results are shown. A complete set of forward dynamic equations for a dual-arm system is derived, and dual-arm dynamic operational quality measures for various modes of dual-arm cooperation allowing sliding contacts are established.

Lee, Sukhan↗

TomoGAN: low-dose synchrotron x-ray tomography with generative adversarial networks: discussion

Synchrotron-based x-ray tomography is a noninvasive imaging technique that allows for reconstructing the internal structure of materials at high spatial resolutions from tens of micrometers to a few nanometers. In order to resolve sample features at smaller length scales, however, a higher radiation dose is required. Therefore, the limitation on the achievable resolution is set primarily by noise at these length scales. We present TomoGAN, a denoising technique based on generative adversarial networks, for improving the quality of reconstructed images for low-dose imaging conditions. We evaluate our approach in two photon-budget-limited experimental conditions: (1) sufficient number of low-dose projections (based on Nyquist sampling), and (2) insufficient or limited number of high-dose projections. In both cases, the angular sampling is assumed to be isotropic, and the photon budget throughout the experiment is fixed based on the maximum allowable radiation dose on the sample. Evaluation with both simulated and experimental datasets shows that our approach can significantly reduce noise in reconstructed images, improving the structural similarity score of simulation and experimental data from 0.18 to 0.9 and from 0.18 to 0.41, respectively. Furthermore, the quality of the reconstructed images with filtered back projection followed by our denoising approach exceeds that of reconstructions with the simultaneous iterative reconstruction technique, showing the computational superiority of our approach.

42 ENGINEERING↗

A Method of Moments Wide Band Adaptive Rational Interpolation Method for High-Quality Factor Resonant Cavities

A new adaptive rational interpolation method is proposed to obtain the wideband frequency response of a resonant cavity simulated with the method of moments (MoM). This interpolation method uses both the Loewner matrix to construct a rational expression for the solution vector of MoM’s matrix system and an error estimator generated by the solution vectors and their derivatives. This error estimator is implemented in the adaptive procedure to gain a minimum set of frequencies and solution vectors required in the interpolation. The resulting set of frequencies and solution vectors is applied to interpolate other system variables, such as shielding effectiveness and input impedance. Here numerical results of a slotted cylindrical cavity supporting high-quality factor resonances are presented, showing that the new rational interpolation method is accurate and efficient in interpolating the complicated resonant response of the solution vector functions.

42 ENGINEERING↗

Large aircraft handling qualities

The short period frequency requirements of MIL-F-8785C are applicable to the very large transport aircraft simulated. The large aircraft simulated meet the requirements of NASA CR-159236 for effective time delay and pitch transient peak ratio. However, the requirements of this reference for the effective rise time parameter are believed to be too convservative for very large transport aircraft. Data obtained to date as well as other data indicate that MIL-SPEC requirements for the parameter t(phi=30 degrees) are too conservative for very large transport aircraft. The results of the present study indicate that t(phi=30 degrees) of less than 6 sec should result in ""acceptable'' roll response characteristics, and t(phi=30 degrees) of less than 4.0 sec should result in ""satisfactory'' roll response.

Grantham, W. D.↗

CME Flux Rope and Shock Identifications and Locations: Comparison of White Light Data, Graduated Cylindrical Shell Model, and MHD Simulations

Coronal mass ejections (CMEs) are major transient phenomena in the solar corona that are observed with ground-based and spacecraft-based coronagraphs in white light or with in situ measurements by spacecraft. CMEs transport mass and momentum and often drive shocks. In order to derive the CME and shock trajectories with high precision, we apply the graduated cylindrical shell (GCS) model to fit a flux rope to the CME directed toward STEREO A after about 19:00 UT on 29 November 2013 and check the quality of the heliocentric distance-time evaluations by carrying out a three-dimensional magnetohydrodynamic (MHD) simulation of the same CME with the Block Adaptive Tree Solar-Wind Roe Upwind Scheme (BATS-R-US) code. Heliocentric distances of the CME and shock leading edges are determined from the simulated white light images and magnetic field strength data. We find very good agreement between the predicted and observed heliocentric distances, showing that the GCS model and the BATS-R-US simulation approach work very well and are consistent. In order to assess the validity of CME and shock identification criteria in coronagraph images, we also compute synthetic white light images of the CME and shock. We find that the outer edge of a cloud-like illuminated area in the observed and predicted images in fact coincides with the leading edge of the CME flux rope and that the outer edge of a faint illuminated band in front of the CME leading edge coincides with the CME-driven shock front.

Schmidt, J. M.↗

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗

Streamflow simulation in data-scarce basins using Bayesian and physics-informed machine learning models

Hydrologic predictions at rural watersheds are important but also challenging due to data shortage. Long short-term memory (LSTM) networks are a promising machine learning approach and have demonstrated good performance in streamflow predictions. However, due to its data-hungry nature, most LSTM applications focus on well-monitored catchments with abundant and high-quality observations. In this work, we investigate predictive capabilities of LSTM in poorly monitored watersheds with short observation records. To address three main challenges of LSTM applications in data-scarce locations, i.e., overfitting, uncertainty quantification (UQ), and out-of-distribution prediction, we evaluate different regularization techniques to prevent overfitting, apply a Bayesian LSTM for UQ, and introduce a physics-informed hybrid LSTM to enhance out-of-distribution prediction. Through case studies in two diverse sets of catchments with and without snow influence, we demonstrate that 1) when hydrologic variability in the prediction period is similar to the calibration period, LSTM models can reasonably predict daily streamflow with Nash–Sutcliffe efficiency above 0.8, even with only 2 years of calibration data; 2) when the hydrologic variability in the prediction and calibration periods is dramatically different, LSTM alone does not predict well, but the hybrid model can improve the out-of-distribution prediction with acceptable generalization accuracy; 3) L2 norm penalty and dropout can mitigate overfitting, and Bayesian and hybrid LSTM have no overfitting; and 4) Bayesian LSTM provides useful uncertainty information to improve prediction understanding and credibility. In conclusion, these insights have vital implications for streamflow simulation in watersheds where data quality and availability are a critical issue.

54 ENVIRONMENTAL SCIENCES↗

Space Shuttle flying qualities and flight control system assessment study, phase 2

A program of flying qualities experiments as part of the Orbiter Experiments Program (OEX) is defined. Phase 1, published as CR-170391, reviewed flying qualities criteria and shuttle data. The review of applicable experimental and shuttle data to further define the OEX plan is continued. An unconventional feature of this approach is the use of pilot strategy model identification to relate flight and simulator results. Instrumentation, software, and data analysis techniques for pilot model measurements are examined. The relationship between shuttle characteristics and superaugmented aircraft is established. STS flights 1 through 4 are reviewed from the point of view of flying qualities. A preliminary plan for a coordinated program of inflight and simulator research is presented.

Myers, T. T.↗

Optimization of a High Bunch Charge ERL Injection Merger for PERLE

Delivery of high charge electron bunches into the main loop of an ERL (energy recovery linac) while preserving the emittance is challenging. This is because at the typical injection momentum, space charge forces still have a significant effect on the beam dynamics. In this work we consider the design of the merger for PERLE, an ERL test facility to be based at IJCLab in France. Previous simulations have shown that the baseline DC gun based injector can achieve the required emittance at the booster linac exit. The quality of the 500 pC bunches must then be preserved with space charge through the merger at total beam energy of 7 MeV keeping the emittance below 6 mm mrad. The beam dynamics in the merger were simulated using the code OPAL and optimised using a genetic algorithm. Three possible merger schemes were investigated. The goal of the optimisation was to minimise the emittance growth while also achieving the required Twiss parameters to match onto the spreader at the main linac exit. A three dipole solution is then examined in more detail.

Hounsell, B.↗

Model-Based Reconstruction for Multi-Frequency Collimated Beam Ultrasound Systems

Collimated beam ultrasound systems are a technology for imaging inside multi-layered structures such as geothermal wells. These systems work by using a collimated narrow-band ultrasound transmitter that can penetrate through multiple layers of heterogeneous material. A series of measurements can then be made at multiple transmit frequencies. However, commonly used reconstruction algorithms such as Synthetic Aperture Focusing Technique (SAFT) tend to produce poor quality reconstructions for these systems both because they do not model collimated beam systems and they do not jointly reconstruct the multiple frequencies. Here, in this article, we propose a multi-frequency ultrasound model-based iterative reconstruction (UMBIR) algorithm designed for multi-frequency collimated beam ultrasound systems. The combined system targets reflective imaging of heterogeneous, multi-layered structures. For each transmitted frequency band, we introduce a physics-based forward model to accurately account for the propagation of the collimated narrow-band ultrasonic beam through the multi-layered media. We then show how the joint multi-frequency UMBIR reconstruction can be computed by modeling the direct arrival signals, detector noise, and incorporating a spatially varying image prior. Results using both simulated and experimental data indicate that multi-frequency UMBIR reconstruction yields much higher reconstruction quality than either single frequency UMBIR or SAFT.

47 OTHER INSTRUMENTATION↗

Flash Graphene Morphologies

Flash Joule heating (FJH) can convert almost any carbon-based precursor into bulk quantities of graphene. This work explores the morphologies and properties of flash graphene (FG) generated from carbon black. It is shown that FG is partially comprised of sheets of turbostratic FG (tFG) that have a rotational mismatch between neighboring layers. The remainder of the FG is wrinkled graphene sheets that resemble non-graphitizing carbon. To generate high quality tFG sheets, a FJH duration of 30 – 100 ms is employed. Beyond 100 ms, the turbostratic sheets have time to AB-stack and form bulk graphite. Atomistic simulations reveal that generic thermal annealing yields predominantly wrinkled graphene which displays minimal to no alignment of graphitic planes, as opposed to the high-quality tFG that might be formed under the direct influence of current conducted through the material. As a result, the tFG was easily exfoliated via shear, hence the FJH process has the potential for bulk production of tFG without the need for pre-exfoliation using chemicals or high energy mechanical shear.

01 COAL, LIGNITE, AND PEAT↗

DOE Hanford Tank Side Cesium Removal System (TSCR) Update: Fabrication and Factory Acceptance Testing - 20300

The Tank Side Cesium Removal (TSCR) System is a modular, quickly-deployable, and passively safe system for efficiently removing solids and cesium from alkaline liquid waste raffinate created by fuel processing on the Hanford site. This paper provides an overview of the project with particular emphasis on completion of fabrication and factory acceptance testing (FAT). All TSCR System components are inside three enclosures known as the Process Enclosure, Ancillary Enclosure, and Control Enclosure. The Process Enclosure is seismic design category II (SDC-II) and built to withstand the rigors of transportation and lifting, high winds, volcanic ash fall, and heavy precipitation (snow) over a temperature range of -25 to 115 deg. F. The Process Enclosure and Control Enclosure were fabricated in AVANTech's Columbia, SC facility, while the Ancillary Enclosure was fabricated in AVANTech's Richland, WA facility. The design and preliminary documented safety analysis were completed in December 2019, with all work being performed in accordance with ASME NQA-1 quality assurance requirements. The Process Enclosure, Ion Exchange Columns (IXCs), and several other TSCR components have a nuclear safety function and are designated as safety significant structures, systems, and components (SS-SSCs). The NQA-1 Commercial Grade Dedication (CGD) process was used to identify critical characteristics in order to accept materials and items used for fabrication of several safety significant TSCR components. This paper also addresses unique aspects and challenges of implementing CGD in DOE facilities and systems. After fabrication and preliminary testing, the Process and Control Enclosures were shipped to AVANTech's Richland facility where the complete TSCR System was setup and connected in a represented layout and configuration that will be used on the Hanford site. Test tanks and pumps were put in place and connected to the TSCR Process Enclosure in a manner that mimicked interfacing Double Shell Tanks in the 241-AP Tank Farm. The TSCR System was then tested/operated under conditions equivalent to those planned for onsite operation. Testing was separated into five separate phases as defined below: - Phase 1: Confirmation of TSCR System readiness for the FAT; - Phase 2: Verification of TSCR System operation under standard conditions; - Phase 3: Demonstration of TSCR System operation under non-standard conditions; - Phase 4: Demonstration of ion exchange media removal from an IXC; and - Phase 5: Completion of TSCR System lift plans and site mobilization. The Hanford Tank Operations Contractor, Washington River Protection Solutions (WRPS), participated in all phases of testing. Testing at AVANTech's Richland Facility allowed WRPS personnel to familiarize themselves with TSCR operations in a radiologically clean facility prior to onsite mobilization. Engineering and operations personnel gained valuable experience that allowed them to better refine onsite operational protocols and procedures - including handling and movement of IXCs with a forklift. This paper describes the technologies used by the TSCR System as well as how these technologies are deployed and operated to advance the Hanford mission of safely, efficiently, and effectively treating tank waste. It also provides detailed information on fabrication and FAT, including operational data from simulated tank waste runs. A project schedule for WRPS' planned construction acceptance testing, operational acceptance testing, and actual tank waste treatment is provided. This presentation should be of great interest to parties responsible for the design, testing, and mobilization of modular NQA-1 technologies for the treatment of raffinate and similar liquid wastes containing high concentrations of radioactive cesium. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Optimizing Error-Bounded Lossy Compression for Scientific Data by Dynamic Spline Interpolation

Today's scientific simulations are producing vast volumes of data that cannot be stored and transferred efficiently because of limited storage capacity, parallel I/O bandwidth, and network bandwidth. The situation is getting worse over time because of the ever-increasing gap between relatively slow data transfer speed and fast-growing computation power in modern supercomputers. Error-bounded lossy compression is becoming one of the most critical techniques for resolving the big scientific data issue, in that it can significantly reduce the scientific data volume while guaranteeing that the reconstructed data is valid for users because of its compression-error-bounding feature. In this paper, we present a novel error-bounded lossy compressor based on a state-of-the-art prediction-based compression framework. Our solution exhibits substantially better compression quality than all of the existing error-bounded lossy compressors, with comparable compression speed. Specifically, our contribution is threefold. (1) We provide an in-depth analysis of why the best-existing prediction-based lossy compressor can only minimally improve the compression quality. (2) We propose a dynamic spline interpolation approach with a series of optimization strategies that can significantly improve the data prediction accuracy, substantially improving the compression quality in turn. (3) We perform a thorough evaluation using six real-world scientific simulation datasets across different science domains to evaluate our solution vs. all other related works. Experiments show that the compression ratio of our solution is higher than that of the second-best lossy compressor by 20%similar to 460% with the same error bound in most of the cases.

Zhao, Kai↗

Effect of surface properties and polymer chain length on polymer adsorption in solution

In polymer nanoparticle composites (PNCs) with attractive interactions between nanoparticles (NPs) and polymers, a bound layer of the polymer forms on the NP surface, with significant effects on the macroscopic properties of the PNCs. The adsorption and wetting behaviors of polymer solutions in the presence of a solid surface are critical to the fabrication process of PNCs. In this study, we use both classical density functional theory (cDFT) and molecular dynamics (MD) simulations to study dilute and semi-dilute solutions of short polymer chains near a solid surface. Using cDFT, we calculate the equilibrium properties of polymer solutions near a flat surface while varying the solvent quality, surface–fluid interactions, and the polymer chain lengths to investigate their effects on the polymer adsorption and wetting transitions. Using MD simulations, we simulate polymer solutions near solid surfaces with three different curvatures (a flat surface and NPs with two radii) to study the static conformation of the polymer bound layer near the surface and the dynamic chain adsorption process. Additionally, we find that the bulk polymer concentration at which the wetting transition in the poor solvent system occurs is not affected by the difference in surface–fluid interactions; however, a threshold value of surface–fluid interaction is needed to observe the wetting transition. We also find that with good solvent, increasing the chain length or the difference in the surface–polymer interaction relative to the surface–solvent interaction increases the surface coverage of polymer segments and independent chains for all surface curvatures. Finally, we demonstrate that the polymer segmental adsorption times are heavily influenced only by the surface–fluid interactions, although polymers desorb more quickly from highly curved surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Propagation of Input Uncertainties in Numerical Simulations of Laser Powder Bed Fusion

Laser powder bed fusion has the potential of redefining state-of-the-art processing and production methods, but defect formation and inconsistent build quality have limited the implementation of this process. Numerical models are widely used to study this process and predict the formation of these defects. Presently, the uncertainties of model input parameters and thermophysical properties used by these numerical simulations have not been investigated. In the present study, the uncertainty in these input parameters and material properties are quantified for laser powder bed fusion, with and without a simulated powder bed, to determine their influence on the predictive accuracy of an experimentally validated numerical model. Accounting for all possible sources of uncertainty quickly becomes computationally expensive on account of the curse of dimensionality. Uncertainty in laser absorption, solid, and liquid specific heat of the metal were found to have the largest effect on model prediction reliability with or without the use of a powder bed. Results also illustrate that accounting for these three uncertain parameters still captures the majority of model prediction uncertainty. Furthermore, the methodology of this study may be used to understand the uncertainty in as-built microstructure through propagation to microstructure prediction models, or applied under processing conditions where high Péclet numbers are observed and the thermal convection and fluid flow within the molten pool are substantial.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The N/Rev phenomenon in simulating a blade-element rotor system

When a simulation model produces frequencies that are beyond the bandwidth of a discrete implementation, anomalous frequencies appear within the bandwidth. Such is the case with blade element models of rotor systems, which are used in the real time, man in the loop simulation environment. Steady state, high frequency harmonics generated by these models, whether aliased or not, obscure piloted helicopter simulation responses. Since these harmonics are attenuated in actual rotorcraft (e.g., because of structural damping), a faithful environment representation for handling qualities purposes may be created from the original model by using certain filtering techniques, as outlined here. These include harmonic consideration, conventional filtering, and decontamination. The process of decontamination is of special interest because frequencies of importance to simulation operation are not attenuated, whereas superimposed aliased harmonics are.

Mcfarland, R. E.↗

The N/rev phenomenon in simulating a blade-element rotor system

When a simulation model produces frequencies that are beyond the bandwidth of a discrete implementation, anomalous frequencies appear within the bandwidth. Such is the case with blade element models of rotor systems, which are used in the real time, man in the loop simulation environment. Steady state, high frequency harmonics generated by these models, whether aliased or not, obscure piloted helicopter simulation responses. Since these harmonics are attenuated in actual rotorcraft (e.g., because of structural damping), a faithful environment representation for handling qualities purposes may be created from the original model by using certain filtering techniques, as outlined here. These include harmonic consideration, conventional filtering, and decontamination. The process of decontamination is of special interest because frequencies of importance to simulation operation are not attenuated, whereas superimposed aliased harmonics are.

Mcfarland, R. E.↗