Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical solution”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamics of electrified liquid metal surface using shallow water model

A shallow water model that incorporates surface tension and electric field effects is developed to investigate the dynamics of an electrified liquid surface. The computational model is verified against the Zakharov–Kuznetsov equation and is applied to study the growth and damping of the electrified liquid surface. A linear wave analysis is performed under a shallow water theory assuming an analytic solution of the electric field, similar to the Tonks–Frenkel instability. The electrified liquid surface grows or dampens based on the balance of the electric field, surface tension, and gravitational forces. As a result, the numerical results obtained from the electrified shallow water solver are in good agreement with the theoretical analysis.

Computational models↗
DOE OSTI2023

Steady-state junction current distribution in p-n GaN diodes measured using low-energy electron microscopy (LEEM)

Here, we report on the measurement of the lateral distribution of the junction current of an electrical biased p-n GaN diode by electron emission microscopy using a low-energy electron microscope. The vacuum level at the surface of the diode was lowered by deposition of cesium to achieve negative electron affinity, allowing overflow electrons at the surface of the biased diodes to be emitted and their spatial distribution imaged. The results were compared to the literature, and a good match with analytical solutions by Joyce and Wemple was obtained.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Thickness dependence of piezo-bimorph adaptive mirror bending

A new generation of adaptive x-ray optics (AXO) is being installed on high-coherent-flux x-ray beamlines worldwide to correct and control the optical wavefront with sub-nm precision. These ultra-smooth mirrors achieve high reflectivities at glancing angles of incidence and can be hundreds of mm long. One type of adaptive x-ray mirror relies on piezoelectric ceramic strips which are segmented into channels and actuated to induce local, longitudinal bending, generating one-dimensional shape changes in the mirror substrate. A recently described mirror model uses a three-layer geometry with parallel actuators on the front and back surfaces of a thicker mirror substrate. By analogy to a solved problem in the thermal actuation of a tri-metal strip, we show that the achievable bending radius varies approximately as the square of the substrate thickness. We provide an analytic solution and simulate bending using a finite-element model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Investigating the origin of the far-field reflection interference fringe (RIF) of microdroplets

We show that the reflection interference fringe (RIF) is formed on a screen far away from the microdroplets placed on a prism-based substrate, which have low contact angles and thin droplet heights, caused by the dual convex–concave profile of the droplet, not a pure convex profile. The geometric formulation shows that the interference fringes are caused by the optical path difference when the reflected rays from the upper convex profile at the droplet–air interface interfere with reflection from the lower concave profile at oblique angles lower than the critical angle. Analytic solutions are obtained for the droplet height and the contact angle out of the fringe number and the fringe radius in RIF from the geometric formulation. Furthermore, the ray tracing simulation is conducted using the custom-designed code. The geometric formulation and the ray tracing show excellent agreement with the experimental observation in the relation between the droplet height and the fringe number and the relation between the contact angle and the fringe radius. This study is remarkable as the droplet's dual profile cannot be easily observed with the existing techniques. However, the RIF technique can effectively verify the existence of a dual profile of the microdroplets in a simple setup. In this work, the RIF technique is successfully developed as a new optical diagnostic technique to determine the microdroplet features, such as the dual profile, the height, the contact angle, the inflection point, and the precursor film thickness, by simply measuring the RIF patterns on the far-field screen.

42 ENGINEERING↗
DOE OSTI2025

FreeMHD: Validation and verification of the open-source, multi-domain, multi-phase solver for electrically conductive flows

The extreme heat fluxes in the divertor region of tokamaks may require an alternative to solid plasma-facing components, for the extraction of heat and the protection of the surrounding walls. Flowing liquid metals are proposed as an alternative, but raise additional challenges that require investigation and numerical simulations. Free surface designs are desirable for plasma-facing components, but steady flow profiles and surface stability must be ensured to limit undesirable interactions with the plasma. Previous studies have mainly used steady-state, 2D, or simplified models for internal flows and have not been able to adequately model free-surface liquid metal (LM) experiments. Therefore, FreeMHD has been recently developed as an open-source magnetohydrodynamics (MHD) solver for free-surface electrically conductive flows subject to a strong external magnetic field. The FreeMHD solver computes incompressible free-surface flows with multi-region coupling for the investigation of MHD phenomena involving fluid and solid domains. The model utilizes the finite-volume OpenFOAM framework under the low magnetic Reynolds number approximation. FreeMHD is validated using analytical solutions for the velocity profiles of closed channel flows with various Hartmann numbers and wall conductance ratios. Next, experimental measurements are then used to verify FreeMHD, through a series of cases involving dam breaking, 3D magnetic fields, and free-surface LM flows. These results demonstrate that FreeMHD is a reliable tool for the design of LM systems under free surface conditions at the reactor scale. Furthermore, it is flexible, computationally inexpensive, and can be used to solve fully 3D transient MHD flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Guiding center equations for the magnetic dipole

Since the discovery of Van Allen radiation belts in the 1960s, observations of energetic ions trapped in the Earth's dipole magnetic field have illustrated the remarkable confinement properties of this configuration. As such, it has been used for confining a hot plasma for nuclear fusion studies, starting from the pioneering work of Bo Lehnert and Akira Hasegawa, in the Levitated Dipole Experiment (LDX) at MIT until 2011 and in the RT-1 experiment at the University of Tokyo. More recently, the dipole has been subject to a renewed interest for fusion studies by a couple of startups and for smaller applications as a cold plasma source. While the equilibrium and magneto-hydrodynamic stability of the dipole have been investigated quite in detail, neoclassical properties of the dipole are comparatively much less known: the dipole is more known in geophysics than in fusion science. For this reason, in this paper, we propose a set of Hamiltonian, guiding-center equations to describe the motion of electrons and ions in a magnetic dipole configuration. We also developed a code, and we show the main features of particle motion, benchmarking our results with the analytical solutions for the bounce and precession motion, which are well documented in the literature. We also draw some general conclusions for the neoclassical transport in usual toroidal confinement schemes, such as the tokamak and the stellarator, pointing out the unique advantages of the dipole in confining energetic particles.

Hamiltonian mechanics↗
DOE OSTI2025

Simple analytic fusion hot spot models for fusion reaction history

The measured fusion reaction history is a combination of the temporal evolution of the fusion hot spot temperature, mass, and volume. Depending on the mechanism of evolution in inertial confinement fusion implosions—shocks, compression, convergence, mass ablation, ignition—the evolution of the reaction history varies. Here, we derive and catalog a set of simplified inertial confinement fusion hot spot models with analytic solutions to infer the evolution of the fusion reaction history for each mechanism. The models give valuable insight into the meaning and cause of fusion reaction history and nuclear burnwidth measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗
DOE OSTI2025

Deconstructing dynamics of symmetry breaking

The Kibble–Zurek mechanism (KZM) successfully predicts the density of topological defects deposited by the phase transitions, but it is not clear why. Its key conjecture is that, near the critical point of the second-order phase transition, critical slowing down will result in a period when the system is too sluggish to follow the potential that is changing faster than its reaction time. The correlation length at the freeze-out instant $\hat{t}$ when the order parameter catches up with the posttransition broken symmetry configuration is then decisive, determining when the mosaic of broken symmetry domains locks in topological defects. To understand why the KZM works so well, we analyze the Landau–Ginzburg model and show why temporal evolution of the order parameter plays such a key role. In conclusion, the analytical solutions we obtain suggest experimentally accessible observables that can shed light on symmetry-breaking dynamics while testing the conjecture on which the KZM is based.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sockeye: A One-Dimensional, Two-Phase, Compressible Flow Heat Pipe Application

Sockeye is a heat pipe analysis application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework. The primary purpose of Sockeye is to provide a transient heat pipe simulation tool to be used in the analysis of nuclear microreactor designs. Sockeye provides the capability to perform one-dimensional, two-phase, compressible flow simulation of a heat pipe working fluid and two-dimensional, axisymmetric heat conduction for the heat pipe cladding and its surroundings. Sockeye is demonstrated against analytical solutions and experimental data from the SAFE-30 heat pipe module test.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

A Characteristics Approach to the Finite Element Method

Herein, we present a new method for solving the linear Boltzmann transport equation. Two commonly used and well-understood methods for solving partial differential equations are the method of characteristics (MOC) and the finite element method (FEM). We propose a new method that combines the fundamental concept of the FEM with the analytic solution from the MOC to obtain coefficients for the FEM basis function expansion. Traditionally, coefficients for the FEM basis function expansion are obtained via matrix inversion. Instead, we solve for the coefficients with the MOC and represent the underlying fields with the basis function expansion using these coefficients. We provide a convergence study for our method with results from two sets of FEM basis functions: Gauss-Legendre and Gauss-Lobatto sets. We also compare two different variations of our method categorized as short characteristics and intermediate characteristics.

42 ENGINEERING↗
DOE OSTI2026

A Monte Carlo Laplace Transform Estimator for Radiation Transport

This work formulates and implements a Laplace transform estimator in a simple Monte Carlo radiation transport code. The estimator maps flux-based quantities of interest, like reaction rates, from a desired phase-space dimension to the complex Laplace domain. This on-the-fly Monte Carlo integration technique enables the spectral analysis of arbitrary nuclear systems via the Laplace transform. A simple code tests the estimator in neutron slowing-down problems across various infinite media, and the results compare well with Ganapol’s uninverted analytical solution of the neutron slowing-down equation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Clarifying the definition of ‘transonic’ screw dislocations

We report a number of recent Molecular Dynamics (MD) simulations have demonstrated that screw dislocations in face centred cubic (fcc) metals can achieve stable steady state motion above the lowest shear wave speed (v shear ) which is parallel to their direction of motion (often referred to as transonic motion). This is in direct contrast to classical continuum analyses which predict a divergence in the elastic energy of the host material at a crystal geometry dependent ‘critical’ velocity v crit . Within this work, we first demonstrate through analytic analyses that the elastic energy of the host material diverges at a dislocation velocity (v crit ) which is greater than v shear , i.e. v crit >v shear . We argue that it is this latter derived velocity (v crit ) which separates ‘subsonic’ and ‘supersonic’ regimes of dislocation motion in the analytic solution. In addition to our analyses, we also present a comprehensive suite of MD simulation results of steady state screw dislocation motion for a range of stresses and several cubic metals at both cryogenic and room temperatures. At room temperature, both our independent MD simulations and the earlier works find stable screw dislocation motion only below our derived v crit . Nonetheless, in real-world polycrystalline materials v crit cannot be interpreted as a hard limit for subsonic dislocation motion. In fact, at very low temperatures our MD simulations of Cu at 10 Kelvin confirm a recent claim in the literature that true ‘supersonic’ screw dislocations with dislocation velocities v>v crit are possible at very low temperatures.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Gravity of two photon decay and its quantum coherence

A linear analytical solution is derived for the gravitational shock wave produced by a particle of mass M that decays into a pair of null particles. The resulting space-time is shown to be unperturbed and isotropic, except for a discontinuous perturbation on a spherical null shell. Formulae are derived for the perturbation as a function of polar angle, as measured by an observer at the origin observing clocks on a sphere at distance R. The effect of the shock is interpreted physically as an instantaneous displacement in time and velocity when the shock passes the clocks. The time displacement is shown to be anisotropic, dominated by a quadrupole harmonic aligned with the particle-decay axis, with a magnitude δτ ~ GM/c 3 , independent of R. The velocity displacement is isotropic. The solution is used to derive the gravitational effect of a quantum state with a superposition of a large number of randomly oriented, statistically isotropic particle decays. This approach is shown to provide a well-controlled approximation to estimate the magnitude of gravitational fluctuations in systems composed of null point particles up to the Planck energy in a causal diamond of duration τ = 2R/c, as well as quantum-gravitational fluctuations of black holes and cosmological horizons. Coherent large-angle quantum distortions of macroscopic geometry from fluctuations up to the Planck scale are shown to grow linearly with the duration, with a variance $\langle$δτ 2 $\rangle$ ~ τtP much larger than that produced in models without causal quantum coherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Development of coated heat flux gauges for fast responding measurements

Thermal systems often exhibit transient behaviors that have important implications for the operation of the system and can be difficult to predict. For these reasons, experimental testing is often required to ensure system durability requirements are achieved. One important parameter governing the survivability of components in hot, high-stress environments is the heat flux into the part that dictates the temperature distribution for the component. However, sensors required to experimentally characterize heat fluxes in extreme environments must also be resilient. Here, this study presents the development of coated heat transfer gauges capable of robust, high-frequency measurements in turbine research facilities. The addition of a protective coating increases the durability of the gauge, but inherent of that coating is the attenuation of high-frequency temperature penetrations. As a result, this study first outlines the use of analytical solutions to define a gauge design for a specific frequency range and heat transfer, ensuring that subsurface signals can be rectified to surface conditions through inverse methods. Then, the fabrication of polyimide substrate sensors with a parylene-F coating is described. Micro surface heaters added to the custom sensors were used to determine important geometric and thermal properties necessary to calculate accurate surface heat flux. Ultimately, this work shows increased sensor robustness in a turbine test bed and experimentally validates that the frequency response of the fabricated sensors meet the design intent.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Magnetic tunnel junction random number generators applied to dynamically tuned probability trees driven by spin orbit torque

Abstract Perpendicular magnetic tunnel junction (pMTJ)-based true-random number generators (RNGs) can consume orders of magnitude less energy per bit than CMOS pseudo-RNGs. Here, we numerically investigate with a macrospin Landau–Lifshitz-Gilbert equation solver the use of pMTJs driven by spin–orbit torque to directly sample numbers from arbitrary probability distributions with the help of a tunable probability tree. The tree operates by dynamically biasing sequences of pMTJ relaxation events, called ‘coinflips’, via an additional applied spin-transfer-torque current. Specifically, using a single, ideal pMTJ device we successfully draw integer samples on the interval [0, 255] from an exponential distribution based on p -value distribution analysis. In order to investigate device-to-device variations, the thermal stability of the pMTJs are varied based on manufactured device data. It is found that while repeatedly using a varied device inhibits ability to recover the probability distribution, the device variations average out when considering the entire set of devices as a ‘bucket’ to agnostically draw random numbers from. Further, it is noted that the device variations most significantly impact the highest level of the probability tree, with diminishing errors at lower levels. The devices are then used to draw both uniformly and exponentially distributed numbers for the Monte Carlo computation of a problem from particle transport, showing excellent data fit with the analytical solution. Finally, the devices are benchmarked against CMOS and memristor RNGs, showing faster bit generation and significantly lower energy use.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators

Abstract The stability of the ion-temperature gradient mode in quasisymmetric stellarators is assessed. This is performed using a set of analytical estimates together with linear gyrokinetic simulations. The peak growth rates, their corresponding real frequencies and wave-vectors are identified. A comparison is made between a first-order near-axis expansion model and eleven realistic designs obtained using numerical optimization methods. It is found that while the near-axis expansion is able to replicate the growth rates, real frequencies and perpendicular wave-vector at the inner core (both using simplified dispersion relations and first-principle gyrokinetic simulations), it leads to an overestimation of the growth rate at larger radii. An approximate analytic solution of the ITG dispersion relation for the non-resonant limit suggests growth rates could be systematically higher in quasi-axisymmetric (QA) configurations compared to quasi-helically (QH) symmetric ones. However except for very close to the axis, linear gyrokinetic simulations do not show systematic differences between QA and QH configurations.

Physics↗
DOE OSTI2023

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗