Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice calculation”

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 127 records · Page 7

Upper Limit on the Photoproduction Cross Section of the Spin-Exotic $π$ 1 (1600)

Here, the spin-exotic hybrid meson π 1 (1600) is predicted to have a large decay rate to the ωππ final state. Using 76.6 pb -1 of data collected with the GlueX detector, we measure the cross sections for the reactions γp → ωπ + π - p, γp → ωπ 0 π 0 p, and γp → ωπ - π 0 Δ ++ in the range E γ = 8 - 10 GeV. Using isospin conservation, we set the first upper limits on the photoproduction cross sections of the π$^0_1$(1600) and $π^-_1$(1600). We combine these limits with lattice calculations of decay widths and find that photoproduction of η' π is the most sensitive two-body system to search for the π 1 (1600).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗

Pion and Kaon decays from lattice QCD and |Vud|, |Vus|

We review lattice calculations of the form factors for Kaon semileptonic decays and the leptonic decay constants of Kaons and Pions. We discuss the determinations of the CKM matrix elements $|V_{ud}|$ and $|V_{us}|$ using the form factors and experimental results. We use the CKM matrix elements to test first-row unitarity of the CKM matrix.

Simone, J. N. [Fermilab] (ORCID:0000000185153337)↗

The Three-Dimensional Structure of the Proton

This project addresses three key questions about the “glue that binds us all” (Nuclear Science Advisory Council’s Long-Range Plan). First, what is the three-dimensional momentum distribution of gluons within the nucleon? Second, how is the motion of quarks and gluons correlated within hadrons and how do those correlations manifest in nucleon structure? Third, how can lattice calculations guide the experimental study of hadron structure at the future electron-ion collider (EIC)?

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Baryon $\sigma$ terms in SU(3) BChPT x 1/N c

ChPT and the 1/N c expansion provide systematic frameworks for the strong interactions at low energy. A combined framework of both expansions has been developed and applied for baryons with three light-quark-flavors. The small scale expansion of the combined approach is identified as the ξ-expansion, in which the power counting of the expansions is linked according to O(p) = O(1/N c ) = O(ξ). The physical baryon masses as well as lattice QCD baryon masses for different quark mass masses are analyzed to O(ξ 3 ) in that framework. σ terms are addressed using the Feynman Hellmann theorem. For the nucleon, a useful connection between the deviation of the Gell-Mann-Okubo relation and the σ term σ 8N associated with the scalar density $\bar{u}u+\bar{d}d-2\bar{s}s$ is identified. In particular, the deviation from the tree level relation σ 8N = $\frac{1}{3}$(2m N -m Σ -m Ξ ), which gives rise to the so called σ-term puzzle, is studied in the ξ-expansion. A large correction non-analytic in ξ results for that relation, making plausible the resolution of the puzzle. Issues with the determination of the strangeness σ terms are discussed, emphasizing the need for lattice calculations at smaller m s for better understanding the range of validity of the effective theory. The analysis presented here leads to σ πN = 69(10) MeV and σ πΔ = 60(10) MeV.

Fernando, Ishara P.↗

Comparison of analytical and experimental steadyand unsteady-pressure distributions at Mach number 0.78 for a high-aspect-ratio supercritical wing model with oscillating control surfaces

The unsteady aerodynamic lifting surface theory, the Doublet Lattice method, with experimental steady and unsteady pressure measurements of a high aspect ratio supercritical wing model at a Mach number of 0.78 were compared. The steady pressure data comparisons were made for incremental changes in angle of attack and control surface deflection. The unsteady pressure data comparisons were made at set angle of attack positions with oscillating control surface deflections. Significant viscous and transonic effects in the experimental aerodynamics which cannot be predicted by the Doublet Lattice method are shown. This study should assist development of empirical correction methods that may be applied to improve Doublet Lattice calculations of lifting surface aerodynamics.

Mccain, W. E.↗

A vortex-lattice method for calculating lifting-surface interference

The interference effect of closely coupled lifting surfaces was investigated. A nonsteady vortex-lattice method was applied to calculate the steady-state and nonsteady lift characteristics of configurations having interference effects. The configurations consist of low-aspect-ratio delta wings, at low and high angles of attack. Leading-edge separation and wake roll-up are simulated by a time-dependent wake-shedding procedure. For steady-state flow conditions, this numerical procedure saves a considerable amount of computer time, compared with iterative methods, and yields the same results. A better understanding of the interference effect can be gained by using the method to study the transient behavior. The nonsteady approach also offers the capability of calculating various nonsteady motions, as is demonstrated in the calculation of the longitudinal damping.

Levin, D.↗

A vortex-lattice method for calculating longitudinal dynamic stability derivatives of oscillating delta wings

A nonsteady vortex-lattice method is introduced for predicting the dynamic stability derivatives of a delta wing undergoing an oscillatory motion. The analysis is applied to several types of small oscillations in pitch. The angle of attack varied between + or - 1 deg, with the mean held at 0 deg when the flow was assumed to be attached and between + or - 1 deg and the mean held at 15 deg when both leading-edge separation and wake roll-up were included. The computed results for damping in pitch are compared with several other methods and with experiments, and are found to be consistent and in good agreement.

Levin, D.↗

General method for calculating derivatives of the lattice electrostatic energy.

A method for calculating the derivatives of lattice electrostatic strain energy is proposed. It offers a computation procedure that is more general, concise, and systematic than any of the procedures previously used by Fuchs (1936), Cousins (1967), and Suzuki et al. (1968). The method can also easily be extended to fourth- and higher-order derivatives without undue difficulty.

Macdonald, D. E.↗

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Diagrams for calculation of airfoil lattices

The field for curved blades is represented by a vortex series with a vortex removed at the blade point. Further, an example of calculation of a curved blade from this series is given, whereby the necessary accuracy required of the different methods in practice is shown according to the case considered.

Betz, Albert↗

A multi-element vortex lattice method for calculating the geometry and effects of a helicopter rotor wake in forward flight

A method is described for the analysis of the unsteady, incompressible potential flow associated with a helicopter rotor and it's wake in forward flight. This method is particularly useful in low advance ratio flight due to the major contribution, in the near field, of the deformed wake. The rotor geometry is prescribed and the unsteady wake geometry is computed from the local flow perturbation velocities. The wake is modeled as a full vortex lattice. The rotor geometry is arbitrary and several rotor blades can be represented. The unsteady airloads on the rotor blades are computed in the presence of the deformed rotor wake by a time-stepping technique. Solution for the load distribution on the blade surfaces is found by prescribing boundary conditions in a reference system which rotates with the blade tips. Transformation tensors are used to describe the contribution of the wake in the inertial system to the rotor in the rotating reference system. The effects of blade cyclic pitch variation are computed using a rotation tensor. The deformation of the wake is computed in the inertial frame. The wake is started impulsively from rest, allowing a natural convection of the wake with time.

Berry, John D.↗

The coverage-dependent effect on the barrier for oxygen penetration into a metal lattice

Model calculations for O/Cu(111) show that the rigid lattice penetration barrier is reduced from 8.2 to 5.3 eV as the O coverage is increased from 1/4 to 1/2 monolayer; it is expected that this large reduction will also apply to the relaxed lattice case. An increase of the distance between the first and second Cu layers, without changing the intralayer spacing, does not affect the barrier.

Bauschlicher, C. W., Jr.↗

Calculation of spin-lattice relaxation during pulsed spin locking in solids

The spin-lattice relaxation time has been calculated for dipolar solids in the case where the spins are locked by an RF pulse sequence with pulses of arbitrary angle and finite width. Expressions are given for the homonuclear case in general and for the heteronuclear case in the delta-function limit. The results for the homonuclear case are experimentally confirmed using solid C6F12. The analysis shows that for small pulse angles, at which the direct spin heating effect is known to be small, the relaxation behavior will be identical to the CW irradiation case.

Rhim, W.-K.↗

Modified lattice-statics approach to dislocation calculations. I - Formalism

A modified lattice-statics method to calculate the atomic displacements associated with a screw dislocation is outlined. The model incorporates an anharmonic region wherein the forces are derived from a pair potential. Appropriate energy and force expressions are derived. The modifications necessary for the implementation of the conjugate-gradient function minimization method are also derived.

Esterling, D. M.↗

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗