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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Multipactor threshold calculation of coaxial transmission lines in microwave applications with nonstationary statistical theory

This paper presents a statistical theory for the initial onset of multipactor breakdown in coaxial transmission lines, taking both the nonuniform electric field and random electron emission velocity into account. A general numerical method is first developed to construct the joint probability density function based on the approximate equation of the electron trajectory. The nonstationary dynamics of the multipactor process on both surfaces of coaxial lines are modelled based on the probability of various impacts and their corresponding secondary emission. The resonant assumption of the classical theory on the independent double-sided and single-sided impacts is replaced by the consideration of their interaction. As a result, the time evolutions of the electron population for exponential growth and absorption on both inner and outer conductor, in response to the applied voltage above and below the multipactor breakdown level, are obtained to investigate the exact mechanism of multipactor discharge in coaxial lines. Furthermore, the multipactor threshold predictions of the presented model are compared with experimental results using measured secondary emission yield of the tested samples which shows reasonable agreement. Finally, the detailed impact scenario reveals that single-surface multipactor is more likely to occur with a higher outer to inner conductor radius ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Mathematical model of massive dynamics in the neighborhood of disturbance focus

We consider a model problem for the research of Earth’s surface state caused by discharge of tectonic stress on the deep cracks in the earth’s crust. Using the numerical method of bicharacteristics, algorithm is developed and solution of the problem of non-stationary dynamics of a homogeneous isotropic elastic half-space is built under plane strain when one resets the vertical stress on crack. Diffraction of waves on the surface is studied, pictures of wave fields and stress-strain state of the medium and the surface for different moments of time are built.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Landau damping of Langmuir twisted waves with kappa distributed electrons

The kinetic theory of Landau damping of Langmuir twisted modes is investigated in the presence of orbital angular momentum of the helical (twisted) electric field in plasmas with kappa distributed electrons. The perturbed distribution function and helical electric field are considered to be decomposed by Laguerre-Gaussian mode function defined in cylindrical geometry. The Vlasov-Poisson equation is obtained and solved analytically to obtain the weak damping rates of the Langmuir twisted waves in a nonthermal plasma. The strong damping effects of the Langmuir twisted waves at wavelengths approaching Debye length are also obtained by using an exact numerical method and are illustrated graphically. The damping rates of the planar Langmuir waves are found to be larger than the twisted Langmuir waves in plasmas which shows opposite behavior as depicted in Fig. 3 by J. T. Mendoça [Phys. Plasmas 19, 112113 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Oxide-cladding aluminum nitride photonic crystal slab: Design and investigation of material dispersion and fabrication induced disorder

Photonic crystal slabs with a lower-index material surrounding the core layer are an attractive choice to circumvent the drawbacks in the fabrication of membranes suspended in air. In this work we propose a photonic crystal (PhC) slab structure composed of a triangular pattern of air holes in a multilayer thin film of aluminum nitride embedded in silicon dioxide layers designed for operating around 450 nm wavelengths. We show the design of an ideal structure and analyze the effects of material dispersion based on a first-order correction perturbation theory approach using dielectric functions obtained by experimental measurements of the thin film materials. Numerical methods were used to investigate the effects of fabrication induced disorder of typical nanofabrication processes on the bandgap size and spectral response of the proposed device. Deviation in holes radii and positions were introduced in the proposed PhC slab model with a Gaussian distribution profile. Impacts of slope in holes sidewalls that might result from the dry etching of AlN were also evaluated. The results show that for operation at the midgap frequency, slope in holes sidewalls is more critical than displacements in holes sizes and positions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A simulation-optimization model for Stone column-supported embankment stability considering rainfall effect

Landslide due to rainfall has been and continues to be one of the most important concerns of geotechnical engineering. The paper presents the variation of factor of safety of stone column-supported embankment constructed over soft soil due to change in water level for an incessant period of rainfall. A combined simulation-optimization based methodology has been proposed to predict the critical surface of failure of the embankment and to optimize the corresponding factor of safety under rainfall conditions using an evolutionary genetic algorithm NSGA-II (Non-Dominated Sorted Genetic Algorithm-II). It has been observed that the position of water table can be reliably estimated with varying periods of infiltration using developed numerical method. The parametric study is presented to study the optimum factor of safety of the embankment and its corresponding critical failure surface under the steady-state infiltration condition. Results show that in case of floating stone columns, period of infiltration has no effect on factor of safety. Even critical failure surfaces for a particular floating column length remain same irrespective of rainfall duration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Shock interaction with deformable particles using a constrained interface reinitialization scheme

In this paper, we present axisymmetric numerical simulations of shock propagation in nitromethane over an aluminum particle for post-shock pressures up to 10 GPa. We use the Mie-Gruneisen equation of state to describe both the medium and the particle. The numerical method is a finite-volume based solver on a Cartesian grid, that allows for multi-material interfaces and shocks, and uses a novel constrained reinitialization scheme to precisely preserve particle mass and volume. We compute the unsteady inviscid drag coefficient as a function of time, and show that when normalized by post-shock conditions, the maximum drag coefficient decreases with increasing post-shock pressure. We also compute the mass-averaged particle pressure and show that the observed oscillations inside the particle are on the particle-acoustic time scale. Finally, we present simplified point-particle models that can be used for macroscale simulations. In the Appendix, we extend the isothermal or isentropic assumption concerning the point-force models to non-ideal equations of state, thus justifying their use for the current problem.

Sridharan, P. (ORCID:0000000287902423)↗
DOE OSTI2016

Explicit accounting of electronic effects on the Hugoniot of porous materials

A generalized enthalpy based equation of state, which includes thermal electron excitations explicitly, is formulated from simple considerations. Its application to obtain Hugoniot of materials needs simultaneous evaluation of pressure-volume curve and temperature, the latter requiring solution of a differential equation. The errors involved in two recent papers [Huayun et al., J. Appl. Phys. 92, 5917 (2002); 92, 5924 (2002)], which employed this approach, are brought out and discussed. In addition to developing the correct set of equations, the present work also provides a numerical method to implement this approach. Constant pressure specific heat of ions and electrons and ionic enthalpy parameter, needed for applications, are calculated using a three component equation of state. The method is applied to porous Cu with different initial porosities. Comparison of results with experimental data shows good agreement. It is found that temperatures along the Hugoniot of porous materials are significantly modified due to electronic effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Magnetization of InAs parabolic quantum dot: An exact diagonalization approach

The magnetization of two electron InAs quantum dot has been studied as a function of magnetic field. The electron-electron interaction has been taken into account by using exact diagonalization method numerically. The magnetization at zero external magnetic field is zero and increases in the negative direction. There is also a paramagnetic peak where the energy levels cross from singlet state to triplet state. Finally, the magnetization falls again to even negative values and saturates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Design of binary SnO{sub 2}-CuO nanocomposite for efficient photocatalytic degradation of malachite green dye

Semiconductor mediated photocatalysis has got enormous consideration as it has shown immense potential in addressing the overall energy and environmental issues. To overcome the earlier drawbacks concerning quick charge recombination and limited visible-light absorption of semiconductor photocatalysts, numerous methods have been produced in the past couple of decades and the most broadly utilized one is to develop the photocatalytic heterojunctions. In our work, a series of SnO{sub 2}-CuO nanocomposites of different compositions were synthesized by a combustion method and have been investigated in detail by various characterization techniques, such as wide angle X-ray diffraction (XRD), UV-vis spectroscopy, transmission electron microscopy (TEM), and field emission scanning electron microscopy (FE-SEM). The results revealed that the crystal structure and optical properties of the nanocomposites were almost same for all the compositions. FE-SEM images showed that the shape of SnO{sub 2}-CuO was spherical in nature and the 1: 1 Sn/Cu sample had a well-proportioned morphology. The malachite green dye was used for the photocatalytic studies in a photoreactor and monitored with a UV-visible spectrometer for different composition ratio of metal (Sn: Cu) such as 1:1, 1:2, 2:1, 1:0.5 and 0.5:1. The 1:1 ratio nanocomposite showed excellent photocatalytic degradation of 96 % compared to pure SnO{sub 2} and CuO. The mechanism of degradation and charge separation ability of the nanocomposite are also explored using photocurrent measurement study.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Kinetic study of ion acoustic twisted waves with kappa distributed electrons

The kinetic theory of Landau damping of ion acoustic twisted modes is developed in the presence of orbital angular momentum of the helical (twisted) electric field in plasmas with kappa distributed electrons and Maxwellian ions. The perturbed distribution function and helical electric field are considered to be decomposed by Laguerre-Gaussian mode function defined in cylindrical geometry. The Vlasov-Poisson equation is obtained and solved analytically to obtain the weak damping rates of the ion acoustic twisted waves in a non-thermal plasma. The strong damping effects of ion acoustic twisted waves at low values of temperature ratio of electrons and ions are also obtained by using exact numerical method and illustrated graphically, where the weak damping wave theory fails to explain the phenomenon properly. The obtained results of Landau damping rates of the twisted ion acoustic wave are discussed at different values of azimuthal wave number and non-thermal parameter kappa for electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Singularity and Bohm criterion in hot positive ion species in the electronegative ion sources

The structure of the discharge for a magnetized electronegative ion source with two species of positive ions is investigated. The thermal motion of hot positive ions and the singularities involved with it are taken into account. By analytical solution of the neutral region, the location of the singular point and also the values of the plasma parameter such as electric potential and ion density at the singular point are obtained. A generalized Bohm criterion is recovered and discussed. In addition, for the non-neutral solution, the numerical method is used. In contrast with cold ion plasma, qualitative changes are observed. The parameter space region within which oscillations in the density and potential can be observed has been scanned and discussed. The space charge behavior in the vicinity of edge of the ion sources has also been discussed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Modulated electromagnetic fields in inhomogeneous media, hyperbolic pseudoanalytic functions, and transmutations

The time-dependent Maxwell system describing electromagnetic wave propagation in inhomogeneous isotropic media in the one-dimensional case reduces to a Vekua-type equation for bicomplex-valued functions of a hyperbolic variable, see Kravchenko and Ramirez [Adv. Appl. Cliord Algebr. 21(3), 547–559 (2011)]. Using this relation, we solve the problem of the transmission through an inhomogeneous layer of a normally incident electromagnetic time-dependent plane wave. The solution is written in terms of a pair of Darboux-associated transmutation operators [Kravchenko, V. V. and Torba, S. M., J. Phys. A: Math. Theor. 45, 075201 (2012)], and combined with the recent results on their construction [Kravchenko, V. V. and Torba, S. M., Complex Anal. Oper. Theory 9, 379-429 (2015); Kravchenko, V. V. and Torba, S. M., J. Comput. Appl. Math. 275, 1–26 (2015)] can be used for efficient computation of the transmitted modulated signals. We develop the corresponding numerical method and illustrate its performance with examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Frequency spirals

We study the dynamics of coupled phase oscillators on a two-dimensional Kuramoto lattice with periodic boundary conditions. For coupling strengths just below the transition to global phase-locking, we find localized spatiotemporal patterns that we call “frequency spirals.” These patterns cannot be seen under time averaging; they become visible only when we examine the spatial variation of the oscillators' instantaneous frequencies, where they manifest themselves as two-armed rotating spirals. In the more familiar phase representation, they appear as wobbly periodic patterns surrounding a phase vortex. Unlike the stationary phase vortices seen in magnetic spin systems, or the rotating spiral waves seen in reaction-diffusion systems, frequency spirals librate: the phases of the oscillators surrounding the central vortex move forward and then backward, executing a periodic motion with zero winding number. We construct the simplest frequency spiral and characterize its properties using analytical and numerical methods. Simulations show that frequency spirals in large lattices behave much like this simple prototype.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Fast Numerical Solution of the Plasma Response Matrix for Real-time Ideal MHD Control

To help effectuate near real-time feedback control of ideal MHD instabilities in tokamak geometries, a parallelized version of A.H. Glasser’s DCON (Direct Criterion of Newcomb) code is developed. To motivate the numerical implementation, we first solve DCON’s δW formulation with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD Riccati equation. We then describe our adaptation of DCON with numerical methods natural to solutions of the Riccati equation, parallelizing it to enable its operation in near real-time. We replace DCON’s serial integration of perturbed modes—which satisfy a singular Euler- Lagrange equation—with a domain-decomposed integration of state transition matrices. Output is shown to match results from DCON with high accuracy, and with computation time < 1s. Such computational speed may enable active feedback ideal MHD stability control, especially in plasmas whose ideal MHD equilibria evolve with inductive timescale $\tau$ ≳ 1s—as in ITER. Further potential applications of this theory are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Benchmarking semiclassical approaches to strong-field QED: Nonlinear Compton scattering in intense laser pulses

The recoil associated with photon emission is key to the dynamics of ultrarelativistic electrons in strong electromagnetic fields, as found in high-intensity laser-matter interactions and astrophysical environments such as neutron star magnetospheres. When the energy of the photon becomes comparable to that of the electron, it is necessary to use quantum electrodynamics (QED) to describe the dynamics accurately. However, computing the appropriate scattering matrix element from strong-field QED is not generally possible due to multiparticle effects and the complex structure of the electromagnetic fields. Therefore, these interactions are treated semiclassically, coupling probabilistic emission events to classical electrodynamics using rates calculated in the locally constant field approximation. Here, we provide comprehensive benchmarking of this approach against the exact QED calculation for nonlinear Compton scattering of electrons in an intense laser pulse. We find agreement at the percentage level between the photon spectra, as well as between the models' predictions of absorption from the background field, for normalized amplitudes a0 > 5. We discuss possible routes towards improved numerical methods and the implications of our results for the study of QED cascades.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A time-to-frequency converter for measuring the shape of short optical pulses

A time-to-frequency converter was constructed using an electro-optic phase modulator as a time lens, allowing the pulse shape in time to be transferred to the frequency domain. We used such a device to record the temporal shape of infrared pulses at a wavelength of 1053 nm (width about 7 ps) and compared these measurements to those made by using both a streak camera and an autocorrelator. This side-by-side comparison illustrates the benefits and limitations of each of the measurement methods. Numerical simulations were used to establish that our time-lens-based system can accurately measure the shape of infrared pulses between 3 ps and 12 ps. We also use our numerical model to determine how such a system can be modified to measure pulses whose width lies in the range of 1–30 ps, a range of interest for the OMEGA-EP laser at the Laboratory for Laser Energetics.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Dynamics of reconfigurable artificial spin ice: Toward magnonic functional materials

Over the past few years, the study of magnetization dynamics in artificial spin ices has become a vibrant field of study. Artificial spin ices are ensembles of geometrically arranged, interacting magnetic nanoislands, which display frustration by design. These were initially created to mimic the behavior in rare earth pyrochlore materials and to study emergent behavior and frustration using two-dimensional magnetic measurement techniques. Recently, it has become clear that it is possible to create artificial spin ices, which can potentially be used as functional materials. In this perspective, we review the resonant behavior of spin ices in the GHz frequency range, focusing on their potential application as magnonic crystals. In magnonic crystals, spin waves are functionalized for logic applications by means of band structure engineering. While it has been established that artificial spin ices can possess rich mode spectra, the applicability of spin ices to create magnonic crystals hinges upon their reconfigurability. Consequently, we describe recent work aiming to develop techniques and create geometries allowing full reconfigurability of the spin ice magnetic state. We also discuss experimental, theoretical, and numerical methods for determining the spectral response of artificial spin ices and give an outlook on new directions for reconfigurable spin ices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗