Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear 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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗
DOE OSTI2022

Local correlation-based transition models for high-Reynolds-number wind-turbine airfoils

Abstract. Modern wind-turbine airfoil design requires robust performance predictions for varying thicknesses, shapes, and appropriate Reynolds numbers. The airfoils of current large offshore wind turbines operate with chord-based Reynolds numbers in the range of 3–15 million. Turbulence transition in the airfoil boundary layer is known to play an important role in the aerodynamics of these airfoils near the design operating point. While the lack of prediction of lift stall through Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) is well known, airfoil design using CFD requires the accurate prediction of the glide ratio (L/D) in the linear portion of the lift polar. The prediction of the drag bucket and the glide ratio is greatly affected by the choice of the transition model in RANS CFD of airfoils. We present the performance of two existing local correlation-based transition models – one-equation model (γ− SA) and two-equation model (γ-Reθt‾- SA) coupled with the Spalart–Allmaras (SA) RANS turbulence model – for offshore wind-turbine airfoils operating at a high Reynolds number. We compare the predictions of the two transition models with available experimental and CFD data in the literature in the Reynolds number range of 3–15 million including the AVATAR project measurements of the DU00-W-212 airfoil. Both transition models predict a larger L/D compared to fully turbulent results at all Reynolds numbers. The two models exhibit similar behavior at Reynolds numbers around 3 million. However, at higher Reynolds numbers, the one-equation model fails to predict the natural transition behavior due to early transition onset. The two-equation transition model predicts the aerodynamic coefficients for airfoils of various thickness at higher Reynolds numbers up to 15 million more accurately compared to the one-equation model. As a result, the two-equation model predictions are more comparable to the predictions from eN transition model. However, a limitation of this model is observed at very high Reynolds numbers of around 12–15 million where the predictions are very sensitive to the inflow turbulent intensity. The combination of the two-equation transition model coupled with the Spalart–Allmaras (SA) RANS turbulence model is a good method for performance prediction of modern wind-turbine airfoils using CFD.

17 WIND ENERGY↗
DOE OSTI2023

CFL Optimized Forward–Backward Runge–Kutta Schemes for the Shallow-Water Equations

Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Interference in nonlinear Compton scattering using a Schrödinger-equation approach

The interference between Compton scattering and nonlinear Compton scattering from a two-color field in the X-ray regime is theoretically examined for bound electrons. Here, the underlying phase shifts are analyzed using a perturbative approach in the incoming classical field. The perturbative approach is bench marked with a non-perturbative approach in the classical field. The interference for different combinations of linear polarization of the two fields is examined when the Compton and the nonlinear Compton scattered waves have the same wave vector and polarization. Only two cases exhibit interference. When there is interference, the calculations reveal an intrinsic phase difference between the Compton scattered wave function and the nonlinear Compton scattered wave function of either 0 or pi depending on the scattering angle.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Estimating basis functions in massive fields under the spatial mixed effects model

Abstract Spatial prediction is commonly achieved under the assumption of a Gaussian random field by obtaining maximum likelihood estimates of parameters, and then using the kriging equations to arrive at predicted values. For massive datasets, fixed rank kriging using the expectation–maximization algorithm for estimation has been proposed as an alternative to the usual but computationally prohibitive kriging method. The method reduces computation cost of estimation by redefining the spatial process as a linear combination of basis functions and spatial random effects. A disadvantage of this method is that it imposes constraints on the relationship between the observed locations and the knots. We develop an alternative method that utilizes the spatial mixed effects model, but allows for additional flexibility by estimating the range of the spatial dependence between the observations and the knots via an alternating expectation conditional maximization algorithm. Experiments show that our methodology improves estimation without sacrificing prediction accuracy while also minimizing the additional computational burden of extra parameter estimation. The methodology is applied to a temperature dataset archived by the United States National Climate Data Center, with improved results over previous methodology.

Pazdernik, Karl↗
DOE OSTI2021

Quantifying the anomalous water absorption behavior of cement mortar in view of its physical sensitivity to water

The one-dimensional water absorption into cement-based materials, possibly the simplest case of unsaturated water flow playing a fundamental role in predicting mass transport and durability, is still far from being completely understood due to several puzzling anomalies. Considering the water sensitivity of C-S-H gels (swell upon wetting and contract upon drying), a modified Richards equation (MRE) is proposed to quantify the anomalous water absorption into cement-based materials of changing porosity. Consistent experimental tests on a cement mortar and calculations indicate that the MRE with time-decaying water permeability can quantitatively capture the long-term absorption with two linear stages and a nonlinear transition stage between them. The water sensitivity is validated to be responsible for the anomalous deviation of water absorption from the classical square root of time law. It also partially contributes to the anomalous dependence of sorptivity on the temperature-dependent surface tension and viscosity of water and many organic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Tensor Network-Based Quantum Algorithm for the Nonlinear 1D Burgers' Equation

In this work, we implement a tensor network-based quantum algorithm to solve unsteady, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the compressible 1-dimensional (1D) Burgers' equation as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts to solve nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. Our framework is based on matrix product states (MPSs) and matrix product operators (MPOs). For example, the velocity field is represented by MPS, whereas the linear and nonlinear spatial differential terms of the velocity field are processed by MPOs. Our primary focus herein is to verify and validate the various tensor network components of the algorithm using solutions obtained by the classical algorithms on high performance computing (HPC) architectures. We use a classical time marching method to demonstrate the functionality of the tensor network operations to model the PDE and their robustness with the time evolution of the system. Our classical simulation results demonstrate the utility of tensor network-based operations in modeling nonlinear PDEs and highlight the necessity as well as potential advantages of using quantum simulations for these techniques.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
DOE OSTI2022

Dimensionless parameters to categorize the failure modes of ductile plate perforation

In this work, two dimensionless parameters were derived organically from impact mechanics equations to remove the need for empirical constants. These two dimensionless parameters effectively collapse the ballistic performance of ductile metal plates along two distinct loci, depending on the projectile geometry and the corresponding failure mode. Projectiles that initiate ductile hole enlargement via cavity expansion fall along a linear curve, while shear-plugging failure modes tend to fall along a cubic curve. Results were demonstrated for various projectiles of different calibers impacting aluminum and Ti-6Al-4V alloys. The main advantage of the approach in this work is to drastically reduce the number of features for input into machine learning algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dimensionless parameters to categorize the failure modes of ductile plate perforation

In this work, two dimensionless parameters were derived organically from impact mechanics equations to remove the need for empirical constants. These two dimensionless parameters effectively collapse the ballistic performance of ductile metal plates along two distinct loci, depending on the projectile geometry and the corresponding failure mode. Projectiles that initiate ductile hole enlargement via cavity expansion fall along a linear curve, while shear-plugging failure modes tend to fall along a cubic curve. Results were demonstrated for various projectiles of different calibers impacting aluminum and Ti-6Al-4V alloys, and further explored for more complex target combinations. The main advantage of the approach in this work is to drastically reduce the number of features for input into machine learning algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗
DOE OSTI2021

Exact results for nonlinear Drude weights in the spin-1/2 XXZ chain

Nonlinear Drude weight (NLDW) is a generalization of the linear Drude weight, which characterizes the nonlinear transport in quantum many-body systems. We investigate these weights for the spin-1/2 XXZ chain in the critical regime. The effects of the Dzyaloshinskii-Moriya interaction and an external magnetic field are also studied. Solving the Bethe equations numerically, we obtain these weights for very large system sizes and identify parameter regimes where the weights diverge in the thermodynamic limit. These divergences appear in all the orders studied in this Letter and can be regarded as a generic feature of the NLDWs. We study the origin of these divergences and reveal that they result from nonanalytic finite-size corrections to the ground-state energy. As a result, we compute closed-form expressions for several weights in the thermodynamic limit and find excellent agreement with the numerical results.

1-dimensional systems↗
DOE OSTI2026

Large-scale white noise and cosmology

The generation of white noise on large scales is a generic property of the dynamics of physical systems described by local nonlinear partial differential equations. Nonlinearities prevent the small-scale dynamics from being erased by smoothing. Unresolved small-scale dynamics act as an uncorrelated (white or Poissonian) noise (seemingly stochastic but actually deterministic) contribution to large-scale dynamics. This white noise exists even when the dynamics is very nearly linear. In cases where the power spectrum is sub-Poissonian on large scales, this noise will dominate on the largest scale power no matter the amplitude of the inhomogeneities. Such is the case in the standard model of cosmology, where the primordial density power spectrum is expected to have an almost Harrison-Zel’dovich, P [ k ] ∼ k , spectrum on a much broader range of scales than can be observed. Even though linear gravitational evolution dominates nonlinear corrections by a factor of ∼ 10 5 , the nonobservation of white noise on the Hubble scale precludes the extrapolation of this power law below the comoving 1 pc scale. More generally, observation or nonobservation of large-scale white noise provides a powerful probe of the Universe on very small scales in the early Universe. Gravitational radiation, phase transitions, vorticity, and running of the spectral index are all phenomena that can be probed with large-scale white noise. Large-scale white noise is a nonoptional feature of all cosmological models but one which has not heretofore been appreciated.

Barenboim, Gabriela [Valencia U., IFIC; Valencia U↗
DOE OSTI2021

ExaWind: Predictive Wind Energy Simulations

This presentation describes the ExaWind project and the team's progress in creating a suite of performance-portable codes designed for predictive simulations of wind farms on next-generation exascale-class supercomputers. Such simulations will require the resolution of scales spanning many orders of magnitude, from blade boundary layers to wind farm flow structures. In the U.S., the first exascale systems will be GPU accelerated, and different GPU manufacturers have been chosen for the different systems. At the heart of the ExaWind software is a hybrid-solver approach based on the codes Nalu-Wind and AMR-Wind, which are computational fluid dynamics solvers for the incompressible Navier-Stokes equations. Nalu-Wind is an unstructured-grid code used to resolve wind turbine geometry and blade boundary layers, whereas AMR-Wind is a structured-grid background solver for atmospheric turbulent flow and turbine wake propagation. The models are coupled with overset meshes and global linear systems are approximated through a loose-coupling algorithm. Results will include validation-quality high-fidelity simulations and strong/weak scaling results from the Summit supercomputer.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method (MPM)

Access to fresh drinking water in the future requires effective management of industrial wastewater and water purification from available resources. In this study, we present the simulation methodology and the analysis of a Reverse Osmosis (RO) membrane under various pressure conditions using the Material Point Method (MPM). In contrast to other numerical methods, MPM solves the continuum governing equations on material points in a Lagrangian framework. The method does not require a grid connecting the material points hence making it suitable to simulate large deformations during membrane compaction. The time integration is carried out using the explicit Euler method, while the spatial discretization is performed using linear or cubic-spline shape functions. A series of planar images containing detailed pore structures obtained from X-ray tomography experiments is converted to a three-dimensional collection of material points to simulate the membrane. Compressive loads are applied to the top layer of the membrane to simulate the experiments. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compaction in detail.

ENGINEERING,MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗
DOE OSTI2020

Hot electrons between cold walls

We consider electron cooling in a collisionless plasma slab between two cold and freely emitting walls. Numerical calculations suggest a counterintuitive behavior of this system: the cooling rate slows down and eventually stops, leaving the system with a significant fraction of its initial thermal energy. Analytical treatment within the Vlasov–Poisson model reveals a set of steady-states with a two-component distribution of electrons: the primary electrons trapped within the potential wells and the secondary electrons forming the counterstreaming beams. We show that such steady-states are linearly stable with respect to one-dimensional perturbations. Establishment of a particular steady-state depends on initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Constraining the release of Sn to the ambient melting point following shock loading using time-resolved x-ray diffraction

The formation of liquid following release from a shocked state governs the transition from spall to cavitation and the formation of ejecta in metals. In order to build physics-based models of these processes, it is necessary to critically evaluate the relative importance of kinetics and entropy generation during the release along with the accuracy of multiphase equations of state. Tin (Sn) has served as a testbed for a variety of experiments examining strength and ejecta due to its accessible melt boundary and solid–solid phase transitions. This work presents experiments examining the phase evolution of high purity Sn following the shock and release to ambient stress near the melting point. Sn is found to release to states between its ambient solidus and liquidus from approximately 19 to 33 GPa under uniaxial loading, with the two-phase region being characterized by a reduction in the intensity of the (220), (211) β -Sn doublet. Jetting experiments performed at 27–28 GPa exhibit comparable diffraction patterns with what is observed following the uniaxial release. The solid fractions of β-Sn in the ambient mixed phase region are found to decrease linearly with increasing shock stress as increasing liquid Sn is formed. Here, the results provide much needed information for interpreting measurements of dynamic strength at a high strain rate and experiments examining cavitation and shallow bubble collapse in Sn.

, spall↗