Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Models”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Issue Update to LWR SiC/SiC Cladding Handbook of Properties

This report aims to provide properties data of silicon carbide (SiC) fiber–reinforced SiC matrix (SiC/SiC) composites for use in modeling for their application as accident-tolerant fuel cladding for light water reactors (LWRs). The properties data presented are specifically for geometries relevant to fuel cladding tubes, although many properties can be derived or predicted based on data from planar geometry specimen data. This report presents various properties, including mechanical properties, thermal properties, chemical stability under normal and off-normal operation conditions, hermeticity, and irradiation resistance. Table S.1 provides a high level summary of those properties mainly for nuclear-grade SiC/SiC composites fabricated via a chemical vapor infiltration (CVI) route, which has been updated from the previous report (ORNL/TM - 2017/385). The newly added data include statistical tensile properties, hydrothermal corrosion behavior of environmental barrier coating on SiC, and accurate temperature- and dose-dependent swelling of SiC/SiC composites.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Evaluation and Validation of the n+ 63,65 Cu Cross Sections [Slides]

This presentation discusses the objective of the research which is to update n+ 63,65 Cu cross section evaluations with recently measured data to resolve discrepancies in benchmark performance. It also discusses the models used which are the R -matrix analysis and the analysis of angular distribution coefficients. Additionally, the validation methods that were used are discussed, including the Rez shielding benchmark and the ICSBEP criticality benchmarks. In conclusion, he n+ 63,65 Cu cross sections have been updated via R-matrix analysis up to 100 keV. An increased average capture cross section and the adoption of experimentally based Legendre coefficients lead to improved performance in reactivity benchmarks. In the fast region, the adoption of the JENDL-4.0 cross sections above 4.0 MeV improves the performance in shielding benchmarks. Ultimately, the n 63,65 Cu ENDF files will be submitted to ENDF/B-VIII.1.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Automated Calculation of Bulk-Stamping Matrices for Circular Blanks

The bulk stamping of circular blanks is considered. There is a clear need for automatic stamping. The input and output data for calculating the working part of the bulk-stamping matrix are analyzed. Software for automated calculation of bulk-stamping matrices is developed. The program code and a possible interface for the calculation of matrix parameters and construction of 3D models by means of the Object Pascal language are presented.

42 ENGINEERING↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2016

JASMINE: a fuel rod thermal-mechanical performance's code

JASMINE, which is developed by China General Nuclear Power Corporation (CGNPC), is a code for calculating thermal-mechanical response of PWR fuel rods. The code models many physical phenomena, including heat conduction through the fuel and cladding to the coolant, cladding elastic and plastic deformation, fuel elastic deformation, fuel swelling and densification,fuel-cladding mechanical interaction, fission gas release from the fuel, rod internal pressure, cladding oxidation and hydriding. Just as most fuel rod codes, JASMINE, which is a quasi-two dimensional code, includes a thermal module, a mechanical module, a fission gas release (FGR) module, and an internal pressure module. For thermal module, the time dependent thermal conduction equation of axisymmetric cylindrical geometry, Halden's fuel thermal conductivity, thermal conductivity integration method and multi-dimension Newton-Raphson procedure are used; for the mechanical module, generalized plane strain hypothesis, cladding-fuel axial lockup model, cladding thermal stress and transfer matrix method are applied; for fission gas release module, the time dependent diffusion equation for a spherical grain, Speight-Turnbull grain surface boundary condition, and Forsberg-Massih's method are used; for the internal pressure module, the ideal gas state equation is used. The code, which is on validation stage, will be used for fuel rod design of nuclear power plants in China in a near future. So far, the preliminary validation process for pellet centerline temperature and fission gas release, which is based on fuel rod test data in FUMEX-III project, has finished. According to the validation results, JASMINE has the ability to reasonably predicts fuel behavior and performance in-pile. JASMINE has good prediction in fuel temperature, and good prediction on FGR as a whole. However, for FGR in the range 1% to 15%, JASMINE over-predicts FGR; for other range, JASMINE underestimates FGR. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Complementarity and the unitarity of the black hole S -matrix

Recently, Akers et al. proposed a non-isometric holographic map from the interior of a black hole to its exterior. Within this model, we study properties of the black hole S-matrix, which are in principle accessible to observers who stay outside the black hole. Specifically, we investigate a scenario in which an infalling agent interacts with radiation both outside and inside the black hole. Because the holographic map involves postselection, the unitarity of the S-matrix is not guaranteed in this scenario, but we find that unitarity is satisfied to very high precision if suitable conditions are met. If the internal black hole dynamics is described by a pseudorandom unitary transformation, and if the operations performed by the infaller have computational complexity scaling polynomially with the black hole entropy, then the S-matrix is unitary up to corrections that are superpolynomially small in the black hole entropy. Furthermore, while in principle quantum computation assisted by postselection can be very powerful, we find under similar assumptions that the S-matrix of an evaporating black hole has polynomial computational complexity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Neutrinoless Double-Beta Decay Investigations of 82Se Using Three Shell Model Hamiltonians

Neutrinoless double-beta decay is considered one of the most promising processes that would help clarify some of the symmetry-breaking problems in our understanding of the observable universe. Recent studies of neutrinoless double-beta decay matrix elements have employed statistical approaches based on modified shell model effective Hamiltonians for 48Ca (Phys. Rev. C 106, 054302 (2022)) and 136Xe (Phys. Rev. C 107, 045501 (2023)). The analyses rely on inducing perturbations in the starting effective Hamiltonians to observe the behavior of a wide range of observables, besides the 0νββ) NME, that are compared with experimental data. Following a Bayesian model averaging approach, the range of probable values for the neutrinoless double-beta decay matrix elements is presented. In this paper, we present a similar study for 82Se, which is described in the same model space as 76Ge that is under experimental observation. Due to its faster calculation time compared to 76Ge, 82Se can be used as an appropriate substitute in our complex statistical study. Using the calculations performed for the statistical analysis of the neutrinoless double-beta decay matrix elements, we also search for the correlations between the observables that we can compare to experimental data.

Science & Technology - Other Topics↗
DOE OSTI2014

Trigonometrical sums connected with the chiral Potts model, Verlinde dimension formula, two-dimensional resistor network, and number theory

We have recently developed methods for obtaining exact two-point resistance of the complete graph minus N edges. We use these methods to obtain closed formulas of certain trigonometrical sums that arise in connection with one-dimensional lattice, in proving Scott’s conjecture on permanent of Cauchy matrix, and in the perturbative chiral Potts model. The generalized trigonometrical sums of the chiral Potts model are shown to satisfy recursion formulas that are transparent and direct, and differ from those of Gervois and Mehta. By making a change of variables in these recursion formulas, the dimension of the space of conformal blocks of SU(2) and SO(3) WZW models may be computed recursively. Our methods are then extended to compute the corner-to-corner resistance, and the Kirchhoff index of the first non-trivial two-dimensional resistor network, 2×N. Finally, we obtain new closed formulas for variant of trigonometrical sums, some of which appear in connection with number theory. -- Highlights: • Alternative derivation of certain trigonometrical sums of the chiral Potts model are given. • Generalization of these trigonometrical sums satisfy recursion formulas. • The dimension of the space of conformal blocks may be computed from these recursions. • Exact corner-to-corner resistance, the Kirchhoff index of 2×N are given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Neutrinoless double beta decay via light neutralinos in R-parity violating supersymmetry

We perform a study of neutrinoless double beta (0νββ) decay mediated by the lightest neutralino of arbitrary mass in the minimal supersymmetric Standard Model (MSSM) under the presence of R-parity violating trilinear interactions. In this scenario, the exchange of the lightest neutralino can result in 0νββ decay of either long-range or short-range behaviour, depending on the neutralino mass. Using nuclear matrix elements calculated in the interacting boson model, we use an interpolation between the long- and short-range behaviours with an approximate formula. The non-observation of 0νββ decay is then used to place constraints on the supersymmetry parameter space, compatible with constraints from collider experiments. We compare these constraints to bounds from pion decays, CKM unitarity and big bang nucleosynthesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Updated constraints on the light-neutrino exchange mechanisms of the 0νββ-decay

The neutrinoless double-beta (0νββ) decay associated with light neutrino exchange mechanisms, which are due to both left-handed V-A and right-handed V+A leptonic and hadronic currents, is discussed by using the recent progress achieved by the GERDA, EXO and KamlandZen experiments. The upper limits for effective neutrino mass m{sub ββ} and the parameters 〈λ〉 and 〈η〉 characterizing the right handed current mechanisms are deduced from the data on the 0νββ-decay of {sup 76}Ge and {sup 136}Xe using nuclear matrix elements calculated within the nuclear shell model and quasiparticle random phase approximation and phase-space factors calculated with exact Dirac wave functions with finite nuclear size and electron screening. The careful analysis of upper constraints on effective lepton number violating parameters assumes a competition of the above mechanisms and arbitrary values of involved CP violating phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗
DOE OSTI2017

In silico investigation of blast-induced intracranial fluid cavitation as it potentially leads to traumatic brain injury

In this paper, we conducted computational macroscale simulations predicting blast-induced intracranial fluid cavitation possibly leading to brain injury. To further understanding of this problem, we developed microscale models investigating the effects of blast-induced cavitation bubble collapse within white matter axonal fiber bundles of the brain. We model fiber tracks of myelinated axons whose diameters are statistically representative of white matter. Nodes of Ranvier are modeled as unmyelinated sections of axon. Extracellular matrix envelops the axon fiber bundle, and gray matter is placed adjacent to the bundle. Cavitation bubbles are initially placed assuming an intracranial wave has already produced them. Pressure pulses, of varied strengths, are applied to the upper boundary of the gray matter and propagate through the model, inducing bubble collapse. Simulations, conducted using the shock wave physics code CTH, predict an increase in pressure and von Mises stress in axons downstream of the bubbles after collapse. This appears to be the result of hydrodynamic jetting produced during bubble collapse. Interestingly, results predict axon cores suffer significantly lower shear stresses from proximal bubble collapse than does their myelin sheathing. Finally, simulations also predict damage to myelin sheathing, which, if true, degrades axonal electrical transmissibility and general health of the white matter structures in the brain.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Multistate Landau-Zener models with all levels crossing at one point

Within this paper, we discuss common properties and reasons for integrability in the class of multistate Landau-Zener models with all diabatic levels crossing at one point. Exploring the Stokes phenomenon, we show that each previously solved model has a dual one, whose scattering matrix can be also obtained analytically. For applications, we demonstrate how our results can be used to study conversion of molecular into atomic Bose condensates during passage through the Feshbach resonance, and provide purely algebraic solutions of the bowtie and special cases of the driven Tavis-Cummings model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Increasing the power density when using inert matrix fuels to reduce production of transuranics

Reducing the production of transuranics is a goal of most advanced nuclear fuel cycles. One way to do this is to recycle the transuranics into the same reactors that are currently producing them using an inert matrix fuel. In previous work we have modeled such a reactor where 72%, of the core is comprised of standard enriched uranium fuel pins, with the remaining 28% fuel made from Yttria stabilized zirconium, in which transuranics are loaded. A key feature of this core is that all of the transuranics produced by the uranium fuel assemblies are later burned in inert matrix fuel assemblies. It has been shown that this system can achieve reductions in transuranic waste of more than 86%. The disadvantage of such a system is that the core power rating must be significantly lower than a standard pressurized water reactor. One reason for the lower power is that high burnup of the uranium fuel precludes a critical level of reactivity at the end of the campaign. Increasing the uranium enrichment and changing the pin pitch are two ways to increase burnup while maintaining criticality. In this paper we use MCNPX and a linear reactivity model to quantify the effect of these two parameters on the end of campaign reactivity. Importantly, we show that in the region of our proposed reactor, enrichment increases core reactivity by 0.02 per percent uranium 235 and pin pitch increases reactivity by 0.02 per mm. Reactivity is lost at a rate of 0.005 per MWd/kgIHM uranium burnup. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Distributions of off-diagonal scattering matrix elements: Exact results

Scattering is a ubiquitous phenomenon which is observed in a variety of physical systems which span a wide range of length scales. The scattering matrix is the key quantity which provides a complete description of the scattering process. The universal features of scattering in chaotic systems is most generally modeled by the Heidelberg approach which introduces stochasticity to the scattering matrix at the level of the Hamiltonian describing the scattering center. The statistics of the scattering matrix is obtained by averaging over the ensemble of random Hamiltonians of appropriate symmetry. We derive exact results for the distributions of the real and imaginary parts of the off-diagonal scattering matrix elements applicable to orthogonally-invariant and unitarily-invariant Hamiltonians, thereby solving a long standing problem. -- Highlights: •Scattering problem in complex or chaotic systems. •Heidelberg approach to model the chaotic nature of the scattering center. •A novel route to the nonlinear sigma model based on the characteristic function. •Exact results for the distributions of off-diagonal scattering-matrix elements. •Universal aspects of the scattering-matrix fluctuations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Exact diagonalization library for quantum electron models

We present an exact diagonalization C++ template library (EDLib) for solving quantum electron models,including the single-band finite Hubbard cluster and the multi-orbital impurity Anderson model. Theobservables that can be computed using EDLib are single particle Green’s functions and spin–spincorrelation functions. This code provides three different types of Hamiltonian matrix storage that canbe chosen based on the model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Mixed quantum/classical theory of rotationally and vibrationally inelastic scattering in space-fixed and body-fixed reference frames

We formulated the mixed quantum/classical theory for rotationally and vibrationally inelastic scattering process in the diatomic molecule + atom system. Two versions of theory are presented, first in the space-fixed and second in the body-fixed reference frame. First version is easy to derive and the resultant equations of motion are transparent, but the state-to-state transition matrix is complex-valued and dense. Such calculations may be computationally demanding for heavier molecules and/or higher temperatures, when the number of accessible channels becomes large. In contrast, the second version of theory requires some tedious derivations and the final equations of motion are rather complicated (not particularly intuitive). However, the state-to-state transitions are driven by real-valued sparse matrixes of much smaller size. Thus, this formulation is the method of choice from the computational point of view, while the space-fixed formulation can serve as a test of the body-fixed equations of motion, and the code. Rigorous numerical tests were carried out for a model system to ensure that all equations, matrixes, and computer codes in both formulations are correct.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Analysis of transverse shear strains in pre-twisted thick beams using variational asymptotic method

The cross-sectional stiffness matrix is derived for a pre-twisted, moderately thick beam made of transversely isotropic materials and having rectangular cross sections. An asymptotically-exact methodology is used to model the anisotropic beam from 3-D elasticity, without any further assumptions. The beam is allowed to have large displacements and rotations, but small strain is assumed. The strain energy is computed making use of the beam constitutive law and kinematical relations derived with the inclusion of geometrical nonlinearities and an initial twist. The energy functional is minimized making use of the Variational Asymptotic Method (VAM), thereby reducing the cross section to a point on the beam reference line with appropriate properties, forming a 1-D constitutive law. VAM is a mathematical technique employed in the current problem to rigorously split the 3-D analysis of beams into two: a 2-D analysis over the beam cross-sectional domain, which provides a compact semi-analytical form of the properties of the cross sections, and a nonlinear 1-D analysis of the beam ref-erence curve. In this method, as applied herein, the cross-sectional analysis is performed asymptotically by taking advantage of a material small parameter and two geometric small parameters. 3-D strain components are derived using kinematics and arranged in orders of the small parameters. Closed-form expressions are derived for the 3-D non-linear warping and stress fields. Warping functions are obtained by the minimization of strain energy subject to certain set of constraints that render the 1-D strain measures well-defined. The zeroth-order 3-D warping field thus yielded is then used to integrate the 3-D strain energy density over the cross section, resulting in the 1-D strain energy density, which in turn helps identify the corresponding cross-sectional stiffness matrix. The model is capable of predicting interlaminar and transverse shear stresses accurately up to first order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗