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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Toward inflation models compatible with the no-boundary proposal

In this paper, we investigate various inflation models in the context of the no-boundary proposal. We propose that a good inflation model should satisfy three conditions: observational constraints, plausible initial conditions, and naturalness of the model. For various inflation models, we assign the probability to each initial condition using the no-boundary proposal and define a quantitative standard, typicality, to check whether the model satisfies the observational constraints with probable initial conditions. There are three possible ways to satisfy the typicality criterion: there was pre-inflation near the high energy scale, the potential is finely tuned or the inflationary field space is unbounded, or there are sufficient number of fields that contribute to inflation. The no-boundary proposal rejects some of naive inflation models, explains some of traditional doubts on inflation, and possibly, can have observational consequences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A model for trapping and re-solution regarding intra-granular bubbles in UO 2 , linked to atomic-scale simulations

In the literature, a clear definition of the irradiation re-solution frequency of gas from bubbles in the UO 2 fuel is absent. Moreover, for intra-granular bubbles, a detailed calculation of the cumulated displaced gas quantities in function of the distance from the radius of the bubble after a re-solution event has never been published. The assessment of these two elements is very useful if we want to increase the adherence of fission gas release codes to our present knowledge of the behavior of fission gases. Hence, we suggest to link the definition of the re-solution frequency to atomic-scale simulations. Furthermore, we present the cumulated displaced gas quantities obtained from Molecular Dynamics calculations, from which we have derived a re-solution profile that can be exploited to better consider the irradiation re-solution phenomenon inside Fission Gas Release codes. On top of that, we have built a new trapping/re-solution model for intra-granular bubbles linked to Molecular Dynamics simulations that can be easily incorporated into Fission Gas Release codes. In conclusion, we also check that the model is properly built through the comparison of the new model against a reference.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Capabilities of multivariate Bayesian inference toward seismic hazard assessment

Multivariate Bayesian analysis can bring significant benefits to seismic hazard analysis: Its multivariate feature enables computing scalar and vector hazard without making any approximations; Correlations between intensity measures are implicitly modeled, permitting direct simulation of ground motion selection tools such as the conditional mean spectrum and the generalized conditioning intensity measure; and Its updating feature enables a seamless integration of new ground motion data into the hazard results. Here, we first develop a multivariate Bayesian ground motion model through the NGA-West2 database. The model functional form considers fault-type, magnitude, and distance dependencies, and also the linear and the rock intensity dependent site response. We use a hybrid Markov Chain Monte Carlo sampling to perform Bayesian inference consisting of Gibbs step and a multilevel Metropolis-Hastings step. We then perform several checks on the model and note that its performance is satisfactory. Finally, we illustrate the merits of this multivariate Bayesian analysis, which include: ground motion model updating with ground motion data recorded in the last four years not part of the NGA-West2 database; computation of scalar and vector seismic hazard using the un-updated and updated ground motion models for Los Angeles, CA; and simulation of the conditional mean spectrum under scalar and vector IM conditioning while accounting for different sources of aleatoric and epistemic uncertainties.

58 GEOSCIENCES↗
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2018

Stability and bifurcation analysis of Japanese encephalitis model with/without effects of some control parameters

In this paper, a mathematical model on transmission of Japanese Encephalitis disease has been developed considering some control parameters and time-dependent environmental-carrying capacity. Here, the total vector population is divided into two subpopulations such as susceptible mosquito and infected mosquito. Here also, total reservoir population (i.e., the population in which the encephalitis virus grows) such as pig, horse, etc., has been considered which is divided into three subpopulations such as susceptible reservoir, infected reservoir and recovered reservoir. Total human population is also divided into three subpopulations such as susceptible human, infected human and recovered human. The dynamical behaviors of the system have been investigated. Here, the basic reproduction number associated with the system has been analyzed with respect to control parameters both theoretically and numerically. The biological feasible equilibria and their stability properties have been discussed and the existence condition of the disease has been illustrated numerically. For a certain set of parametric values, effectiveness of control parameters of our proposed model has been checked numerically. At last, Hopf bifurcations have been made numerically without considering control parameters for the case of constant environmental-carrying capacity of mosquito.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Infrared dynamics of a large N QCD model, the massless string sector and mesonic spectra

A consistency check for any UV complete model for large N QCD should be, among other things, the existence of a well-defined vector and scalar mesonic spectra. In this paper, we use our UV complete model in type IIB string theory to study the IR dynamics and use this to predict the mesonic spectra in the dual type IIA side. The advantage of this approach is two-fold: not only will this justify the consistency of the supergravity approach, but it will also give us a way to compare the IR spectra and the model with the ones proposed earlier by Sakai and Sugimoto. Interestingly, the spectra coming from the massless stringy sector are independent of the UV physics, although the massive string sector may pose certain subtleties regarding the UV contributions as well as the mappings to actual QCD. Additionally, we find that a component of the string landscape enters the picture: there are points in the landscape where the spectra can be improved somewhat over the existing results in the literature. These points in the landscape in-turn also determine certain background supergravity components and fix various pathologies that eventually lead to a consistent low energy description of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Application of MMC model on simulation of shearing process of thick hot-rolled high strength steel plate

Shear operation is widely used as the first step in sheet metal forming to cut the sheet or plate into the required size. The shear of thick hot-rolled High Strength Steel (HSS) requires large shearing force and the sheared edge quality is relatively poor because of the large thickness and high strength compared with the traditional low carbon steel. Bad sheared edge quality will easily lead to edge cracking during the post-forming process. This study investigates the shearing process of thick hot-rolled HSS plate metal, which is generally exploited as the beam of heavy trucks. The Modified Mohr-Coulomb fracture criterion (MMC) is employed in numerical simulation to calculate the initiation and propagation of cracks during the process evolution. Tensile specimens are designed to obtain various stress states in tension. Equivalent fracture strains are measured with Digital Image Correlation (DIC) equipment to constitute the fracture locus. Simulation of the tension test is carried out to check the fracture model. Then the MMC model is applied to the simulation of the shearing process, and the simulation results show that the MMC model predicts the ductile fracture successfully.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Implementing the three-particle quantization condition for π + π + K + and related systems

Recently, the formalism needed to relate the finite-volume spectrum of systems of nondegenerate spinless particles has been derived. In this work we discuss a range of issues that arise when implementing this formalism in practice, provide further theoretical results that can be used to check the implementation, and make available codes for implementing the three-particle quantization condition. Specifically, we discuss the need to modify the upper limit of the cutoff function due to the fact that the left-hand cut in the scattering amplitudes for two nondegenerate particles moves closer to threshold; we describe the decomposition of the three-particle amplitude K df,3 into the matrix basis used in the quantization condition, including both s and p waves, with the latter arising in the amplitude for two nondegenerate particles; we derive the threshold expansion for the lightest three-particle state in the rest frame up to O(1/L 5 ); and we calculate the leading-order predictions in chiral perturbation theory for K df,3 in the π + π + K + and π + K + K + systems. We focus mainly on systems with two identical particles plus a third that is different (“2+1” systems). We describe the formalism in full detail, and present numerical explorations in toy models, in particular checking that the results agree with the threshold expansion, and making a prediction for the spectrum of π + π + K + levels using the two- and three-particle interactions predicted by chiral perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Subgrid-scale model for radiative transfer in turbulent participating media

The simulation of turbulent flows of radiating gases, taking into account all turbulence length scales with an accurate radiation transport solver, is computationally prohibitive for high Reynolds or Rayleigh numbers. This is particularly the case when the small structures are not optically thin. We develop in this paper a radiative transfer subgrid model suitable for the coupling with direct numerical simulations of turbulent radiating fluid flows. Owing to the linearity of the Radiative Transfer Equation (RTE), the emission source term is spatially filtered to define large-scale and subgrid-scale radiation intensities. The large-scale or filtered intensity is computed with a standard ray tracing method on a coarse grid, and the subgrid intensity is obtained analytically (in Fourier space) from the Fourier transform of the subgrid emission source term. A huge saving of computational time is obtained in comparison with direct ray tracing applied on the fine mesh. Model accuracy is checked for three 3D fluctuating temperature fields. The first field is stochastically generated and allows us to discuss the effects of the filtering level and of the optical thicknesses of the whole medium, of the integral length scale, and of the cutoff wave length. The second and third cases correspond respectively to turbulent natural convection of humid air in a cubical box, and to the flow of hot combustion products inside a channel. In all cases, the achieved accuracy on radiative powers and wall fluxes is about a few percents.

42 ENGINEERING↗
DOE OSTI2025

Missing components in ΛCDM from DESI Y1 baryonic acoustic oscillation measurements: Insights from redshift remapping

We explore transformations of the Friedman-Lemaître-Robertson-Walker (FLRW) metric and cosmological parameters that align with observational data while aiming to gain insights into potential extensions of standard cosmological models. We modified the FLRW metric by introducing a scaling factor, e 2Θ(a) –the cosmological scaling function (CSF), which alters the standard relationship between cosmological redshift and the cosmic scale factor without affecting angular measurements or cosmic microwave background (CMB) anisotropies. Using data from DESI Year 1, Pantheon+ supernovae, and the Planck CMB temperature power spectrum, we constrained both the CSF and cosmological parameters through a Markov chain Monte Carlo approach. Our results indicate that the CSF model fits observational data with a lower Hubble constant (although it is compatible with the value given by Planck 2018 within 1σ) and is predominantly dark matter dominated. Additionally, the CSF model produces temperature and lensing power spectra similar to those predicted by the standard model, though with lower values in the CSF model at large scales. We also checked that when fitting a CSF model without dark energy to the data, we obtain a more negative conformal function. This suggests that the CSF model may offer hints about missing elements and opens up a new avenue for exploring physical interpretations of cosmic acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Possible detection of a cyclotron resonance scattering feature in the X-ray pulsar 4U 1909+07

We present timing and broad band spectral studies of the high-mass X-ray binary pulsar 4U 1909+07 using data from Suzaku observations during 2010 November 2-3. The pulse period of the pulsar is estimated to be 604.11 ± 0.14 s. Pulsations are seen in the X-ray light curve up to ∼70 keV. The pulse profile is found to be strongly energy-dependent: a complex, multi-peaked structure at low energy becomes a simple single peak at higher energy. We found that the 1-70 keV pulse-averaged continuum can be fit by the sum of a blackbody and a partial covering Negative and Positive power law with Exponential cutoff model. A weak iron fluorescence emission line at 6.4 keV was detected in the spectrum. An absorption-like feature at ∼44 keV was clearly seen in the residuals of the spectral fitting, independent of the continuum model adopted. To check the possible presence of a cyclotron resonance scattering feature (CRSF) in the spectrum, we normalized the pulsar spectrum with the spectrum of the Crab Nebula. The resulting Crab ratio also showed a clear dip centered at ∼44 keV. We performed statistical tests on the residuals of the spectral fitting and also on the Crab spectral ratio to determine the significance of the absorption-like feature and identified it as a CRSF of the pulsar. We estimated the corresponding surface magnetic field of the pulsar to be 3.8 × 10{sup 12} G.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Higgs-flavon mixing and LHC phenomenology in a simplified model of broken flavor symmetry

Here, the LHC phenomenology of a low-scale gauged flavor symmetry model with inverted hierarchy is studied, through introduction of a simplified model of broken flavor symmetry. A new scalar (a flavon) and a new neutral top-philic massive gauge boson emerge with mass in the TeV range, along with a new heavy fermion associated with the standard model top quark. After checking constraints from electroweak precision observables, we investigate the influence of the model on Higgs boson physics, notably on its production cross section and decay branching fractions. Limits on the flavon φ from heavy Higgs boson searches at the LHC at 7 and 8 TeV are presented. The branching fractions of the flavon are computed as a function of the flavon mass and the Higgs-flavon mixing angle. We also explore possible discovery of the flavon at 14 TeV, particularly via the φ → Z 0 Z 0 decay channel in the 2ℓ2ℓ' final state, and through standard model Higgs boson pair production φ → hh in the b¯bγγ final state. We conclude that the flavon mass range up to 500 GeV could be probed down to quite small values of the Higgs-flavon mixing angle with 100 fb –1 of integrated luminosity at 14 TeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

toughio v1.2.2

TOUGH is a widely used general purpose numerical simulator designed for fluid and heat flows of multiphase, multicomponent fluid mixtures in porous and fractured media, which has been applied to many real-world problems such as underground geological storage, geothermal reservoir engineering or nuclear waste disposal, to name a few. `toughio` is a lightweight, object-oriented and vectorized Python library that aims to provide user-friendly routines to facilitate pre- and post-processing of a TOUGH simulation. Although its primary target is the latest version TOUGH3, `toughio` can also be used to read/write TOUGH2 input/output files. Currently, to the best of our knowledge, only PyTOUGH offers an exhaustive number of features to carry out a complete TOUGH simulation using a scripting language. `toughio` and PyTOUGH share the same objectives, yet with different approaches and philosophies. On the one hand, a PyTOUGH mesh is represented as a MULGRAPH geometry where elements can be unstructured horizontally but only layered vertically. It also checks that the model parameters are valid to ensure that a TOUGH simulation runs smoothly. On the other hand, although it provides basic meshing features, `toughio` mostly relies on common third-party softwares (e.g. Abaqus, FLAC3D, Gmsh, LaGriT) to generate the mesh by importing and converting it to a TOUGH mesh, which also conveniently facilitates the coupling of TOUGH with any other simulator that also supports the same mesh formats. Besides, `toughio` is designed to be just a TOUGH input and output file reader/writer. It assumes that the user knows what he is doing and does not verify the validity of the model parameters (`toughio` does verify the functional arguments).

Luu, Keurfon↗
DOE OSTI2016

SU-F-T-158: Experimental Characterization of Field Size Dependence of Dose and Lateral Beam Profiles of Scanning Proton and Carbon Ion Beams for Empirical Model in Air

Purpose: The Gaussian model for the lateral profiles in air is crucial for an accurate treatment planning system. The field size dependence of dose and the lateral beam profiles of scanning proton and carbon ion beams are due mainly to particles undergoing multiple Coulomb scattering in the beam line components and secondary particles produced by nuclear interactions in the target, both of which depend upon the energy and species of the beam. In this work, lateral profile shape parameters were fitted to measurements of field size dependence dose at the center of field size in air. Methods: Previous studies have employed empirical fits to measured profile data to significantly reduce the QA time required for measurements. From this approach to derive the weight and sigma of lateral profiles in air, empirical model formulations were simulated for three selected energies for both proton and carbon beams. Results: The 20%–80% lateral penumbras predicted by the double model for proton and single model for carbon with the error functions agreed with the measurements within 1 mm. The standard deviation between measured and fitted field size dependence of dose for empirical model in air has a maximum accuracy of 0.74% for proton with double Gaussian, and of 0.57% for carbon with single Gaussian. Conclusion: We have demonstrated that the double Gaussian model of lateral beam profiles is significantly better than the single Gaussian model for proton while a single Gaussian model is sufficient for carbon. The empirical equation may be used to double check the separately obtained model that is currently used by the planning system. The empirical model in air for dose of spot scanning proton and carbon ion beams cannot be directly used for irregular shaped patient fields, but can be to provide reference values for clinical use and quality assurance.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Constitutive Modeling for Hot Deformation Behavior of Al-5083 + SiC Composite

In the present research, isothermal hot deformation behavior of aluminum 5083 alloy + 15 (wt.%) SiC composite was obtained through compression test on INSTRON 8801 universal tensile testing machine (UTM) under a wide temperature range of 473-773 K and strain rate range of 0.01-10 s{sup −1}. The experimental true stress–strain data were employed to establish constitutive equations based on modified Johnson–Cook (JC) model and modified Zerilli–Armstrong (ZA) model to predict the hot flow behavior of the composite. The flow stress values obtained from these two models were plotted against the experimental flow curves to check the accuracy of these models. Suitability of the models was evaluated by comparing correlation coefficient (R), average absolute relative error and relative errors of prediction. The results show that the hot flow stresses of the present material depend on temperature and strain rate significantly. Both the models give good description of the hot deformation behavior of the composite. The prediction accuracy is found to be higher for modified ZA model compared to modified JC model, though the number of materials constants involved and time needed to evaluate them to establish the model are lower for modified JC model.

36 MATERIALS SCIENCE↗