Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Recent progress in microscale modeling of RF sheaths

The microscale properties of RF sheaths in the ion cyclotron range of frequencies (ICRF) are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) code modeling. Previous work that parametrized RF sheath properties, specifically the RF sheath impedance and the rectified (DC) sheath potential, is generalized to include the effect of net DC current flow through the sheath. Analytical results are presented in the low frequency limit where the displacement current is negligible, and tested against results from a fluid numerical model. Here, it is shown that when the sheath draws DC electron current, the voltage rectification is reduced from the zero current case, and the electron admittance is increased. In separate but related work on the microscale model, selected cases have been simulated with PIC codes to validate, further illuminate and extend fluid model results and their parametrizations. Quantitative agreement in trends for voltage rectification and sheath admittance vs. RF driving voltage is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Roles of bandgrading, lifetime, band alignment, and carrier concentration in high-efficiency CdSeTe solar cells

CdSeTe alloying has significantly increased the efficiency of CdTe-based solar technology. Here, computational modeling compares how different CdSeTe bandgrading, carrier lifetimes, band alignment, and carrier concentrations contribute to transport, recombination, and performance. We find that the gain in photocurrent caused by bandgap narrowing alone is insufficient to describe experimental efficiency gains. Performance can be increased by adjusting CdSeTe compositions and bandgrading depths. However, these performance gains are small relative to the contributions of enhanced lifetime by Se alloying, which can explain record cell efficiency gains with minimal open-circuit voltage loss despite significant bandgap narrowing. Similarly, CdSeTe band alignment shifts can significantly increase performance if front interface recombination is prevalent. For a wide range of CdSeTe grading profiles, the hole density is a critical component to achieve efficiencies exceeding 25%.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ionization wave propagation in a He plasma jet in a controlled gas environment

Characterizing ionization wave propagation in low temperature plasma jets is critical to predicting production of reactive species and plasma–surface interactions for biomedical applications and surface functionalization. In this work, results from optical emission and laser induced fluorescence measurements of the ionization wave in a He plasma jet operating in a controlled gas environment are discussed and used for comparison with numerical modeling. Additionally, the ionization wave was observed using ICCD (Intensified Charge Coupled Device) imaging and characterized by time and spatially resolved electron density measurements using laser-collision-induced fluorescence. The plasma jet was initially characterized using pure He (nominally at 200 Torr), while varying pressure and voltage. When operating in pure He, the ionization wave broadly expands exiting the plasma tube. Increasing the operating pressure reduces the speed and isotropic expansion of the ionization wave. Furthermore, the jet operated with a humid He shroud was also studied. The humid He shroud results in the electron density increasing and having an annular profile due to the lower ionization potential of H2O compared to He and localized photoionization in the mixing region. Numerical modeling highlighted the importance of resonance radiation emitted by excited states of He, photoelectron emission from the quartz tube, and the kinetic behavior of the electrons produced by photoionization ahead of the ionization front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A fluid-structure interaction solver for investigating torsional galloping in solar-tracking photovoltaic panel arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. Furthermore, these simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects. These results are supplemented by a discussion of stability trends noted during testing and possible implications when considering multi-row array interactions.

14 SOLAR ENERGY↗
DOE OSTI2021

Tokamak ITG-KBM transition benchmarking with the mixed variables/pullback transformation electromagnetic gyrokinetic scheme

Electromagnetic gyrokinetic simulation of high temperature plasma is required to predict confinement in magnetic fusion devices and has posed challenges for existing codes. In this paper, we demonstrate successful global gyrokinetic simulation of the ion temperature gradient-driven mode-kinetic ballooning mode transition in a toroidal fusion plasma test case using the mixed variables/pullback transformation (MV/PT) scheme with the particle-in-cell codes XGC and ORB5, and compare to results from a conventional continuum code from the literature. Furthermore, the MV/PT scheme combines explicit time integration with mitigation of the well-known electromagnetic gyrokinetic “cancelation problem.” We calculate eigenmodes in the electrostatic and parallel vector potentials, and find good agreement in growth rate, real frequency, and the normalized plasma pressure of mode transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of decelerating streamers in inhomogeneous atmosphere with implications for runaway electron generation

The dynamics of positive and negative streamers is numerically simulated in atmospheric pressure air in the range of parameters corresponding to the streamer deceleration and termination in the middle of the discharge gap. A detailed comparison with experiments in air at constant and variable density demonstrates good agreement between the 2D simulation results and the observations. It is shown that positive and negative streamers behave in radically different ways when decelerating and stopping. When the head potential drops, the negative streamer transits to the mode in which the propagation is due to the forward electron drift. In this case, the radius of the ionization wave front increases, whereas the electric field at the streamer head decreases further and the streamer stops. Its head diameter continues to increase due to the slow drift of free electrons in the residual under-breakdown field. On the contrary, the only advancement mechanism for a positive streamer with a decreasing head potential is a decrease in the effective radius of the ionization wave, leading to a local increase in the electric field. This mechanism makes it possible to compensate for the reduction in the efficiency of gas photoionization at small head diameters. A qualitative 1D model is suggested to describe streamer deceleration and stopping for different discharge polarities. Estimates show that, during positive streamer stopping, the local electric field at the streamer head can exceed the threshold corresponding to the transition of electrons to the runaway mode when the head potential (relative to the surrounding space) decreases to~1.2 kV in atmospheric pressure air. In this case, pulsed generation of a beam of runaway electrons directed into the channel of a stopping positive streamer can occur. The energy of the formed pulsed electron beam depends on the intensity of photoionization in front of the streamer head. This energy can vary from 700 V (when increasing the photoionization rate by a factor of 10 with respect to the value in atmospheric pressure air) to 2.6 kV (when decreasing the photoionization rate by a factor of 1000). It is possible that this behavior of decelerating positive streamers can explain the observed bursts of x-ray radiation during the streamer propagation in long air gaps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Formation of edge pressure pedestal and reversed magnetic shear due to toroidal rotation in a tokamak equilibrium

Toroidal rotation is well known to play significant roles in the edge transport and L-H transition dynamics of tokamaks. Our recent calculation finds that a sufficiently strong localized toroidal rotation can directly bring out the formation of edge pressure pedestal with reversed magnetic shear that is reminiscent of an H-mode plasma, purely through the effects of toroidal rotation on the tokamak MHD equilibrium itself. In particular, the enhanced edge toroidal rotation enables a substantial peaking of the parallel current profile near edge in higher β regimes, which leads to the flattening or reversal of the local q (safety factor) profile. Here the formation of pressure pedestal along with the reversed magnetic shear region is shown to be the natural outcome of the MHD tokamak equilibrium in a self-consistent response to the presence of a localized toroidal rotation typically observed in H-mode or QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Filamentation and focusing of electron beams due to interactions with plasma waves

Results of numerical modeling of the interaction of an electron beam propagating across relativistic plasma waves indicate that electron beam filamentation and focusing may occur under certain conditions. The model is based on solving the relativistic equation of motion in three dimensions for the individual electrons in a tenuous Gaussian beam, as they pass through a relativistic plasma wave. Several electron beam and plasma wave parameters were varied, and the results are summarized. One of the results is that the spacing of the electron beam filaments correlates with the wavelength of the plasma wave. The electron beam filaments appear as vertical slabs after the beam exits the plasma waves. The electron beam also compresses to a focus after it exits the plasma, and the focal distance depends on several parameters including the electron beam energy and phase velocity of the relativistic plasma wave. It is suggested that these focusing and filamentation phenomena may be the basis for diagnostics schemes for laser plasma interactions. The parameters used in the model were electron beam energies in the 5–50 keV range and plasma wave properties typical for the beat-wave produced by CO2 lasers, which correspond to the facilities available in our laboratory. The limitations of these results to lower energy density beam and plasma regimes and to higher energy density regimes will be discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Modeling of tungsten filament in gas discharge in H - ion source

Sources of negative ions such as H - are essential elements of proton accelerators and tokamaks. They have limited lifetime. The replacement of an ion source is a costly process causing delays in the operation of the entire machine. The hot tungsten filament is the key element that limits the lifetime of the H - ion source at the Los Alamos Neutron Science Center (LANSCE) facility. An accurate model, which describes the filament physics and provides a reliable estimate on its lifetime, may significantly improve the operation of the facility. In this report we describe such a model, including a comprehensive list of relevant physical processes from first principles. The model can be used to describe different regimes of operation using diagnostics data as input parameters. It has been benchmarked against the data collected during the production cycle at LANSCE and shows good agreement with experimental data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Computation of the Biot–Savart line integral with higher-order convergence using straight segments

One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot–Savart field from each straight segment is analytically known. However, if the endpoints of the straight segments are chosen to lie on the coil, then the accuracy of the Biot–Savart computation is generally only the second order in the number of endpoints. In this work, we propose a simple modification: shift each end point of the coil in the outward normal direction by an amount proportional to the local curvature. With this modification, the Biot–Savart accuracy increases to the fourth order and the numerical error is dramatically reduced for a given number of discretization points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Diverse simulations of time-resolved photoluminescence in thin-film solar cells: A SnO 2 /CdSe y Te 1-y case study

Time-resolved photoluminescence (TRPL) is widely used to measure carrier lifetime in thin-film solar cell absorbers. However, the injection dependence of data and frequent non-exponential decay shapes complicate the interpretation. Here, we develop a numerical model to simulate injection-dependent TRPL measurements in a SnO 2 /CdSe y Te 1-y case study solar cell structure, considering parameters of interest to researchers in industry and academia. Previous simulations have shown that, in low injection, excess electrons and holes injected by the laser pulse are rapidly separated in the electric field formed by the pn junction. As a result, at early times, the PL signal can decay faster than the Shockley-Read-Hall lifetime in the absorber bulk (t_bulk). Prior simulations have shown that the charge stored in the junction can slowly leak out to affect decays at late times. However, it has not been clear if and to what degree charge storage can affect the slopes extracted from TRPL decays - t_2 - commonly cited as the TRPL-measured lifetime. Here, we show that charge storage can, in some cases, result in t_2 values that substantially overestimate t_bulk. Previous simulations indicate that high injection conditions can screen the junction field and minimize charge separation. Here, we show that continued injection increases can drive down t_2 below t_bulk as radiative recombination becomes dominant. We catalog charge storage and radiative recombination impacts for a diverse set of material parameters and compare results to double-heterostructure models.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,S↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

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

Quantum computing without quantum computers: Database search and data processing using classical wave superposition

Quantum computers are proven to be more efficient at solving a specific class of problems compared to traditional digital computers. Superposition of states and quantum entanglement are the two key ingredients that make quantum computing so powerful. However, not all quantum algorithms require quantum entanglement (e.g., search through an unsorted database). Is it possible to utilize classical wave superposition to speed up database searching as much as by using quantum computers? There were several attempts to mimic quantum computers using classical waves. It was concluded that the use of classical wave superposition comes with the cost of an exponential increase in resources. In this work, we consider the feasibility of building classical wave-based devices able to provide fundamental speedup over digital counterparts without the exponential overhead. We present experimental data on database searching through a magnetic database using spin wave superposition. The results demonstrate the same speedup as expected for quantum computers. Also, we present examples of numerical modeling demonstrating classical wave interference for period finding. This approach may not compete with quantum computers with efficiency but outperform classical digital computers. We argue that classical wave-based devices can perform some of the quantum algorithms with the same efficiency as quantum computers as long as quantum entanglement is not required.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Hybrid simulation of Alfvén wave parametric decay instability in a laboratory relevant plasma

Large-amplitude Alfvén waves are subject to parametric decays, which can have important consequences in space, astrophysical, and fusion plasmas. Although the Alfvén wave parametric decay instability (PDI) was predicted decades ago, observational evidence is limited, stimulating considerable interest in laboratory demonstration of the instability and associated numerical modeling. Here, we report an important step toward direct hybrid simulation of the Alfvén wave PDI in a laboratory plasma, using antenna-like wave injection of a circularly polarized wave and realistic wave-plasma parameters. Considering collisionless damping, we identify the threshold Alfvén wave amplitudes and frequencies required for triggering the instability. These threshold behaviors are corroborated by simple theoretical analysis, considering the balance between PDI growth and ion Landau damping and PDI development in a bounded plasma. Other effects not included in the present model such as finite transverse wave scale and ion–neutral collisions are briefly discussed. These hybrid simulations promise to be a useful tool for investigating laboratory Alfvén wave dynamics and may provide guidance for future laboratory demonstration of the PDI process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Constrained model calibration of grain structure dependent spall dynamics in shock-loaded tantalum

Here, we perform a gas gun experiment by shock loading tantalum samples of varying grain structures to assess the suitability of a numerical model for simulating spall behavior. The observed differences in spall strength, as well spallation and re-compression history, are not captured in uncalibrated hydrodynamic simulations. An optimization is performed on the Johnson spall model to determine the best parameters that fit the observed trends. Linear stability analysis is employed to motivate bounds on those parameters. Herein, optimized simulations agree well with the experimental results, reproducing pullback depth and recompression timescales across the different samples tested. Further, the observed pullback time of the single crystal sample was found to imply, via the stability analysis, a percolation threshold in good agreement with the theoretical value for a body centered cubic lattice. Therefore, the combined linear stability and percolation analysis shows promise and may be applied to other materials with diverse microstructures. Collectively, the findings demonstrate that the model is suitable for reproducing spall-induced free surface behavior across various microstructures, but also points to caution in using model coefficients for uncalibrated microstructures.

58 GEOSCIENCES↗
DOE OSTI2023

Non-ideal MHD growth of current interchange tearing modes at plasma edge and response to externally-imposed flow

Current interchange tearing modes (CITMs) have been proposed to explain an intermittent eruption of electric current toward a divertor at an edge region of a tokamak device. Extended magnetohydrodynamic (MHD) simulations combined with the Scrape-Off-Layer diffusion model, which has been developed by the authors [Miura et al., Phys. Plasmas 24, 092111 (2017)], are carried out in order to study the growth of a CITM under influences of two-fluid and gyro-viscous effects, or non-ideal MHD effects. Numerical simulations both with and without an externally imposed azimuthal flow show that the growth of a CITM is observed for a relatively small flow velocity, and that a CITM can be a candidate mechanism for the intermittent eruption of the current in a tokamak, whether the growth is under influences of non-ideal MHD effects or not. Here, the growth can be suppressed when an externally imposed azimuthal flow is sufficiently large to cause a finite radial displacement because such a displacement prevents a transition from an interchange to a tearing mode that is the key process of a CITM. Furthermore, it is also found that a stripe pattern similar to so-called streamers is formed in the course of the CITM growth. The basic nature of this streamer-like structure is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Numerical study of magnetic island coalescence using magnetohydrodynamics with adaptively embedded particle-in-cell model

Collisionless magnetic reconnection typically requires kinetic treatment that is, in general, computationally expensive compared to fluid-based models. In this study, we use the magnetohydrodynamics with an adaptively embedded particle-in-cell (MHD-AEPIC) model to study the interaction of two magnetic flux ropes. This innovative model embeds one or more adaptive PIC regions into a global MHD simulation domain such that the kinetic treatment is only applied in regions where the kinetic physics is prominent. We compare the simulation results among three cases: (1) MHD with adaptively embedded PIC regions, (2) MHD with statically (or fixed) embedded PIC regions, and (3) a full PIC simulation. The comparison yields good agreement when analyzing their reconnection rates and magnetic island separations as well as the ion pressure tensor elements and ion agyrotropy. In order to reach good agreement among the three cases, large adaptive PIC regions are needed within the MHD domain, which indicates that the magnetic island coalescence problem is highly kinetic in nature, where the coupling between the macro-scale MHD and micro-scale kinetic physics is important.

79 ASTRONOMY AND ASTROPHYSICS↗