Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Physics”

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

Computational framework for modeling of physical process

Techniques, systems, and devices are described for providing a computational frame for estimating high-dimensional stochastic behaviors. In one exemplary aspect, a method for performing numerical estimation includes receiving a set of measurements of a stochastic behavior. The set of correlated measurements follows a non-standard probability distribution and is non-linearly correlated. Also, a non-linear relationship exists between a set of system variables that describes the stochastic behavior and a corresponding set of measurements. The method includes determining, based on the set of measurements, a numerical model of the stochastic behavior. The numerical model comprises a feature space comprising non-correlated features corresponding to the stochastic behavior. The non-correlated features have a dimensionality of M and the set of measurements has a dimensionality of N, M being smaller than N. The method includes generating a set of approximated system variables corresponding to the set of measurements based on the numerical model.

Chen, Xiao↗
DOE OSTI2024

Computational framework for modeling of physical process

Techniques, systems, and devices are described for providing a computational frame for estimating high-dimensional stochastic behaviors. In one exemplary aspect, a method for performing numerical estimation includes receiving a set of measurements of a stochastic behavior. The set of correlated measurements follows a non-standard probability distribution and is non-linearly correlated. Also, a non-linear relationship exists between a set of system variables that describes the stochastic behavior and a corresponding set of measurements. The method includes determining, based on the set of measurements, a numerical model of the stochastic behavior. The numerical model comprises a feature space comprising non-correlated features corresponding to the stochastic behavior. The non-correlated features have a dimensionality of M and the set of measurements has a dimensionality of N, M being smaller than N. The method includes generating a set of approximated system variables corresponding to the set of measurements based on the numerical model.

Chen, Xiao↗
DOE OSTI2024

Chiral limit of 2d QCD revisited with lightcone conformal truncation

We study the chiral limit of 2d QCD with a single quark flavor at finite Nc using LCT. By modifying the LCT basis according to the quark mass in a manner motivated by ’t Hooft’s analysis, we are able to restore convergence for quark masses much smaller than the QCD strong coupling scale. For such small quark masses, the IR of the theory is expected to be well described by the Sine-Gordon model. We verify that LCT numerics are able to capture in detail the spectrum and correlation functions of the Sine-Gordon model. This opens up the possibility for studying deformations of various integrable CFTs using LCT by considering the chiral limit of QCD like theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Emergent geometry and duality in the carbon nucleus

Abstract The carbon atom provides the backbone for the complex organic chemistry composing the building blocks of life. The physics of the carbon nucleus in its predominant isotope, 12 C, is similarly full of multifaceted complexity. Here we provide a model-independent density map of the geometry of the nuclear states of 12 C using the ab initio framework of nuclear lattice effective field theory. We find that the well-known but enigmatic Hoyle state is composed of a “bent-arm” or obtuse triangular arrangement of alpha clusters. We identify all of the low-lying nuclear states of 12 C as having an intrinsic shape composed of three alpha clusters forming either an equilateral triangle or an obtuse triangle. The states with the equilateral triangle formation also have a dual description in terms of particle-hole excitations in the mean-field picture.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Understanding Peelle’s Pertinent Puzzle bias in generalized least squares regression through eigenspectrum analysis

Certain correlation structures in the data covariance matrix (DCM) used for generalized least squares (GLS) regression can result in biased estimates, commonly known in the field of nuclear data evaluation as Peele’s Pertinent Puzzle (PPP). This article introduces a generative, forward modeling framework within which the PPP bias is characterized through an eigenspectrum analysis of the DCM. This analysis highlights the root cause of the bias, generalizes the problem beyond the nuclear data field, and provides insight to the problem regimes where it can occur. What follows is an understanding that the bias can show up for any experimental neutron time-of-flight data for which systematic uncertainties have been quantified. Lastly, a discussion of the adaptation of cross validation approaches that require pre-whitening to incorporate the known ‘fix’ to the PPP bias in the GLS estimator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Target normal sheath acceleration with a large laser focal diameter

In this work, the dependence of the laser-driven ion acceleration from thin titanium foils in the Target Normal Sheath Acceleration (TNSA) regime on target and laser parameters is explored using two dimensional particle-in-cell simulations. The oblique incidence (θ L = 45°) and large focal spot size (w 0 = 40μm) are chosen to take an advantage of quasi one-dimensional geometry of sheath fields and effective electron heating. This interaction setup also reveals low and achromatic angular divergence of a proton beam. It is shown that the hot electron temperature deviates from the ponderomotive scaling for short laser pulses and small pre-plasmas. This deviation is mainly due to the laser sweeping, as the short duration laser pulse each moment in time effectively heats only a fraction of a focal spot on the foil. This instantaneous partial heating results in an electron temperature deviation from the ponderomotive scaling and, thus, lower maximum proton energies than it could have been expected from the TNSA theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Finding the order in complexity: The electronic structure of 14-1-11 zintl compounds

Yb 14 MnSb 11 and Yb 14 MgSb 11 have rapidly risen to prominence as high-performing p-type thermoelectric materials. However, the fairly complex crystal structure of A 14 MX 11 Zintl compounds renders the interpretation of the electronic band structure obscure, making it difficult to chemically guide band engineering and optimization efforts. In this work, we delineate the valence-balanced Zintl chemistry of A 14 MX 11 compounds using the molecular orbital theory. By analyzing the electronic band structures of Yb 14 MgSb 11 and Yb 14 AlSb 11 , we show that the conduction band minimum is composed of either an antibonding molecular orbital originating from the (Sb 3 ) 7– trimer or a mix of atomic orbitals of A, M, and X. The singly degenerate valence band is comprised of non-bonding Sb pz orbitals primarily from the Sb atoms in the (MSb 4 ) m– tetrahedra and of isolated Sb atoms distributed throughout the unit cell. Such a chemical understanding of the electronic structure enables strategies to engineer electronic properties (e.g., the bandgap) of A 14 MX 11 compounds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Excitation of long electrostatic solitary waves in ion beam neutralization process

Unusually long electrostatic solitary waves (ESWs) are discovered in a particle-in-cell simulation study of the process of ion beam neutralization by electron emission from a filament. These ESWs are long because the density perturbation responses to the potential wells created by the ESWs are very small. Notably, the density perturbation is small because the trapped (positive) and untrapped (negative) electron density perturbations nearly compensate each other because of a non-Maxwellian electron velocity distribution in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Optimal Protocols in Quantum Annealing and Quantum Approximate Optimization Algorithm Problems

Quantum annealing (QA) and the quantum approximate optimization algorithm (QAOA) are two special cases of the following control problem: apply a combination of two Hamiltonians to minimize the energy of a quantum state. Which is more effective has remained unclear. Here we analytically apply the framework of optimal control theory to show that generically, given a fixed amount of time, the optimal procedure has the pulsed (or “bang-bang”) structure of QAOA at the beginning and end but can have a smooth annealing structure in between. This is in contrast to previous works which have suggested that bang-bang (i.e., QAOA) protocols are ideal. To support this theoretical work, we carry out simulations of various transverse field Ising models, demonstrating that bang-anneal-bang protocols are more common. Futher, the general features identified here provide guideposts for the nascent experimental implementations of quantum optimization algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Emerging applications: Neuromorphic computing and reservoir computing

The emergence of doped hafnium oxide (HfO 2 )-based ferroelectric films has enabled highly scalable and silicon-compatible ferroelectric devices, opening new frontiers in neuromorphic and reservoir computing. Among these, ferroelectric field-effect transistors (FeFETs) are particularly promising due to their analog memory characteristics and unique polarization dynamics. These properties make FeFETs ideal candidates for artificial synapses in neuromorphic architectures, supporting deep neural networks and spiking neural networks based on leaky-integrate-and-fire (LIF) mechanisms. Beyond neuromorphic computing, FeFETs also play a crucial role in physical reservoir computing, leveraging their intrinsic nonlinear and history-dependent behavior for efficient real-time learning. This approach offers significant advantages for time-series processing and edge artificial intelligence (AI) applications, addressing the growing need for energy-efficient computing. As a result, this article explores the principles, key demonstrations, and future potential of FeFET-based neuromorphic and reservoir computing, highlighting their impact on next-generation AI hardware.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimization of Geometric Perturbations on a Rod Moving Through a High Explosive Target

After completing a study to ensure the simulation results were converged, several high resolution 3D Smoothed Particle Hydrodynamic (SPH) simulations of copper rods impacting a high explosive (LX14) target were performed. This was then formulated into an optimization problem: I wanted to find the optimum shape and location of a perturbation on the rod that would maximize its erosion after it left the target. The shape of the perturbation was modeled as a 2D Gaussian bump and parameterized by its location along the rod axis (z 0 ) and amplitude (A). The final mass of the coherent part of the rod as it leaves the target was used as a metric to represent the erosion of the rod, and the optimization was formulated to maximize this metric with respect to the aforementioned design variables. Due to the expensive nature of the high-fidelity 3D SPH simulations, a surrogate model needed to be chosen so that many function calls to the optimizer would be feasible. Thus, a strategic full factorial sampling plan was chosen to build a dataset, which consisted of 24 high-fidelity simulations. Two surrogate models, a third order polynomial regression model and a Gaussian Process Model, were analyzed using a 14%/86% test/train holdout technique. The root mean square and R2 score of the test set was used to determine the best model, and the third order polynomial regression model was chosen as the surrogate model. Finally, the Nelder-Mead Simplex and Basin-hopping optimization algorithms were implemented, and it was found that the two algorithms gave slightly different optimum values. Nelder-Mead gave an optimum point of [z* 0 ;A*] = [9:9;0:4] and Basin-Hopping gave an optimum value of x* = [z* 0 ;A*] = [9:2;0:1].

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗