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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A new approach to the evaluation and solution of the relativistic kinetic dispersion relation and verification with continuum kinetic simulation

Here, the present work describes a new approach to evaluation and root finding for the kinetic dispersion relation of Langmuir waves, which is central to the analytical understanding of collisionless damping in plasmas. The plasma dispersion function is solved to machine precision using direct integration in the complex plane in combination with an analytic evaluation of the residue to account for the deformation along the Landau contour. To efficiently attain machine precision, the contour is displaced in the complex plane prior to integration, and numerical subtleties related to the placement of the contour are discussed. The approach is generic in that it applies to arbitrary distribution functions, with the present manuscript focused on relativistic cases. Detailed verification of results via direct kinetic simulation in a variety of configuration space dimensions is also presented. Finally, the technique is applied to the challenging case of highly relativistic (i.e. extremely hot) plasmas. Here we show both qualitative agreement with prior work, as well as the disappearance of the Landau root which would have significant implication for real-life observation or experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Protection Against Graph-Based False Data Injection Attacks on Power Systems

Graph signal processing (GSP) has emerged as a powerful tool for practical network applications, including power system monitoring. By representing power system voltages as smooth graph signals, recent research has focused on developing GSP-based methods for state estimation, attack detection, and topology identification. Included, efficient methods have been developed for detecting false data injection (FDI) attacks, which until now were perceived as non-smooth with respect to the graph Laplacian matrix. Consequently, these methods may not be effective against smooth FDI attacks. In this paper, we propose a graph FDI (GFDI) attack that minimizes the Laplacian-based graph total variation (TV) under practical constraints. In addition, we develop a low-complexity algorithm that solves the non-convex GDFI attack optimization problem using ell_1-norm relaxation, the projected gradient descent (PGD) algorithm, and the alternating direction method of multipliers (ADMM). We then propose a protection scheme that identifies the minimal set of measurements necessary to constrain the GFDI output to high graph TV, thereby enabling its detection by existing GSP-based detectors. Our numerical simulations on the IEEE-57 bus test case reveal the potential threat posed by well-designed GSP-based FDI attacks. Moreover, we demonstrate that integrating the proposed protection design with GSP-based detection can lead to significant hardware cost savings compared to previous designs of protection methods against FDI attacks.

Morgenstern, Gal↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗
DOE OSTI2020

Numerical investigation of microbial quorum sensing under various flow conditions

Microorganisms efficiently coordinate phenotype expressions through a decision-making process known as quorum sensing (QS). We investigated QS amongst distinct, spatially distributed microbial aggregates under various flow conditions using a process-driven numerical model. Model simulations assess the conditions suitable for QS induction and quantify the importance of advective transport of signaling molecules. In addition, advection dilutes signaling molecules so that faster flow conditions require higher microbial densities, faster signal production rates, or higher sensitivities to signaling molecules to induce QS. However, autoinduction of signal production can substantially increase the transport distance of signaling molecules in both upstream and downstream directions. We present empirical approximations to the solutions of the advection–diffusion–reaction equation that describe the concentration profiles of signaling molecules for a wide range of flow and reaction rates. These empirical relationships, which predict the distribution of dissolved solutes along pore channels, allow to quantitatively estimate the effective communication distances amongst multiple microbial aggregates without further numerical simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Numerical simulations of onset and growth of Rayleigh–Taylor instability involving solids in converging geometry

Numerical simulation of Rayleigh–Taylor instabilities involving solids has been carried out in a spherically converging geometry. A solid shell is driven toward the center by high-pressure gas outside. Low-pressure gas inside the shell is then compressed, causing the solid shell to slow down and eventually bounce back. Here, simulation results show the absence of the onset condition unlike in the planar geometry equivalent. Results also show that disturbance growth seems to be strongly associated with the deformation direction instead of the acceleration direction. For example, when the interface is compressed in the cross-stream direction, the disturbance at the gas–solid interface grows even under the Rayleigh–Taylor stable situation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Continuum gyrokinetic simulations of edge plasmas in single-null geometries

The first continuum gyrokinetic calculations of electrostatic ion scale turbulence are presented for the case of a diverted tokamak geometry. The simulation model solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to the quasi-neutrality equation for electrostatic potential variations, where a fluid model is used for an electron response. Furthermore, in order to facilitate simulations of highly-anisotropic microturbulence in the presence of strong magnetic shear and a magnetic X-point, a numerical algorithm utilizing a locally field-aligned multiblock coordinate system has been developed and implemented in the 5D finite-volume code COGENT. In this approach, the toroidal direction is divided into blocks, such that within each block, the cells are field-aligned and a non-matching grid interface is allowed at block boundaries. The toroidal angle corresponds to the “coarse” field-aligned coordinate, whereas the poloidal cross section, comprised of the radial and poloidal directions, is finely gridded to resolve short-scale perpendicular turbulence structures and to support accurate re-mapping (interpolation) at block boundaries. The 5D simulations explore cross-separatrix ion scale turbulence in the presence of a self-consistent radial electric field and address the effects of magnetic-shear stabilization in the X-point region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Trace impurity transport in multi-species plasmas with large particle fluxes

A burning fusion plasma has a large inwards flux of fuel and outwards flux of ash. Existing impurity transport theories do not account for a steady-state environment of such large fluxes. In this paper, we extend classical transport theory to account for a background of such fluxes. In a mainly two-ion species magnetized plasma, with oppositely directed density gradients maintained by oppositely directed sources of the two species, trace impurity ions are shown to align with the larger mass ions. This unexpected behavior is derived analytically and simulated numerically using the MITNS (Multiple-Ion Transport Numerical Solver) code. Overall, the result suggests a potential advantage of edge-fueled p-B11 fusion in extracting high-Z impurities, and possible applications in plasma-based separation methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Preferential turbulence enhancement in two-dimensional compressions

When initially isotropic three-dimensional (3D) turbulence is compressed along two dimensions, the compression supplies energy directly to the flow components in the compressed directions, while the flow component in the noncompressed direction experiences the effects of compression only indirectly through the nonlinearity of the hydrodynamic equations. Here we present our study of such 2D compressions using numerical simulations. For initially isotropic turbulence, we find that the nonlinearity can be insufficient to maintain isotropy, with the energy components parallel to the compression coming to dominate the turbulent energy, with a number of consequences. Among these are the possibilities for stronger and more easily sustained growth of turbulent energy than in 3D compressions and for an increasing turbulent Mach number even in a compression without thermal losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Monte-Carlo Neutrino Transport in Neutron Star Merger Simulations

Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage observations of these systems, we need numerical simulations that provide reliable predictions for the properties of the matter unbound in these mergers. An important limitation of current simulations is the use of approximate methods for neutrino transport that do not converge to a solution of the transport equations as numerical resolution increases, and thus have errors that are impossible to quantify. Here, we report on a first simulation of a binary neutron star merger that directly solves the transport equations using Monte-Carlo techniques. We demonstrate that a computationally cheap Monte-Carlo transport can be implemented without ill effects from sampling noise. We apply Monte-Carlo transport to an unequal mass neutron star binary merger up to 5ms past merger, and report on the properties of the matter and neutrino outflows. Finally, we compare our results to the output of our best approximate `M1' transport scheme, demonstrating that an M1 scheme that carefully approximates the neutrino energy spectrum only leads to ~10% uncertainty in the composition and velocity of the ejecta, and ~20% uncertainty in the electron (anti)neutrino luminosities and energies. Here, the most significant disagreement found between M1 and Monte-Carlo results is a factor of ~2 difference in the luminosity of heavy-lepton neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Imaging surface topography with coherent x-ray reflectivity: Theory, kinematics, and simulations

A theoretical formalism is described for understanding coherent x-ray reflectivity (CXR) from the surface of a semi-infinite crystal having a variable surface topography, described by the height profile ℎ(𝑥,𝑦). The surface topography is imaged as a complex “effective density,” obtained from the phasing and inversion of the coherent x-ray reflectivity data, measured through a rocking scan centered at a vertical momentum transfer 𝑄$^{0}_{𝑧}$ and a vertical range Δ⁢𝑄 𝑧 . The formalism predicts that the effective density has an amplitude with a maximum located at the surface height for each position within the surface plane. The phase of the effective density has a lateral variation that is controlled by the surface height and a vertical variation that reflects a combination of the interfacial structure and specific choice of measurement conditions. This understanding enables direct observation of nanometer-scale interfacial topography, i.e., ℎ⁡(𝑥,𝑦)⁢𝑐 𝑠 (where 𝑐 𝑠 is the vertical substrate lattice parameter) with Å-scale sensitivity to surface height. Numerical simulations illustrate and confirm the theoretical results. These results show how the interpretation of the interfacial density phase obtained by CXR data inversion (i.e., surface topography with respect to a flat surface) is conceptually similar to that previously known for Bragg coherent diffraction imaging (BCDI) measurements of isolated nanoparticles (i.e., lattice displacements with respect to an ideal crystal lattice). This suggests that CXR can be thought of as a form of dark field imaging with respect to the bright field BCDI approach. An implication of these results is that interfacial imaging may bypass some of the significant challenges associated with BCDI imaging of multiple particles having different orientations.

X-ray imaging↗
DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine-learned closure of URANS for stably stratified turbulence: connecting physical timescales & data hyperparameters of deep time-series models

Stably stratified turbulence (SST), a model that is representative of the turbulence found in the oceans and atmosphere, is strongly affected by fine balances between forces and becomes more anisotropic in time for decaying scenarios. Moreover, there is a limited understanding of the physical phenomena described by some of the terms in the Unsteady Reynolds-Averaged Navier–Stokes (URANS) equations—used to numerically simulate approximate solutions for such turbulent flows. Rather than attempting to model each term in URANS separately, it is attractive to explore the capability of machine learning (ML) to model groups of terms, i.e. to directly model the force balances. We develop deep time-series ML for closure modeling of the URANS equations applied to SST. We consider decaying SST which are homogeneous and stably stratified by a uniform density gradient, enabling dimensionality reduction. We consider two time-series ML models: long short-term memory and neural ordinary differential equation. Both models perform accurately and are numerically stable in a posteriori (online) tests. Furthermore, we explore the data requirements of the time-series ML models by extracting physically relevant timescales of the complex system. We find that the ratio of the timescales of the minimum information required by the ML models to accurately capture the dynamics of the SST corresponds to the Reynolds number of the flow. The current framework provides the backbone to explore the capability of such models to capture the dynamics of high-dimensional complex dynamical system like SST flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

On the subject of large-scale pool fires and turbulent boundary layer interactions

The role to which a realistic inflow turbulent boundary layer (TBL) influences transient and mean large-scale pool fire quantities of interest (QoIs) is numerically investigated. High-fidelity, low-Mach large-eddy simulations that activate low-dissipation, unstructured numerics are conducted using an unsteady flamelet combustion modeling approach with mutiphysics coupling to soot and participating media radiation transport. Three inlet profile configurations are exercised for a large-scale, high-aspect rectangular pool that is oriented perpendicular to the flow direction: a time-varying, TBL inflow profile obtained from a periodic precursor simulation, the time-mean of the transient TBL, and a steady power-law inflow profile that replicates the mean TBL crosswind velocity of 10.0 m/s at a vertical height of 10 m. Results include both qualitative transient flame evolution and quantitative flame shape with ground-level temperature and convective/radiative heat flux profiles. While transient fire events, which are driven by burst-sweep TBL coupling, such as blow-off and reattachment are vastly different in the TBL case (contributing to increased root mean square QoI fluctuation prediction and disparate flame lengths), mean surface QoI magnitudes are similar. Quadrant analysis demonstrates that the TBL configuration modifies burst-sweep phenomena at windward pool locations, while leeward recovery is found. In conclusion, positive fluctuations of convective heat flux correlate with fast moving fluid away from the pool surface due to intermittent combustion events.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Investigation of Best-Practices and Computationally Inexpensive Radiative Exchange Models for Discrete Element Method Modeling of Aluminosilicate Particles in Concentrating Solar Power Environments

Chemically inert, aluminosilicate based particles have been investigated as both a thermal transport and sensible energy storage medium for concentrating solar power facilities. These particles will experience a wide range of operating temperatures (300-1000 K) and handling conditions (dense to dilute falling particle curtains, dense granular flows, or dense structures), requiring specially-designed and optimized infrastructures. The relative influence of collisional and frictional interactions between particles varies based on temperature-dependent particulate properties and greatly impacts the bulk, granular flow behavior. These underlying physics are captured using discrete element method modeling tools. However, this modeling method is computationally expensive as each particle position and interaction is tracked during the simulation. These modeling methods are further complicated by introducing temperature-dependent particle properties, high-temperature radiative exchange, and directional irradiation sources experienced by granular flows in concentrating solar power environments. In this study, coupled experimental and numerical slump testing of aluminosilicate particles was performed and computationally efficient radiative exchange models were evaluated to establish best-practices for discrete element method models for concentrating solar power environments. The three particle types investigated included Carbobead HSP 30 /60, Carbobead CP 30/60, and Granusil 4030. Existing modeling limitations and computationally-efficient multi-modal heat transfer models were evaluated using Aspherix®, a commercial discrete element method software. High-temperature (< 1073 K) slump testing of aluminosilicate particles was performed to investigate the deviation between experimentally-observed and numerically-predicted angles of repose introduced by computation-time reduction practices including the relaxation of the particle elastic modulus and coarse-graining. Coarse-graining is used to use a single modeled particle that is representative of a collection of smaller particles, decreasing the computational cost at the expense of geometric accuracy. Additionally, relaxation of the elastic modulus is used to reduce computational time at the expense of an increased, modeled particle overlap. Prior studies have determined that aluminosilicate particles retain a high elastic modulus at high temperatures (< 1073 K), requiring small simulation timesteps to ensure resolved contact forces resemble appropriate solid mechanics. A parametric study was performed to evaluate the influence of computation time improvements on the deviation between experimental and modeled angle of repose across high temperatures < 1073 K. Additionally, numerical case studies were performed on candidate particle systems at varying porosities and temperatures. These studies were performed to investigate the influence of computationally-efficient radiative-exchange modeling methods coupled to Aspherix® on modeled accuracy and computation time. The recently-developed distance-based approximation was evaluated in estimating radiative exchange between particles and participating surfaces located in close proximity. The distance based approximation was developed to use tabulated estimates of the radiative distribution factor between individual particles and surfaces in close proximity (< 40 particle radii). These methods were expanded to the aluminosilicate particles of interest, including the influence of particle size distributions. To capture radiative exchange between particles and surfaces not in close proximity (> 40 particle radii) and to capture the absorption of directional irradiation from concentrating solar resources, a volumetrically-averaged radiative distribution factor was calculated between the modeled granular flow and surfaces using Monte Carlo ray-tracing for participating media. Volume-averaged absorption and scattering coefficients were predicted using a volumetric discretization of the modeled domain with monodisperse approximations based on geometric optics and experimentally-determined scattering phase functions for aluminosilicate particles.

14 SOLAR ENERGY↗
DOE OSTI2020

Programmable topotaxis of magnetic rollers in time-varying fields

We discuss how spatially uniform, time-periodic magnetic fields can be designed to power and direct the migration of ferromagnetic spheres up (or down) local gradients in the topography of a solid substrate. Our results are based on a dynamical model that considers the time-varying magnetic torques on the particle and its motion through the fluid at low Reynolds number. We use both analytical theory and numerical simulation to design magnetic fields that maximize the migration velocity up (or down) an inclined plane. We show how “topotaxis” of spherical particles relies on differences in the hydrodynamic resistance to rotation about axes parallel and perpendicular to the plane. Importantly, the designed fields can drive multiple independent particles to move simultaneously in different directions as determined by gradients in their respective environments. Experiments on ferromagnetic spheres provide evidence for topotactic motions up inclined substrates. The ability to program the autonomous navigation of driven particles within anisotropic environments is relevant to the design of colloidal robots.

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

Multicycle large-eddy simulations of a direct-injection hydrogen-fueled optical engine

Hydrogen (H 2 ) is a carbon-free chemical energy carrier and one promising solution for achieving effective decarbonization of the transportation sector, particularly for internal combustion engines (ICEs). With a focus on ICEs, and compared to port-fuel injection, direct injection (DI) of gaseous H 2 during the compression stroke offers potential advantages, which include backfire avoidance and reduction of preignition occurrence. In these last two decades, much research, experimental and numerical, has been devoted to understanding H 2 's mixing and combustion processes in ICEs. Computational fluid dynamics modeling efforts commonly rely on unsteady Reynolds-averaged Navier Stokes (URANS) turbulence frameworks, mostly due to their computational affordability. However, many authors have pointed out the opportunity to perform large-eddy simulations (LESs) to investigate the cyclic variability of H 2 engines and assess potential advantages of using LES in place of URANS, especially for lean operation. This study addresses this knowledge gap and presents a computational fluid dynamics (CFD) study of the H 2 DI process in an optical engine operating at relatively low tumble conditions, using multicycle LESs. In conclusion, the manuscript presents a thorough validation of the results against experimental data available from the literature as well as direct comparison with URANS, demonstrating the feasibility of multicycle LESs for CFD modeling of DI H 2 -fueled ICEs.

Direct injection↗
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
DOE OSTI2026

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗