Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Benchmarking of hydrodynamic plasma waveguides for multi-GeV laser-driven electron acceleration

Hydrodynamic plasma waveguides initiated by optical field ionization have recently become a key component of multi-GeV laser wakefield accelerators. Here, we present the most complete and accurate experimental and simulation-based characterization to date, applicable to current multi-GeV experiments and future 100 GeV-scale laser plasma accelerators. Crucial to the simulations is the correct modeling of intense Bessel beam interaction with meter-scale gas targets, the results of which are used as initial conditions for hydrodynamic simulations. The simulations are in good agreement with our experiments measuring evolving plasma and neutral hydrogen density profiles using two-color short pulse interferometry, enabling realistic determination of the guided mode structure for application to laser-driven plasma accelerator design. Published by the American Physical Society 2024

Physics↗
NASA NTRS2004

Investigation of Volumetric Sources in Airframe Noise Simulations

Hybrid methods for the prediction of airframe noise involve a simulation of the near field flow that is used as input to an acoustic propagation formula. The acoustic formulations discussed herein are those based on the Ffowcs Williams and Hawkings equation. Some questions have arisen in the published literature in regard to an apparently significant dependence of radiated noise predictions on the location of the integration surface used in the solution of the Ffowcs Williams and Hawkings equation. These differences in radiated noise levels are most pronounced between solid-body surface integrals and off-body, permeable surface integrals. Such differences suggest that either a non-negligible volumetric source is contributing to the total radiation or the input flow simulation is suspect. The focus of the current work is the issue of internal consistency of the flow calculations that are currently used as input to airframe noise predictions. The case study for this research is a computer simulation for a three-element, high-lift wing profile during landing conditions. The noise radiated from this flow is predicted by a two-dimensional, frequency-domain formulation of the Ffowcs Williams and Hawkings equation. Radiated sound from volumetric sources is assessed by comparison of a permeable surface integration with the sum of a solid-body surface integral and a volume integral. The separate noise predictions are found in good agreement.

Casper, Jay H.↗
DOE OSTI2020

Analytic model of electron transport through and over non-linear barriers

Tunneling barriers are an essential component of electron sources, sensors, detectors, and vacuum nanoelectronics, and a pivotal factor in their performance, but the barriers themselves routinely depart from the analytic models used to model their behavior. Here, a new formalism is developed to analytically and accurately model emission through and over barriers associated with depletion layers, nanotip barriers, and MIM structures. The transmission probability for depletion layer and metal-insulator-metal (MIM) and metaloxide-semiconductor (MOS) barriers is accurately modeled as the electron energy exceeds the barrier height using approaches designed for rapid implementation demanded by simulation codes, and extensible to general barriers. The models supersede conventional thermal and field models in depletion and MIM/MOS barrier studies. Thermal-field methods are used to treat the transmission probability, and shape factor methods to treat the tunneling factor. Analytic formulae for current density are obtained. The methods ease device simulation and characterization of current-voltage relations for emerging technologically interesting barriers with better accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Insights into the Properties of MXenes and MXene Analogs from Atomistic Simulation

We review the properties that have been predicted for MXenes and MXene analogs from computational simulation with a focus on structural and electronic properties, energy storage, and ion transport. Methods considered range from quantum mechanical approaches to classical molecular dynamics. We conclude by reviewing current limitations and outstanding questions for simulation and MXene properties that have been little explored to-date.

Muraleedharan, Murali Gopal↗
NASA NTRS2019

Simulation of Non-Markovian Dynamics on IBM QX

Currently available quantum computers allow us to run proof of principle algorithms that are unitary in their nature. Therefore, this architecture is unoptimized for simulation of an open quantum system. Here we present a method that helps us to overcome unitarity. We show how to run a non-Markovian evolution of a qubit system. We discuss all the discrepancies from theoretical predictions.

Wudarski, Filip A.↗
DOE OSTI2019

Beam Current Effect on Microstructure and Properties of Electron-Beam-Melted Ti-6Al-4V Alloy

In this study, a noncommercial 3D printing machine was used to fabricate Ti-6Al-4V alloy by electron-beam melting (EBM). The influence of beam current on the microstructure, phase composition and mechanical properties of electron-beam-melted Ti-6Al-4V alloy was investigated. Numerical simulation is implemented to evaluate thermal fields during electron-beam melting of Ti-6Al-4V powder. The decrease in beam current from 3.5 to 2.5 mA leads to refinement of microstructure: The average width of α plates decreases from 10 down to 6 μm. The formation of finer microstructure is attributed to higher cooling rate at lower beam current confirmed by simulation. The phase composition of EBM Ti-6Al-4V indirectly depends on the beam current. High content of β phase (7%) was achieved at the beam current of 3 mA. The produced Ti-6Al-4V samples are characterized by high microhardness (470-520 HV).

36 MATERIALS SCIENCE↗
DOE OSTI2018

Finding Maximum Cliques on the D-Wave Quantum Annealer

This work assesses the performance of the D-Wave 2X (DW) quantum annealer for finding a maximum clique in a graph, one of the most fundamental and important NP-hard problems. Because the size of the largest graphs DW can directly solve is quite small (usually around 45 vertices), we also consider decomposition algorithms intended for larger graphs and analyze their performance. For smaller graphs that fit DW, we provide formulations of the maximum clique problem as a quadratic unconstrained binary optimization (QUBO) problem, which is one of the two input types (together with the Ising model) acceptable by the machine, and compare several quantum implementations to current classical algorithms such as simulated annealing, Gurobi, and third-party clique finding heuristics. We further estimate the contributions of the quantum phase of the quantum annealer and the classical post-processing phase typically used to enhance each solution returned by DW. We demonstrate that on random graphs that fit DW, no quantum speedup can be observed compared with the classical algorithms. On the other hand, for instances specifically designed to fit well the DW qubit interconnection network, we observe substantial speed-ups in computing time over classical approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effect of Diffusion Length in Modeling of Equiaxed Dendritic Solidification under Buoyancy Flow in a Configuration of Hebditch–Hunt Experiment

Modeling of equiaxed solidification is vital for understanding the solidification process of metallic alloys. In this work, an extended literature review is given for the models currently used for equiaxed solidification simulations. Based on this analysis, we present a three-phase multiscale equiaxed solidification model in which some approximations regarding solute transport at microscopic scale are put together in a new way and incorporated into macroscopic transport equations. For the latter, a term relating to the momentum exchange between the two phases is revised, and a modification for the grain packing algorithm is proposed. A modernized model for equiaxed dendrite growth is tested using a case of solidification of Sn-5 wt pct Pb alloy in a parallelepiped cavity that mimics the Hebditch–Hunt experiment. The results obtained using two approaches to calculate diffusion length are presented and compared both with each other and with numerical results from elsewhere. It is demonstrated that diffusion length has a crucial effect on the final segregation pattern.

36 MATERIALS SCIENCE↗
DOE OSTI2020

PERC silicon PV infrared to ultraviolet optical model

Here, optical modeling of aluminum back surface field (Al-BSF) and passivated emitter and rear contact (PERC) silicon wafer solar cells are extended into the infrared wavelength range (300–2500 nm). The model incorporates ray tracing to account for pyramidal texture of the front silicon surface, free carrier absorption in the front emitter; bulk silicon wafer, BSF; and Lambertian scattering at the silicon back contact. Total reflectance of pyramidal textured cells with different contact geometries with and without encapsulation are simulated. Simulated reflectance and short circuit current density match well with experimental values. These physics-based optical models serve as input for evaluating optical and thermal management strategies for silicon wafer solar cells.

14 SOLAR ENERGY↗
DOE OSTI2021

Mechanisms of Dynamic Near-Wake Modulation of a Utility-Scale Wind Turbine

The current study uses large eddy simulations to investigate the transient response of a utility-scale wind turbine wake to dynamic changes in atmospheric and operational conditions, as observed in previous field-scale measurements. Most wind turbine wake investigations assume quasi-steady conditions, but real wind turbines operate in a highly stochastic atmosphere, and their operation (e.g. blade pitch, yaw angle) changes constantly in response. Furthermore, dynamic control strategies have been recently proposed to optimize wind farm power generation and longevity. Therefore, improved understanding of dynamic wake behaviours is essential. First, changes in blade pitch are investigated and the wake expansion response is found to display hysteresis as a result of flow inertia. The time scales of the wake response to different pitch rates are quantified. Next, changes in wind direction with different time scales are explored. Under short time scales, the wake deflection is in the opposite direction of that observed under quasi-steady conditions. Finally, yaw changes are implemented at different rates, and the maximum inverse wake deflection and time scale are quantified, showing a clear dependence on yaw rate. To gain further physical understanding of the mechanism behind the inverse wake deflection, the streamwise vorticity in different parts of the wake is quantified. The results of this study provide guidance for the design of advanced wake flow control algorithms. The lag in wake response observed for both blade pitch and yaw changes shows that proposed dynamic control strategies must implement turbine operational changes with a time scale of the order of the rotor time scale or slower.

flow-structure interactions↗
DOE OSTI2024

Chemistry Informed Machine Learning-Based Heat Capacity Prediction of Solid Mixed Oxides

Knowing heat capacity is crucial for modeling temperature changes with the absorption and release of heat and for calculating the thermal energy storage capacity of oxide mixtures with energy applications. The current prediction methods (ab initio simulations, computational thermodynamics, and the Neumann–Kopp rule) are computationally expensive, not fully generalizable, or inaccurate. Machine learning has the potential of being fast, accurate, and generalizable, but it has been scarcely used to predict mixture properties, particularly for mixed oxides. Here, we demonstrate a method for the generalizable prediction of heat capacity of solid oxide pseudobinary mixtures using heat capacity data obtained from computational thermodynamics and descriptors from ab initio databases. Further, models trained through this workflow achieved an error (mean absolute error of 0.43 J mol –1 K –1 ) lower than the uncertainty in differential scanning calorimetry measurements, and the workflow can be extended to predict other properties derived from the Gibbs free energy and for higher-order oxide mixtures.

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

Polymer Topology Effects on Dynamics of Comb Polymer Melts

Controlling polymer viscosity and flow is key to their many applications through strength and processability. The topology of the polymer i.e., linear, stars, and branched, affects the macroscopic flow characteristics of melts, where introducing one branch is sufficient to increase the viscosity significantly. While a number of studies have probed the effects of polymer topology on their rheology, the molecular understanding that underlies the macroscopic behavior remains an open question. The current study uses molecular dynamics simulations to resolve the effects of topology of polymer melts on chain mobility and viscosity in the comb regime using polyethylene as a model system. A coarse-grained model where four methylene groups constitute one bead is used, and the results are transposed to the atomistic level. Here, we find that while the number of branches only slightly affects the chain mobility and viscosity, their length strongly impacts their behavior. Finally, the results are discussed in terms of interplay between the relaxation of the branches and reptation of the backbone where the topology of the polymer affects the tube dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure of iso-density sets in supersonic isothermal turbulence

Context: The gas density structure of the cold molecular phase of the interstellar medium is the main controller of star formation. Aims: A theoretical framework is proposed to describe the structural content of the density field in isothermal supersonic turbulence. Methods: It makes use of correlation and structure functions of the phase indicator field defined for different iso-density values. The relations between these two-point statistics and the geometrical features of iso-density sets such as the volume fraction, the surface density, the curvature, and fractal characteristics are provided. An exact scale-by-scale budget equation is further derived revealing the role of the turbulent cascade and dilation on the structural evolution of the density field. Although applicable to many flow situations, this tool is here first invoked for characterising supersonic isothermal turbulence, using data from the currently best-resolved numerical simulation. Results: We show that iso-density sets are surface fractals rather than mass fractals, with dimensions that markedly differ between dilute, neutral, and dense regions. The surface–size relation is established for different iso-density values. We further find that the turbulent cascade of iso-density sets is directed from large towards smaller scales, in agreement with the classical picture that turbulence acts to concentrate more surface into smaller volumes. Intriguingly, there is no range of scales that complies with a constant transfer rate in the cascade, challenging our fundamental understanding of interstellar turbulence. Finally, we recast the virial theorem in a new formulation drawing an explicit relation between the aforementioned geometrical measures and the dynamics of iso-density sets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Automated and Distributed Monte Carlo Generation for GlueX

MCwrapper is a set of systems that manages the entire Monte Carlo production workflow for GlueX and provides standards for how that Monte Carlo is produced. MCwrapper was designed to be able to utilize a variety of batch systems in a way that is relatively transparent to the user, thus enabling users to quickly and easily produce valid simulated data at home institutions worldwide. Additionally, MCwrapper supports an autonomous system that takes user’s project submissions via a custom web application. The system then atomizes the project into individual jobs, matches these jobs to resources, and monitors the jobs status. The entire system is managed by a database which tracks almost all facets of the systems from user submissions to the individual jobs themselves. Users can interact with their submitted projects online via a dashboard or, in the case of testing failure, can modify their project requests from a link contained in an automated email. Beginning in 2018 the GlueX Collaboration began to utilize the Open Science Grid (OSG) to handle a bulk of simulation tasks; these tasks are currently being performed on the OSG automatically via MCwrapper. This talk will outline the entire system of MCwrapper, its use cases, and the unique challenges facing the system.

Britton, Thomas↗
DOE OSTI2025

Re-evaluating the prompt fission neutron spectrum of spontaneously fissioning 252 Cf

The prompt fission neutron spectrum (PFNS) of spontaneously fissioning 252 Cf is a Neutron Data Standards observable. Nearly all fission spectra of actinides were measured relative to it, using efficiencies derived from it, or analyzed with simulations validated by it. The current Standards evaluation was published by W. Mannhart in 1987. It could not be updated because the evaluation input, experimental mean values and covariances, were lost. First, we attempt to reproduce it. However, Mannhart’s evaluation can only be reproduced within its one-σ uncertainties as some of its aspects (e.g., experimental covariances, rejected data points) remain unknown. Therefore, a new evaluation is presented: We revisit all existing experimental 252 Cf(sf) PFNS data, including those published after the release of the current Standards evaluation, and re-estimate associated covariances. The newly evaluated 252 Cf(sf) PFNS differs distinctly from Mannhart’s below 300 keV and extends it to lower and higher outgoing neutron energies (500 eV–25 MeV). The new evaluated uncertainties are larger from 3–9 MeV and smaller otherwise. Spectrum averaged cross sections of importance to the International Reactor Dosimetry and Fusion File community calculated with the new spectrum are close to those calculated with Mannhart’s evaluation and agree with experimental values well within their uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Radio-frequency wave interactions with a plasma sheath in oblique-angle magnetic fields using a sheath impedance model

Here, the physics of interactions between waves in plasmas and sheaths for background magnetic fields which make oblique angles with the sheath surfaces is studied with the use of the self-consistent finite element code rfSOL incorporating the recently developed sheath impedance model. The calculation based on this model employs the generalized sheath boundary condition (sheath BC), which surpasses the previously used capacitive sheath BC in reliability by taking into account the contributions of the ion and electron currents in the sheath as well as the displacement current. A series of numerical simulations is carried out in two-dimensional slab geometry with a flat or curved sheath surface as part of the boundary. It is shown that the sheath–plasma wave appears when the equilibrium magnetic field line angle with respect to the sheath surface is small, the absolute value of the radio-frequency (RF) sheath voltage is large, and the plasma density is slightly higher than the lower hybrid resonance density (LHR density), all of which bring the sheath property closer to being capacitive. It is also shown that the sharp variation of the magnetic field line angle along the sheath surface can sensitively affect the maximum absolute value of the RF sheath voltage at a plasma density slightly lower than the LHR density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. We examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. We give suggestions on future experimental and modeling work on the liquid–ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Sliding dynamics of skyrmion molecular crystals

Using both atomistic and particle-based simulations, we investigate the current-driven dynamics of skyrmions on two-dimensional periodic substrates when there are multiple skyrmions per substrate minimum. At zero drive, the system forms pinned skyrmion molecular crystal states consisting of dimers, trimers, or dimer-trimer mixtures that have both positional and orientational order. On a square substrate lattice, the motion above depinning occurs via a running soliton that travels completely transverse to the applied current. This motion is generated by a torque from the Magnus force, which rotates the n-mer states perpendicular to the applied current. At higher drives, the flow becomes disordered while the Hall angle diminishes and gradually approaches the intrinsic value. In some cases, we also find directional locking where the Hall angle becomes locked to certain symmetry directions of the substrate over a range of currents. The transitions into and out of directionally locked states are accompanied by negative differential mobility in which the net velocity decreases as the drive increases. In conclusion, on a triangular substrate, we find no transverse mobility effects, but still observe directionally locked motion.

36 MATERIALS SCIENCE↗