Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Partial differential equation”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Asymptotic vacuum solution at tokamak X-point tip

In the H-mode regime of diverted tokamaks, the presence of strong pressure gradients in the pedestal gives rise to a sizable bootstrap current, together with the Ohmic and Pfirsch–Schlueter currents, close to the separatrix. For such equilibria, the presence of finite current density close to the separatrix requires the reexamination of equilibrium properties. It is almost universally assumed that the two branches of the separatrix (the stable and unstable manifolds) are straight as they cross at the X-point. However, the opposite angles of the plasma-filled segment and vacuum one cannot be equal if the current density does not vanish at the separatrix on the plasma side. We solve this difficulty by chipping off a thin layer of plasma edge so that the sharp corner of the plasma-filled segment becomes a hyperbola. Using the conformal transformation, we found that in the assumption of a hyperbolic boundary, the X point moves beyond the plasma boundary to fall in the vacuum region. An acute angle of the plasma-filled segment leads to an obtuse opposite angle of vacuum segment and vice versa. In the case of an acute angle of the plasma-filled segment, the new X point shifts inside the X point formed by the asymptotes of a hyperbolic boundary; in the case of an obtuse angle of the plasma-filled segment, the new X point shifts outside the X point formed by the asymptotes of a hyperbolic plasma boundary. Furthermore, the results are important for understanding the X point features, which affect the tokamak edge stability and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Theory and simulations of cross-beam energy transfer between speckled laser beams

In this paper, we introduce a realistic modeling of laser beams smoothed by random phase plates to investigate cross-beam energy transfer (CBET). We compute the total field of the smoothed laser beams, which allows us to perform accurate simulations and to build a model accounting for the realistic structure of the speckles. This lets us confirm the results already discussed in a previous work given by Oudin et al. [Phys. Plasmas 29, 112112 (2022)]. When CBET is induced by a plasma flow, plane wave models can predict the exchange between smoothed laser beams. By contrast, when the exchange is caused by a wavelength shift between the laser fields of each beam, plane wave models overestimate the energy transfer, even for strongly Landau-damped ion acoustic waves. We also show that the laser angular spread, which results from its focusing, affects the resonance between the ponderomotive force induced by the laser fields and the ion acoustic wave. It increases the resonance width, which weakly depends on the Landau damping rate when there is a wavelength shift and no flow. In addition, the power transfer obtained with or without a phase plate is essentially the same in the perturbative regime and for small crossing angles.

Laser plasma interactions↗
DOE OSTI2020

2D microspatial distribution uniformity of photon detection efficiency and crosstalk probability of multi-pixel photon counters

Two-dimensional (2D) microspatial distribution uniformity of photon detection efficiency (PDE) and optical crosstalk probability P{sub ct} of multi-pixel photon counters (MPPCs) is studied. The experimental results show that the 2D spatial distribution of P{sub ct} is obviously uneven, i.e. P{sub ct} is larger at the corners and edges of a single pixel in MPPCs, which suggest a higher electrical field in the depletion region of the pixel at the corners and edges. The nonuniformity of the 2D spatial distribution of PDE also become evident when the size of the pixels of MPPCs is small, which signifies higher nonuniformity of the electric field distribution in MPPCs with small pixel size. A method is proposed for characterization of the 2D electrical field spatial distribution uniformity in a single pixel of MPPCs, which can be used for guiding the optimisation of the fabrication process of MPPCs and their properties. This promising method can naturally be extended to any Geiger avalanche photodiodes (G-APDs) and their arrays. (laser applications and other topics in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Autocrine and paracrine interferon signalling as ‘ring vaccination’ and ‘contact tracing’ strategies to suppress virus infection in a host

The innate immune response, particularly the interferon response, represents a first line of defence against viral infections. The interferon molecules produced from infected cells act through autocrine and paracrine signalling to turn host cells into an antiviral state. Although the molecular mechanisms of IFN signalling have been well characterized, how the interferon response collectively contribute to the regulation of host cells to stop or suppress viral infection during early infection remain unclear. Here, we use mathematical models to delineate the roles of the autocrine and the paracrine signalling, and show that their impacts on viral spread are dependent on how infection proceeds. In particular, we found that when infection is well-mixed, the paracrine signalling is not as effective; by contrast, when infection spreads in a spatial manner, a likely scenario during initial infection in tissue, the paracrine signalling can impede the spread of infection by decreasing the number of susceptible cells close to the site of infection. Furthermore, we argue that the interferon response can be seen as a parallel to population-level epidemic prevention strategies such as ‘contact tracing’ or ‘ring vaccination’. Thus, our results here may have implications for the outbreak control at the population scale more broadly.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Generalized Stress Inversion Approach With Application to Residual Stress Estimation

In this work develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

42 ENGINEERING↗
DOE OSTI2023

Design of InGaN-ZnSnGa 2 N 4 quantum wells for high-efficiency amber light emitting diodes

A novel type-II InGaN-ZnSnGa 2 N 4 quantum well (QW) structure is proposed based on recent experimental achievements for the successful epitaxy of ZnSnN 2 -GaN alloys and the determination of their band offsets with GaN. The simulation results indicate that this structure is promising as the active region for high-efficiency InGaN-based amber (λ ~ 590 nm) light-emitting diodes (LEDs). The hole wavefunction in the valence band is better confined with the insertion of a monolayer scale of ZnSnGa 2 N 4 into the InGaN QW while the electron wavefunction in the conduction band is better confined with the incorporation of an AlGaN layer in the GaN quantum barrier. The band structure of the InGaN-ZnSnGa 2 N 4 QW is numerically simulated based on the experimentally measured band offsets between ZnSnGa 2 N 4 and GaN. With the InGaN-ZnSnGa 2 N 4 QW design, a low In content (20%) is required in the InGaN layer to reach a peak emission wavelength of ~590 nm, yet an In composition of 25% is needed to reach the same emission wavelength for a conventional InGaN QW with the same layer thicknesses. Moreover, the electron-hole wavefunction overlap (Гe1-hh1) for the InGaN-ZnSnGa 2 N 4 QW design reaches 18% for an emission wavelength at ~590 nm. This result is much improved over the conventional InGaN QW overlap of 5% emitting at the same wavelength. The increase in electron-hole wavefunction overlap results in an approximately 14 times enhancement in the predicted spontaneous emission radiative recombination rate of the InGaN-ZnSnGa 2 N 4 QW as compared to that of the conventional InGaN QW. This InGaN-ZnSnGa 2 N 4 QW structure design can be promising to pave a new way to achieve high efficiency amber LEDs.

42 ENGINEERING↗
DOE OSTI2022

Voronoi cell analysis: The shapes of particle systems

Many physical systems can be studied as collections of particles embedded in space, often evolving in time. Natural questions arise concerning how to characterize these arrangements—are they ordered or disordered? If they are ordered, how are they ordered and what kinds of defects do they possess? Voronoi tessellations, originally introduced to study problems in pure mathematics, have become a powerful and versatile tool for analyzing countless problems in pure and applied physics. Here we explain the basics of Voronoi tessellations and the shapes that they produce and describe how they can be used to characterize many physical systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Computing Edge States without Hard Truncation

We present a numerical method which accurately computes the discrete spectrum and associated bound states of semi-infinite Hamiltonians which model electronic “edge” states localized at boundaries of one- and two-dimensional crystalline materials. The problem is nontrivial since arbitrarily large finite “hard” (Dirichlet) truncations of the Hamiltonian in the infinite bulk direction tend to produce spurious bound states partially supported at the truncation. Our method, which overcomes this difficulty, is to compute the Green's function of the semi-infinite Hamiltonian by imposing an appropriate boundary condition in the bulk direction; then, the spectral data is recovered via Riesz projection. We demonstrate our method's effectiveness by studies of edge states at a graphene zig-zag edge in the presence of defects modeled both by a discrete tight-binding model and a continuum PDE model under finite difference discretization. Our method may also be used to study states localized at domain wall-type edges in one- and two-dimensional materials where the edge Hamiltonian is infinite in both directions; we demonstrate this for the case of a tight-binding model of distinct honeycomb structures joined along a zig-zag edge. Here, we expect our method to be useful for designing novel devices based on precise wave-guiding by edge states.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Extended Galerkin Neural Network Approximation of Singular Variational Problems with Error Control

We present extended Galerkin neural networks, a variational framework for approximating general boundary value problems (BVPs) with error control. The main contributions of this work are (1) a rigorous theory guiding the construction of new weighted least squares variational formulations suitable for use in neural network approximation of general BVPs, and (2) an “extended” feedforward network architecture which incorporates and is even capable of learning singular solution structures, thus greatly improving approximability of singular solutions. Furthermore, numerical results are presented for several problems, including steady Stokes flow around reentrant corners and in convex corners with Moffatt eddies in order to demonstrate efficacy of the method.

a posteriori error estimate↗
DOE OSTI2022

Bringing Trimmed Serendipity Methods to Computational Practice in Firedrake

We present an implementation of the trimmed serendipity finite element family, using the open-source finite element package Firedrake. The new elements can be used seamlessly within the software suite for problems requiring H 1 , H (curl), or H (div)-conforming elements on meshes of squares or cubes. To test how well trimmed serendipity elements perform in comparison to traditional tensor product elements, we perform a sequence of numerical experiments including the primal Poisson, mixed Poisson, and Maxwell cavity eigenvalue problems. Overall, we find that the trimmed serendipity elements converge, as expected, at the same rate as the respective tensor product elements, while being able to offer significant savings in the time or memory required to solve certain problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗
DOE OSTI2026

Intrusive Uncertainty Quantification and Optimal Experiment Design in the Open-Source Pyomo Ecosystem

This contribution describes ParmEst and Pyomo.DoE, two pillars of the open-source Python-based Pyomo ecosystem for computational optimization with (partial differential) algebraic equation mathematical models. Specifically, ParmEst facilitates intrusive frequentist parameter estimation (PE) and uncertainty quantification (UQ) through built-in features, such as covariance matrix estimation, bootstrapping, and likelihood ratio tests. Complementary, Pyomo.DoE enables optimal experiment design by maximizing various metrics of the Fisher information matrix, such as A-optimality (trace), D-optimality (determinant), E-optimality (minimum eigenvalue), and ME-optimality (condition number). ParmEst and Pyomo.DoE can solve high-dimensional optimization problems by leveraging the model structure and exact derivative information. Finally, we will discuss future opportunities to integrate PE and UQ capabilities with optimization under uncertainty, including robust optimization with non-convex models via PyROS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Pyomo.DOE: An open-source package for model-based design of experiments in Python

Predictive mathematical models are a cornerstone of science and engineering. Yet selecting, calibrating, and validating said science-based models often remains an art in practice. Model-based design of experiments (MBDoE) provides a systematic framework to maximize information gain from experiments while minimizing time and resource costs. But MBDoE remains limited to niche application areas, in part because practitioners must integrate expertise in statistics, computational optimization, and modeling. To help reduce this barrier, we introduce Pyomo.DOE, an open-source package for MBDoE. Pyomo.DOE uses a nonlinear sensitivity analysis code k_aug to quickly approximate the Fisher information matrix and leverages a new stochastic programming abstraction. We demonstrate Pyomo.DOE with the first application of MBDoE to fixed-bed breakthrough experiments, which highlights the power of Pyomo.DOE to quantify the value of experimental modifications a priori for large-scale partial differential-algebraic equation (PDAE) models. Here we also provide a mathematical primer on MBDoE targeted at general chemical engineers.

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

BMINN: Learning chemical potentials and parameters from voltage data for multi-phase battery modeling

Free-energy landscapes and chemical potentials govern the dynamics of phase transitions, transport, and stability in functional materials, yet they remain experimentally inaccessible under realistic operating conditions. Here we introduce a Bayesian model-integrated neural network (BMINN) that embeds physics-based formulations of non-autonomous partial differential-algebraic equations into probabilistic learning. This approach reconstructs hidden thermodynamics directly from macroscopic current-voltage data, providing quantitative access to metastable states, staging transitions, and energy barriers without synchrotron probes. Demonstrated on lithium-graphite electrodes, BMINN recovers full Gibbs free-energy landscapes with fidelity validated against operando X-ray diffraction. The framework generalizes across dynamical regimes, enabling accurate voltage prediction, internal state estimation, and inference of governing parameters. Beyond batteries, BMINN exemplifies a broadly applicable strategy for learning missing physics in multiphase, non-equilibrium systems, offering a new pathway to uncover hidden thermodynamic functions across condensed matter and materials physics.

25 ENERGY STORAGE↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗