Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “THREE-DIMENSIONAL CALCULATIONS”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Three-Dimensional Heat Flux and Thermal Analysis of Angled Tungsten Samples on DIII-D

ITER-grade tungsten and dispersoid-strengthened tungsten samples with the top surface angled at ~15° towards the incident plasma flux were exposed to 9 H-mode discharges with edge-localized modes (ELMs) in the lower divertor of DIII-D tokamak using the Divertor Material Evaluation System (DiMES). Surface damage included cracking and flaking of material on the two samples farthest away from the plasma strike point, and significant melting of the two samples closest to the strike point. Heat flux and thermal analysis tools new to DIII-D have been applied to better understand this material response and to help optimize the exposure conditions for future experiments. SMITER field-line tracing simulations based on IRTV data and EFIT equilibria estimate an average inter-ELM perpendicular heat flux, 𝑞⊥,𝑖nter−𝐸LM , on the angled surfaces of 10.1 – 19.6 MW/m² for a majority of the 9 discharges, increasing to 15.6 – 24.5 MW/m² for the single, higher-power shot where samples melted. Fast camera data showed shallow intra-ELM melting and re-solidification, which transitioned to bulk inter-ELM melting with melt motion in the 𝐽⃗ 𝑥 𝐵⃗ direction. About 50% of the protruding volume of the most affected sample was displaced via melt-motion. SIERRA thermal modeling software was able to reproduce an onset time of melting consistent with fast camera data and final sample conditions, within < 200 ms. Maximum surface temperatures of 3122 K and 2787 K are estimated for the samples farthest away from the strike point, while the closest samples achieve melting at 4067 ms and 4750 ms into the ~5000 ms plasma exposure. A +10% increase in both the SMITER 𝑞⊥,𝑖nter−𝐸LM calculations and the estimated ELM heat loads 𝑞⊥, 𝐸LM was required to achieve this result, which is within the uncertainty of the diagnostic data but likely accounts for non-ideal geometry effects plus other physics uncertainties not included in this first iteration of modeling. This work provided valuable estimates of the 3D temperature evolution to help better understand the observed surface morphology and internal recrystallization of samples, which are discussed in detail in a complementary manuscript [1]. Benchmarking efforts with more diagnosed DIII-D experiments are underway to further refine the SMITER and SIERRA models for DiMES. Future use of these tools will enable researchers to precisely target heat flux exposure conditions in DIII-D to test, but not exceed, the thermomechanical limitations of novel plasma-facing materials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2025

Lattice QCD Benchmark of Proton Helicity and Flavor-Dependent Unpolarized Transverse Momentum-Dependent Parton Distribution Functions at Physical Quark Masses

We present the first lattice QCD calculations of the isovector helicity transverse momentum-dependent parton distribution function (TMDPDF) and the flavor-dependent unpolarized TMDPDFs for up and down quarks in the proton. Our computations utilize domain-wall fermion discretization with physical quark masses. Employing Coulomb-gauge-fixed quark correlation functions within the large-momentum effective theory framework, we access nonperturbative transverse quark separations 𝑏 𝑇 up to approximately 1 fm, corresponding to transverse momenta as low as 200 MeV. Based on the quasi-TMD factorization theorem, we construct renormalization-group–invariant ratios that are equal to the corresponding light-cone TMDPDF ratios. At moderate 𝑥, our results reveal that the isovector helicity and unpolarized TMDPDFs exhibit nearly identical transverse structure up to a normalization factor, and the unpolarized distributions display only mild flavor dependence. These findings not only support key trends observed in recent global analyses but also provide robust, nonperturbative constraints that can distinguish between different parameterizations. This Letter establishes a first-principles benchmark for TMDPDFs, offering valuable input for ongoing and future experimental efforts to map the proton’s three-dimensional structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Three-dimensional higher-order saddle-point-induced flatbands in Co-based kagome metals

The saddle point (Van Hove singularity) exhibits a divergent density of states in two-dimensional systems, leading to fascinating phenomena such as strong correlations and unconventional superconductivity, yet it is seldom observed in three-dimensional (3D) systems. In this work we find two types of 3D higher-order saddle points (HOSPs) in emerging 3D kagome metals YbCo 6⁢ Ge 6 and MgCo 6 ⁢Ge 6 . Both HOSPs exhibit a singularity in their density of states, which is significantly enhanced compared to the ordinary saddle point. The HOSP near the Fermi energy generates a flatband extending a large area in the Brillouin zone, potentially amplifying the correlation effect and fostering electronic instabilities. Two types of HOSPs exhibit distinct robustness upon element substitution and lattice distortions in these kagome compounds. Our work paves the way for engineering exotic band structures, such as saddle points and flatbands, and exploring interesting phenomena in Co-based kagome materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Determining the Conformation of Supported Complexes Using an 17 O TEDOR-like NMR Experiment

Dynamic nuclear polarization surface-enhanced nuclear magnetic resonance (NMR) spectroscopy has enabled the determination of the three-dimensional configuration of surface sites, in particular supported metal complexes of relevance to single-site heterogeneous catalysis. These approaches have chiefly leveraged the application of NMR double-resonance experiments that either reveal the complex conformation via point-to-point intramolecular distances between spin-labeled atoms or the complex-surface orientation via distances between the spins and the surface plane. Either method typically requires expensive isotope labeling and each reports on different structural features. The application of an experiment that simultaneously reveals both types of distances with chemical resolution would be ideal. Here, in this article, we describe an 17 O{ 1 H} pseudo-3D correlation experiment that achieves this goal. Specifically, Si–O–Si and Si–O–M oxygens are well-resolved by 17 O NMR; therefore, distances can be simultaneously measured radially, between Si– 17 O–M and the 1 H’s of the ligands, and vertically to the Si– 17 O–Si linkages of the silica support. We demonstrate the experiment using supported yttrium and zirconium complexes. Good agreement is obtained when comparing the experimental results to theoretical predictions from density functional theory calculations, highlighting the reliability of this relatively simple experiment.

alkyls↗
DOE OSTI2025

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗
DOE OSTI2024

Effect of insulator end cap thickness on time-dependent Hartmann flow in a rotating mirror

We present a framework for analysing plasma flow in a rotating mirror. By making a series of physical assumptions, we reduce the magnetohydrodynamic (MHD) equations in a three-dimensional cylindrical system to a one-dimensional system in a shallow, cuboidal channel within a transverse magnetic field, similar to the Hartmann flow in ducts. We then solve the system both numerically and analytically for a range of values of the Hartmann number and calculate the dependence of the plasma flow speed on the thickness of the insulating end cap. We observe that the mean flow overshoots and decelerates before achieving a steady-state value, a phenomenon that the analytical model cannot capture. This overshoot is directly proportional to the thickness of the insulating end cap and the external electric field, with a weak dependence on the external magnetic field. Our simplified model can act as a benchmark for future simulations of the supersonic mirror device CMFX (centrifugal magnetic fusion experiment), which will employ more sophisticated physics and realistic magnetic field geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Pb:Sn Ratio-Driven Polytypism and Band Modulation in Photoresponsive Hexagonal Perovskitoids

Three-dimensional (3D) perovskites of the formula AMX 3 are known for their excellent optoelectronic properties, but their design is limited by the narrow range of A-site cations that can template the 3D corner-sharing structure, and many of the viable options have already been explored. These materials also face structural instability under environmental conditions. In contrast, 3D hexagonal perovskitoids, with the same chemical formula, may offer enhanced stability and richer structural diversity through a range of corner- and face-sharing octahedral configuration options, providing greater opportunities for structural design; however, challenges such as synthesis complexity and wide band gaps have hindered optoelectronic performance achieved to date. Herein, we synthesized a structural homologous series of mixed-metal hexagonal perovskitoids with the formula APb 1-x Sn x I 3 (x = 0, 0.25, 0.50, 0.75, 1; A = ethylammonium, guanidinium) and identified three polytypes (9R, 12R, 6H) using single-crystal X-ray diffraction (SCXRD). These structures exhibited increasing corner-sharing connectivity with higher Sn content, revealing a previously unobserved relationship between metal composition and structural evolution in perovskitoid materials. The incorporation of Sn reduced the band gap (tunable from 2.51 eV to 1.87 eV) and drove structural transformations, a trend seen also in density functional theory (DFT) calculations, which suggest a thermodynamic preference for Pb at face-sharing sites and Sn at corner-sharing sites. DFT band structure calculations and optical spectroscopy also reveal an anomalous behavior in these materials, which we term polytypic band modulation. This phenomenon combines conventional band bowing with the structural transformations that occur as Sn content increases, as observed in the band gap and in photoluminescence spectra. Photodiodes fabricated from thin films of these materials exhibited stable and pronounced photoresponses across various light intensities over time. Furthermore, the combination of templating 3D perovskitoids with multiple cations and alloying Pb and Sn, suggests a vastly underexplored phase space that offers new parameters to tune perovskitoids.

Cations↗
DOE OSTI2025

Dimensional Evolution Guides Property Control in the A n Cu 4– n TiS 4 Semiconductor Series

Through progressive reduction of the three-dimensional (3D) covalent network of Cu 4 TiS 4 , we isolate seven new members of the A n Cu 4–n TiS 4 family (A = alkali metal; n = 0–4), spanning 3D, 2D, 1D, and 0D structural fragments. The dimensional reduction is rational, as it preserves the edge-sharing connectivity between [CuS 4 ] 7– and [TiS 4 ] 4– tetrahedra across the series. This structural evolution is driven by the stepwise substitution of Cu with alkali metals, guiding the formation of fragments with reduced dimensionality. The effects of “n” and “A” on the crystal structures, stabilities, electronic structures, and optoelectronic properties are profound, demonstrating that the manipulation of alkali metal size and A n Cu 4–n TiS 4 stoichiometry enables predictable variations in structure and properties. For example, the n = 0 and n = 4 end members of the A n Cu 4–n TiS 4 family set the range of achievable band gaps with 2.00 eV for Cu 4 TiS 4 , 2.60 eV for Na 4 TiS 4 , and intermediate values for the n = 1–3 members. Notably, CsCu 3 TiS 4 exhibits exceptional air stability and congruent melting, with density functional theory (DFT) calculating moderate hole and electron effective masses in specific crystallographic directions (mh = 1.24m 0 , me = 0.87m 0 ). Additionally, A 3 CuTiS 4 (A = Na, K, Rb) displays direct band gap behavior and long photoluminescence lifetimes of 2.3–8.6 μs, and K 3 CuTiS 4 has a PLQY of 5.19%. These findings underscore the potential of the A n Cu 4–n TiS 4 family for applications in optoelectronics and demonstrate widely applicable design concepts that unveil rational stoichiometries within a given composition space to generate a series of crystal structures related through an evolving covalent dimensionality that corresponds to a predictable electronic structure and property progression.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Single nuclear spin detection and control in a van der Waals material

Optically active spin defects in solids are leading candidates for quantum sensing and quantum networking. Recently, single spin defects were discovered in hexagonal boron nitride (hBN), a layered van der Waals (vdW) material. Owing to its two-dimensional structure, hBN allows spin defects to be positioned closer to target samples than in three-dimensional crystals, making it ideal for atomic-scale quantum sensing, including nuclear magnetic resonance (NMR) of single molecules. However, the chemical structures of these defects remain unknown and detecting a single nuclear spin with a hBN spin defect has been elusive. Here we report the creation of single spin defects in hBN using 13 C ion implantation and the identification of three distinct defect types based on hyperfine interactions. We observed both S = 1/2 and S = 1 spin states within a single hBN spin defect. We demonstrated atomic-scale NMR and coherent control of individual nuclear spins in a vdW material, with a π-gate fidelity up to 99.75% at room temperature. By comparing experimental results with density functional theory (DFT) calculations, we propose chemical structures for these spin defects. Our work advances the understanding of single spin defects in hBN and provides a pathway to enhance quantum sensing using hBN spin defects with nuclear spins as quantum memories.

Quantum metrology↗
DOE OSTI2025

Quantum Spin Ice in Three-Dimensional Rydberg Atom Arrays

Quantum spin liquids are exotic phases of matter whose low-energy physics is described as the deconfined phase of an emergent gauge theory. With recent theory proposals and an experiment showing preliminary signs of Z 2 topological order [G. Semeghini , ], Rydberg atom arrays have emerged as a promising platform to realize a quantum spin liquid. In this work, we propose a way to realize a U(1) quantum spin liquid in three spatial dimensions, described by the deconfined phase of U(1) gauge theory in a pyrochlore lattice Rydberg atom array. We study the ground state phase diagram of the proposed Rydberg system as a function of experimentally relevant parameters. Within our calculation, we find that by tuning the Rabi frequency, one can access both the confinement-deconfinement transition driven by a proliferation of “magnetic” monopoles and the Higgs transition driven by a proliferation of “electric” charges of the emergent gauge theory. We suggest experimental probes for distinguishing the deconfined phase from ordered phases. This work serves as a proposal to access a confinement-deconfinement transition in three spatial dimensions on a Rydberg-based quantum simulator. Published by the American Physical Society 2025

Shah, Jeet (ORCID:0000000158738129)↗
NASA NTRS1971

Contribution to methods for calculating the flow about thin lifting wings at transonic speeds: Analytic expressions for the far field

The problem of determining the small-disturbance flow about two-dimensional airfoils at transonic speeds has been successfully treated by the process of matching a numerical solution of the near field to analytic expressions for the far field. The three-dimensional problem, it would appear, can be treated in a similar way with the aid of algorithms adapted to high-speed and high-capacity computers. The far-field potential for both lifting and nonlifting three-dimensional wings at transonic speeds is developed herein for a subsonic free stream. This potential could be used for a three-dimensional-wing computation similar to the computation made for the two-dimensional wing.

Klunker, E. B.↗
NASA NTRS1974

Three-dimensional crack problems

Methods for solving three dimensional crack problems are discussed. The stress intensity factor for the case of uniform pressure on the surface crack is analysis. An analysis of concentrated forces acting on the crack faces parallel to the crack front of a half-plane crack was made. Two dimensional calculations were made using the exact solution for the stress field around an elliptical hole in a uniformly stretched plate. The strain energy density was computed in the neighborhood of the ends of narrow, but sharp, ellipses.

Hartranft, R. J.↗
NASA NTRS1974

Comparison of approximate and numerical analyses of nonlinear combustion instability

At the present time, there are three general analytical techniques available to study problems of unsteady motions in rocket motors: linear stability analysis; approximate nonlinear analysis, founded on examining the behavior of coupled normal modes; and numerical calculations based on the conservation equations for one-dimensional flows. The last two yield the linear results as a limit. It is the main purpose of this paper to check the accuracy of the approximate analysis against the numerical analysis for some special cases. The results provide some justification for using the approximate analysis to study three-dimensional problems.

Culick, F. E. C.↗
NASA NTRS1975

Particle motion in atmospheric boundary layers of Mars and Earth

To study the eolian mechanics of saltating particles, both an experimental investigation of the flow field around a model crater in an atmospheric boundary layer wind tunnel and numerical solutions of the two- and three-dimensional equations of motion of a single particle under the influence of a turbulent boundary layer were conducted. Two-dimensional particle motion was calculated for flow near the surfaces of both Earth and Mars. For the case of Earth both a turbulent boundary layer with a viscous sublayer and one without were calculated. For the case of Mars it was only necessary to calculate turbulent boundary layer flow with a laminar sublayer because of the low values of friction Reynolds number; however, it was necessary to include the effects of slip flow on a particle caused by the rarefied Martian atmosphere. In the equations of motion the lift force functions were developed to act on a single particle only in the laminar sublayer or a corresponding small region of high shear near the surface for a fully turbulent boundary layer. The lift force functions were developed from the analytical work by Saffman concerning the lift force acting on a particle in simple shear flow.

White, B. R.↗
NASA NTRS1975

Further numerical calculations of the circulation of the atmosphere of Venus

Results of two-dimensional simulations of the deep circulation of Venus are presented. They prove that the high surface temperature can only be explained by the greenhouse effect, and that Goody and Robinson's (1966) dynamical model is not valid. Analytical radiative equilibrium solutions for a semi-grey atmosphere, both with and without an internal heat source, are presented. It is shown that the greenhouse effect is sufficient to produce the high surface temperature if the thermal optical depth is much greater than 100 and the ratio of solar to thermal optical depths does not exceed 0.005. A two-dimensional version of a three-dimensional model is used to test the validity of the new mechanism proposed by Gierasch (1975) to explain the 4-day circulation. Numerical experiments with large horizontal viscosities failed to show strong zonal velocities even for the case of large Prandtl numbers. It is observed that the dissipation of angular momentum introduced by the strong horizontal diffusion more than compensates for the upward transport of angular momentum due to the Hadley cell.

De Rivas, E. K.↗
NASA NTRS1974

Computation of steady axisymmetric flow using a one-dimensional time-dependent method

An iterative numerical method for computing steady, three dimensional, viscous, compressible flow fields, about aerodynamic bodies was studied. In order to develop the iterative method economically, the primary emphasis was directed towards supersonic, axisymmetric flow. However, the technique readily extends to three spatial dimensions. The viscous flow field about a cone-cylinder-flare body was calculated and compared to existing experimental data. Numerical predictions of the cone boundary layer and the flow field shock wave structure agreed with corresponding measurements. A separation was calculated at the cylinder-flare junction in six iterations; however, the size of the vortex did not correspond to the measured size. It was estimated that fifty iterations would be required to properly define the vortex. It was concluded that the iteration technique is of limited value for plane two dimensional and axisymmetrix flows, but of great value for three-dimensional flows.

Walitt, L.↗