Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “limiting techniques”

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 19 records

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic Features of Cu-Ceria Interface under CO 2 Hydrogenation to Methanol

It is generally accepted that metal–support interaction is very important for the hydrogenation of CO 2 to methanol, but little has been revealed about the feature of interfacial active sites under real reaction conditions since there are only limited techniques that can be applied under high-pressure conditions. Here, in this work, by combining multiple in situ and operando techniques on a model Cu/ceria catalyst, we have tracked Cu and ceria sites for methanol formation. Under the reaction condition, it is found that upon reaching the reaction temperature, oxidized Cu species in the as-synthesized catalyst immediately change into metallic Cu species. Following this, it is the gradual formation of methanol, the changing rate of which coincides with the formation of a unique Ce 3+ species. The combined experimental results and density functional theory (DFT) calculations have determined that the formed Ce 3+ sites driven by the reaction conditions are bound to hydrides, adsorbed carbonate species, and interfacial active Cu sites. The Cu-ceria interaction in this complex moiety is weak and can be easily disturbed with reaction environment variations, leading to dynamic changes at the interface upon the hydrogenation of active carbonate intermediates, which are precursors for the formation of methanol. The formation of this unique Cu–Ce 3+ interface and its dynamicity lead to an increase of methanol selectivity from less than 20% to 60%. These results suggest that reactant-derived species (H – and carbonate in this work) can be essential components of the active center with the functions of manipulating the metal−oxide interaction and directing reaction pathways.

CO2 hydrogenation to methanol↗

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

Aerosol jet printed 3 omega sensors for thermal conductivity measurement

The 3 omega (3ω) method is a trusted technique for measuring thermal conductivity—a fundamental material property of critical importance in a broad range of applications. However, traditional 3ω sensor processing requires some form of physical vapor deposition, such as metal evaporation or sputtering. These 3ω sensor deposition techniques limit the materials and sample sizes applicable to the 3ω method. This work demonstrates an aerosol jet printing method to directly print silver 3ω sensors that yield accurate temperature-dependent measurement up to 300 °C on materials with thermal conductivity ranging from 1 to 150 W/m K. The interrelationship between printed sensor geometry, sensor thermal stability, and applicability to the 3ω method is examined. Thermal conductivity measurement with 3ω sensors conventionally sintered at 300 °C agrees to independent laser flash measurement within 4% from room temperature to 150 °C. An unconventional rapid high-temperature sintering method is shown to produce sensors that agree within 3% of the laser flash measurements from room temperature to 300 °C. The rapid sintering profiles also reduced the sensor–substrate thermal boundary resistance of the printed sensors by as much as 88%. The direct printing of 3ω sensors creates opportunities for measurement of thermal transport properties in applications previously inapplicable to the 3ω method.

47 OTHER INSTRUMENTATION↗

Single-molecule 3D imaging of HIV cellular entry by liquid-phase electron tomography

Enveloped viruses, including human immunodeficiency virus (HIV) and SARS-CoV-2, target cells through membrane fusion process. The detailed understanding of the process is sought after for vaccine development but remains elusive due to current technique limitations for direct three-dimensional (3D) imaging of an individual virus during its viral entry. Recently, we developed a simple specimen preparation method for real-time imaging of metal dynamic liquid-vaper interface at nanometer resolution by transmission electron microscopy (TEM). Here, we extended this method to study biology sample through snapshot 3D structure of a single HIV (pseudo-typed with the envelope glycoprotein of vesicular stomatitis virus, VSV-G) at its intermediate stage of viral entry to HeLa cells in a liquid-phase environment. By individual-particle electron tomography (IPET), we found the viral surface release excess lipids with unbound viral spike proteins forming ~50-nm nanoparticles instead of merging cell membrane. Moreover, the spherical-shape shell formed by matrix proteins underneath the viral envelope does not disassemble into a cone shape right after fusion. Further, the snapshot 3D imaging of a single virus provides us a direct structure-based understanding of the viral entry mechanism, which can be used to examine other viruses to support the development of vaccines combatting the current ongoing pandemic.

Kong, Lingli↗

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Membrane Profiling by Free Flow Electrophoresis and SWATH-MS to Characterize Subcellular Compartment Proteomes in Mesembryanthemum crystallinum

The study of subcellular membrane structure and function facilitates investigations into how biological processes are divided within the cell. However, work in this area has been hampered by the limited techniques available to fractionate the different membranes. Free Flow Electrophoresis (FFE) allows for the fractionation of membranes based on their different surface charges, a property made up primarily of their varied lipid and protein compositions. In this study, high-resolution plant membrane fractionation by FFE, combined with mass spectrometry-based proteomics, allowed the simultaneous profiling of multiple cellular membranes from the leaf tissue of the plant Mesembryanthemum crystallinum. Comparisons of the fractionated membranes’ protein profile to that of known markers for specific cellular compartments sheds light on the functions of proteins, as well as provides new evidence for multiple subcellular localization of several proteins, including those involved in lipid metabolism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-of-flight detection of terahertz phonon-polariton

A polariton is a fundamental quasiparticle that arises from strong light-matter interaction and as such has attracted wide scientific and practical interest. When light is strongly coupled to the crystal lattice, it gives rise to phonon-polaritons (PPs), which have been proven useful in the dynamical manipulation of quantum materials and the advancement of terahertz technologies. Yet, current detection and characterization methods of polaritons are still limited. Traditional techniques such as Raman or transient grating either rely on fine-tuning of external parameters or complex phase extraction techniques. To overcome these inherent limitations, we propose and demonstrate a technique based on a time-of-flight measurement of PPs. We resonantly launch broadband PPs with intense terahertz fields and measure the time-of-flight of each spectral component with time-resolved second harmonic generation. The time-of-flight information, combined with the PP attenuation, enables us to resolve the real and imaginary parts of the PP dispersion relation. We demonstrate this technique in the van der Waals magnets NiI 2 and MnPS 3 and reveal a hidden magnon-phonon interaction. We believe that this approach will unlock new opportunities for studying polaritons across diverse material systems and enhance our understanding of strong light-matter interaction.

condensed-matter physics↗

Preserving Contextual Awareness during Selection of Moving Targets in Animated Stream Visualizations

In many types of dynamic interactive visualizations, it is often desired to interact with moving objects. Stopping moving objects can make selection easier, but pausing animated content can disrupt perception and understanding of the visualization. To address such problems, we explore selection techniques that only pause a subset of all moving targets in the visualization. We present various designs for controlling pause regions based on cursor trajectory or cursor position. We then report a dual-task experiment that evaluates how different techniques affect both target selection performance and contextual awareness of the visualization. Our findings indicate that all pause techniques significantly improved selection performance as compared to the baseline method without pause, but the results also show that pausing the entire visualization can interfere with contextual awareness. However, the problem with reduced contextual awareness was not observed with our new techniques that only pause a limited region of the visualization. Thus, our research provides evidence that region-limited pause techniques can retain the advantages of selection in dynamic visualizations without imposing a negative effect on contextual awareness.

Ragan, Eric↗

Answering old questions with new techniques: Understanding performance-limiting factors in transition metal dichalcogenide photoelectrochemical solar cells

In the late 1970s and early 1980s, several research groups reported high-efficiency and stable photoelectrochemical solar cells based on transition metal dichalcogenides (TMDS; e.g., MoS 2 , WSe 2 , and MoSe 2 ) immersed in iodide/tri-iodide electrolytes. A consensus emerged that smooth crystals were necessary for high efficiency based on significant evidence that rough crystals with exposed edge sites produced lower photocurrents and fill factors. However, anecdotal observations in the literature hinted at significant performance variation among apparently smooth crystals with the possibility of highly active “hot” edge sites. Furthermore, this mini-review article is a case study on how spatially resolved photoelectrochemical techniques developed in the 2020s are answering old questions regarding the origin of performance variation in high efficiency n-type TMD|I – ,I 3 – |Pt photoelectrochemical solar cells.

14 SOLAR ENERGY↗

Non-target RNA depletion strategy to improve sensitivity of next-generation sequencing for the detection of RNA viruses in poultry

PCR-based assays have become the benchmark for detecting pathogens of poultry and other livestock; however, these techniques are limited in their ability to detect multiple infecting agents, provide limited genetic information on the pathogen, and, for RNA viruses, must be reviewed frequently to assure high sensitivity and specificity. In contrast, untargeted, high-throughput sequencing can rapidly detect all infecting agents in a sample while providing genomic sequence information to allow more in-depth characterization of viruses. Although next-generation sequencing (NGS) offers many advantages, one of its primary limitations is low sensitivity to pathogens given the abundance of host and other non-target sequences in sequencing libraries. We explored methods for improving the sensitivity of NGS to detect respiratory and enteric viruses in poultry from RNA extracts of swab samples. We employed commercial and custom-designed negative enrichment strategies to selectively deplete the most abundant rRNA reads from the host and non-target bacteria; host RNA was diminished from up to 40% of total reads to as low as 3%, and the total number of reads assigned to abundant bacterial classes were reduced greatly. Our treatment resulted in up to a 700-fold increase in the number of viral reads, detection of a greater number of viral agents, and higher average genome coverage for pathogens. Depletion assays added only 2 h to the NGS library preparation workflow. Custom depletion probe design offered significant cost savings (US$7–12 per sample) compared to commercial kits (US$30–50 per sample).

Parris, D. Joshua↗

Conservation Voltage Reduction (CVR) via Two-Timescale Control in Unbalanced Power Distribution Systems

Voltage control devices are employed in power distribution systems to reduce the power consumption by operating the system closer to the lower acceptable voltage limits; this technique is called conservation voltage reduction (CVR). The different modes of operation for system’s legacy devices (with binary control) and new devices (e.g. smart inverters with continuous control) coupled with variable photovoltaic (PV) generation results in voltage fluctuations which makes it challenging to achieve CVR objective. This paper presents a two-timescale control of feeder’s voltage control devices to achieve CVR that includes (1)a centralized controller operating in a slower time-scale to coordinate voltage control devices across the feeder and (2) local controllers operating in a faster timescale to mitigate voltage fluctuations due to PV variability. The centralized controller utilizes a three-phase optimal power flow model to obtain the decision variables for both legacy devices and smart inverters. The local controllers operate smart inverters to minimize voltage fluctuations and restore nodal voltages to their reference values by adjusting the reactive power support. The proposed approach is validated using the IEEE-123 bus (medium-size) and R3-12.47-2 (large-size) feeders. It is demonstrated that the proposed approach is effective in achieving the CVR objective for unbalanced distribution systems.

optimal control, distributed generation, fluctuati↗

In situ and operando force-based atomic force microscopy for probing local functionality in energy storage materials

Electrochemical energy storage is the key enabling component of electric vehicles and solar-/wind-based energy technologies. The enhancement of energy stored requires the detailed understanding of charge storage mechanisms and local electrochemical and electromechanical phenomena over a variety of length scales from atoms to full cells. Classical electrochemical techniques, such as voltammetry, represent the macroscopic electrochemical properties, and consequently do not allow to extract important information about local electrochemical reactions, ions adsorption, intercalations, and transport. Understanding, controlling, and tuning the local electrochemical functionalities in functional energy materials require in situ/operando techniques which limit the use of structural and functional characterization techniques that provide local information. Here, force-based atomic force microscopies (AFMs) have provided novel insights into locally probing electrochemical mechanisms on tens of nanometer and even molecular length scales and provide a viable pathway to probe electrochemical processes in situ/operando. In this review, we highlight the contributions in the development and application of force-based AFM methods to elucidate the local charge storage mechanism in a variety of energy-related materials. We will focus in particular on methods or AFM modalities in a liquid electrolyte environment.

25 ENERGY STORAGE↗

Defect Characterization of Monocrystalline Silicon Solar Cells with Polysilicon Passivated Contact Using Electrically-Detected Magnetic Resonance (EDMR) Spectroscopy

As the c-Si based solar cell efficiencies are approaching over 26%, it is becoming critical to characterize the low concentrations of the defects – as low as 10^10-10^11 cm-3 (for e.g., iron contamination in high-lifetime Ga-doped wafers3 and n-type wafers), and further reduce them. Also, atomistic level understanding of the mechanisms of the low concentration process-induced-defects and reliability limiting defects (such as light and elevated temperature induced degradation, surface passivation degradation) is needed to design the mitigation strategies. The conventional characterization techniques are limited due to their detection limitations. Some of the techniques based on lifetime spectroscopies can still be used for low concentration characterization however, they are based on estimations and theoretical models and hence, indirect and cannot fully reveal information about the microscopic mechanism of the defects. Thus, we present the application of an ultrasensitive magnetic resonance-based technique for the direct spectroscopic detection of the defects in Si PV - electrically detected magnetic resonance (EDMR). In this work, we aim to focus on establishing a process flow for fabrication of minicells with (miniature replica of the larger-area cells) and setting up the routine for EDMR measurements on them with the EDMR instrumentation capability at NREL. For the EDMR measurements, sample size is limited by the dimensions of sample holder tube (width less than 3.2 mm, active area - 20 mm). Thus, we have designed c-Si based minicells with polysilicon (poly-Si) passivated contacts same as the larger-area cells that we fabricate in our group at NREL. We also modified our minicell process flow for fabricating the textured minicells for preserving the texture during processing and taking care of the laser-ablation edge damage which can significantly affect the performance of such small devices. We have achieved comparable performance on these newly fabricated minicells as that of our 4 cm2 devices with same structure (comparable VOC, JSC, FF). We also conducted EDMR measurements on the minicells and observed a distinct EDMR signal at g-value ~2.005 at temperatures 30K and above, as shown in Fig. 2. We associate this signal to the presence of silicon dangling bonds based on the g-value. We also observed an EDMR signal at g-value ~1.998 at temperature ~5K. The origin of this signal is still being investigated. Thus, we show the proof of concept of minicells and EDMR measurements with which we now aim to study some of the unknown defects in silicon solar cell devices.

EDMR↗