Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗
DOE OSTI2024

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2020

Validation of neutron emission and neutron energy spectrum calculations on a Mega Ampere Spherical Tokamak with directional relativistic spectrum simulator

The recently developed directional relativistic spectrum simulator (DRESS) code has been validated for the first time against numerical calculations and experimental measurements performed on the Mega Ampere Spherical Tokamak . In this validation, the neutron emissivities and rates computed by DRESS are benchmarked against TRANSP/NUBEAM predictions while the neutron energy spectra provided by DRESS taking as input TRANSP/NUBEAM and ASCOT/BBNBI in Gyro-Orbit mode fast ion distributions are validated against proton pulse height spectra measured by the neutron flux monitor. Finally, excellent agreement was found between DRESS and TRANSP/NUBEAM predictions of local and total neutron emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Bus Clustering for Distribution Grid Topology Identification

Recovering the distribution grid topology is essential to perform several distribution system operator functions. Many algorithms that address the topology recovery problem have already been proposed in the literature. Most are based on a priori information regarding which buses are fed by which substation; however, this information might not be available because frequent grid reconfigurations change the distribution grid portion connected to each substation. This paper addresses the problem of assigning every substation the set of buses that it is feeding, given field data. First, the aforementioned task is cast as a nonconvex optimization problem. Second, a relaxed version of the optimization problem is solved via the alternating direction method of multipliers. Finally, the performance of our approach is validated through numerical simulations of realistic scenarios using a standard IEEE benchmark feeder.

bus clustering↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers part I: Computational models and workflow

Computational simulations have become central to the seismic analysis and design of major infrastructure over the past several decades. Most major structures are now “proof tested” virtually through representative simulations of earthquake-induced response. More recently, with the advancement of high-performance computing (HPC) platforms and the associated massively parallel computational ecosystems, simulation is beginning to play a role in increased understanding and prediction of ground motions for earthquake hazard assessments. However, the computational requirements for regional-scale geophysics-based ground motion simulations are extreme, which has restricted the frequency resolution of direct simulations and limited the ability to perform the large number of simulations required to numerically explore the problem parametric space. In this article, recent developments toward an integrated, multidisciplinary earth science-engineering computational framework for the regional-scale simulation of both ground motions and resulting structural response are described with a particular emphasis on advancing simulations to frequencies relevant to engineered systems. This multidisciplinary computational development is being carried out as part of the US Department of Energy (DOE) Exascale Computing Project with the goal of achieving a computational framework poised to exploit emerging DOE exaflop computer platforms scheduled for the 2022–2023 timeframe.

58 GEOSCIENCES↗
DOE OSTI2022

Numerical and Experimental Investigation of Nanostructure-Based Asymmetric Light Transmission Interfaces for Solar Concentrator Applications

Research in asymmetric light transmission interfaces has been recently gaining traction. While traditionally considered for optical circuitry applications, there is a new interest to use these interfaces in luminescent solar concentrators. Previous studies have shown that applying them to the top surface of a concentrator could mitigate surface losses. This paper presents experimental results for proof-of-concept asymmetric light transmission interfaces that may have potential applications in luminescent solar concentrators. The interfaces and the underneath substrate were created in a single step from polydimethylsiloxane using silicon molds fabricated on <100> wafers via anisotropic wet etching. The resulting structures were pyramidal in shape. Large surface areas of nanostructures repeating at 800 nm, 900 nm, and 1000 nm were tested for backward and forward transmission using a spectrometer. Results showed a 21%, 10%, and 0% average transmissivity difference between the forward and backward directions for each periodicity, respectively. The trends seen experimentally were confirmed numerically via COMSOL simulations.

Oliveto, Vincent↗
DOE OSTI2023

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Prediction of dynamic behavior of a system undergoing an exchange of components using a transmission simulator method

In systems with modular design, it is of interest to perform component exchanges to upgrade components and explore new system configurations. If the dynamics of the components are known, then dynamic substructuring can be used to obtain a prediction of system performance for the new configuration. The transmission simulator method is a substructuring technique that utilizes an additional fixture to apply an elastic boundary condition to a test article such that an appropriate dynamic response can be obtained for a given substructuring application, and the technique has been used widely for performing coupling and decoupling operations. Here, in this work, a modified transmission simulator method is proposed that allows a direct exchange of components to be performed. In a numerical demonstration, the proposed method is shown to accurately predict the dynamics of an assembly that has undergone a component exchange in which the exchanged components have different material properties and different geometric properties. This approach can be applied to dynamic qualification frameworks in which structures undergo modifications or upgrades rather than complete redesigns.

Component exchange↗
DOE OSTI2021

Multi-Area Model-Free State Estimation via Distributed Tensor Decomposition

This paper proposes a model-free method for distribution system state estimation based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A third-order tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. To demonstrate the efficacy of the developed algorithm, numerical simulations are carried out using an IEEE test system.

alternating direction method of multipliers↗
DOE OSTI2025

Experimental and computational evaluation of alpha particle production from laser-driven proton–boron nuclear reaction in hole-boring scheme

The majority of studies on laser-driven proton–boron nuclear reaction is based on the measurement of α-particles with solid-state nuclear tracks detector (Cr39). However, Cr39's interpretation is difficult due to the presence of several other accelerated particles which can bias the analysis. Furthermore, in some laser irradiation geometries, cross-checking measurements are almost impossible. In this case, numerical simulations can play a very important role in supporting the experimental analysis. In our work, we exploited different laser irradiation schemes (pitcher–catcher and direct irradiation) during the same experimental campaign, and we performed numerical analysis, allowing to obtain conclusive results on laser-driven proton–boron reactions. A direct comparison of the two laser irradiation schemes, using the same laser parameters is presented.

Huault, M. (ORCID:0000000186999166)↗
DOE OSTI2020

Ptychographic phase retrieval by proximal algorithms

We derive a set of ptychography phase-retrieval iterative engines based on proximal algorithms originally developed in convex optimization theory, and discuss their connections with existing ones. The use of proximal operator creates a simple frame work that allows us to incorporate the effect of noise from a maximum-likelihood (ML) principle. We focus on three particular algorithms, namely proximal minimization, alternating direction method of multiplier and accelerated proximal gradient (APG). We benchmark their performance with numerical simulations, and discuss their optimal conditions for convergence and accuracy. An experimental dataset is used to demonstrate their effectiveness as well, in which case an array of cubic Au nanoparticles with a size of 50 nm is imaged. We show that with the presence of Poisson noise, a dataset with photon counts up to 10 4 at one detector pixel already requires ML-based methods to achieve a stable convergence. Among the three algorithms derived in this work, APG method is reported first time for its application in ptychographic reconstruction and shows superior performance in terms of both accuracy and convergence rate with a noisy dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Turbulence-level dependence of cosmic ray parallel diffusion

ABSTRACT Understanding the transport of energetic cosmic rays belongs to the most challenging topics in astrophysics. Diffusion due to scattering by electromagnetic fluctuations is a key process in cosmic ray transport. The transition from a ballistic to a diffusive-propagation regime is presented in direct numerical calculations of diffusion coefficients for homogeneous magnetic field lines subject to turbulent perturbations. Simulation results are compared with theoretical derivations of the parallel diffusion coefficient’s dependences on the energy and the fluctuation amplitudes in the limit of weak turbulence. The present study shows that the widely used extrapolation of the energy scaling for the parallel diffusion coefficient to high turbulence levels predicted by quasi-linear theory does not provide a universally accurate description in the resonant-scattering regime. It is highlighted here that the numerically calculated diffusion coefficients can be polluted for low energies due to missing resonant interaction possibilities of the particles with the turbulence. Five reduced-rigidity regimes are established, which are separated by analytical boundaries derived in this work. Consequently, a proper description of cosmic ray propagation can only be achieved by using a turbulence-level-dependent diffusion coefficient and can contribute to solving the Galactic cosmic ray gradient problem.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Reducing warpage in a hybrid large-scale additive manufacturing and compression molding process

In recent years, a hybrid manufacturing process, developed by combining extrusion-based large-scale additive manufacturing (AM) and compression molding (CM) techniques, has shown promising outcomes for producing structurally functional parts. The process can be used with both short fiber-reinforced composites and neat polymers and hence, even multi-material parts can be manufactured easily. This process offers the advantages of structural enhancement by having a desired fiber orientation using a large-scale AM process, as well as rapid manufacturing capability using a CM process. In the large-scale AM process, the alignment of fibers in the deposition direction enables significant improvement in the mechanical properties of the manufactured parts. However, the anisotropy resulting from the directional arrangement of fibers also introduces challenges related to warpage in the produced parts. This study aims to identify the causes of warpage and propose strategies to mitigate it. The research involves the use of preforms manufactured through the large-scale AM, which are then combined with the neat resin for CM manufacturing. A finite element-based numerical simulation model is developed, employing a sequentially coupled thermomechanical approach. Through a parametric study using the simulation models, optimization of printing direction and preform geometry is performed to minimize warpage. This contributes to the advancement and wider adoption of AM/CM hybrid manufacturing to produce structurally functional parts.

Jo, Eonyeon↗
DOE OSTI2024

Stress interference in multilayer additive friction stir deposition of AA6061 aluminum

Due to the multilayer deposition nature of metal additive manufacturing processes, each layer being printed experiences the state of thermokinetic and thermomechanical stress that in turn interfere with the state of thermokinetics and thermomechanical stress of subsequently deposited layers. Especially, this multilayer interference significantly affects the resultant properties of the component fabricated using solid-state additive friction stir deposition due to evolution of asymmetric state of planar stress. Due to the lack of comprehensive and suitable in situ diagnosis technique, the complex interference of inter- and multi-layer stresses during additive friction stir deposition was studied in an integrated approach of numerical simulation of fluidic state and experimental probing of stress influenced ultrasonic elastography. The uni-directional and bi-directional layer deposition configurations adopted during additive friction stir deposition result in the generation of constructive and destructive interference of the interlayer stress and hence, asymmetric and symmetric dynamic elasticity distribution respectively within the subsequent layers. With subsequent deposition of additional layers, the odd and even numbers of deposited layers generate asymmetric and nearly symmetric dynamic elasticity distributions.

Yang, Teng↗
DOE OSTI2024

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗
DOE OSTI2023

Wigner time delay in photoionization: a 1D model study

Abstract In scattering theory, the Wigner–Smith time delay, calculated through a phaseshift derivative or its multichannel generalization, has been demonstrated to measure the amount of delay or advance experienced by colliding particles during their interaction with the scattering potential. Fetic, Becker, and Milosevic argue that this concept cannot be extended to include photoionization, viewed as a half-scattering experiment. Their argument is based on the lack of information about scattering phaseshifts in the part of the wavefunction (satisfying the ingoing-wave boundary condition) going to the detector. This article aims to test this claim by examining a photodetachment process in a simple 1D model with a short-range symmetrical potential. Using time-dependent perturbation theory with a dipole interaction, the relevant wavepacket of the outgoing particle is analyzed and compared to the free wavepacket as a reference. Our findings confirm that, indeed, a time delay arises in the liberated fragmentation wavepacket, which is expressed as an energy derivative of the scattering phaseshift. Our study highlights that the source of the phaseshift content in the wavepacket arriving at the detector is the dipole matrix element, which is a direct consequence of imposing the ingoing-wave boundary condition. We illustrate our results through numerical simulations of both the non-free and free wavepackets. The amount of the observed time delay is found to be half of that appearing in a typical scattering experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Lattice QED in an external magnetic field: Evidence for dynamical chiral symmetry breaking

We simulate quantum electrodynamics (QED) in a strong constant homogeneous external magnetic field on a Euclidean space-time lattice using the rational hybrid Monte Carlo method, developed for simulating lattice quantum chromodynamics (QCD). Our primary goal is to measure the chiral condensate in the limit when the input electron mass m is zero. We observe a nonzero value, indicating that the external magnetic field catalyzes chiral symmetry breaking as predicted by approximate truncated Schwinger-Dyson methods. Such behavior is associated with dominance by the lowest Landau level which causes the effective dimensional reduction from 3 + 1 dimensions to 1 + 1 dimensions for charged particles (electrons and positrons) where the attractive forces of QED can produce chiral symmetry breaking with a dynamical electron mass and associated chiral condensate. Since our lattice simulations use bare (lattice) parameters, while the Schwinger-Dyson analyses work with renormalized quantities, direct numerical comparison will require renormalization of our lattice results. Published by the American Physical Society 2024

Kogut, J. B. (ORCID:000000030151548X)↗