Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential mode”

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 37 records · Page 2

Using X-ray radiography to study oxygen flow in a proton exchange membrane electrolyzer operating under balanced pressure conditions

Of the various water electrolyzer technologies, the proton exchange membrane electrolyzer (PEMWE) is one of the best solutions for producing clean hydrogen without releasing CO 2 . In order to allow for widespread use of clean hydrogen, it is necessary to decrease its cost, which is intrinsically related to system operation. Current PEMWE plants operate in differential mode, directly pressurizing hydrogen and benefiting from thermodynamic compression, which increases overall system efficiency. However, high differential pressure above 30 bar can cause membrane stress, resulting in membrane creeping and failure. Pressurizing the water and operating at balanced pressure allows hydrogen to be produced at higher pressures while preserving the integrity of the membrane and porous layers. Nevertheless, the impact of pressurizing water on PEMWE performance must be better understood to maximize performance under balanced pressure conditions. Here, this study examined the impact of water pressure on two-phase flow. A high-pressure electrolyzer setup was developed to perform operando X-ray radiography and examine oxygen transport with high temporal resolution. The imaging segmentation process, developed to capture bubble properties in the channel, was applied to a specific experiment. The results clearly showed that as pressure increased up to 30 bars, the initial bubbly flow transitioned to slug flow, which led to channel saturation with oxygen. This work demonstrates that two-phase flow in an electrolyzer can be studied using X-ray radiography, which has the advantages of fast measurements and the ability to probe dense materials, such as those required for pressurized electrolyzers.

Balanced high pressure operation↗

Minor Elements and Solidification Cracking During Laser Powder-Bed Fusion of a High $\gamma ^{\prime }$ CoNi-Base Superalloy

Here, the cracking behavior of a high γ' volume fraction CoNi-base superalloy fabricated via laser powder bed fusion (LPBF) is studied in relation to the content of carbon and boron. Severe cracking occurred with the increase in boron content from 0.08 to 0.16 at. pct (0.015 to 0.029 wt pct), while compositions with 0.1 to 0.36 at. pct C (0.02 to 0.076 wt pct) and 0.08 at. pct B exhibited minimal cracking. Assessment of cracks in the high-boron composition shows a variation in crack density with printing parameters, and alignment of the cracks with the build direction. Scanning electron microscopy (SEM) of the crack surfaces shows evidence of a solidification cracking mode. Differential thermal analysis (DTA) reveals a decreased incipient melting temperature for the high-boron composition, and atom probe tomography (APT) is used to measure the enrichment at grain boundaries, revealing distinct boron segregation. Scheil-Gulliver solidification simulations for the different C and B levels are consistent with the incipient melting behavior observed with DTA. Evaluation of the solidification cracking susceptibility from the simulations allow for comparison of the CoNi alloy behavior to Ni-base superalloys studied for LPBF fabrication and displays how such metrics may aid in the design of new precipitation-strengthened superalloys for additive manufacturing (AM).

36 MATERIALS SCIENCE↗

Inverse design of plasma metamaterial devices with realistic elements

In an expansion of a previous study, we apply inverse design methods to produce two-dimensional plasma metamaterial devices with realistic plasma elements which incorporate quartz envelopes, collisionality (loss), non-uniform density profiles, and resistance to experimental error/perturbation. Finite difference frequency domain simulations are used along with forward-mode differentiation to design waveguides and demultiplexers operating under the transverse magnetic polarization. Additionally, optimal devices with realistic elements are compared to previous devices with idealized elements, and several parameter initialization schemes for the optimization algorithm are explored, yielding a robust procedure for producing such devices. Demultiplexing and waveguiding are demonstrated for microwave-regime devices composed of plasma elements with reasonable space-averaged plasma frequencies ~10GHz and a collision frequency ~1 GHz, allowing for future in-situ training and experimental realization of these designs.

42 ENGINEERING↗

EMI Mitigation of a Ćuk-Based Power-Electronic System Using Switching-Sequence-Based Control

Switching-sequence-based control (SBC) laws when designed based on topological switching behavior can have positive effects on slow- and fast-scale dynamics of a power-electronic system (PES). The slow-scale control can encompass fast PES state regulation and tracking, based on predefined objective, while fastscale control can address differential-mode (DM) and commonmode (CM) spectral-peak energy associated with PES switching operation. Such control laws may offer enhanced programmability to conventional PES design where bulky electromagnetic interference (EMI) filters have been traditionally used to reduce EMI of switching power converters to meet EMI regulatory standards. An EMI filter is always a less programmable solution since it is usually designed for the worst-case EMI mitigation and usually overkill for a PES operating under reduced load condition. The control scheme outlined in this article offers EMI mitigation across wide operating regions without compromising PES regulation. Moreover, it does so by use of switching sequences that guarantee the reachability of the PES dynamics using an advanced Lyapunov-function-based approach. SBC is a powerful tool to generate control actions for a PES based on multivariate PES state constraints. Hence, contemporary EMI regulatory standards are used as constraints in the SBC formulation to operate the PES under wide operating regime while autonomously mitigating the EMI levels. The work may be of paramount importance for operating the ultra-fast-transition recent wide-bandgap semiconductor devices like GaN–FET and SiC MOSFET under higher power with increasing switching frequencies, which is usually desirable for increased power density and reduced switching losses. Here, a hardware Cuk–PES operated ´ with GaN–FETs is fabricated and is used for case illustration. It is shown by experimental results how SBC mitigate DM and CM EMI noise of the PES while maintaining regulation even for the higher order nonminimum phase PES, while reducing sensor requirements using state observer derived from the switching model of the PES.

42 ENGINEERING↗

Instability Issue of Paralleled Dies in an SiC Power Module in Solid-State Circuit Breaker Applications

Paralleled dies in a power module could have instability issues during high current switching transients. Here, the instability is caused by the differential-mode oscillation among paralleled MOSFETs. Conventional analyses of paralleled MOSFETs’ stability are normally limited to a single operating point, which ignores the influences of the switching trajectory and nonlinear device parameters on stability. This article reveals that the switching trajectory can significantly influence parallel stability. The analysis is improved by solving eigenvalues of state-space modeling system matrices of all operating points that the switching trajectory goes through considering nonlinear device parameters. Higher voltage and current stresses result in greater real parts of complex eigenvalues, which explains why the paralleled MOSFETs are more unstable with higher voltage and current stresses. To improve stability in solid-state circuit breaker applications, we propose a method to manipulate the switching trajectory to avoid the unstable region where the conventional hard switching trajectory normally goes through. Experimental results show that the turn- off current capability can be increased from ~five times of rated current with the gate oscillation using the conventional turn- off trajectory to ~ten times of rated current without the gate oscillation using the optimal turn- off trajectory.

42 ENGINEERING↗

Precision measurement of the electron energy-loss function in tritium and deuterium gas for the KATRIN experiment

Abstract The KATRIN experiment is designed for a direct and model-independent determination of the effective electron anti-neutrino mass via a high-precision measurement of the tritium $$\upbeta $$ β -decay endpoint region with a sensitivity on $$m_\nu $$ m ν of 0.2 $$\hbox {eV}/\hbox {c}^2$$ eV / c 2 (90% CL). For this purpose, the $$\upbeta $$ β -electrons from a high-luminosity windowless gaseous tritium source traversing an electrostatic retarding spectrometer are counted to obtain an integral spectrum around the endpoint energy of 18.6 keV. A dominant systematic effect of the response of the experimental setup is the energy loss of $$\upbeta $$ β -electrons from elastic and inelastic scattering off tritium molecules within the source. We determined the energy-loss function in-situ with a pulsed angular-selective and monoenergetic photoelectron source at various tritium-source densities. The data was recorded in integral and differential modes; the latter was achieved by using a novel time-of-flight technique. We developed a semi-empirical parametrization for the energy-loss function for the scattering of 18.6-keV electrons from hydrogen isotopologs. This model was fit to measurement data with a 95% $$\hbox {T}_2$$ T 2 gas mixture at 30 K, as used in the first KATRIN neutrino-mass analyses, as well as a $$\hbox {D}_2$$ D 2 gas mixture of 96% purity used in KATRIN commissioning runs. The achieved precision on the energy-loss function has abated the corresponding uncertainty of $$\sigma (m_\nu ^2)< {{10}^{-2}}{\hbox {eV}^{2}}$$ σ ( m ν 2 ) < 10 - 2 eV 2 [1] in the KATRIN neutrino-mass measurement to a subdominant level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Stem CO2 Efflux measurements from Manaus, Brazil, 2017

This data was collected in order to test the effects of temperatures on the autotrophic respiratory process in the tropics. This data package contains raw real-time stem CO2 efflux files during the night and day from 3 different species near the K34 in Manaus during the day and the night. Raw and derived data files are included in the format of .csv and .xlsx and information about the canopy trees and temperatures recorded can be found in the field event log using microsoft excel. The data was collected from three canopy dominant trees in the central Amazon near the K34 tower. Real-time stem CO2 efflux was determined using a Li7000 gas analyzer configured in differential mode and connected to a dynamic stem chamber secured to the stem at breast height using straps to generate a reasonable seal. Ambient air was continuously pumped into the stem chamber and CO2 concentrations of the air entering and exiting the chamber were continuously monitored. CO2 efflux was calculated based on the CO2 concentration difference between ambient and stem chambers, the flow rate of air through the chamber, and the area of the enclosed stem. Additional auxiliary data, including sap flow and crown temperature were also measured. No data processing or QA/QC was done on the raw data packages.

54 ENVIRONMENTAL SCIENCES↗

High Altitude Electromagnetic Pulse Testing of Photovoltaic Modules

This report details the test setup, process, and results for radiated susceptibility testing of multicrystalline silicon photovoltaic (PV) modules as part of the EMP-Resilient Electric Grid Grand Challenge Laboratory Directed Research and Development (LDRD) project at Sandia National Laboratories. Testing was conducted over October 10-17, 2019, where 8 photovoltaic modules were exposed to E1 transient pulses with peak field levels up to 100 kV/m. Modules were terminated in a resistive load representing connected components. State of health testing conducted via I-V curve tracing of the photovoltaic modules showed no observable loss of device function due to large electric field transients. Differential mode currents were measured on the order of 10's of amps for up to a microsecond following the radiated field pulse. Common mode currents took the form of a damped sinusoid with a maximum peak of 10's to 100's of amps with a resonance near 60 MHz.

14 SOLAR ENERGY↗

Single-cell and spatiotemporal transcriptomic profiling of brain immune infiltration following Venezuelan equine encephalitis virus infection

Neurotropic alphaviruses such as Venezuelan equine encephalitis virus (VEEV) are critical human pathogens that continually expand to naïve populations and for which there are no licensed vaccines or therapeutics. VEEV is highly infectious via the aerosol route and is a recognized weaponizable biothreat that causes neurological disease in humans. The neuropathology of VEEV has been attributed to an inflammatory immune response in the brain yet the underlying mechanisms and specific immune cell populations involved are not fully elucidated. This study uses single-cell RNA sequencing to produce a comprehensive transcriptional profile of immune cells isolated from the brain over a time course of infection in a mouse model of VEEV. Analyses reveal differentially activated subpopulations of microglia, including a distinct type I interferon-expressing subpopulation. This is followed by the sequential infiltration of myeloid cells and cytotoxic lymphocytes, also comprising subpopulations with unique transcriptional signatures. We identify a subpopulation of myeloid cells that form a distinct localization pattern in the hippocampal region whereas lymphocytes are widely distributed, indicating differential modes of recruitment, including that to specific regions of the brain. Altogether, this study provides a high-resolution analysis of the immune response to VEEV in the brain and highlights potential avenues of investigation for therapeutics that target neuroinflammation in the brain.

59 BASIC BIOLOGICAL SCIENCES↗

A general approach to seismic inversion with automatic differentiation

Imaging Earth structure or seismic sources from seismic data involves minimizing a target misfit function, and is commonly solved through gradient-based optimization. The adjoint-state method has been developed to compute the gradient efficiently; however, its implementation can be time-consuming and difficult. We develop a general seismic inversion framework to calculate gradients using reverse-mode automatic differentiation. The central idea is that adjoint-state methods and reverse-mode automatic differentiation are mathematically equivalent. Here, the mapping between numerical PDE simulation and deep learning allows us to build a seismic inverse modeling library, ADSeismic, based on deep learning frameworks, which supports high performance reverse-mode automatic differentiation on CPUs and GPUs. We demonstrate the performance of ADSeismic on inverse problems related to velocity model estimation, rupture imaging, earthquake location, and source time function retrieval. ADSeismic has the potential to solve a wide variety of inverse modeling applications within a unified framework.

58 GEOSCIENCES↗

Differential Credibility of Climate Modes in CMIP6

This work quantitatively evaluates the fidelity with which the northern annular mode (NAM), southern annular mode (SAM), Pacific–North American pattern (PNA), El Niño–Southern Oscillation (ENSO), Pacific decadal oscillation (PDO), Atlantic multidecadal oscillation (AMO), and the first-order mode interactions are represented in Earth system model (ESM) output from the CMIP6 archive. Several skill metrics are used as part of a differential credibility assessment (DCA) of both spatial and temporal characteristics of the modes across ESMs, ESM families, and specific ESM realizations relative to ERA5. The spatial patterns and probability distributions are generally well represented but skill scores that measure the degree to which the frequencies of maximum variance are captured are consistently lower for most ESMs and climate modes. Substantial variability in skill scores manifests across realizations from individual ESMs for the PNA and oceanic modes. Further, the ESMs consistently overestimate the strength of the NAM–PNA first-order interaction and underestimate the NAM–AMO connection. Furthermore, these results suggest that the choice of ESM and ESM realizations will continue to play a critical role in determining climate projections at the global and regional scale at least in the near term.

54 ENVIRONMENTAL SCIENCES↗

Source-to-Source Automatic Differentiation of OpenMP Parallel Loops

This article presents our work toward correct and efficient automatic differentiation of OpenMP parallel worksharing loops in forward and reverse mode. Automatic differentiation is a method to obtain gradients of numerical programs, which are crucial in optimization, uncertainty quantification, and machine learning. The computational cost to compute gradients is a common bottleneck in practice. For applications that are parallelized for multicore CPUs or GPUs using OpenMP, one also wishes to compute the gradients in parallel. Here, we propose a framework to reason about the correctness of the generated derivative code, from which we justify our OpenMP extension to the differentiation model. We implement this model in the automatic differentiation tool Tapenade and present test cases that are differentiated following our extended differentiation procedure. Performance of the generated derivative programs in forward and reverse mode is better than sequential, although our reverse mode often scales worse than the input programs.

97 MATHEMATICS AND COMPUTING↗

Nonideal stability analysis of differentially rotating plasmas with global curvature effects

The linear stability of global nonaxisymmetric modes in differentially rotating, magnetized, nonideal plasma is critical to classifying turbulence and transport phenomena. We investigate the competition between the local magneto-rotational instability (MRI) and the magneto-curvature instability (MCI)—a distinct nonaxisymmetric low-frequency curvature-driven global branch that appears alongside MRI. Here, to accomplish this, we developed a nonideal global spectral method, which is validated against NIMROD code simulations. This spectral approach allows for the direct derivation of an extended effective potential formalism and a resistive Alfvénic resonance condition, providing a framework for direct analysis of energy contributions and confinement mechanisms. Our study reveals that the global, low-frequency MCI persists at low magnetic Reynolds numbers (Rm), whereas the localized, high-frequency MRI is stabilized by diffusive broadening of its structure around its Alfvénic resonances. Consequently, we identify the global MCI branch as the primary onset mechanism for nonaxisymmetric magnetohydrodynamic instability in systems with finite curvature, e.g., astrophysical rotators. We establish distinct parameter regimes for mode dominance: MCI prevails in geometrically moderate-thickness disks with intermediate curvature and radial gaps, while MRI dominates in thin, low-curvature disks with large radial gaps. Mode competition is also highly sensitive to the flow profile, particularly vorticity and its gradient, with nonuniform shear profiles exhibiting more robust instability due to flow curvature (i.e., the second derivative of the flow profile) and shear contributions. A key outcome is the development of spectral diagrams derived from the global spectral method. These diagrams comprehensively map dominant instabilities and their characteristics, offering a predictive tool for critical onset parameters (i.e., flow curvature, magnetic field, and Rm) and facilitating the interpretation of experimental and simulation results. Notably, these diagrams demonstrate that the global MCI is generally the sole unstable mode at the initial onset of nonaxisymmetric instability.

Haywood, Alexander [Princeton Univ., NJ (United St↗

Qualitative and quantitative enhancement of parameter estimation for model-based diagnostics using automatic differentiation with an application to inertial fusion

Parameter estimation using observables is a fundamental concept in the experimental sciences. Mathematical models that represent the physical processes can enable reconstructions of the experimental observables and greatly assist in parameter estimation by turning it into an optimization problem which can be solved by gradient-free or gradient-based methods. In this work, the recent rise in flexible frameworks for developing differentiable scientific computing programs is leveraged in order to dramatically accelerate data analysis of a common experimental diagnostic relevant to laser–plasma and inertial fusion experiments, Thomson scattering. A differentiable Thomson-scattering data analysis tool is developed that uses reverse-mode automatic differentiation (AD) to calculate gradients. By switching from finite differencing to reverse-mode AD, three distinct outcomes are achieved. First, gradient descent is accelerated dramatically to the extent that it enables near real-time usage in laser–plasma experiments. Second, qualitatively novel quantities which require $\mathcal{O}(10^3)$ parameters can now be included in the analysis of data which enables unprecedented measurements of small-scale laser–plasma phenomena. Third, uncertainty estimation approaches that leverage the value of the Hessian become accurate and efficient because reverse-mode AD can be used for calculating the Hessian.

97 MATHEMATICS AND COMPUTING↗

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗

Use of differential plasma rotation to prevent disruptive tearing mode onset from 3-wave coupling

Abstract Plasma differential rotation is found to be capable of preventing disruptive neoclassical tearing modes (NTMs) seeded by nonlinear three-wave coupling. As tearing modes degrade confinement and can lead to disruptions, stabilization strategies are crucial to the successful operation of future devices. In ITER-relevant scenarios on DIII-D, rotationally coupled m / n = 1/1 and 3/2 modes have been observed to drive 2/1 islands through three-wave coupling. The frequency of the driven 2/1 mode is set by matching conditions and the frequencies of the driving modes. When the driven mode frequency matches the local plasma rotation frequency, e.g. at low differential rotation, the driven 2/1 island can grow into a disruptive NTM. Using neutral beam torque as an actuator to scan the differential rotation, these experiments demonstrate that a sufficiently large frequency mismatch prevents destabilization of disruptive 2/1 NTMs by three-wave coupling. This work indicates that differential rotation can be used as an actuator to prevent NTMs seeded by three-wave coupling.

Richner, N. J. (ORCID:0000000155443915)↗

Hypercomplex Automatic Differentiation in the Eulerian Hydrocode PAGOSA

Enabling the computation of partial derivatives or sensitivities in production hydrocodes is beneficial for design, optimization, sensitivity analysis, and uncertainty quantification. Traditional finite difference approximations of these sensitivities are inefficient since convergence studies of the step size is required for each parameter of interest. For these reasons, HYPercomplex Automatic Differentiation (HYPAD) was implemented in the Eulerian hydrocode PAGOSA. HYPAD is analogous to forward-mode automatic differentiation except hypercomplex numbers (numbers with multiple imaginary parts) are used instead of dual numbers. Accurate partial derivatives can be computed of all state variables with respect to multiple input variables in a single run. The method was implemented using operator overloading to handle hypercomplex algebra. HYPAD was demonstrated and verified on Sod’s shock tube problem to compute derivatives of the state variables with respect to a material parameter, initial conditions, and geometry.

97 MATHEMATICS AND COMPUTING↗

CHEMNODE: CHEMICAL KINETICS SOLVER APPROACH BASED ON NEURAL ORDINARY DIFFERENTIAL EQUATIONS

This software contains an algorithm to model and predict the time evolution of chemical kinetics in reacting flow simulations using neural ordinary differential equations (NODEs). The approach works by using artificial neural networks to predict the chemical source terms. It achieves this by using forward-mode automatic differentiation and the Levenberg-Marquardt algorithm to adjust the neural network parameters, such that the discrepancies between the actual and predicted species profiles are minimized.

OWOYELE, OPEOLUWA↗