Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite scale equations”

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

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗
DOE OSTI2022

Modelling of electrostatic ion-scale turbulence in divertor tokamaks with the gyrokinetic code $\mathrm{COGENT}$

Here, continuum gyrokinetic simulations of electrostatic ion scale turbulence are presented for the case of a diverted (single-null) tokamak geometry. The simulation model, implemented in the finite-volume code COGENT, solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to a vorticity equation for electrostatic potential variations, where a fluid model is used for an electron response. The model describes the ion scale ion temperature gradient (ITG) and resistive drift modes as well as neoclassical ion physics effects. Different turbulence regimes are observed depending on the plasma profiles, and the roles of a self-consistent background electric field and an X-point geometry are explored. In particular, increasing the pedestal density gradient and the corresponding radial electric field is demonstrated to suppress the ITG turbulence, whereas the same edge plasma background can still be destabilized by the resistive modes when the plasma resistivity is increased. The effects of X-point geometry are assessed by comparing cross-separatrix simulations with counterpart calculations performed for a toroidal annulus geometry. For the simulation parameters considered, similar global behaviour is observed in both cases, whereas strong local suppression of turbulence fluctuations is demonstrated near the X-point for the case of a single-null geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Advanced Turbulence Models for Large-Scale Atmospheric Boundary Layer Flows

We present high-fidelity large-eddy-simulation (LES) modeling approaches for the turbulent atmospheric boundary layer (ABL) flows. Wind energy is a prime example of an application driven by ABL. Generation of electrical energy from farms of wind turbines at night in the stable ABL is a particularly interesting situation. In this report, we consider the well-known GEWEX (Global Energy and Water Cycle Experiment) Atmospheric Boundary Layer Study (GABLS) stably stratified benchmark LES case. We use a high-order spectral element code Nek5000/RS, which is supported under the DOE's Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms. In our earlier ANL report, we demonstrated our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. In this report, we extend the range of our SGS modeling approaches in the context of the mean-field eddy viscosity (MFEV), to include the solution of an SGS turbulent kinetic energy equation (TKE). We demonstrate the model fidelity of Nek5000/RS in comparison to that of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, with which we studied scaling performance for both codes in comparison on DOE's leadership computing platforms.

17 WIND ENERGY↗
DOE OSTI2023

High-Fidelity Simulations of Plasma-Assisted Oxidation of Hydrocarbon Fuels Using Nanosecond Pulsed Discharges

Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.

adaptive mesh refinement↗
DOE OSTI2026

A sweeping positivity-preserving high-order finite difference WENO scheme for Euler equations

We develop a simple, high-order, conservative and robust positivity-preserving sweeping procedure for the density and the nonlinear pressure function in the compressible Euler equations. Using the scaling limiter in Zhang and Shu (J Comput Phys 229:3091–3120, 2010), we obtain a non-trivial extension of the scalar sweeping technique in Liu et al. (J Sci Comput 73:1028–1071, 2017) for the positivity of pressure. The sweeping procedure developed in this paper is a post-processing technique, which can be applied to any concave functions of the conserved variables in hyperbolic conservation law systems. Thus, it has applications beyond the Euler equations. This procedure preserves positivity and conservation of physical quantities without destroying the accuracy of the underlying scheme. The algorithm works for general schemes including finite difference, finite volume and discontinuous Galerkin methods; however, in this paper we focus on finite difference weighted essentially non-oscillatory (WENO) methods. As a result, we provide numerical tests of the fifth-order finite difference WENO scheme to demonstrate the accuracy and robustness of the technique.

Compressible Euler equations↗
DOE OSTI2024

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Adaptive Space-Time Methods for Large Scale Optimal Design

When modeling complex physical systems with advanced dynamics, such as shocks and singularities, many classic methods for solving partial differential equations can return inaccurate or unusable results. One way to resolve these complex dynamics is through r-adaptive refinement methods, in which a fixed number of mesh points are shifted to areas of high interest. The mesh refinement map can be found through the solution of the Monge-Ampére equation, a highly nonlinear partial differential equation. Due to its nonlinearity, the numerical solution of the Monge-Ampére equation is nontrivial and has previously required computationally expensive methods. In this report, we detail our novel optimization-based, multigrid-enabled solver for a low-order finite element approximation of the Monge-Ampére equation. This fast and scalable solver makes r-adaptive meshing more readily available for problems related to large-scale optimal design. Beyond mesh adaptivity, our report discusses additional applications where our fast solver for the Monge-Ampére equation could be easily applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Fast Solver for the Fractional Helmholtz Equation

The purpose of this paper is to study a Helmholtz problem with a spectral fractional Laplacian, instead of the standard Laplacian. Recently, it has been established that such a fractional Helmholtz problem better captures the underlying behavior in Geophysical Electromagnetics. In this work, we establish the well-posedness and regularity of this problem. We introduce a hybrid spectral-finite element approach to discretize it and show well-posedness of the discrete system. In addition, we derive a priori discretization error estimates. Finally, we introduce an efficient solver that scales aswell as the best possible solver for the classical integer-order Helmholtz equation. We conclude withseveral illustrative examples that confirm our theoretical findings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Butterfly Factorization Via Randomized Matrix-Vector Multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of an m × n matrix with m ≈ n provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of O(log n) sparse factors, each containing O(n) nonzero entries. The factorization can be attained using O(n 3/2 log n) computation and O(n log n) memory resources. Furthermore, the proposed algorithm can be implemented in parallel and can apply to matrices with strong or weak admissibility conditions arising from surface integral equation solvers as well as multi-frontal-based finite-difference, finite-element, or finite-volume solvers. A distributed-memory parallel implementation of the algorithm demonstrates excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Scalability of high-performance PDE solvers

Performance tests and analyses are critical to effective high-performance computing software development and are central components in the design and implementation of computational algorithms for achieving faster simulations on existing and future computing architectures for large-scale application problems. In this article, we explore performance and space-time trade-offs for important compute-intensive kernels of large-scale numerical solvers for partial differential equations (PDEs) that govern a wide range of physical applications. We consider a sequence of PDE-motivated bake-off problems designed to establish best practices for efficient high-order simulations across a variety of codes and platforms. We measure peak performance (degrees of freedom per second) on a fixed number of nodes and identify effective code optimization strategies for each architecture. In addition to peak performance, we identify the minimum time to solution at 80% parallel efficiency. The performance analysis is based on spectral and p-type finite elements but is equally applicable to a broad spectrum of numerical PDE discretizations, including finite difference, finite volume, and h-type finite elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Turbulence-driven flow dynamics in general axisymmetric toroidal geometry

The present work gives the equations governing the generation of toroidally axisymmetric flows by turbulent Reynolds and Maxwell stresses in finite aspect ratio, general cross section tokamak plasmas. Inclusion of the divergence-free flow constraint in lowest order changes the time scale for evolution of the poloidal flows driven by turbulence by substantial factors. In the pedestal region for present day machines, comparing to earlier cylindrical models, the time scale evaluated using a large aspect ratio circular cross section model can be two orders of magnitude longer while the present, general geometry result can be about one order of magnitude longer. Inclusion of gyroviscosity in the calculation shows that the only lowest order radial velocity fluctuations that enter the problem are those due to fluctuating E×B flows. Toroidal and poloidal flow effects on the toroidally axisymmetric flows are inextricably coupled due to the neoclassical poloidal viscosity. Accordingly, the physics is inherently three dimensional and measurements of all three velocity components are required to obtain the information needed to quantitatively test the theory. The parallel and angular momentum equations for the lowest order, toroidally axisymmetric flows look like radial transport equations when the turbulence is included. The turbulence terms provide the radial transport fluxes. In standard neoclassical theory, the parallel flow equation is local on each flux surface; there is no radial derivative term. However, adding turbulence gives a way, in principle, for radial transport to lead to poloidal flows that deviate from the neoclassical prediction. As a result, the inclusion of the Maxwell stress provides a mechanism for MHD fluctuations to alter the toroidally axisymmetric flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Extending the predictive power of perturbative QCD using the principle of maximum conformality and the Bayesian analysis

In addition to the evaluation of high-order loop contributions, the precision and predictive power of perturbative QCD (pQCD) predictions depends on two important issues: (1) how to achieve a reliable, convergent fixed-order series, and (2) how to reliably estimate the contributions of unknown higher-order terms. The recursive use of renormalization group equation, together with the Principle of Maximum Conformality (PMC), eliminates the renormalization scheme-and-scale ambiguities of the conventional pQCD series. The result is a conformal, scale-invariant series of finite order which also satisfies all of the principles of the renormalization group. In this paper we propose a novel Bayesian-based approach to estimate the size of the unknown higher order contributions based on an optimized analysis of probability distributions. We show that by using the PMC conformal series, in combination with the Bayesian analysis, one can consistently achieve high degree of reliability estimates for the unknown high order terms. Thus the predictive power of pQCD can be greatly improved. We illustrate this procedure for two pQCD observables: R e+e– and R τ , which are each known up to four loops in pQCD. Numerical analyses confirm that by using the scale-independent and more convergent PMC conformal series, one can achieve reliable Bayesian probability estimates for the unknown higher-order contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗