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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
NASA NTRS1978

NASTRAN analysis of a wheel-rail loading on its foundation

One type of azimuth bearing for a large ground antenna (100 m) will consist of steel wheels, mounted at four corners of the alidade, rolling on a circular flat rail which provides the vertical restraints; a radial constraining bearing at the center of the alidade provides the horizontal restraints. One important design feature is the compressive stresses in the grout or concrete foundation under the wheel-rail load. A finite element analysis check was made of a particular design that consists of a steel rail resting on a concrete foundation. Symmetry was assumed as much as possible in order to minimize the models, but meaningful element sizes were used. Recently developed isoparametric hexahedron elements available in the NASTRAN computing program, which minimizes the number of elements required while maintaining the accuracy of the computed stresses, were used with two versions of NASTRAN. Test cases to check with the analytical solutions were made. A side loading was also applied to calculate the increase in the concrete stresses.

Katow, M. S.↗
NASA NTRS2016

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.↗
DOE OSTI2014

The abelian confinement mechanism revisited: New aspects of the Georgi–Glashow model

The confinement problem remains one of the most difficult problems in theoretical physics. An important step toward the solution of this problem is Polyakov’s work on abelian confinement. The Georgi–Glashow model is a natural testing ground for this mechanism which has been surprising us by its richness and wide applicability. In this work, we shed light on two new aspects of this model in 2+1 D. First, we develop a many-body description of the effective degrees of freedom. Namely, we consider a non-relativistic gas of W-bosons in the background of monopole–instanton plasma. Many-body treatment is a standard toolkit in condensed matter physics. However, we add a new twist by supplying the monopole–instantons as external background field. Using this construction along with a mean-field approximation, we calculate the form of the potential between two electric probes as a function of their separation. This potential is expressed in terms of the Meijer-G function which interpolates between logarithmic and linear behavior at small and large distances, respectively. Second, we develop a systematic approach to integrate out the effect of the W-bosons at finite temperature in the range 0≤T<M{sub W}, where M{sub W} is the W-boson mass, starting from the full relativistic partition function of the Georgi–Glashow model. Using a heat kernel expansion that takes into account the non-trivial thermal holonomy, we show that the partition function describes a three-dimensional two-component Coulomb gas. We repeat our analysis using the many-body description which yields the same result and provides a check on our formalism. At temperatures close to the deconfinement temperature, the gas becomes essentially two-dimensional recovering the partition function of the dual sine-Gordon model that was considered in a previous work. -- Highlights: • We consider the problem of abelian confinement in the Georgi–Glashow model from a new perspective. • We develop a many-body description of the degrees of freedom of this model in the background of monopole–instanton plasma. • We use the many-body formalism to find an analytic expression of the potential between two probes at zero temperature. • We also use a systematic approach to integrate out the W-bosons starting from the full relativistic partition function. • This results in a three-dimensional two-component Coulomb gas with long range and Aharonov–Bohm phase interaction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Numerical challenges for energy conservation in N -body simulations of collapsing self-interacting dark matter halos

Dark matter (DM) halos can be subject to gravothermal collapse if the DM is not collisionless, but engaged in strong self-interactions instead. When the scattering is able to efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. However, a comparison with observations requires accurate predictions. One method to obtain such predictions is to use N-body simulations. Collapsed halos are extreme systems that pose severe challenges when applying state-of-the-art codes to model self-interacting dark matter (SIDM). In this work, we investigate the root of such problems, with a focus on energy non-conservation. Moreover, we discuss possible strategies to avoid them. We ran N-body simulations, both with and without SIDM, of an isolated DM-only halo and we adjusted the numerical parameters to check the accuracy of the simulation. We find that not only the numerical scheme for SIDM can lead to energy non-conservation, but also the modelling of gravitational interaction and the time integration are problematic. The main issues we find are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; and (c) SIDM velocity kicks breaking the time symmetry. Tuning the parameters of the simulation to achieve a high level of accuracy allows us to conserve energy not only at early stages of the evolution, but also later on. However, the cost of the simulations becomes prohibitively large as a result. Some of the problems that make the simulations of the gravothermal collapse phase inaccurate can be overcome by choosing appropriate numerical schemes. However, other issues still pose a challenge. Our findings motivate further works on addressing the challenges in simulating strong DM self-interactions.

dark matter↗
NASA NTRS2011

Preprocessing Inconsistent Linear System for a Meaningful Least Squares Solution

Mathematical models of many physical/statistical problems are systems of linear equations. Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

Sen, Syamal K.↗
NASA NTRS2011

Preprocessing in Matlab Inconsistent Linear System for a Meaningful Least Squares Solution

Mathematical models of many physical/statistical problems are systems of linear equations~ Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the . minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

Sen, Symal K.↗
DOE OSTI2021

The Atacama Cosmology Telescope: delensed power spectra and parameters

We present ΛCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (TT), 30% (EE), 26% (TE) and 20% (BB); this results in detections of the delensing effect at 8.7σ (TT), 5.1σ (EE), 2.6σ (TE), and 2.4σ (BB) significance. The combination of 150 and 98 GHz TT, EE, and TE delensed spectra is well fit by a standard ΛCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Furthermore, fitting with a ΛCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

The Atacama Cosmology Telescope: delensed power spectra and parameters

We present ΛCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (T T), 30% (EE), 26% (T E) and 20% (BB); this results in detections of the delensing effect at 8.7σ (T T), 5.1σ (EE), 2.6σ (T E), and 2.4σ (BB) significance. The combination of 150 and 98 GHz T T, EE, and T E delensed spectra is well fit by a standard ΛCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Fitting with a ΛCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model.

Dongwon Han↗
DOE OSTI2014

Dynamical model for the toroidal sporadic meteors

More than a decade of radar operations by the Canadian Meteor Orbit Radar have allowed both young and moderately old streams to be distinguished from the dispersed sporadic background component. The latter has been categorized according to broad radiant regions visible to Earth-based observers into three broad classes: the helion and anti-helion source, the north and south apex sources, and the north and south toroidal sources (and a related arc structure). The first two are populated mainly by dust released from Jupiter-family comets and new comets. Proper modeling of the toroidal sources has not to date been accomplished. Here, we develop a steady-state model for the toroidal source of the sporadic meteoroid complex, compare our model with the available radar measurements, and investigate a contribution of dust particles from our model to the whole population of sporadic meteoroids. We find that the long-term stable part of the toroidal particles is mainly fed by dust released by Halley type (long period) comets (HTCs). Our synthetic model reproduces most of the observed features of the toroidal particles, including the most troublesome low-eccentricity component, which is due to a combination of two effects: particles' ability to decouple from Jupiter and circularize by the Poynting-Robertson effect, and large collision probability for orbits similar to that of the Earth. Our calibrated model also allows us to estimate the total mass of the HTC-released dust in space and check the flux necessary to maintain the cloud in a steady state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Multiple soft limits of cosmological correlation functions

We derive novel identities satisfied by inflationary correlation functions in the limit where two external momenta are taken to be small. We derive these statements in two ways: using background-wave arguments and as Ward identities following from the fixed-time path integral. Interestingly, these identities allow us to constrain some of the O(q{sup 2}) components of the soft limit, in contrast to their single-soft analogues. We provide several nontrivial checks of our identities both in the context of resonant non-Gaussianities and in small sound speed models. Additionally, we extend the relation at lowest order in external momenta to arbitrarily many soft legs, and comment on the many-soft extension at higher orders in the soft momentum. Finally, we consider how higher soft limits lead to identities satisfied by correlation functions in large-scale structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Scalar and vector perturbations in a universe with discrete and continuous matter sources

We study a universe filled with dust-like matter in the form of discrete inhomogeneities (e.g., galaxies and their groups and clusters) and two sets of perfect fluids with linear and nonlinear equations of state, respectively. The background spacetime geometry is defined by the FLRW metric. In the weak gravitational field limit, we develop the first-order scalar and vector cosmological perturbation theory. Our approach works at all cosmological scales (i.e. sub-horizon and super-horizon ones) and incorporates linear and nonlinear effects with respect to energy density fluctuations. We demonstrate that the scalar perturbation (i.e. the gravitational potential) as well as the vector perturbation can be split into individual contributions from each matter source. Each of these contributions satisfies its own equation. The velocity-independent parts of the individual gravitational potentials are characterized by a finite time-dependent Yukawa interaction range being the same for each individual contribution. We also obtain the exact form of the gravitational potential and vector perturbation related to the discrete matter sources. The self-consistency of our approach is thoroughly checked. The derived equations can form the theoretical basis for numerical simulations for a wide class of cosmological models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Photon asymmetry measurements of $\overrightarrow{\gamma}p \rightarrow \pi^{0} p$ γ → p → π 0 p for $E_{\gamma}=$ 320-650 MeV

High-statistics measurements of the photon asymmetry Σ for the $\overrightarrow{\gamma}p \rightarrow \pi^{0} p$ reaction have been made in the center-of-mass energy range W = 1214–1450 MeV. The data were measured with the MAMI A2 real photon beam and Crystal Ball/TAPS detector systems in Mainz, Germany. The resulting measurements significantly improve the existing world data and are shown to be in good agreement with previous measurements, and with the MAID, SAID, and Bonn-Gatchina predictions. We have also combined the photon asymmetry results with recent cross-section measurements from Mainz to calculate the profile functions, $\check{Σ}$ (= σ 0 Σ), and perform a moment analysis. Comparison with calculations from the Bonn-Gatchina model shows that the precision of the data is good enough to further constrain the higher partial waves, and there is an indication of interference between the very small F-waves and the N(1520)3/2 - and N(1535)1/2 - resonances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

NVL-C: Static Analysis Techniques for Efficient, Correct Programming of Non-Volatile Main Memory Systems

Computer architecture experts expect that non-volatile memory (NVM) hierarchies will play a more significant role in future systems including mobile, enterprise, and HPC architectures. With this expectation in mind, we present NVL-C: a novel programming system that facilitates the efficient and correct programming of NVM main memory systems. The NVL-C programming abstraction extends C with a small set of intuitive language features that target NVM main memory, and can be combined directly with traditional C memory model features for DRAM. We have designed these new features to enable compiler analyses and run-time checks that can improve performance and guard against a number of subtle programming errors, which, when left uncorrected, can corrupt NVM-stored data. Moreover, to enable recovery of data across application or system failures, these NVL-C features include a flexible directive for specifying NVM transactions. So that our implementation might be extended to other compiler front ends and languages, the majority of our compiler analyses are implemented in an extended version of LLVM's intermediate representation (LLVM IR). We evaluate NVL-C on a number of applications to show its flexibility, performance, and correctness.

Denny, Joel↗
DOE OSTI2021

On the Existence of the Plateau Emission in High-energy Gamma-Ray Burst Light Curves Observed by Fermi-LAT

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (Fermi) shows long-lasting high-energy emission in many gamma-ray bursts (GRBs), similar to X-ray afterglows observed by the Neil Gehrels Swift Observatory (Swift). Some LAT light curves (LCs) show a late-time flattening reminiscent of X-ray plateaus. We explore the presence of plateaus in LAT temporally extended emission analyzing GRBs from the second Fermi-LAT GRB Catalog from 2008 to 2016 May with known redshifts, and check whether they follow closure relations corresponding to four distinct astrophysical environments predicted by the external forward shock model. We find that three LCs can be fit by the same phenomenological model used to fit X-ray plateaus and show tentative evidence for the existence of plateaus in their high-energy extended emission. The most favorable scenario is a slow-cooling regime, whereas the preferred density profile for each GRBs varies from a constant-density interstellar medium to an r –2 wind environment. We also compare the end time of the plateaus in γ-rays and X-rays using a statistical comparison with 222 Swift GRBs with plateaus and known redshifts from 2005 January to 2019 August. Within this comparison, the case of GRB 090510 shows an indication of chromaticity at the end time of the plateau. Finally, we update the 3D fundamental plane relation among the rest-frame end time of the plateau, its correspondent luminosity, and the peak prompt luminosity for 222 GRBs observed by Swift. We find that these three LAT GRBs follow this relation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Study of fault tolerant software technology for dynamic systems

The major aim of this study is to investigate the feasibility of using systems-based failure detection isolation and compensation (FDIC) techniques in building fault-tolerant software and extending them, whenever possible, to the domain of software fault tolerance. First, it is shown that systems-based FDIC methods can be extended to develop software error detection techniques by using system models for software modules. In particular, it is demonstrated that systems-based FDIC techniques can yield consistency checks that are easier to implement than acceptance tests based on software specifications. Next, it is shown that systems-based failure compensation techniques can be generalized to the domain of software fault tolerance in developing software error recovery procedures. Finally, the feasibility of using fault-tolerant software in flight software is investigated. In particular, possible system and version instabilities, and functional performance degradation that may occur in N-Version programming applications to flight software are illustrated. Finally, a comparative analysis of N-Version and recovery block techniques in the context of generic blocks in flight software is presented.

Caglayan, A. K.↗
NASA NTRS1956

Theory of wing-body drag at supersonic speeds

The relation of Whitcomb's "area rule" to the linear formulas for wave drag at lightly supersonic speeds is discussed. By adopting an approximate relation between the source strength and the geometry of a wing-body combination, the wave-drag theory is expressed in terms involving the areas intercepted by oblique planes or Mach planes. The resulting formulas are checked by comparison with the drag measurements obtained in wind-tunnel experiments and in experiments with falling models in free air. Finally, a theory for determining wing-body shapes of minimum drag at supersonic Mach numbers is discussed and some preliminary experiments are reported.

Jones, Robert T↗
NASA NTRS2002

Error-Based Design Space Windowing

Windowing of design space is considered in order to reduce the bias errors due to low-order polynomial response surfaces (RS). Standard design space windowing (DSW) uses a region of interest by setting a requirement on response level and checks it by a global RS predictions over the design space. This approach, however, is vulnerable since RS modeling errors may lead to the wrong region to zoom on. The approach is modified by introducing an eigenvalue error measure based on point-to-point mean squared error criterion. Two examples are presented to demonstrate the benefit of the error-based DSW.

Papila, Melih↗