Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

BLOC Site - NOAA PSL Wind Retrievals WINDoe / Derived Data

This dataset contains daily netcdf files with horizontal wind profiles retrieved with the WINDoe retrieval (Gebauer and Bell 2024) at Block Island (BLOC). WINDoe retrievals datasets are also available at Nantucket Island (NANT, nant.windoe.z01.c1) and Rhode Island (RHOD, rhod.windoe.z01.c1). WINDoe is an optimal estimation algorithm to retrieve wind profiles combining multiple instruments. The code is available in this github repository (https://github.com/OAR-atmospheric-observations/WINDoe/tree/main), and the retrieval is described by Gebauer and Bell (2024). WINDoe allows combining the individual datasets and outputs into one profile taking into account the information and uncertainties of each dataset. The use of WINDoe minimizes data gaps and maximizes data availability, compared to using wind profiles from only one of the instruments. The regular height grid eases comparisons to numerical weather prediction models. Code modifications have been made that include reading in WFIP3 specific instruments, averaging Doppler lidar radial velocities at various azimuth angles to avoid overfitting, and allowing the user to define a height grid by the user in the vipfile. The instruments used as input to the retrieval are a radar wind profiler (low- and high resolution mode) providing data in and above the boundary layer, a scanning Doppler lidar usually providing data throughout the boundary layer, a profiling lidar providing data from 50 to 200 m at BLOC and NANT, and from 10 to 280 m at Rhode Island, and a surface tower (4 m at NANT and RHOD and 10 m at BLOC). From the scanning lidars, we used radial velocity measurements at 60 deg elevation angle at six different azimuth angles with a resolution of approximately 30 m along the line of sight and the lowest range gate at approximately 70 m. The wind profiles are retrieved with WINDoe up to 3.74 km with 10 m vertical resolution. The profiles are retrieved every 15 min at BLOC and NANT and every 60 min at RHOD.

17 WIND ENERGY↗
DOE OSTI2026

NANT Site - NOAA PSL Wind Retrievals WINDoe / Derived Data

This dataset contains daily NetCDF files with horizontal wind profiles retrieved with the WINDoe retrieval (Gebauer and Bell 2024) at Nantucket Island (NANT). WINDoe retrievals datasets are also available at Block Island (BLOC, bloc.windoe.z01.c1) and Rhode Island (RHOD, rhod.windoe.z01.c1). WINDoe is an optimal estimation algorithm to retrieve wind profiles combining multiple instruments. The code is available in this github repository (https://github.com/OAR-atmospheric-observations/WINDoe/tree/main), and the retrieval is described by Gebauer and Bell (2024). WINDoe allows combining the individual datasets and outputs into one profile taking into account the information and uncertainties of each dataset. The use of WINDoe minimizes data gaps and maximizes data availability, compared to using wind profiles from only one of the instruments. The regular height grid eases comparisons to numerical weather prediction models. Code modifications have been made that include reading in WFIP3 specific instruments, averaging Doppler lidar radial velocities at various azimuth angles to avoid overfitting, and allowing the user to define a height grid by the user in the vipfile. The instruments used as input to the retrieval are a radar wind profiler (low- and high resolution mode) providing data in and above the boundary layer, a scanning Doppler lidar usually providing data throughout the boundary layer, a profiling lidar providing data from 50 to 200 m at BLOC and NANT, and from 10 to 280 m at Rhode Island, and a surface tower (4 m at NANT and RHOD and 10 m at BLOC). From the scanning lidars, we used radial velocity measurements at 60 deg elevation angle at six different azimuth angles with a resolution of approximately 30 m along the line of sight and the lowest range gate at approximately 70 m. The wind profiles are retrieved with WINDoe up to 3.74 km with 10 m vertical resolution. The profiles are retrieved every 15 min at BLOC and NANT and every 60 min at RHOD.

17 WIND ENERGY↗
DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗
DOE OSTI2022

Chromosome aberrations among atomic-bomb survivors exposed in utero: updated analysis accounting for revised radiation doses and smoking

A previous study of peripheral blood lymphocyte translocations around age 40 among atomic-bomb survivors exposed in utero revealed no overall association with radiation dose—despite a clear association between translocations and dose among their mothers—but the data suggested an increase at doses below 100 mGy with a definite peak. That analysis of the in utero-exposed survivors did not adjust for their subsequent smoking behavior, an established cause of chromosomal aberrations, or their subsequent exposures to medical irradiation, a potential mediator. In addition, atomic-bomb survivor radiation dose estimates have subsequently been updated and refined. We therefore re-estimated the dose response using the latest DS02R1 dose estimates and adjusting for smoking as well as for city and proximal–distal location at the time of exposure to the atomic bomb. Sex of the survivor, mother’s age around the time of conception, and approximate trimester of gestation at the time of exposure were also considered as explanatory variables and modifiers. Precision of the estimated dose response was slightly lower due to greater variability near zero in the updated dose estimates, but there was little change in evidence of a low-dose increase and still no suggestion of an overall increase across the entire dose range. Adjustment for smoking behavior led to a decline in background number of translocations (the dose–response intercept), but smoking did not interact with dose overall (across the entire dose range). Adjustment for medical irradiation did not alter the association between dose and translocation frequency. Sex, mother’s age, and trimester were not associated with number of translocations, nor did they interact with dose overall. Interactions with dose in the low-dose range could not be evaluated because of numerical instability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Quantum state reduction: Generalized bipartitions from algebras of observables

Reduced density matrices are a powerful tool in the analysis of entanglement structure, approximate or coarse-grained dynamics, decoherence, and the emergence of classicality. It is straightforward to produce a reduced density matrix with the partial-trace map by tracing out part of the quantum state, but in many natural situations this reduction may not be achievable. We investigate the general problem of identifying how the quantum state is reduced given a restriction on the observables. For example, in an experimental setting, the set of observables that can actually be measured is usually modest (compared to the set of all possible observables) and their resolution is limited. In such situations, the appropriate state-reduction map can be defined via a generalized bipartition, which is associated with the structure of irreducible representations of the algebra generated by the restricted set of observables. One of our main technical results is a general, not inherently numeric, algorithm for finding irreducible representations of matrix algebras. In our work, we demonstrate the viability of this approach with two examples of limited-resolution observables. The definition of quantum state reductions can also be extended beyond algebras of observables. To accomplish this task we introduce a more flexible notion of bipartition, the partial bipartition, which describes coarse grainings preserving information about a limited set (not necessarily algebra) of observables. We describe a variational method to choose the coarse grainings most compatible with a specified Hamiltonian, which exhibit emergent classicality in the reduced state space. We apply this construction to the concrete example of the one-dimensional Ising model. Our results have relevance for quantum information, bulk reconstruction in holography, and quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗
DOE OSTI2023

Raman Wavelength Conversion in Ionic Liquids

Here, we explore the use of room-temperature ionic liquids (ILs) as Raman wavelength converters. ILs provide an engineerable framework to design suitable liquids for wavelength conversion over a broad spectral range, through careful selection of the molecular structures of the IL anions and cations so that specific characteristics can be obtained, such as a desirable Raman shift, low Brillouin scattering, and good optical transmission in the pump and Stokes wavelengths. Applying such criteria, we demonstrate that 1-ethyl-3-methylimidazolium dicyanamide (EMIM DCA) is an effective medium for conversion of 532-nm pulses from a Q-switched Nd :YAG laser to 603 nm. This corresponds to an approximate 2200 cm –1 shift, which can be used to generate mid-infrared radiation through subsequent difference frequency generation for optical pumping of CO 2 lasers. Threefold-higher Raman conversion efficiency is obtained in EMIM DCA compared with water under identical conditions in a proof-of-principle single-pass conversion setup, resulting in an efficient generation of multimillijoule, <6 ns duration, high-quality orange laser pulses in a wavelength region that is difficult to access at high energies. Consequently, we examine ILs representing two other classes of Raman-active functional groups and obtain conversion up to the fifth-order Stokes shift and first anti-Stokes shift. Through the tunable selection of their components and their useful dynamical properties, ILs provide a platform for efficient, simple, and alignment-tolerant high-energy Raman shifting with numerous industrial and technological applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Spectral properties and enhanced superconductivity in renormalized Migdal-Eliashberg theory

Migdal-Eliashberg theory describes the properties of the normal and superconducting states of electron-phonon-mediated superconductors based on a perturbative treatment of the electron-phonon interactions. It is necessary to include both electron and phonon self-energies self-consistently in Migdal-Eliashberg theory in order to match numerically exact results from determinantal quantum Monte Carlo in the adiabatic limit. Here in this work we provide a method to obtain the real-axis solutions of the Migdal-Eliashberg equations with electron and phonon self-energies calculated self-consistently. Our method avoids the typical challenge of computing cumbersome singular integrals on the real axis and is numerically stable and exhibits fast convergence. Analyzing the resulting real-frequency spectra and self-energies of the two-dimensional Holstein model, we find that self-consistently including the lowest-order correction to the phonon self-energy significantly affects the solution of the Migdal-Eliashberg equations. The calculation captures the broadness of the spectral function, renormalization of the phonon dispersion, enhanced effective electron-phonon coupling strength, minimal increase in the electron effective mass, and the enhancement of superconductivity which manifests as a superconducting ground state despite strong competition with charge-density-wave order. We discuss surprising differences in two common definitions of the electron-phonon coupling strength derived from the electron mass and the density of states, quantities which are accessible through experiments such as angle-resolved photoemission spectroscopy and electron tunneling. An approximate upper bound on 2Δ/T c for conventional superconductors mediated by retarded electron-phonon interactions is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Triangular lattice Hubbard model physics at intermediate temperatures

Moiré systems offer an exciting playground to study the many-body effects of strongly correlated electrons in regimes that are not easily accessible in conventional material settings. Motivated by a recent experiment on WSe 2 /WS 2 moiré bilayers [Y. Tang et al., Nature (London) 579, 353 (2020)], which realized a triangular superlattice with a small hopping (of approximately 10 Kelvin), with tunable density of holes, we explore the Hubbard model on the triangular lattice for intermediate temperatures t ≲ T < U. Employing finite-temperature Lanczos calculations and closely following the fitting protocols used in the experiment, we recover the observed trends in the reported Curie-Weiss temperature Θ with filling, using the reported interaction strength U/t=20. Here, we focus on the large increase of |Θ| on decreasing the density below half filling and the sign change of Θ at higher fillings, which signals the onset of ferromagnetism. The increase in |Θ| is also seen in the t–J model (the low-energy limit of the Hubbard model) in the intermediate temperature range, which we clarify is opposite to the trend in its high temperature limit. Differences between the low, intermediate, and high temperature regimes are discussed. Our numerical calculations also capture the crossover between short-range antiferromagnetic to ferromagnetic order in the intermediate temperature regime, a result broadly consistent with the experimental findings. We find that this behavior is a finite-temperature remnant of the underlying zero-temperature phase transition, which we explore with ground-state density matrix renormalization group calculations. We provide evidence of ferromagnetism characterized by weak (but robust) correlations that explain the small Θ seen in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Infrared thermography NDT for in-situ defect detection in sandwich composite panel manufacturing

Composite manufacturing presents numerous challenges, as defects can arise from various sources throughout the process. In sandwich composite structures, the integration of a foam core introduces additional complexity and increases the likelihood of defect formation like delamination. To mitigate these issues and reduce the risk of future structural failures, in-situ monitoring during manufacturing is essential. This study investigates infrared (IR) thermography as a non-destructive technique for detecting manufacturing defects in foam-core sandwich composite panels under thermally excited conditions representative of in-situ processing. A stationary FLIR A8590 IR camera (640 × 512 pixels, 30Hz, 17mm lens, 9 ft stand-off distance) was used to monitor prefabricated panels subjected to controlled external heating simulating compression molding and resin cure exotherm. Interlaminar delamination defects with characteristic sizes ranging from 0.25 × 0.25in² to 5 × 5in² produced measurable surface temperature depressions of approximately 4–10°C during transient cooling, exceeding the effective noise floor of the camera by more than two standard deviations. Thicker laminates exhibited prolonged defect detectability windows due to increased thermal diffusion time. In contrast, embedded Teflon inclusions generated weak thermal contrasts of ≤ 3°C, approaching the measurement noise floor, due to limited thermal property contrast with the surrounding glass fiber composite. These results establish quantitative detectability limits for stationary thermographic inspection of sandwich composite panels under manufacturing-representative thermal cycles.

Barakat, Abdallah [ORNL] (ORCID:0000000296141398)↗
DOE OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{&#x221E;}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Stage-local partitioned two-step runge-kutta methods for large systems of ordinary differential equations

We introduce stage-local partitioned two-step Runge-Kutta methods are an extension of standard two-step Runge-Kutta methods, which are an alternative to the standard additive two-step Runge-Kutta methods currently existing in the literature. Furthermore, these new schemes are designed with an eye towards truly N-partitioned systems and leverage local stage approximations to make several computationally interesting approximations viable. Specifically, the focus on local stage approximations makes possible the construction of truly asynchronous schemes, in the parallel sense, possible. In addition, we show that an implicit-explicit approach to these schemes can lead to methods that require the inversion of only local nonlinear systems.

Applied Dynamical Systems↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sensitivity of the Shallow‐To‐Deep Convective Transition to Moisture and Wind Shear in the Amazon

Deep convection is the primary influence on weather and climate in tropical regions. However, understanding and simulating the shallow‐to‐deep (STD) convective transition has long been challenging. Here, we conduct high‐resolution numerical simulations to assess the environmental controls on the evolution of isolated convection in the Amazon during the wet season. The large‐scale forcing derived through a constrained variational analysis approach for the GoAmazon2014/5 Experiment is used in the simulations. Through sensitivity experiments, we examine the relative importance of moisture and wind shear in controlling the shallow‐to‐deep convective transition for isolated convective events. Convection exhibits the greatest sensitivity to humidity within the lowest 1.5 km, where a 4 mm reduction in column water vapor nearly suppresses ice water formation on deep convective days. In contrast, a reduction in column water vapor in the free troposphere by a factor of two or more is necessary to produce a comparable impact on convection. Increasing low‐level wind speed from 6 to 9 m s -1 enhances afternoon deep convection, raising the cloud ice mixing ratio by approximately 25%. Conversely, upper‐level wind shear reveals the weakest correlation with daytime convection in our simulations. Our results help characterize the role of moisture and wind shear on the STD transition and our understanding of the underlying mechanisms.

environmental controls↗
DOE OSTI2024

Towards dynamical low-rank approximation for neutrino kinetic equations. Part I: Analysis of an idealized relaxation model

Dynamical low-rank approximation (DLRA) is an emerging tool for reducing computational costs and provides memory savings when solving high-dimensional problems. Here, in this work, we propose and analyze a semi-implicit dynamical low-rank discontinuous Galerkin (DLR-DG) method for the space homogeneous kinetic equation with a relaxation operator, modeling the emission and absorption of particles by a background medium. Both DLRA and the discontinuous Galerkin (DG) scheme can be formulated as Galerkin equations. To ensure their consistency, a weighted DLRA is introduced so that the resulting DLR-DG solution is a solution to the fully discrete DG scheme in a subspace of the standard DG solution space. Similar to the standard DG method, we show that the proposed DLR-DG method is well-posed. We also identify conditions such that the DLR-DG solution converges to the equilibrium. Numerical results are presented to demonstrate the theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗