Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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 775 records · Page 43

Gate-Tunable Semiconductor Heterojunctions from 2D/3D van der Waals Interfaces

Van der Waals (vdW) semiconductors are attractive for highly scaled devices and heterogeneous integration since they can be isolated into self-passivated, two-dimensional (2D) layers that enable superior electrostatic control. These attributes have led to numerous demonstrations of field-effect devices ranging from transistors to triodes. By exploiting the controlled, substitutional doping schemes in covalently-bonded, three-dimensional (3D) semiconductors and the passivated surfaces of 2D semiconductors, one can construct devices that can exceed performance metrics of “all-2D” vdW heterojunctions. Here, we demonstrate, 2D/3D semiconductor heterojunctions using MoS2 as the prototypical 2D semiconductor laid upon Si and GaN as the 3D semiconductor layers. By tuning the Fermi levels in MoS2, we demonstrate devices that concurrently exhibit over seven orders of magnitude modulation in rectification ratios and conductance. Our results further suggest that the interface quality does not necessarily affect Fermi-level tuning at the junction opening up possibilities for novel 2D/3D heterojunction device architectures

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Unveiling sectoral coupling for resilient electrification of the transportation sector

Abstract Electrifying the transportation sector is crucial for reducing greenhouse gas emissions and offers numerous benefits including increased energy efficiency, lower total ownership costs, enhanced national energy security, and improved air quality. Despite the availability of necessary technologies, fully integrating the transportation and electricity sectors presents challenges in understanding all benefits and risks. Previous studies have not highlighted the role of coupling between these sectors. To better understand this coupling, this work reviews the structure of the current fossil-fuel-based transportation sector (including its dependence on the electricity sector) and case studies of its vulnerabilities to key risks. By adopting a systemic perspective, we uncover the indispensable interplay between the transportation and electricity sectors, shedding light on previously neglected dynamics. Leveraging the principles of grid architecture (GA), we introduce a hierarchical approach to assess vulnerabilities within the prevailing fuel-based transportation system and elucidate pathways for enhancement through electrification.

Mitra, Bhaskar↗

Surface chemistry models for GaAs epitaxial growth and hydride cracking using reacting flow simulations

Hydride vapor phase epitaxy (HVPE) is a promising technology that can aid in the cost reduction of III-V materials and devices manufacturing, particularly high-efficiency solar cells for space and terrestrial applications. However, recent demonstrations of ultrafast growth rates (~500 µm/h) via uncracked hydrides are not well described by present models for the growth. Therefore, it is necessary to understand the kinetics of the growth process and its coupling with transport phenomena, so as to enable fast and uniform epitaxial growth. In this work, we derive a kinetic model using experimental data and integrate it into a computational fluid dynamics simulation of an HVPE growth reactor. We also modify an existing hydride cracking model that we validate against numerical simulations and experimental data. Here, we show that the developed growth model and the improved cracking model are able to reproduce experimental growth measurements of GaAs in an existing HVPE system.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Metallic and porous Ti nanorod arrays for visible-IR light absorption and dendrite-free stable lithium-metal batteries

It is still a challenging task to prepare highly porous nanorod arrays of metals formed on substrates for optical and energy storage applications. Herein, we demonstrate the design and synthesis of black color, metallic and highly porous Ti nanorod arrays as novel current collectors for dendrite-free and highly stable Li-metal anodes. The high porosity of metallic nanorod arrays provides numerous heterogeneous nucleation sites and huge contact area and large space for the accommodation of Li metal. The conductive metallic Ti nanorod arrays enhance electrode integration. Effectively, it eliminates formation Li dendrites and demonstrates superior cycling stability over 300 cycles. Additionally, the unique porous structures of the nanorod arrays can decrease the amplitude of forced vibration in narrow space leading to light absorption. Interestingly, the metal is black instead of metallic color. The black metallic nanorod arrays can absorb more than 96% of both visible and infra-red lights. Furthermore, this black color metallic porous nanorod arrays may find additional applications in aerospace, energy, biomedical, defence, and chemical industries.

25 ENERGY STORAGE↗

Spectral function for 4 He using the Chebyshev expansion in coupled-cluster theory

Here, we compute spectral function for 4 He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows us to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark our predictions with other calculations in the literature and with electron-scattering data in the quasi-elastic peak. The spectral function formalism allows one to extend ab initio lepton-nucleus cross sections into the relativistic regime. This makes it a promising tool for modeling this process at higher-energy transfers. The results we present open the door for studies of heavier nuclei, important for the neutrino oscillation programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Two transitions in complex eigenvalue statistics: Hermiticity and integrability breaking

Open quantum systems have complex energy eigenvalues which are expected to follow non-Hermitian random matrix statistics, when chaotic, or two-dimensional (2d) Poisson statistics, when integrable. We investigate the spectral properties of a many-body quantum spin chain, i.e., the Hermitian Heisenberg model with imaginary disorder. Its rich complex eigenvalue statistics is found to separately break both Hermiticity and integrability at different scales of the disorder strength. With no disorder, the system is integrable and Hermitian, with spectral statistics corresponding to the 1d Poisson point process. At very small disorder, we find a transition from 1d Poisson statistics to an effective D -dimensional Poisson point process, showing Hermiticity breaking. At intermediate disorder, we find integrability breaking, as inferred from the statistics matching that of non-Hermitian complex symmetric random matrices in class AI † . For large disorder, as the spins align, we recover the expected integrability (now in the non-Hermitian setup), indicated by 2d Poisson statistics. These conclusions are based on fitting the spin-chain data of numerically generated nearest- and next-to-nearest-neighbor spacing distributions to an effective 2d Coulomb gas description at inverse temperature β . We confirm that such an effective description of random matrices also applies in classes AI † and AII † up to next-to-nearest-neighbor spacings. Published by the American Physical Society 2025

Akemann, Gernot (ORCID:0000000217104258)↗

Friction Stir Scribe Joining of CFRP to Aluminum (Abstract)

The objective of this project is to develop and demonstrate friction stir scribe (FSS) joining of carbon fiber reinforced polymer (CFRP) to aluminum on light-duty vehicle body-in-white (BIW) joints. The project will also overcome the major obstacles of implementing FSS technology in a fully 3-D production robotic work cell and demonstrate the required CFRP-aluminum joint strengths in industrially relevant light-duty, BIW components. Specifically the project will develop the critical process technology, models and tools necessary to advance the FSS method through experimentation, validation at the laboratory scale, and integration into a production-like robotic environment. The ultimate goal is to demonstrate the technology on industrially relevant components based on process development and numerical simulation performed at scale.

42 ENGINEERING↗

Investigation of powertrain system decarbonization using electrically assisted turbocharging and hybridization in off-road vehicles

In recent years, the causes, effects, and existential threat of a global anthropogenic climate shift have drawn significant attention and stimulated mitigation efforts from all genres of scientific, political, and industrial bodies. The greenhouse effect and the detrimental environmental impact of excess greenhouse gas (GHG) emissions, like carbon dioxide, is well studied and targeted as the primary culprit for aversive action. However, some fields exist where reducing GHG emissions is met by considerable challenges. One such field is the production and operations of off-road vehicles. Applications of these vehicles are highly diverse and are often characterized by rugged, especially transient, and power intensive duty cycles that make one-fits-all vehicle configuration solutions impractical. This research surveys existing literature to identify the modern technologies packages and challenges that face development and configuration of powertrain systems which currently navigate the regulated emissions and unique duty cycle requirements of this market space. A novel powertrain concept is proposed and evaluated with respect to the pinnacle objectives of load performance improvement, criteria pollutant reduction, and reducing the life cycle GHG intensity of its operation. A sure-fire path to vehicular GHG reduction is through improving the fuel efficiency of internal combustion engines (ICEs), an entrenched component of the off-road vehicle sector. Unfortunately, this is more easily said than done. An approach that has proven successful in this endeavor is engine downsizing and turbocharging, where a larger engine is replaced with a smaller one with added air system boosting via turbocharger to reduce frictional and pumping losses while also enabling access to additional fuel energy. However, practical realization of these potential benefits is often impeded by the transient response capability of the smaller engine across the operating space. For this reason, downsized ICE powertrains have turned to electrified forced induction systems (EFISs) for a decoupling of exhaust energy and engine speed from boost capability. Platformed on 48V hybrid technology, these systems introduce the need for more sophisticated controls around the engine gas-exchange process for the management of boost performance and exhaust gas emissions. On this account, simulation studies utilizing a GT-POWER model of a turbocharged 4.5 L engine outfit with an EFIS using an electrically driven compressor (eBooster®) are conducted to provide insight into the performance of this air handling architecture on an off-road engine. The results show that improvement in transient torque response time in sync with reduction in engine-out soot and NOx emissions are possible with an engine recalibration that leverages the transient air-fuel ratio authority of the eBooster®. Benefits are further demonstrated when duty cycle simulations of the 48V mild-hybrid engine concept are exercised, showing an acute decrease in cumulative fuel usage and soot production. Powertrain hybridization is another technological pathway achieving pronounced GHG reduction successes in modern on-road vehicle applications through integration of Li-ion battery technology. In the off-road vehicle segment, a review of available literature concludes that hybridized architectures are present but generally lack the depth of technologies that have both high specific energy and power capabilities, and thus are limited in their inclusion of Li-ion batteries for ICE assistance and enhanced energy storage capability. Therefore, building on the mild-hybrid engine results, the downsized and eBoosted engine concept was integrated into a larger high-voltage battery-hybrid series-electric powertrain system. Hybrid powertrain parameter sensitivity studies were carried out in a numerical charge-sustaining framework, providing novel insights into power flows between the battery and the engine and how their respective capabilities and operation contribute to GHG and criteria pollutant emissions of diverse duty cycles. Application of supervisory power management introduced robustness into the power sourcing and battery SOC control process and showed that optimum specification of battery properties can yield synergies between GHG emissions and battery energy capacity. Furthermore, examination of recent literature on Li-ion batteries has shown that pack manufacturing is a highly energy intensive process, thereby producing considerable quantities of GHGs that scale with energy storage capacity. Also scaling with a battery’s energy storage capacity is its investment cost. In consideration of these factors, an inclusive technoeconomic and GHG life cycle analysis is conducted. This analysis systematically compares the carbon footprint and total cost of ownership associated with the proposed hybrid powertrain concept to reference and alternative powertrain configurations, facilitating a thorough evaluation of the decarbonization effectiveness and economic viability.

99 GENERAL AND MISCELLANEOUS↗

GPU acceleration of all-electron electronic structure theory using localized numeric atom-centered basis functions

We present an implementation of all-electron density-functional theory for massively parallel GPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3–4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi2Se3 bilayer show that the present GPU strategy scales for large-scale distributed-parallel simulations.

42 ENGINEERING↗

Study of recirculating liquid fuel in a 1D critical stationary system

Several studies have been conducted to investigate the physics of liquid fuel reactors, showing also applications with molten salts. A liquid nuclear fuel implies changes in the neutron balance equation to take into account the precursors' displacement and the emission of delayed neutrons in a different position than at the original fission. This requires to upgrade the computer codes normally used to calculate nuclear reactors using only solid fuel. In this work, we revisit a simple problem with liquid fuel, which is proposed for the verification of the numerical solutions obtained by advanced computer codes. This problem studies criticality with constant coefficients, thus neglecting thermal feedback. We elaborate on the analytical solution of the problem, deriving also a generalized eigenvalue problem by finite-volume integration over the cells of a discretized mesh to study the evolution of the dominance ratio with fuel velocity. Finally, we investigate the influence of the fuel velocity on the reactivity of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and Integration of Predictive Models for Manufacturing and Structural Performance of Carbon Fiber Composites in Automotive Applications (Final Report)

The goal of this project is to develop integrated, state-of-the-art, computational modeling tools based on an integrated computational materials engineering (ICME) methodology that are critically needed to enable structural carbon fiber (CF) applications in automobiles. These tools are designed to predict the manufacturing and structural performance of CF composites, including stochastic effects. During the first phase of the project, the manufacturing and structural performance tools, including a stochastic driver, were developed, calibrated, and validated against coupon and component level test results. For this project, from the manufacturing side, the development efforts were focused on high-pressure resin transfer molding (HP-RTM), a potentially game-changing manufacturing technology that is capable of reaching the 3-5-minute cycle time that is necessary for volume manufacturing of composites. On the structural performance side, the crashworthiness of CF structures was studied. Accordingly, computational tools were developed and validated. The resulting differences between the numerical predictions and the experimental results in the first phase of the project were required to be less than 15%. During the second phase of the project, the manufacturing and performance tools were integrated by mapping the manufacturing outcomes (e.g., fiber angles, residual stresses, degree of cure, and defects) into the structural models. Also, using the tools developed in this project, an automotive assembly currently manufactured in steel (2016 GM Malibu) was redesigned to be manufactured using advanced CF composites with the objective of comparing the performance of the CF assembly with the developed ICME model predictions. The weight, performance and cost of the CF assembly were compared with the baseline steel assembly as well. The CF underbody design was shown to be 30% lighter than the corresponding steel design. The CF underbody assembly withstood the side pole impact with less than half the intrusion for the steel assembly. Excellent correlations within 10% for both peak impact load and intrusion were observed between the numerical predictions and experimental results. Based on the cost models developed in the project, the cost increase per Kg saved was determined to be $22. The project team has prepared a total of 28 publications in various journals and conferences and broadly communicated the project results for the benefit of the composite industry. The project team has also prepared a total of 10 patent applications that were submitted to the US Patent and Trademark office.

36 MATERIALS SCIENCE↗

Bound-preserving finite element approximations of the Keller–Segel equations

We report this paper aims to develop numerical approximations of the Keller–Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a non-negative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear and can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and nonblowup phenomena. In the blowup setting, we identify a locking phenomenon that relates the L ∞ (Ω)-norm to the L 1 (Ω)-norm limiting the growth of the singularity when supported on a macroelement.

97 MATHEMATICS AND COMPUTING↗

Decoupling Carrier Dynamics and Energy Transport in Ultrafast Near-Field Nanoscopy

Ultrafast near-field optical nanoscopy has emerged as a powerful platform to characterize low-dimensional materials. While analytical and numerical models have been established to account for photoexcited carrier dynamics, quantitative evaluation of the associated pulsed laser heating remains elusive. Here, we decouple the photocarrier density and temperature increase in near-field nanoscopy by integrating the two-temperature model (TTM) with finite-difference time-domain (FDTD) simulations. These results reveal that the electron–phonon coupling in a silicon film after femtosecond laser excitation is most pronounced within approximately 3 ps–substantially shorter than the photocarrier decay time scale at tens of picoseconds. Moreover, the coupled TTM-FDTD method indicates that ultrafast laser heating can cause up to a 14% variation in the near-field signal at a 220 μJ/cm 2 pump pulse fluence. Our numerical results are further validated by transient experiments, highlighting the potential of this method for investigations of carrier and thermal phenomena in emerging nanomaterials and nanodevices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗

Conservative DG method for the micro-macro decomposition of the Vlasov–Poisson–Lenard–Bernstein model

The micro-macro (mM) decomposition approach is considered for the numerical solution of the Vlasov– Poisson–Lenard–Bernstein (VPLB) system, which is relevant for plasma physics applications. In the mM approach, the kinetic distribution function is decomposed as $f$ = $\mathscr{E}$[$ρ_f$] + g , where $\mathscr{E}$ is a local equilibrium distribution, depending on the macroscopic moments $ρ_f$ = $∫_{\mathbb{R}}$ $efdv$ = $\langle$$ef$$\rangle$$_{\mathbb{R}}$, where e = (1, $v$, $\frac{1}{2}$ $v^2$) T , and $g$, the microscopic distribution, is defined such that $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0. We aim to design numerical methods for the mM decomposition of the VPLB system, which consists of coupled equations for $ρ_f$ and $g$. To this end, we use the discontinuous Galerkin (DG) method for phase-space discretization, and implicit-explicit (IMEX) time integration, where the phase-space advection terms are integrated explicitly and the collision operator is integrated implicitly. We give special consideration to ensure that the resulting mM method maintains the $\langle$$eg$$\rangle$$_{\mathbb{R}}$ = 0 constraint, which may be necessary for obtaining (i) satisfactory results in the collision dominated regime with coarse velocity resolution, and (ii) unambiguous conservation properties. The constraint-preserving property is achieved through a consistent discretization of the equations governing the micro and macro components. Here, we present numerical results that demonstrate the performance of the mM method. The mM method is also compared against a corresponding DG-IMEX method solving directly for $f$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗

Access Patterns and Performance Behaviors of Multi-layer Supercomputer I/O Subsystems under Production Load

Scientific computing workloads at HPC facilities have been shifting from traditional numerical simulations to AI/ML applications for training and inference while processing and producing ever-increasing amounts of scientific data. To address the growing need for increased storage capacity, lower access latency, and higher bandwidth, emerging technologies such as non-volatile memory are integrated into supercomputer I/O subsystems. With these emerging trends, we need a better understanding of the multilayer supercomputer I/O systems and ways to use these subsystems efficiently. In this work, we study the I/O access patterns and performance characteristics of two representative supercomputer I/O subsystems. Through an extensive analysis of year-long I/O logs on each system, we report new observations in I/O reads and writes, unbalanced use of storage system layers, and new trends in user behaviors at the HPC I/O middleware stack.

Bez, JL↗

Marine ice sheet experiments with the Community Ice Sheet Model

Ice sheet models differ in their numerical treatment of dynamical processes. Simulations of marine-based ice are sensitive to the choice of Stokes flow approximation and basal friction law and to the treatment of stresses and melt rates near the grounding line. We study the effects of these numerical choices on marine ice sheet dynamics in the Community Ice Sheet Model (CISM). In the framework of the Marine Ice Sheet Model Intercomparison Project 3d (MISMIP3d), we show that a depth-integrated, higher-order solver gives results similar to a 3D (Blatter–Pattyn) solver. We confirm that using a grounding line parameterization to approximate stresses in the grounding zone leads to accurate representation of ice sheet flow with a resolution of ~2 km, as opposed to ~0.5 km without the parameterization. In the MISMIP+ experimental framework, we compare different treatments of sub-shelf melting near the grounding line. In contrast to recent studies arguing that melting should not be applied in partly grounded cells, it is usually beneficial in CISM simulations to apply some melting in these cells. This suggests that the optimal treatment of melting near the grounding line can depend on ice sheet geometry, forcing, or model numerics. In both experimental frameworks, ice flow is sensitive to the choice of basal friction law. To study this sensitivity, we evaluate friction laws that vary the connectivity between the basal hydrological system and the ocean near the grounding line. CISM yields accurate results in steady-state and perturbation experiments at a resolution of ~2 km (arguably 4 km) when the connectivity is low or moderate and ~1 km (arguably 2 km) when the connectivity is strong.

54 ENVIRONMENTAL SCIENCES↗