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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,747 records · Page 97

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A light weight regularization for wave function parameter gradients in quantum Monte Carlo

The parameter derivative of the expectation value of the energy, ∂E/∂p, is a key ingredient in variational Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve estimate of this derivative suffers from an infinite variance, which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator, which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias, which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of ∂E/∂p for a multi-Slater–Jastrow trial wave function on the LiH molecule and in the optimization of a multi-Slater–Jastrow trial wave function on the CuO molecule. This regularized estimator is a simple and efficient estimator of ∂E/∂p for VMC optimization techniques.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Super-resolution surface slope metrology of x-ray mirrors

We present experimental, analytical, and numerical methods developed for reconstruction (deconvolution) of one-dimensional (1D) surface slope profiles over the spatial frequency range where the raw data are significantly perturbed due to the limited resolution of the measurement instrument. We characterize the spatial resolution properties of a profiler with the instrument's transfer function (ITF). To precisely measure the ITF, we apply a recently developed method utilizing test surfaces with 1D linear chirped height profiles of constant slope amplitude. Based on the results of the ITF calibration, we determine parameters of an analytical model for the ITF that is used in the original reconstruction software. Here, we treat surface slope metrology data obtained with the Optical Surface Measuring System (OSMS), using as a sensor an electronic autocollimator (AC) ELCOMAT-3000. The spatial resolution of the OSMS is limited by the size of the AC light-beam-collimating aperture. For the purposes of this investigation, the OSMS is equipped with a circular aperture with a diameter of 2.5 mm. This is a typical arrangement of most AC-based slope profilers developed for surface slope metrology of state-of-the-art x-ray mirrors. Using the example of surface slope metrology of two state-of-the-art elliptically shaped x-ray focusing mirrors, we demonstrate that the developed data reconstruction procedure allows us to significantly improve the accuracy of surface slope metrology with the OSMS over the spatial wavelength range from ~1.6 mm to 7 mm. Thus, the amplitude of the quasi-periodic error characteristic of the deterministic polishing process used appears to be higher by a factor of~2 than is apparent from the rough metrology data. Underestimation of the surface slope errors in this spatial wavelength range can lead to serious errors in the expected performance of x-ray mirrors in synchrotron beamlines, especially at modern light sources utilizing coherent x rays, where the perturbations can lead to increased speckle-like intensity variation.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Facile and quantitative estimation of strain in nanobubbles with arbitrary symmetry in 2D semiconductors verified using hyperspectral nano-optical imaging

When layers of van der Waals materials are deposited via exfoliation or viscoelastic stamping, nanobubbles are sometimes created from aggregated trapped fluids. Though they can be considered a nuisance, nanobubbles have attracted scientific interest in their own right owing to their ability to generate large in-plane strain gradients that lead to rich optoelectronic phenomena, especially in the semiconducting transition metal dichalcogenides. Determination of the strain within the nanobubbles, which is crucial to understanding these effects, can be approximated using elasticity theory. However, the Föppl–von Kármán equations that describe strain in a distorted thin plate are highly nonlinear and often necessitate assuming circular symmetry to achieve an analytical solution. Here, we present an easily implemented numerical method to solve for strain tensors of nanobubbles with arbitrary symmetry in 2D crystals. The method only requires topographic information from atomic force microscopy and the Poisson ratio of the 2D material. We verify that this method reproduces the strain for circularly symmetric nanobubbles that have known analytical solutions. Finally, we use the method to reproduce the Grüneisen parameter of the E' mode for 1L-WS2 nanobubbles on template-stripped Au by comparing the derived strain with measured Raman shifts from tip-enhanced Raman spectroscopy, demonstrating the utility of our method for estimating localized strain in 2D crystals.

2D materials↗

DOE OSTI2020

Sneezing and asymptomatic virus transmission

The novel coronavirus disease (COVID-19) spread pattern continues to show that geographical barriers alone cannot contain a virus. Asymptomatic carriers play a critical role in the nature of this virus quickly escalating into a global pandemic. Asymptomatic carriers may transmit the virus unintentionally through sporadic sneezing. A novel Computational Fluid Dynamics (CFD) approach has been proposed with a realistic modeling of a human sneeze achieved by the combination of state-of-the-art experimental and numerical methods. This modeling approach may be suitable for future engineering analyses aimed at reshaping public spaces and common areas, with the main objective to accurately predict the spread of aerosol and droplets that may contain pathogens. This study shows that the biomechanics of a human sneeze, including complex muscle contractions and relaxations, can be accurately modeled by the angular head motion and the dynamic pressure response during sneezing. These have been considered as the human factors and were implemented in the CFD simulation by imposing a momentum source term to the coupled Eulerian–Lagrangian momentum equations. The momentum source was modeled by the measured dynamic pressure response in conjunction with the angular head motion. This approach eliminated the need to create an ad hoc set of inlet boundary conditions. With this proposed technique, it is easier to add multiple fixed and/or moving sources of sneezes in complex computational domains. Additionally, extensive sensitivity analyses based on different environmental conditions were performed, and their impact was described in terms of potential virus spread.

42 ENGINEERING↗

DOE OSTI2020

Frequency-dependent sensitivity of AC elastocaloric effect measurements explored through analytical and numerical models

We present a comprehensive study of the frequency-dependent sensitivity for measurements of the AC elastocaloric effect by applying both exactly soluble models and numerical methods to the oscillating heat flow problem. These models reproduce the finer details of the thermal transfer functions observed in experiments, considering here representative data for single-crystal Ba(Fe1−xCox)2As2. Based on our results, we propose a set of practical guidelines for experimentalists using this technique. This work establishes a baseline against which the frequency response of the AC elastocaloric technique can be compared and provides intuitive explanations of the detailed structure observed in experiments.

Straquadine, J. A. W. (ORCID:0000000165309541)↗

DOE OSTI2021

First coupled GENE–XGC microturbulence simulations

Covering the core and the edge region of a tokamak, respectively, the two gyrokinetic turbulence codes Gyrokinetic Electromagnetic Numerical Experiment (GENE) and X-point Gyrokinetic Code (XGC) have been successfully coupled by exchanging three-dimensional charge density data needed to solve the gyrokinetic Poisson equation over the entire spatial domain. Certain challenges for the coupling procedure arise from the fact that the two codes employ completely different numerical methods. This includes, in particular, the necessity to introduce mapping procedures for the transfer of data between the unstructured triangular mesh of XGC and the logically rectangular grid (in a combination of real and Fourier space) used by GENE. Constraints on the coupling scheme are also imposed by the use of different time integrators. First, coupled simulations are presented. We have considered collisionless ion temperature gradient turbulence, in both circular and fully shaped plasmas. Coupled simulations successfully reproduce both GENE and XGC reference results, confirming the validity of the code coupling approach toward a whole device model. Here, many lessons learned in the present context, in particular, the need for a coupling procedure as flexible as possible, should be valuable to our and other efforts to couple different kinds of codes in pursuit of a more comprehensive description of complex real-world systems and will drive our further developments of a whole device model for fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Designing high-power, octave spanning entangled photon sources for quantum spectroscopy

Entangled photon spectroscopy is a nascent field that has important implications for measurement and imaging across chemical, biology, and materials fields. Entangled photon spectroscopy potentially offers improved spatial and temporal-frequency resolutions, increased cross sections for multiphoton and nonlinear measurements, and new abilities in inducing or measuring quantum correlations. A critical step in enabling entangled photon spectroscopies is the creation of high-flux entangled sources that can use conventional detectors as well as provide redundancy for the losses in realistic samples. Here, we report a periodically poled, chirped, lithium tantalate platform that generates entangled photon pairs with ~10 -7 efficiency. For a near watt level diode laser, this results in a near μW-level flux. The single photon per mode limit that is necessary to maintain non-classical photon behavior is still satisfied by distributing this power over up to an octave-spanning bandwidth. The spectral–temporal photon correlations are observed via a Michelson-type interferometer that measures the broadband Hong–Ou–Mandel two-photon interference. A coherence time of 245 fs for a 10 nm bandwidth in the collinear case and a coherence time of 62 fs for a 125 nm bandwidth in the non-collinear case are measured using a CW pump laser and, essentially, collecting the full photon cone. We outline in detail the numerical methods used for designing and tailoring the entangled photons source, such as changing center wavelength or bandwidth, with the ultimate aim of increasing the availability of high-flux UV–Vis entangled photon sources in the optical spectroscopy community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Learning the structure of wind: A data-driven nonlocal turbulence model for the atmospheric boundary layer

In this work, we develop a novel data-driven approach to modeling the atmospheric boundary layer. This approach leads to a nonlocal, anisotropic synthetic turbulence model which we refer to as the deep rapid distortion (DRD) model. Our approach relies on an operator regression problem that characterizes the best fitting candidate in a general family of nonlocal covariance kernels parameterized in part by a neural network. This family of covariance kernels is expressed in Fourier space and is obtained from approximate solutions to the Navier–Stokes equations at very high Reynolds numbers. Each member of the family incorporates important physical properties such as mass conservation and a realistic energy cascade. The DRD model can be calibrated with noisy data from field experiments. After calibration, the model can be used to generate synthetic turbulent velocity fields. To this end, we provide a new numerical method based on domain decomposition which delivers scalable, memory-efficient turbulence generation with the DRD model as well as others. We demonstrate the robustness of our approach with both filtered and noisy data coming from the 1968 Air Force Cambridge Research Laboratory Kansas experiments. Using these data, we witness exceptional accuracy with the DRD model, especially when compared to the International Electrotechnical Commission standard.

17 WIND ENERGY↗

DOE OSTI2022

Simulating a pulsed-power-driven plasma with ideal MHD

We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the magnetohydrodynamics (MHD) domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This much more affordable ideal MHD representation is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Anion diffusion in two-dimensional halide perovskites

Commercialization of halide perovskites in the semiconductor industry is hindered by their short-term stability. The instability of perovskites is closely interlinked with ionic diffusion. Historically, attempts to study diffusion in 2D perovskites mostly utilized electrical characterizations, but these characterizations pose a challenge in deconvoluting the impact of device architecture, interlayers, and ionic species. In this Perspective, we focus our attention on simple optical characterizations employed in the literature to investigate halide diffusion in 2D perovskites using lateral and vertical heterostructure platforms. We review the various synthesis techniques used for fabrication of halide perovskite heterostructures and discuss the qualitative and quantitative diffusion studies performed using these platforms. We discuss the numerical methods used to validate and supplement the experimental halide diffusion kinetics. Finally, we highlight the need to conduct further research on the impact of device operating conditions, lattice structure, and vacancy concentration on halide diffusion. Through this Perspective, we aim to emphasize the need of developing a comprehensive understanding of halide diffusion in perovskites for their successful deployment in optoelectronics.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Δ, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. In this work, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Δ. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Δ variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfvén wave resonances as rotation slows down. These are important features that can alter the non-axisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Preface for the Sherwood Fusion Theory 2021–2022 special collection

The Sherwood Fusion Theory conference has a 60-year history and has been held nearly annually for at least five decades. With theoretical and computational research directed toward fusion energy as an organizing principle, the typical Sherwood program encompasses a wide range of plasma physics topics pertaining to a variety of confinement concepts along with advancements in algorithms and numerical methods for fusion plasma simulations. Additionally, the typical annual meeting comprises roughly a dozen invited talks plus a small number of plenary speakers and several poster sessions. The program committee presents several awards each year to graduate (or undergraduate) students with outstanding poster presentations. Beginning in 2021, both invited speakers and student poster award winners were asked to submit manuscripts for this special collection, which combines contributions from the 2021 and 2022 Sherwood meetings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Superconductivity found in meteorites

Meteorites can contain a wide range of material phases due to the extreme environments found in space and are ideal candidates to search for natural superconductivity. However, meteorites are chemically inhomogeneous, and superconducting phases in them could potentially be minute, rendering detection of these phases difficult. To alleviate this difficulty, we have studied meteorite samples with the ultrasensitive magnetic field modulated microwave spectroscopy (MFMMS) technique [J. G. Ramírez, A. C. Basaran, J. de la Venta, J. Pereiro, I. K. Schuller, Rep. Prog. Phys. 77, 093902 (2014)]. Here, we report the identification of superconducting phases in two meteorites, Mundrabilla, a group IAB iron meteorite [R. Wilson, A. Cooney, Nature 213, 274–275 (1967)] and GRA 95205, a ureilite [J. N. Grossman, Meteorit. Planet. Sci. 33, A221–A239 (1998)]. MFMMS measurements detected superconducting transitions in samples from each, above 5 K. By subdividing and remeasuring individual samples, grains containing the largest superconducting fraction were isolated. The superconducting grains were then characterized with a series of complementary techniques, including vibrating-sample magnetometry (VSM), energy-dispersive X-ray spectroscopy (EDX), and numerical methods. Finally, these measurements and analysis identified the likely phases as alloys of lead, indium, and tin.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

SAM: A Modern System Code for Advanced Non-LWR Safety Analysis

The System Analysis Module (SAM), developed at Argonne National Laboratory and by collaborators at other organizations, is for advanced non–light water reactor safety analysis. SAM aims to provide fast-running, modest-fidelity, whole-plant transient analysis capabilities that are essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. To facilitate code development, SAM utilizes the MOOSE object-oriented application framework, its underlying finite element library, and linear and nonlinear solvers to leverage modern advanced software environments and numerical methods. SAM aims to solve tightly coupled physical phenomena, including fission reaction, heat transfer, fluid dynamics, and thermal-mechanical responses in advanced reactor structures, systems, and components with high accuracy and efficiency. Finally, this paper gives an overview of the SAM code development, including goals and functional requirements, physical models, current capabilities, verification and validation, software quality assurance, and examples of simulations for advanced nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗