SEARCH · Engineering Papers

Results for “Grid singularity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Interacting tilt and kink instabilities in repelling current channels

We present a numerical study in resistive magnetohydrodynamics (MHD) where the initial equilibrium configuration contains adjacent, oppositely directed, parallel current channels. Since oppositely directed current channels repel, the equilibrium is liable to an ideal magnetohydrodynamic tilt instability. This tilt evolution, previously studied in planar settings, involves two magnetic islands or flux ropes, which on Alfvénic timescales undergo a combined rotation and separation. This in turn leads to the creation of (near) singular current layers, posing severe challenges to numerical approaches. Using our open-source grid-adaptive MPI-AMRVAC software, we revisit the planar evolution case in compressible MHD, as well as its extension to two-and-a-half-dimensional (2.5D) and full three-dimensional (3D) scenarios. As long as the third dimension can be ignored, pure tilt evolutions result that are hardly affected by out of plane magnetic field components. In all 2.5D runs, our simulations do show secondary tearing type disruptions throughout the near singular current sheets in the far nonlinear saturation regime. In full 3D runs, both current channels can be liable to additional ideal kink deformations. We discuss the effects of having both tilt and kink instabilities acting simultaneously in the violent, reconnection-dominated evolution. In 3D, both the tilt and the kink instabilities can be stabilized by tension forces. As a concrete space plasma application, we argue that interacting tilt-kink instabilities in repelling current channels provide a novel route to initiate solar coronal mass ejections, distinctly different from the currently favored pure kink or torus instability routes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Integration of microseismic monitoring data into coupled flow and geomechanical models with ensemble Kalman filter

Abstract Hydraulic stimulation of low‐permeability rocks in enhanced geothermal systems, shale resources, and CO 2 storage aquifers can trigger microseismic events, also known as microearthquakes (MEQs). The distribution of microseismic source locations in the reservoir may reveal important information about the distribution of hydraulic and geomechanical rock properties. In this paper, we present a framework for conditioning heterogeneous rock permeability and geomechanical property distributions on microseismic data. To simulate the multiphysics processes in these systems, we combine a fully coupled flow and geomechanical model with the Mohr‐Coulomb type rock failure criterion. The resulting multiphysics simulation constitutes the forecast model that relates microseismic source locations to reservoir rock properties. We adopt this forward model in an ensemble Kalman filter (EnKF) data assimilation framework to jointly estimate reservoir permeability and geomechanical property distributions from injection‐induced microseismic response measurements. We show that integration of a large number of spatially correlated microseismic data with practical ensemble sizes can lead to severe underestimation of ensemble spread, and eventually ensemble collapse. To mitigate the variance underestimation issue, two low‐rank data representation schemes are presented and discussed. In the first approach, microseismic data are projected onto a low‐dimensional subspace defined by the left singular vectors of the perturbed observation matrix. The second method uses a coarser grid for representing the microseismic data. A series of numerical experiments is presented to evaluate the performance of the proposed methods and to illustrate their applicability for assimilating microseismic data into coupled flow and geomechanical forward models to estimate multiphysics rock properties.

Tarrahi, Mohammadali↗

DOE OSTI2015

Euler–Lagrange Simulations of Bubble Cloud Dynamics Near a Wall

We present in this paper a two-way coupled Eulerian–Lagrangian model to study the dynamics of clouds of microbubbles subjected to pressure variations and the resulting pressures on a nearby rigid wall. The model simulates the two-phase medium as a continuum and solves the Navier–Stokes equations using Eulerian grids with a time and space varying density. The microbubbles are modeled as interacting singularities representing moving and oscillating spherical bubbles, following a modified Rayleigh–Plesset–Keller–Herring equation and are tracked in a Lagrangian fashion. A two-way coupling between the Euler and Lagrange components is realized through the local mixture density determined by the bubbles' volume change and motion. Using this numerical framework, simulations involving a large number of bubbles were conducted under driving pressures at different frequencies. The results show that the frequency of the driving pressure is critical in determining the overall dynamics: either a collective strongly coupled cluster behavior or nonsynchronized weaker multiple bubble oscillations. The former creates extremely high pressures with peak values orders of magnitudes higher than that of the excitation pressure. This occurs when the driving frequency matches the natural frequency of the bubble cloud. The initial distance between the bubble cloud and the wall also affects significantly the resulting pressures. A bubble cloud collapsing very close to the wall exhibits a cascading collapse, with the bubbles farthest from the wall collapsing first and the nearest ones collapsing last, thus the energy accumulates and this results in very high pressure peaks at the wall. Finally, at farther cloud distances from the wall, the bubble cloud collapses quasi-spherically with the cloud center collapsing last.

42 ENGINEERING↗

DOE OSTI2020

Editors’ Choice—Perspective—Challenges in Moving to Multiscale Battery Models: Where Electrochemistry Meets and Demands More from Math

There has been significant recent interest in studying multiscale characteristics of current and next-generation batteries, including lithium-metal and lithium-sulfur batteries. Advances in computing power make researchers believe that the detailed multiscale models can be efficiently simulated to arrive at the insights for the degradation and performance loss; however, this is not true and special attention needs to be paid to local singularities, boundary layers, moving boundaries, etc. This article presents 2D examples that illustrate the importance of grid convergence studies, provides well-defined detailed models to test the efficiency of numerical schemes, and discusses the associated simulation challenges.

Shah, Krishna (ORCID:0000000278026361)↗

DOE OSTI2014

Hybrid energy systems (HESs) using small modular reactors (SMRs)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2014

SMR Handbook: Hybrid Energy Systems Involving SMRs

Large-scale nuclear reactors are traditionally operated for a singular purpose: steady-state production of dispatchable baseload electricity that is distributed broadly on the electric grid. While this implementation is key to a sustainable, reliable energy grid, SMRs offer new opportunities for increased use of clean nuclear energy for both electric and thermal applications in more locations – while still accommodating the desire to support renewable production sources. This chapter considers a scenario in which renewable generation would be tightly coupled with the nuclear generation source – behind the grid – to meet the grid demand as an integrated energy system while simultaneously producing other commodities with the available thermal energy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

Effect of impurity doping on tunneling conductance in AB-stacked bi-layer graphene: A tight-binding study

We report here a microscopic tight-binding model calculation for AB-stacked bilayer graphene in presence of biasing potential between the two layers and the impurity effects to study the evolution of the total density of states with special emphasis on opening of band gap near Dirac point. We have calculated the electron Green’s functions for both the A and B sub-lattices by Zubarev technique. The imaginary part of the Green’s function gives the partial and total density of states of electrons. The density of states are computed numerically for 1000 × 1000 grid points of the electron momentum. The evolution of the opening of band gap near van-Hove singularities as well as near Dirac point is investigated by varying the different interlayer hoppings and the biasing potentials. The inter layer hopping splits the density of states at van-Hove singularities and produces a V-shaped gap near Dirac point. Further the biasing potential introduces a U shaped gap near Dirac point with a density minimum at the applied potential(i.e. at V/2).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Reliability and Geographic Trends of 50,000 Photovoltaic Systems in the USA: Preprint

This paper presents performance and reliability data from nearly 50,000 photovoltaic (PV) systems totaling 1.7 gigawatts installed capacity in the USA from 2009 to 2012 and their geographic trends. About 90% of the normal systems and about 85% of all systems, including systems with known issues, performed to within 10% or better of expected performance. Although considerable uncertainty may exist due to the nature of the data, hotter climates appear to exhibit some degradation not seen in the more moderate climates. Special causes of underperformance and their impacts are delineated by reliability category. Hardware-related issues are dominated by inverter problems (totaling less than 0.5%) and underperforming modules (totaling less than 0.1%). Furthermore, many reliability categories show a significant decrease in occurrence from year 1 to subsequent years, emphasizing the need for higher-quality installations but also the need for improved standards development. The probability of PV system damage because of hail is below 0.05%. Singular weather events can have a significant impact such as a single lightning strike to a transformer or the impact of a hurricane. However, grid outages are more likely to have a significant impact than PV system damage when extreme weather events occur.

14 SOLAR ENERGY↗

DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Singularity-EOS XCAP Report

Solving the Euler equations is a fundamental component of simulating many physical phenomena ranging from high explosives to astrophysics. An equation of state (EOS) is a required piece that relates any two thermodynamic quantities to all other thermodynamic values. The presence of multiple materials within a control volume further complicates the solution requiring additional equations to describe the interaction of materials at a sub-grid level. Equations of state themselves can also come in many forms ranging from simple algebraic relations to more complicated differential equation models that describe material interactions over a broad range of physical conditions. In the latter case, the EOS is often pre-computed at a given set of grid points and provided in a tabular form where additional properties can be derived from the interpolation functions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗

DOE OSTI2024

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

DOE OSTI2018

A purely Lagrangian method for simulating the shallow water equations on a sphere using smooth particle hydrodynamics

It has long been suggested that a purely Lagrangian solution to global-scale atmospheric/oceanic flows can potentially outperform tradition Eulerian schemes. Meanwhile, a demonstration of a scalable and practical framework remains elusive. Motivated by recent progress in particle-based methods when applied to convection dominated flows, this work presents a fully Lagrangian method for solving the inviscid shallow water equations on a rotating sphere in a smooth particle hydrodynamics framework. To avoid singularities at the poles, the governing equations are solved in Cartesian coordinates, augmented with a Lagrange multiplier to ensure that fluid particles are constrained to the surface of the sphere. An underlying grid in spherical coordinates is used to facilitate efficient neighbor detection and parallelization. The method is applied to a suite of canonical test cases, and conservation, accuracy, and parallel performance are assessed..

Computer Science↗

DOE OSTI2019

Adaptive Grid Lens Modeling of the Cosmic Horseshoe Using Hubble Space Telescope Imaging

The Cosmic Horseshoe gravitational lensing system is an extraordinary example of strong gravitational lensing both due to the nearly complete Einstein ring formed in this system and due to the star-forming nature of the lensed $z$ = 2.38 Lyman-break galaxy. In this paper, we describe the development of a new lens modeling package and the lens models produced using the Hubble Space Telescope imaging. Our new lens modeling package uses adaptive grid methods, based especially on the algorithm described by Vegetti & Koopmans. The new lens modeling package introduces the K-means method to deal with multiple background sources. We utilize two parameterized models, the singular isothermal ellipsoid and the Navarro–Frenk–White model in order to optimize the Bayesian penalty function. Color–color diagrams show two distinct colors exist in the Einstein ring, which suggests that the Cosmic Horseshoe is formed from two background galaxies or from a highly irregular galaxy. We find that the best lens model includes two components and the total mass within the Einstein ring is $(5.43\pm 0.08)\times {10}^{12}\,{M}_{\odot }$. Furthermore, the background source reconstruction shows two peaks in the source plane and the most probable lens profile is consistent with the light profile of the foreground galaxy when comparing axis ratios and position angles.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

What do we mean by the word “Shock”?

From one vantage point, a shock is a continuous but drastic change in state variables that occurs over very small time and length scales. These scales and associated changes in state variables can be measured experimentally. From another vantage point, a shock is a mathematical singularity consisting of instantaneous changes in state variables. This more mathematical view gives rise to analytical solutions to idealized problems. And from a third vantage point, a shock is a structure in a hydrocode prediction. Its width depends on the simulation’s grid resolution and artificial viscosity. These three vantage points can be in conflict when ideas from the associated fields are combined, and yet combining them is an important goal of an integrated modeling program. This presentation explores an example of how models for real materials in the presence of real shocks react to a hydrocode’s numerical shocks of finite width. The presentation will include an introduction to plasticity for the novice, an historical view of plasticity algorithms, a demonstration of how pursuing the meaning of “shock” has resulted in hydrocode improvements, and will conclude by answering some of the questions that arise from that pursuit. After the technical part of the presentation, a few slides advertising LANL’s Computational Physics Student Summer Workshop will be shown.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Implicitly solving phase appearance and disappearance problems using two-fluid six-equation model

Phase appearance and disappearance issue presents serious numerical challenges in two-phase flow simulations using the two-fluid six-equation model. Numerical challenges arise from the singular equation system when one phase is absent, as well as from the discontinuity in the solution space when one phase appears or disappears. In this work, a high-resolution spatial discretization scheme on staggered grids and fully implicit methods were applied for the simulation of two-phase flow problems using the two-fluid six-equation model. A Jacobian-free Newton-Krylov (JFNK) method was used to solve the discretized nonlinear problem. An improved numerical treatment was proposed and proved to be effective to handle the numerical challenges. The treatment scheme is conceptually simple, easy to implement, and does not require explicit truncations on solutions, which is essential to conserve mass and energy. Various types of phase appearance and disappearance problems relevant to thermal-hydraulics analysis have been investigated, including a sedimentation problem, an oscillating manometer problem, a non-condensable gas injection problem, a single-phase flow with heat addition problem and a subcooled flow boiling problem. Successful simulations of these problems demonstrate the capability and robustness of the proposed numerical methods and numerical treatments. As a result, volume fraction of the absent phase can be calculated effectively as zero.

42 ENGINEERING↗

DOE OSTI2016

Discretizing singular point sources in hyperbolic wave propagation problems

Here, we develop high order accurate source discretizations for hyperbolic wave propagation problems in first order formulation that are discretized by finite difference schemes. By studying the Fourier series expansions of the source discretization and the finite difference operator, we derive sufficient conditions for achieving design accuracy in the numerical solution. Only half of the conditions in Fourier space can be satisfied through moment conditions on the source discretization, and we develop smoothness conditions for satisfying the remaining accuracy conditions. The resulting source discretization has compact support in physical space, and is spread over as many grid points as the number of moment and smoothness conditions. In numerical experiments we demonstrate high order of accuracy in the numerical solution of the 1-D advection equation (both in the interior and near a boundary), the 3-D elastic wave equation, and the 3-D linearized Euler equations.

58 GEOSCIENCES↗

DOE OSTI2018

Steep Cliffs and Saturated Exponents in Three-Dimensional Scalar Turbulence

Herein the intermittency of a passive scalar advected by three-dimensional Navier-Stokes turbulence at a Taylor-scale Reynolds number of 650 is studied using direct numerical simulations on a 4096 3 grid; the Schmidt number is unity. By measuring scalar increment moments of high orders, while ensuring statistical convergence, we provide unambiguous evidence that the scaling exponents saturate to 1.2 for moment orders beyond about 12, indicating that scalar intermittency is dominated by the most singular shocklike cliffs in the scalar field. We show that the fractal dimension of the spatial support of steep cliffs is about 1.8, whose sum with the saturation exponent value of 1.2 adds up to the space dimension of 3, thus demonstrating a deep connection between the geometry and statistics in turbulent scalar mixing. The anomaly for the fourth and sixth order moments is comparable to that in the Kraichnan model for the roughness exponent of 4/3.

42 ENGINEERING↗