Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Near-axis expansion of stellarator equilibrium at arbitrary order in the distance to the axis

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight into the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect ratio of the ideal magnetohydrodynamics equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet–Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren–Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest-order solution and the W7-X magnetic field is performed.

Physics↗
DOE OSTI2019

From the Ternary Phase Ca 14 Zn 1+δ Sb 11 (δ ≈ 0.4) to the Quaternary Solid Solutions Ca 14– x RE x ZnSb 11 (RE = La–Nd, Sm, Gd, x ≈ 0.9). A Tale of Electron Doping via Rare-Earth Metal Substitutions and the Concomitant Structural Transformations

In this work, the ternary compound Ca 14 Zn 1.37(1) Sb 11 and its six rare-earth metal substituted derivatives Ca 14–x RE x ZnSb 11 (RE = La–Nd, Sm, Gd; x ≈ 0.90±0.06) have been synthesized and structurally characterized by single-crystal X-ray diffraction methods. All com-pounds formally crystallize in the tetragonal Ca 14 AlSb 11 structure type (space group $I$4 1 $/ acd$, No. 142, $Z$ = 8). The crystal structure of Ca 14 Zn 1.37(1) Sb 11 subtly differs from the structure of the remaining six, as well as from the structure of the archetype, due to the presence of a partially occupied interstitial Zn position. The extra zinc atom is needed in this structure to alleviate the unfavorable number of valence electrons in the imaginary Ca 14 ZnSb 11 . Electron doping, via substitution of RE 3+ ions on Ca 2+ sites is shown as an alternative route to achieve electron balance in these Zn-based analogs of the Ca 14 AlSb 11 structure, which does not require the in-corporation of interstitial atoms. Electrical resistivity measurements done on single-crystalline samples are in agreement with the notion that Ca 14–x RE x ZnSb 11 behave as either bad metals or heavily-doped semi-conductors. Magnetization measurements show Curie-Weiss paramagnetic behavior related to the local-moment magnetism of the RE 3+ ions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

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

Distinguishing cause from correlation in tokamak experiments to trigger edge-localised plasma instabilities

The generic question is considered: How can we determine the probability of an otherwise quasi-random event, having been triggered by an external influence? A specific problem is the quantification of the success of techniques to trigger, and hence control, edge-localised plasma instabilities (ELMs) in magnetically confined fusion (MCF) experiments. The development of such techniques is essential to ensure tolerable heat loads on components in large MCF fusion devices, and is necessary for their development into economically successful power plants. Bayesian probability theory is used to rigorously formulate the problem and to provide a formal solution. Accurate but pragmatic methods are developed to estimate triggering probabilities, and are illustrated with experimental data. These allow results from experiments to be quantitatively assessed, and rigorously quantified conclusions to be formed. Example applications include assessing whether triggering of ELMs is a statistical or deterministic process, and the establishment of thresholds to ensure that ELMs are reliably triggered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Confinement effects of magnetic field on two-dimensional hydrogen atom in plasmas

In this study, for the first time, the Schrödinger equation with more general exponential cosine screened Coulomb (MGECSC) potential is solved numerically in the presence and in the absence of an external magnetic field within two-dimensional formalism using the asymptotic iteration method. The MGECSC potential includes four different potential forms when considering different sets of the parameters in the potential. The plasma screening effects in the weak and strong magnetic field regimes as well as the confinement effects of magnetic field on the two-dimensional hydrogen atom in Debye and quantum plasmas are investigated by solving the corresponding equations. It is found that applying a uniform magnetic field on the hydrogen atom embedded in a plasma leads to change in the profile of the total interaction potential. Thus, confinement effects of magnetic field on hydrogen atom embedded in Debye and quantum plasmas modeled by a MGECSC potential lead to shift bound state energies. This effect would be important to isolate the plasma from the external environment in the experimental applications of plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Full-zone spectral envelope function formalism for the optimization of line and point tunnel field-effect transistors

Efficient quantum mechanical simulation of tunnel field-effect transistors (TFETs) is indispensable to allow for an optimal configuration identification. We therefore present a full-zone 15-band quantum mechanical solver based on the envelope function formalism and employing a spectral method to reduce computational complexity and handle spurious solutions. We demonstrate the versatility of the solver by simulating a 40 nm wide In{sub 0.53}Ga{sub 0.47}As lineTFET and comparing it to p-n-i-n configurations with various pocket and body thicknesses. We find that the lineTFET performance is not degraded compared to semi-classical simulations. Furthermore, we show that a suitably optimized p-n-i-n TFET can obtain similar performance to the lineTFET.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum algorithm for polaritonic chemistry based on an exact ansatz

Abstract Cavity-modified chemistry uses strong light-matter interactions to modify the electronic properties of molecules in order to enable new physical phenomena such as novel reaction pathways. As cavity chemistry often involves critical regions where configurations become nearly degenerate, the ability to treat multireference problems is crucial to understanding polaritonic systems. In this Letter, we show through the use of a unitary ansatz derived from the anti-Hermitian contracted Schrödinger equation that cavity-modified systems with strong correlation, such as the deformation of rectangular H 4 coupled to a cavity mode, can be solved efficiently and accurately on a quantum device. In contrast, while our quantum algorithm can be made formally exact, classical-computing methods as well as other quantum-computing algorithms often yield answers that are both quantitatively and qualitatively incorrect. Additionally, we demonstrate the current feasibility of the algorithm on near intermediate-scale quantum hardware by computing the dissociation curve of H 2 strongly coupled to a bosonic bath.

Warren, Samuel (ORCID:0000000157134454)↗
DOE OSTI2018

Rare behavior of growth processes via umbrella sampling of trajectories

We compute probability distributions of trajectory observables for reversible and irreversible growth processes. These results reveal a correspondence between reversible and irreversible processes, at particular points in parameter space, in terms of their typical and atypical trajectories. Thus key features of growth processes can be insensitive to the precise form of the rate constants used to generate them, recalling the insensitivity to microscopic details of certain equilibrium behavior. We obtained these results using a sampling method, inspired by the "s-ensemble" large-deviation formalism, that amounts to umbrella sampling in trajectory space. The method is a simple variant of existing approaches, and applies to ensembles of trajectories controlled by the total number of events. It can be used to determine large-deviation rate functions for trajectory observables in or out of equilibrium.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Effect of interfacial structures on phonon transport across atomically precise Si/Al heterojunctions

Phonons are important carriers of energy and information in many cryogenic devices used for quantum information science and in fundamental physics experiments such as dark matter detectors. In these systems phonon behavior can be dominated by interfaces, hence there is increasing demand for a more detailed understanding of interfacial phonon transport in relevant material systems. Previous studies have focused on understanding thermal transport over the entire phonon spectrum at and above room temperature. At ultra-low temperatures, however, knowledge is missing regarding athermal phonon behavior due to the challenge in modeling the extreme conditions in micro-scale, heterogeneous cryogenic systems, as well as extracting single-phonon information from a large ensemble. In this work, we delineate interfacial phonon transport using a combination of classical molecular dynamics (MD) and phonon wave packet simulations, to illustrate the consistency and differences between the ensemble- and single-phonon dynamics. We consider three single-crystal Si surface reconstructions --- (1×1), (v3×v3) and (7×7) --- and model both experimentally observed Si(1×1)/Al interfaces and hypothesized Si(v3×v3)/Al and Si(7×7)/Al interfaces. The overall interfacial thermal conductance calculated from non-equilibrium MD shows that for the Si(1×1)/Al system, the presence of Al twin grain boundaries can hinder phonon transport and reduce thermal conductance by 2--12% relative to single-crystal Al; whereas the Si (v3×v3) and (7×7) reconstructions can enhance it by 6--19%. Normal mode decomposition reveals that both the increase and decrease in conductance are related to inelastic phonon scattering. Single-phonon wave packet simulations show results consistent with NEMD, while further suggesting that phonon polarization conversion is significant even when elastic transmission dominates, and that the interfacial structures have anisotropic impacts on atomic vibrations along different lattice directions. Our findings suggest alternative perspectives in achieving selective phonon transport through controlling interfacial structures of materials with atomically-precise fabrication techniques, and that the phonon wave packet formalism is a potentially powerful method for developing a detailed understanding of non-equilibrium phenomena in the low-temperature limit.

Lu, Zexi↗
DOE OSTI2022

Light-front wavefunctions of mesons by design

Abstract We develop a mechanism to build the light-front wavefunctions (LFWFs) of meson bound states on a small-sized basis function representation. Unlike in a standard Hamiltonian formalism, the Hamiltonian in this method is implicit, and the information of the system is carried directly by the functional form and adjustable parameters of the LFWFs. In this work, we model the LFWFs for four charmonium states, $$\eta _c$$ η c , $$J/\psi $$ J / ψ , $$\psi '$$ ψ ′ , and $$\psi (3770)$$ ψ ( 3770 ) as superpositions of orthonormal basis functions. We choose the basis functions as eigenfunctions of an effective Hamiltonian, which has a longitudinal confining potential in addition to the transverse confining potential from light-front holographic QCD. We determine the basis function parameters and superposition coefficients by employing both guidance from the nonrelativistic description of the meson states and the experimental measurements of the meson decay widths. With the obtained wavefunctions, we study the features of those meson states, including charge radii and parton distribution functions. We use the $$J/\psi $$ J / ψ LFWF to calculate the meson production in diffractive deep inelastic scattering and ultra-peripheral heavy-ion collisions, and the $$\eta _c$$ η c LFWF to calculate its diphoton transition form factor. Both results show good agreement with experiments. The obtained LFWFs have simple-functional forms and can be readily used to predict additional experimental observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Connecting quasinormal modes and heat kernels in 1-loop determinants

We connect two different approaches for calculating functional determinants on quotients of hyperbolic spacetime: the heat kernel method and the quasinormal mode method. For the example of a rotating BTZ background, we show how the image sum in the heat kernel method builds up the logarithms in the quasinormal mode method, while the thermal sum in the quasinormal mode method builds up the integrand of the heat kernel. More formally, we demonstrate how the heat kernel and quasinormal mode methods are linked via the Selberg zeta function. We show that a 1-loop partition function computed using the heat kernel method may be cast as a Selberg zeta function whose zeros encode quasinormal modes. We discuss how our work may be used to predict quasinormal modes on more complicated spacetimes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

Random motion analysis of flexible satellite structures

A singular perturbation formulation is used to study the responses of a flexible satellite when random measurement errors can occur. The random variables, at different instants of time, are assumed to be uncorrelated. Procedures for obtaining maxima and minima are described, and a variation of the linear method is developed for the formal solution of the two-point boundary-value problems represented by the variational equations. Random and deterministic solutions for the structural position coordinates are studied, and an analytic algorithm for treating the force equation of motion is developed. Since the random system indicated by the variational equation will always be asymptotically unstable, any analysis of stability must be based on the deterministic system.

Huang, T. C.↗
NASA NTRS1981

Langley's Space Shuttle Technology: A bibliography

This bibliography documents most of the major publications, research reports, journal articles, presentations, and contractor reports, which have been published since the inception of the Space Shuttle Technology Task Group at the NASA Langley Reseach Center on July 11, 1969. This research work was performed in house by the Center staff or under contract, monitored by the Center staff. The report is arranged according to method of publication: (1) NASA Formal Reports; (2) Contractor Reports; and (3) Articles and Conferences. Disciplines covered are in the areas of aerothermodynamics, structures, dynamics and aeroelasticity, environmental, and materials. The publications are listed without abstracts for quick reference and planning.

Champine, G. R.↗
NASA NTRS1990

An analysis of scatter decomposition

A formal analysis of a mapping method known as scatter decomposition (SD) is presented. SD divides an irregular domain into many equal-size pieces and distributes them modularly among processors. It is shown that, if a correlation in workload is a convex function of distance, then scattering a more finely decomposed domain yields a lower average processor workload variance; if the workload process is stationary Gaussian and the correlation function decreases linearly in distance to zero and then remains zero, scattering a more finely decomposed domain yields a lower expected maximum processor workload. Finally, if the correlation function decreases linearly across the entire domain, then (among all mappings that assign an equal number of domain pieces to each processor) SD minimizes the average processor workload variance. The dependence of these results on the assumption of decreasing correlation is illustrated with cases where a coarser granularity actually achieves better load balance.

Nicol, David M.↗
NASA NTRS1994

Incorporating interfacial phenomena in solidification models

A general methodology is available for the incorporation of microscopic interfacial phenomena in macroscopic solidification models that include diffusion and convection. The method is derived from a formal averaging procedure and a multiphase approach, and relies on the presence of interfacial integrals in the macroscopic transport equations. In a wider engineering context, these techniques are not new, but their application in the analysis and modeling of solidification processes has largely been overlooked. This article describes the techniques and demonstrates their utility in two examples in which microscopic interfacial phenomena are of great importance.

Beckermann, Christoph↗
NASA NTRS1994

Abstraction Planning in Real Time

When a planning agent works in a complex, real-world domain, it is unable to plan for and store all possible contingencies and problem situations ahead of time. This thesis presents a method for planning a run time that incrementally builds up plans at multiple levels of abstraction. The plans are continually updated by information from the world, allowing the planner to adjust its plan to a changing world during the planning process. All the information is represented over intervals of time, allowing the planner to reason about durations, deadlines, and delays within its plan. In addition to the method, the thesis presents a formal model of the planning process and uses the model to investigate planning strategies.

COMPUTER PROGRAMMING↗