Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model checking”

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 289 records · Page 16

On implementing self-checking microprocessors

A simple and general model of the interfaces and check circuits used for comparing and detecting faults in a pair of 16-bit processors is described, and problems encountered in the application of TI 9900 processors are discussed. The greatest incompatibility is found to lie between the rollback structures of the CPUs and the interface and check logic (ICL) model. The ICL model generates a reset when an error is detected, and a rollback is expected to occur when it is released. The TI 9900 requires a reset of minimum duration, and after release goes through an initialization cycle, obtains rollback parameters from fixed memory locations, and executes the rollback, consistent with the ICL. The ICL is relatively simple, having a complexity equivalent to fewer than 1000 gates.

Rennels, D. A.↗

Thermal System Verification and Model Validation for NASA's Cryogenic Passively Cooled James Webb Space Telescope

A thorough and unique thermal verification and model validation plan has been developed for NASA s James Webb Space Telescope. The JWST observatory consists of a large deployed aperture optical telescope passively cooled to below 50 Kelvin along with a suite of several instruments passively and actively cooled to below 37 Kelvin and 7 Kelvin, respectively. Passive cooling to these extremely low temperatures is made feasible by the use of a large deployed high efficiency sunshield and an orbit location at the L2 Lagrange point. Another enabling feature is the scale or size of the observatory that allows for large radiator sizes that are compatible with the expected power dissipation of the instruments and large format Mercury Cadmium Telluride (HgCdTe) detector arrays. This passive cooling concept is simple, reliable, and mission enabling when compared to the alternatives of mechanical coolers and stored cryogens. However, these same large scale observatory features, which make passive cooling viable, also prevent the typical flight configuration fully-deployed thermal balance test that is the keystone to most space missions thermal verification plan. JWST is simply too large in its deployed configuration to be properly thermal balance tested in the facilities that currently exist. This reality, when combined with a mission thermal concept with little to no flight heritage, has necessitated the need for a unique and alternative approach to thermal system verification and model validation. This paper describes the thermal verification and model validation plan that has been developed for JWST. The plan relies on judicious use of cryogenic and thermal design margin, a completely independent thermal modeling cross check utilizing different analysis teams and software packages, and finally, a comprehensive set of thermal tests that occur at different levels of JWST assembly. After a brief description of the JWST mission and thermal architecture, a detailed description of the three aspects of the thermal verification and model validation plan is presented.

Cleveland, Paul E.↗

Kate's Model Verification Tools

Kennedy Space Center's Knowledge-based Autonomous Test Engineer (KATE) is capable of monitoring electromechanical systems, diagnosing their errors, and even repairing them when they crash. A survey of KATE's developer/modelers revealed that they were already using a sophisticated set of productivity enhancing tools. They did request five more, however, and those make up the body of the information presented here: (1) a transfer function code fitter; (2) a FORTRAN-Lisp translator; (3) three existing structural consistency checkers to aid in syntax checking their modeled device frames; (4) an automated procedure for calibrating knowledge base admittances to protect KATE's hardware mockups from inadvertent hand valve twiddling; and (5) three alternatives for the 'pseudo object', a programming patch that currently apprises KATE's modeling devices of their operational environments.

Morgan, Steve↗

Bridging the Gap Between Requirements and Model Analysis : Evaluation on Ten Cyber-Physical Challenge Problems

Formal verfication and simulation are powerful tools to validate requirements against complex systems. [Problem] Requirements are developed in early stages of the software lifecycle and are typically written in ambiguous natural language. There is a gap between such requirements and formal notations that can be used by verification tools, and lack of support for proper association of requirements with software artifacts for verification. [Principal idea] We propose to write requirements in an intuitive, structured natural language with formal semantics, and to support formalization and model/code verification as a smooth, well-integrated process. [Contribution] We have developed an end-to-end, open source requirements analysis framework that checks Simulink models against requirements written in structured natural language. Our framework is built in the Formal Requirements Elicitation Tool (fret); we use fret's requirements language named fretish, and formalization of fretish requirements in temporal logics. Our proposed framework contributes the following features: 1) automatic extraction of Simulink model information and association of fretish requirements with target model signals and components; 2) translation of temporal logic formulas into synchronous dataflow cocospec specifications as well as Simulink monitors, to be used by verification tools; we establish correctness of our translation through extensive automated testing; 3) interpretation of counterexamples produced by verification tools back at requirements level. These features support a tight integration and feedback loop between high level requirements and their analysis. We demonstrate our approach on a major case study: the Ten Lockheed Martin Cyber-Physical, aerospace-inspired challenge problems.

Mavridou, Anastasia↗

Verification of RESRAD-BUILD Code Version 4

This report documents the verification of the RESRAD-BUILD code, Version 4.0, which was released on December 22, 2022. Two earlier reports verifying Versions 3.0 and 3.1, respectively, were published in 2001 (Kamboj, et al. 2001) and 2003 (Tetra Tech NUS 2003). Version 4.0 of the RESRAD-BUILD code has many new features and modeling enhancements over the earlier versions, including the previously released Version 3.5. Chapter 2 of this report focuses on verifying the external dose and risk modeling for point, line, area, and volume sources, as well as for floor deposition. Besides verification, the external radiation doses calculated by RESRAD-BUILD were also benchmarked with those calculated by the MCNP code (Briemeister 1993). Section J.3 of the RESRAD-BUILD User’s Manual Vol. 1 (Yu et al. 2022) documents the results of that benchmarking effort. Chapter 3 of this report focuses on verifying the ventilation modeling, from checking the remaining source inventory, releases of radionuclides to the air, air concentrations and deposited floor concentrations over time, to the radiation dose and risk associated with inhalation, ingestion, and air submersion, with and without vacuuming. The verification efforts involve designing spreadsheets to perform calculations the same as or like those performed by the RESRAD-BUILD code and then comparing the spreadsheet results with those produced by the code. When the results agree or the differences are within acceptable range, the accuracy of model implementation in the code is verified. In addition to model implementation, the implementation of key functions and features that facilitate the modeling or the use of the code were also verified during the release testing of the code. Appendix A presents the test cases developed for these verification testing, and Appendix B presents the testing results that verify proper implementation of key functions and features.

54 ENVIRONMENTAL SCIENCES↗

Recent Improvements and Verification of A Full Body Model in Opensim

BACKGROUND: The dynamic feasibility [1,2] criterion, that the subject’s Center of Pressure (COP) be located within the Base of Support (BOS) which outlines the feet, has aided in assessing the stability of human motion recorded on earth while performing the recorded tasks in lunar gravity or as countermeasures exercises on a vibration isolation and stabilization system in microgravity. The convex hull of the BOS on the platform under the subject’s feet was estimated using virtual markers on the feet of the scaled subject model. The COP was calculated using the ground reaction forces and moments determined from motion capture data with biomechanical modeling tools [3]. Occasionally, large-amplitude oscillatory spikes or “artifacts” were observed in the subject’s linear and angular momentum derivatives, affecting some COP data derived from motion capture. The purpose of this investigation was to assess and improve the accuracy of model scaling and BOS estimation as well as to determine the efficacy of model adjustments in mitigating artifacts influencing motion capture-derived ground reaction force and COP results. METHODS: To aid evaluation of proposed process and model updates, motion capture data were collected for two subjects during unit test and range of motion trials, lunar tasks, and countermeasures exercise motions. Markers were added to the full body Plug-in Gait marker set [4] during data collection. New markers were placed on the front, back, sides, and top of the head to improve scaling using distances between marker pairs. Medial elbow markers were added to stabilize the upper arm during OpenSim Inverse Kinematics (IK) [5]. Finally, markers were added on the outer edge of the heels and on the outside edges of the first and last toes on each foot. These additional foot markers were made available to test new automated foot scaling techniques and to calculate the error between the subject’s estimated and recorded BOS. The modified unscaled OpenSim Full Body Rajagopal Model [6,7] was adjusted using some previously investigated techniques [8] to mitigate rapid shifts in joint angles occurring during IK, as these were found to cause the spike artifacts observed in subsequent stages of analysis. Since OpenSim models use Euler angles and rotation sequences, the arm axes of rotation were adjusted, and the pelvis order of rotation was changed to minimize the likelihood of encountering “gimbal lock” during common human motion. The model clavicle, arm, elbow, wrist, pelvis, and ankle angle limits were adjusted to better accommodate the full human range of motion seen in exercise and lunar data. The shoulder joint center was calculated using a “pivoting” algorithm [9], and both shoulder joint center and upper arm markers were included during IK to provide additional shoulder stability on a case-by-case basis. The quality of IK results was assessed by three criteria: minimizing error between recorded motion capture markers and model markers, checking for reasonable rates of change in joint angles between fames (i.e., no IK artifacts), and ensuring the absence of spikes in the inertial forces and angular momentum derivatives calculated using a custom OpenSim plugin [10]. RESULTS: The additional markers placed on the subject during data collection allowed the head and feet to be scaled more accurately using distances between new marker pairs. Scaling with BOS markers placed on the subject and removing the limit on subtalar angle resulted in more accurate BOS determination. Unrealistically large changes in joint angles between frames could be reduced by including clavicle, sternum, and medial elbow markers during IK. In cases with large arm ranges of motion, results could be further improved by running IK using medial elbow and virtual shoulder joint center markers. Model adjustments significantly improved the IK results affecting COP calculation and increased the accuracy of BOS estimation.

C A Bell↗

Connection between plasma response and resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression in DIII-D [Connection between plasma response and RMP ELM suppression in DIII-D]

Calculations of the plasma response to applied non-axisymmetric fields in several DIII-D discharges show that predicted displacements depend strongly on the edge current density. This result is found using both a linear two-fluid-MHD model (M3D-C1) and a nonlinear ideal-MHD model (VMEC). Furthermore, it is observed that the probability of a discharge being edge localized mode (ELM)-suppressed is most closely related to the edge current density, as opposed to the pressure gradient. It is found that discharges with a stronger kink response are closer to the peeling–ballooning stability limit in ELITE simulations and eventually cross into the unstable region, causing ELMs to reappear. Thus for effective ELM suppression, the RMP has to prevent the plasma from generating a large kink response, associated with ELM instability. Experimental observations are in agreement with the finding; discharges which have a strong kink response in the MHD simulations show ELMs or ELM mitigation during the RMP phase of the experiment, while discharges with a small kink response in the MHD simulations are fully ELM suppressed in the experiment by the applied resonant magnetic perturbation. The results are cross-checked against modeled 3D ideal MHD equilibria using the VMEC code. The procedure of constructing optimal 3D equilibria for diverted H-mode discharges using VMEC is presented. As a result, kink displacements in VMEC are found to scale with the edge current density, similar to M3D-C1, but the displacements are smaller. A direct correlation in the flux surface displacements to the bootstrap current is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optoelectronic Properties of C60 and C70 Fullerene Derivatives: Designing and Evaluating Novel Candidates for Efficient P3HT Polymer Solar Cells

Ten novel fullerene-derivatives (FDs) of C60 and C70 had been designed as acceptor for polymer solar cell (PSC) by employing the quantitative structure-property relationship (QSPR) model, which was developed strategically with a reasonably big pool of experimental power conversion efficiency (PCE) data. The QSPR model was checked and validated with stringent parameter and reliability of predicted PCE values of all designed FDs. They were assessed by the applicability domain (AD) and process randomization test. The predicted PCE of FDs range from 7.96 to 23.01. The obtained encouraging results led us to the additional theoretical analysis of the energetics and UV-Vis spectra of isolated dyes employing Density functional theory (DFT) and Time-dependent-DFT (TD-DFT) calculations using PBE/6-31G(d,p) and CAM-B3LYP/6-311G(d,p) level calculations, respectively. The FD4 is the best C60-derivatives candidates for PSCs as it has the lowest exciton binding energy, up-shifted lowest unoccupied molecular orbital (LUMO) energy level to increase open-circuit voltage (VOC) and strong absorption in the UV region. In case of C70-derivatives, FD7 is potential candidate for future PSCs due to its strong absorption in UV-Vis region and lower exciton binding energy with higher VOC. Our optoelectronic results strongly support the developed QSPR model equation. Analyzing QSPR model and optoelectronic parameters, we concluded that the FD1, FD2, FD4, and FD10 are the most potential candidates for acceptor fragment of fullerene-based PSC. The outcomes of tactical molecular design followed by the investigation of optoelectronic features are suggested to be employed as a significant resource for the synthesis of FDs as an acceptor of PSCs.

36 MATERIALS SCIENCE↗

Combining static analysis and state transition graphs for verification of event-condition-action systems in the RERS 2012 and 2013 challenges

In this work, we present a combination of approaches for the verification of event-condition-action (ECA) systems. The analyzed ECA systems range from structurally simple to structurally complex systems. We address the verification of reachability properties and behavioral properties. Reachability properties are represented by assertions in the program and we determine statically whether an assertion holds for all execution paths. Behavioral properties are represented as linear temporal logic formulas specifying the input/output behavior of the program. Our approach assumes a finite state space. We compare a symbolic analysis with an exhaustive state space exploration and discuss the trade-offs between the approaches in terms of the number of computed states and run-time behavior. All variants compute a state transition graph which can also be passed to an LTL verifier. The variants have a different impact on the number of computed states in the state transition graph which in turn impacts the run-time and memory consumption of subsequent phases. Finally, we evaluate the different analysis variants with the RERS benchmarks.

97 MATHEMATICS AND COMPUTING↗

Extended convex hull-based distributed optimal energy flow of integrated electricity-gas systems

Integrated electricity and gas systems are constructed to facilitate the gas-fired generation, and the distributed operation of these integrated systems have received much attention due to the increased emphasis on data security and privacy between different agencies. This paper proposes an extended convex hull based method to address optimal energy flow problems for the integrated electricity and gas systems in a distributed manner. First, a multi-block electricity-gas system model is constructed by dividing the whole system into N blocks considering both physical and regional differences. This multi-block model is then convexified by replacing the nonconvex gas transmission equation with the extended convex hull-based constraints. The Jacobi-Proximal alternating direction method of multipliers algorithm is adopted to solve the convexified model and minimize its operation cost. Finally, the feasibility of the optimal solution for the convexified model is checked, and a sufficient condition is developed. If the sufficient condition is satisfied, the optimal solution for the original nonconvex problem can be recovered from that for the convexified problem. Simulation results demonstrate that the proposed method is tractable and effective in obtaining feasible optimal solutions for multi-block optimal energy flow problems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

New approach for precise computation of Lyman-α forest power spectrum with hydrodynamical simulations

Current experiments are providing measurements of the flux power spectrum from the Lyman-α forests observed in quasar spectra with unprecedented accuracy. Their interpretation in terms of cosmological constraints requires specific simulations of at least equivalent precision. In this paper, we present a suite of cosmological N-body simulations with cold dark matter and baryons, specifically aiming at modeling the low-density regions of the inter-galactic medium as probed by the Lyman-α forests at high redshift. The simulations were run using the GADGET-3 code and were designed to match the requirements imposed by the quality of the current SDSS-III/BOSS or forthcoming SDSS-IV/eBOSS data. They are made using either 2 × 768{sup 3} ≅ 1 billion or 2 × 192{sup 3} ≅ 14 million particles, spanning volumes ranging from (25 Mpc h{sup −1}){sup 3} for high-resolution simulations to (100 Mpc h{sup −1}){sup 3} for large-volume ones. Using a splicing technique, the resolution is further enhanced to reach the equivalent of simulations with 2 × 3072{sup 3} ≅ 58 billion particles in a (100 Mpc h{sup −1}){sup 3} box size, i.e. a mean mass per gas particle of 1.2 × 10{sup 5}M{sub ⊙} h{sup −1}. We show that the resulting power spectrum is accurate at the 2% level over the full range from a few Mpc to several tens of Mpc. We explore the effect on the one-dimensional transmitted-flux power spectrum of four cosmological parameters (n{sub s}, σ{sub 8}, Ω{sub m} and H{sub 0}) and two astrophysical parameters (T{sub 0} and γ) that are related to the heating rate of the intergalactic medium. By varying the input parameters around a central model chosen to be in agreement with the latest Planck results, we built a grid of simulations that allows the study of the impact on the flux power spectrum of these six relevant parameters. We improve upon previous studies by not only measuring the effect of each parameter individually, but also probing the impact of the simultaneous variation of each pair of parameters. We thus provide a full second-order expansion, including cross-terms, around our central model. We check the validity of the second-order expansion with independent simulations obtained either with different cosmological parameters or different seeds. Finally, a comparison to the one-dimensional Lyman-α forest power spectrum obtained with BOSS by [1] shows an excellent agreement.

79 ASTRONOMY AND ASTROPHYSICS↗

A 14 h {sup −3} Gpc{sup 3} study of cosmic homogeneity using BOSS DR12 quasar sample

The BOSS quasar sample is used to study cosmic homogeneity with a 3D survey in the redshift range 2.2 < z < 2.8. We measure the count-in-sphere, N (< r ), i.e. the average number of objects around a given object, and its logarithmic derivative, the fractal correlation dimension, D {sub 2}( r ). For a homogeneous distribution N (< r ) ∝ r {sup 3} and D {sub 2}( r ) = 3. Due to the uncertainty on tracer density evolution, 3D surveys can only probe homogeneity up to a redshift dependence, i.e. they probe so-called ''spatial isotropy'. Our data demonstrate spatial isotropy of the quasar distribution in the redshift range 2.2 < z < 2.8 in a model-independent way, independent of any FLRW fiducial cosmology, resulting in 3 − ( D {sub 2}) < 1.7 × 10{sup −3} (2 σ) over the range 250 < r < 1200 h {sup −1} Mpc for the quasar distribution. If we assume that quasars do not have a bias much less than unity, this implies spatial isotropy of the matter distribution on large scales. Then, combining with the Copernican principle, we finally get homogeneity of the matter distribution on large scales. Alternatively, using a flat ΛCDM fiducial cosmology with CMB-derived parameters, and measuring the quasar bias relative to this ΛCDM model, our data provide a consistency check of the model, in terms of how homogeneous the Universe is on different scales. D {sub 2}( r ) is found to be compatible with our ΛCDM model on the whole 10 < r < 1200 h {sup −1} Mpc range. For the matter distribution we obtain 3 − ( D {sub 2}) < 5 × 10{sup −5} (2 σ) over the range 250 < r < 1200 h {sup −1} Mpc, consistent with homogeneity on large scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Program Verification for Extreme-Scale Applications (Final Scientific/Technical Report)

The project seeks to develop tools and techniques to help software developers verify the correctness and accuracy of their code. The target domain is scientific software of the kind widely used and developed in the Department of Energy research community, with a particular focus on "extreme scale" programs - those that are expected to involve possibly millions of parallel threads of execution. The report covers the University of Delaware contribution to the collaborative project. The project had a number of successful outcomes, especially regarding the development and extension of the CIVL software verification framework. CIVL is a verification tool for C or Fortran programs that use MPI, OpenMP, CUDA, and/or Pthreads for parallelization. CIVL went through vast improvements and extensions, and was successfully applied to a number of challenging codes. It found a subtle bug in the Devito PDE framework. It was able to verify the functional correctness of a conjugate gradient solver using a novel probabilistic technique with vanishingly small chance of error. CIVL was used very successfully in verification competitions, and the CIVL solutions to the competition challenges were presented and published. CIVL was also used successfully by other researchers on a computational chemistry kernel.

97 MATHEMATICS AND COMPUTING↗

Gamma and Neutron Measurement and Modeling of Irradiated TRISO Fuel

Given the unique characteristics of the PBR fuel cycle, both gamma and neutron measurements are expected to play important roles in performing and maintaining nuclear material control and accounting for spent pebbles to safeguard the fuel cycle. Given the lack of irradiated pebbles in the US, a variety of irradiated TRISO fuel samples with wide ranges of burnups and cooling times available at ORNL were used in this work. A large number of gamma and neutron measurements have been performed on these samples to collect data to test the various detectors and to benchmark the computer models to simulate the depletion and decay of the fuel and the measurements themselves. Two neutron detectors, including a custom-made detector and the Very High-Performance Neutron Multiplicity Counting, were used to measure the neutrons emitted by these TRISO samples. Three gamma spectrometry detectors, including an HPGe and the M400 CZT detector, were used to measure gamma-ray emissions from these samples. The M400 was recently adopted by the IAEA for fresh uranium measurements, but it was tested for spent fuel measurements prior to this project. Detailed MCNP models were developed to simulate these neutron and gamma measurements. Some GADRAS models were also developed to cross check the MCNP models for the gamma measurements. It was found challenging to perform neutron measurements in the hot cell due to the high background counts. Close agreements were observed between the simulated and measured neutron count rates in both detectors’ measurements of californium calibration sources. Both the HPGe and M400 detectors were able to measure the 604 and 662 keV peaks from these samples, which are the two most important peaks used to infer fuel burnup. Although the M400 detector did not have nearly good energy resolution and did not detect some of the minor peaks as the HPGe detector, it was found to be capable of handling significantly higher dose rates than HPGe. Given the complexities in the TRISO samples (e.g., different samples sizes) and uncertainties in the alignments between the detector and the TRISO fuel inside the containers, large scatters were found between the peak area rates and the samples’ burnups. However, the 604/662 peak ratios were found to trend well with the samples’ burnups among most samples in both measured and simulated results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Force testing manual for the Langley 20-inch Mach 6 tunnel

Data reduction and procedures for conducting force tests in a 20 inch Mach 6 tunnel are described. A discussion of pretest and testing phases are included. Items that are to be checked during model design and construction are outlined as well as safety requirements, starting loads tests, instructions for data acquisition and model installation. Measurement of balance and model misalignment and instructions for calibrating the angle of attack screen are covered. Procedures for making reference pressure, attitude tare, and data runs are included. The 20 inch tunnel force program is examined, and a description of data recording system input and load contrast sheets is given. An appendix presents a description, operating characteristics, and Mach number calibration of the tunnel, as well as tunnel characteristics.

Keyes, J. W.↗

Detection of SO2 in the UV spectrum of Venus

The broad absorption feature below 3300 A in the Venus UV spectrum is identified as primarily due to SO2 absorption based on new higher resolution spectra of the 3000-3400 A region showing broad (10 A), unresolved absorptions in the regions at all SO2 band origins between 3000 and 3300 A. SO2 mixing ratios vary from 5 x 10 to the -7th down to an upper limit of 2 x 10 to the -8th at a phase angle of 138 deg. Previous observational determinations of the SO2 mixing ratio were biased toward large phase angles, and consequently did not detect any SO2 absorption at the 10 to the -8th level. The upper limit derived from the CS2 band head at 3206 A is not greater than 5 x 10 to the -8th. The observed range of SO2 mixing ratios is consistent with model predictions based on the sulfur photochemistry at the cloud tops. Ground-based observations of SO2 mixing ratio will provide constraints on models and check on the Venera and Pioneer Venus measurements of the mixing ratios of SO2 and other sulfur-bearing gases with altitude.

Barker, E. S.↗

The isotopic composition of galactic cosmic ray lithium, beryllium and boron

The isotopic composition of galactic-cosmic-ray Li, Be, and B has been measured near 100 MeV/nucleon by using the University of Chicago IMP 7 and IMP 8 cosmic-ray telescopes during 1973-1975. The measured abundances allow detailed checks of models of interstellar propagation and solar modulation to be made and conclusions to be drawn concerning the spectral forms at the source and the minimum solar modulation level. For example, comparing these results with local interstellar spectra calculated by using a 'leaky box' model, it is found that if solar modulation is ignored, there is no unique leakage mean free path consistent with all the observations. However, by taking account of a sizable level of residual solar modulation, excellent agreement is obtained between the calculated and measured abundances. Thus, these isotopic abundances confirm the old hypothesis that cosmic-ray Li, Be, and B are produced as secondaries in interstellar space.

Garcia-Munoz, M.↗