Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational physics applications”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Overview of NASA Langley's Piezoelectric Ceramic Packaging Technology and Applications

Over the past decade, NASA Langley Research Center (LaRC) has developed several actuator packaging concepts designed to enhance the performance of commercial electroactive ceramics. NASA LaRC focused on properly designed actuator and sensor packaging for the following reasons, increased durability, protect the working material from the environment, allow for proper mechanical and electrical contact, afford "ready to use" mechanisms that are scalable, and develop fabrication methodology applicable to any active material of the same physical class. It is more cost effective to enhance or tailor the performance of existing systems, through innovative packaging, than to develop, test and manufacture new materials. This approach led to the development of several solid state actuators that include THUNDER, the Macrofiber Composite or (MFC) and the Radial Field Diaphragm or (RFD). All these actuators are fabricated using standard materials and processes derived from earlier concepts. NASA s fabrication and packaging technology as yielded, piezoelectric actuators and sensors that are easy to implement, reliable, consistent in properties, and of lower cost to manufacture in quantity, than their predecessors (as evidenced by their continued commercial availability.) These piezoelectric actuators have helped foster new research and development in areas involving computational modeling, actuator specific refinements, and engineering system redesign which led to new applications for piezo-based devices that replace traditional systems currently in use.

Bryant, Robert G.↗
DOE OSTI2020

Demonstration of a new unstructured mesh IMC x-ray transport capability in LAP codes

The Advanced Simulation and Computing (ASC) Transport project’s Jayenne Implicit Monte Carlo (IMC) transport library now includes an unstructured mesh capability and is available in a Lagrangian Applications Project (LAP) code. In this presentation, we discuss recent work by the LAP and Transport projects that provides an IMC transport capability for radiation hydrodynamics in the Lagrangian frame. Verification problems testing the new capabilities have been simulated and analyzed, i.e. Marshak wave, Mach 45, Su-Olsen, picket fence, and crooked pipe, both in one and two dimensions. We also present results on two stretch goal problems: a simplified COAX high energy density physics experiment and a supernova shock simulation. Finally, we identify current limitations and future work needed to bring a full capability to the user community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Analyzing Production Rate and Carbon Utilization Trade-offs in CO 2 RR Electrolyzers

Carbon utilization is a crucially important performance characteristic of CO 2 electrolyzers; however, applications of simple and accurate descriptive models to experimental data on this topic are currently lacking. Here, we apply a simple analytical reactor model to parameterize single-pass conversion as a function of feed-gas flow rate and show that it captures a wide body of experimental data in the literature exceptionally well. In doing so, we demonstrate that this simple conceptual approach can characterize progress and capture recent innovations in reactor design. To validate the extracted parameters, we employ a well-established comprehensive model to confirm that the results agree with physical expectations. Lastly, we explore the implications of the descriptive model by examining the inherent trade-off between single-pass conversion efficiency and reactor productivity. Understanding such trade-offs is crucially important for advancing this electrolyzer technology toward commercial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Summer Internship Project: Set Topology [Slides]

Why a Set Topology? 1. Expand the array of tools given to physicists, educators, and students. 2. Through usages of PIC, MPM, and Monte Carlo. 3. My work at LANL was specified to handle requirements for the usage of MPM with FleCSI. 4. What we want it to do - phase change, simulating plasma, and other Multiphysics applications. 5. In FleCSI we already have mesh structures, we needed a way to represent those particles and certain behavior's of these particles that clients desire.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Computationally designed coiled coil ‘bundlemers’ as model colloidal nanoparticles for solution assembly and materials design (Final Report)

As a collaborative team at the University of Delaware and the University of Pennsylvania, Kloxin, Pochan and Saven designed new biomimetic nanomaterials de novo, leveraging a variety of complementary areas of expertise: computational design of biopolymers (Saven at the University of Pennsylvania), and synthesis and characterization (Kloxin and Pochan at the University of Delaware). Overall activities included: sequence-specific peptide synthesis; covalent crosslinking; noncovalent assembly; site-specific functionalization; and nanostructural characterization using electron microscopy and solution-phase (x-ray and neutron) scattering. Using natural and non-natural amino acids, the team created modular, functional peptide building blocks for elaboration of new nanostructured materials. Ultimately, the development of robust peptide-based, building blocks provides tools for researchers to readily produce complex nanomaterial structures in a wide range of applications. The project had three, interconnecting goals in an effort to provide the broader scientific community with a new peptide-based paradigm for materials design and characterization. First, we further developed the coiled-coil bundle-based toolbox (otherwise known as the ‘bundlemer’ toolbox) via computational design with experimental bundle assembly verification. Second, we developed new uses of covalent interactions, in addition to desired physical (noncovalent) interactions, to assemble bundlemers into 1-D polymer chains with targeted chain rigidity, length, and dispersity. Thirds, we used the above designs to experimentally realize (physical or covalent) polymers to target the creation of liquid crystals or to realize interparticle assembly into nanoporous lattices. The close integration of the three groups was instrumental in success of the biomolecular materials design, formation, and understanding for future designs.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Solutions of some problems in applied mathematics using MACSYMA

Various Symbolic Manipulation Programs (SMP) were tested to check the functioning of their commands and suitability under various operating systems. Support systems for SMP were found to be relatively better than the one for MACSYMA. The graphics facilities for MACSYMA do not work as expected under the UNIX operating system. Not all commands for MACSYMA function as described in the manuals. Shape representation is a central issue in computer graphics and computer-aided design. Aside from appearance, there are other application dependent, desirable properties like continuity to certain order, symmetry, axis-independence, and variation-diminishing properties. Several shape representations are studied, which include the Osculatory Method, a Piecewise Cubic Polynomial Method using two different slope estimates, Piecewise Cubic Hermite Form, a method by Harry McLaughlin, and a Piecewise Bezier Method. They are applied to collected physical and chemical data. Relative merits and demerits of these methods are examined. Kinematics of a single link, non-dissipative robot arm is studied using MACSYMA. Lagranian is set-up and Lagrange's equations are derived. From there, Hamiltonian equations of motion are obtained. Equations suggest that bifurcation of solutions can occur, depending upon the value of a single parameter. Using the characteristic function W, the Hamilton-Jacobi equation is derived. It is shown that the H-J equation can be solved in closed form. Analytical solutions to the H-J equation are obtained.

Punjabi, Alkesh↗
DOE OSTI2022

Floquet engineering of tilted and gapped Dirac bandstructure in 1T$$^\prime$$-MoS$$_2$$

Abstract We have developed a rigorous theoretical formalism for Floquet engineering, investigating, and subsequently tailoring most crucial electronic properties of 1T $$^\prime$$ ′ -MoS $$_2$$ 2 by applying an external high-frequency dressing field within the off-resonance regime. It was recently demonstrated that monolayer semiconducting 1T $$^\prime$$ ′ -MoS $$_2$$ 2 exhibits tunable and gapped spin- and valley-polarized tilted Dirac bands. The electron-photon dressed states depend strongly on the polarization of the applied irradiation and reflect a full complexity of the low-energy Hamiltonian for non-irradiated material. We have calculated and analyzed the properties of the electron dressed states corresponding to linear and circular polarization of a dressing field by focusing on their symmetry, anisotropy, tilting, direct and indirect band gaps. Circularly polarized dressing field is known for transition into a new electronic state with broken time-reversal symmetry and a non-zero Chern number, and therefore, the combination of these topologically non-trivial phases and transitions between them could reveal some truly unique and previously unknown phenomena and applications. We have also computed and discussed the density of states for various types of 1T $$^\prime$$ ′ -MoS $$_2$$ 2 materials and its modification in the presence of a dressing field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Inherent stochasticity during insulator–metal transition in VO 2

Vanadium dioxide (VO 2 ), which exhibits a near-room-temperature insulator–metal transition, has great potential in applications of neuromorphic computing devices. Although its volatile switching property, which could emulate neuron spiking, has been studied widely, nanoscale studies of the structural stochasticity across the phase transition are still lacking. In this study, using in situ transmission electron microscopy and ex situ resistive switching measurement, we successfully characterized the structural phase transition between monoclinic and rutile VO 2 at local areas in planar VO 2 /TiO 2 device configuration under external biasing. After each resistive switching, different VO 2 monoclinic crystal orientations are observed, forming different equilibrium states. We have evaluated a statistical cycle-to-cycle variation, demonstrated a stochastic nature of the volatile resistive switching, and presented an approach to study in-plane structural anisotropy. Furthermore, our microscopic studies move a big step forward toward understanding the volatile switching mechanisms and the related applications of VO 2 as the key material of neuromorphic computing.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗
NASA NTRS2018

Application of a First-Principles Anomalous Transport Model for Electrons to Multiple Hall Thrusters and Operating Conditions

We have developed a physics-based model based on a pseudo-particle description of the electron cyclotron drift instability. A key improvement of the model with respect to previous work is that linear theory is not applied in the event of wave saturation and deviations of electrons or ions from a Maxwellian distribution function. In the acceleration region, the anomalous collision frequency is computed as the minimum value necessary to prevent the electron drift velocity from exceeding the thermal velocity. A functional based on the electron equilibration time is defined to control the transition from high to low resistivity regions. The model was previously applied to a single Hall thruster at its nominal operating condition, showing promising results that captured accurately the location of the thruster’s acceleration region. In this paper, we extend the use of this firstprinciples models to two additional thrusters, also considering multiple operating conditions for each of them. Numerical results are compared to experimental measurements obtained with non-invasive laser induced fluorescence. In general, the agreement between experiments and simulations is good. The model is able to predict the location of the acceleration region for all cases. We observe however that fine details, such as changes in the plasma potential gradient within the acceleration regions, are not captured. The model is also insensitive to changes in the magnetic field strength while experiments show that small shifts in location (of less than 5% of the acceleration channel length) occur. We plan to address the weaknesses of our method with the help of physical insight gained from kinetic simulations of the acceleration region.

Chaplin, Vernon H.↗
DOE OSTI2023

Accurate real space iterative reconstruction (RESIRE) algorithm for tomography

Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Cross-property deep transfer learning framework for enhanced predictive analytics on small materials data

Abstract Artificial intelligence (AI) and machine learning (ML) have been increasingly used in materials science to build predictive models and accelerate discovery. For selected properties, availability of large databases has also facilitated application of deep learning (DL) and transfer learning (TL). However, unavailability of large datasets for a majority of properties prohibits widespread application of DL/TL. We present a cross-property deep-transfer-learning framework that leverages models trained on large datasets to build models on small datasets of different properties. We test the proposed framework on 39 computational and two experimental datasets and find that the TL models with only elemental fractions as input outperform ML/DL models trained from scratch even when they are allowed to use physical attributes as input, for 27/39 (≈ 69%) computational and both the experimental datasets. We believe that the proposed framework can be widely useful to tackle the small data challenge in applying AI/ML in materials science.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Broadband coherent XUV light from e − / e + microbunching in an intense laser pulse

Attosecond pulses of coherent extreme ultraviolet (XUV) light are instrumental for investigating subatomic dynamics and can be produced using a free-electron laser (FEL). It has been suggested that an optical FEL, which employs a laser pulse in place of a conventional undulator, could enable a dramatically more compact implementation of such a light source. Yet, the high electron density and subsequent high emittance implied by an optical FEL makes this concept challenging to realize with an electron beam. There has been impressive progress in recent years producing collimated dense and relativistic beams of electrons and positrons in the laboratory. As we demonstrate here, the inherent stability of a quasi-neutral electron-positron beam mitigates Coulomb expansion, and renders it a promising alternative source of coherent light. Specifically, we show via computer simulations that broadband coherent light in the XUV domain, which takes the form of 8-as pulses at 92-as intervals, can be generated by microbunching of relativistic electrons and positrons in a laser pulse. This process occurs over a sub-millimeter length scale, enabling the development of light sources which are orders-of-magnitude more compact than existing sources, with potential applications in physics, chemistry, biology, and industry.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Physics of E × B discharges relevant to plasma propulsion and similar technologies

This paper provides perspectives on recent progress in understanding the physics of devices in which the external magnetic field is applied perpendicular to the discharge current. This configuration generates a strong electric field that acts to accelerate ions. The many applications of this set up include generation of thrust for spacecraft propulsion and separation of species in plasma mass separation devices. These “E × B” plasmas are subject to plasma–wall interaction effects and to various micro- and macroinstabilities. In many devices we also observe the emergence of anomalous transport. This perspective presents the current understanding of the physics of these phenomena and state-of-the-art computational results, identifies critical questions, and suggests directions for future research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Basis light-front quantization for a chiral nucleon-pion Lagrangian

We present the first application of the basis light-front quantization method to a simple chiral model of the nucleon-pion system as a relativistic bound state for the physical proton. In this work, the light-front mass-squared matrix of the nucleon-pion system is obtained within a truncated basis. The mass and the corresponding light-front wave function (LFWF) of the proton are computed by numerical diagonalization of the resulting mass-squared matrix. With the boost-invariant LFWF, we calculate the probability density distribution of the pion's longitudinal momentum fraction and the Dirac form factor of the proton.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A massively parallel time-domain coupled electrodynamics–micromagnetics solver

We present a high-performance coupled electrodynamics–micromagnetics solver for full physical modeling of signals in microelectronic circuitry. The overall strategy couples a finite-difference time-domain approach for Maxwell’s equations to a magnetization model described by the Landau–Lifshitz–Gilbert equation. The algorithm is implemented in the Exascale Computing Project software framework, AMReX, which provides effective scalability on manycore and GPU-based supercomputing architectures. Furthermore, the code leverages ongoing developments of the Exascale Application Code, WarpX, which is primarily being developed for plasma wakefield accelerator modeling. Our temporal coupling scheme provides second-order accuracy in space and time by combining the integration steps for the magnetic field and magnetization into an iterative sub-step that includes a trapezoidal temporal discretization for the magnetization. The performance of the algorithm is demonstrated by the excellent scaling results on NERSC multicore and GPU systems, with a significant (59×) speedup on the GPU using a node-by-node comparison. We demonstrate the utility of our code by performing simulations of an electromagnetic waveguide and a magnetically tunable filter.

97 MATHEMATICS AND COMPUTING↗