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.

SourceResearch connectionFrom yearThrough year

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

At least 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Instability of an electron-plasma shear layer in an externally imposed strain flow

The E x B shear instability of a two-dimensional (2D) filament (i.e., a thin, rectangular strip perpendicular to the magnetic field) of magnetized pure electron plasma is investigated experimentally in the presence of an externally imposed strain flow. Data are acquired using a specialized Penning–Malmberg trap in which strain flows can be applied in 2D by biasing segmented electrodes surrounding the plasma. The E x B drift dynamics are well-described by the Drift-Poisson equations, which are isomorphic to the 2D Euler equations describing ideal fluids. Thus, the experimental results correspond to the Rayleigh instability of a shear layer in a 2D ideal fluid, where the electron density is analogous to the fluid vorticity. Shear layers are prepared by stretching initially axisymmetric electron vortices using a strong, applied strain flow. The data at early times are in quantitative agreement with a linear model which extends Rayleigh’s work to account for the influence of an external strain flow. In the presence of weak strain, the system approximately maintains a phase relationship that corresponds to an instantaneous Rayleigh eigenmode. The instability develops into the nonlinear regime later in time and at smaller spatial scales as the strain rate is increased. A secondary vortex pairing instability is observed, but it is suppressed when the strain-to-vorticity ratio exceeds roughly 0.025. In this way, vorticity transport perpendicular to the filament is diminished due to the applied strain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases

MoorDyn, an open-source mooring dynamics model, is being expanded with capabilities for additional mooring system features and load cases. As floating wind turbine technology matures, mooring systems are becoming more sophisticated and more complex scenarios need to be considered in the design process. Mooring systems may have synthetic line materials, ballast/buoyancy bodies along the lines, or interconnections between platforms. Failure modes may involve multiple cascading line failures that depend on mooring system dynamics. Features recently added to MoorDyn aim to address these emerging needs. MoorDyn’s linear elasticity model has been supplemented to support user-defined stress-strain curves, which can be adjusted to represent synthetic mooring materials. Rigid six-degree-of-freedom bodies in the mooring system can now be modeled using two new model objects. “Rod” objects provide an option for rigid cylindrical bodies. They use the existing Morison equation-based hydrodynamics model and can be connected to mooring lines at either end. “Body” objects provide a generic six-degree-of-freedom rigid-body representation based on a lumped-parameter model of translational and rotational properties. Rod objects can be added to Body objects and mooring lines can be attached at any location, allowing a wide variety of submerged structures to be integrated into the mooring system. Lastly, a means of dynamically simulating mooring line failures has been implemented. These new features, currently in the C++ version of MoorDyn, are described and then demonstrated on a two-turbine shared-mooring array. A qualitative view of the results suggests the new features are functioning as expected.

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

DOE OSTI2024

Interpolating the ‘t Hooft model between Instant and Light-Front dynamics in the Coulomb Gauge

The 1+1D model of quantum chromodynamics (QCD) in the infinite number of colors, or 't Hooft model, can be interpolated between the instant form dynamics (IFD) and the light-front dynamics (LFD) using an interpolation parameter 0 (IFD) ≤ δ ≤ π/4 (LFD). This was realized in the interpolating axial gauge which links the axial gauge (A 1 = 0) in IFD and the light-front gauge (A + = 0) [1]. In this presentation, we discuss the corresponding realization in the interpolating Coulomb gauge which links the temporal gauge (A 0 = 0) in IFD and the light-front gauge (A + = 0) and its benefit of resolving the issue associated with the absence of the conjugate field to the gauge field A 0 in the axial gauge. In both gauges, all degrees of freedom are physical making these gauge choices ideal for finding the bound-state equations and for renormalizability. Although the gauge independence of the physical observables such as the meson mass spectra following Regge trajectories may be guaranteed due to the gauge symmetry of QCD, the realization and interpretation of the identical physical results may depend on the gauge choices. Here, we discuss such difference in the realization of the confinement phenomena ala linear potential in the two different gauges, Coulomb vs. Axial, and highlight the gauge independent physical results expected. We also comment on the utility of the interpolation which leads to an alternative quasi-PDF that can be implemented in the lattice QCD without suffering from the large momentum boost.

Duggin, Hunter↗

DOE OSTI2020

MoorDyn V2: New Capabilities in Mooring System Components and Load Cases: Preprint

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

DOE OSTI2021

Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities

We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

The Cauchy Problem for Boltzmann Bi-linear Systems: The Mixing of Monatomic and Polyatomic Gases

Abstract From a unified vision of vector valued solutions in weighted Banach spaces, this paper establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unified approach for vector valued solutions in weighted Banach spaces. Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities, angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for p -binomial forms producing sharper estimates for the k -moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only a $$2^+$$ 2 + moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.

Physics↗

DOE OSTI2021

Differentiating densities on smooth manifolds

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Uncertainty quantification for Joule heating processes in fibrous pore-resolved media

Joule heating (JH) is an energy-efficient and sustainable technique for heating materials. Its application for industrial heating, particularly, has been gaining attention due to its potential for increasing the yield of various chemical products. The process involves the use of heating elements (materials that are highly conductive electrically and thermally) to heat up other materials or substances. These conductors, however, can exhbit varying degrees of uncertainty due to non-linearities in their temperature-dependent properties, which could result in variable material behavior. In this work, we carry out uncertainty quantification (UQ) at the pore scale to describe the uncertainty of such materials. In so doing, we applied the non-intrusive polynomial chaos expansion (PCE) technique to quantify the uncertainty within the system. The steady state Joule heating equation was solved numerically at the pore scale mimicking conditions within a heating chamber for propane dehydrogenation, and various electro-thermal profiles were obtained. We also examined the effect of the number of sampling points (20 – 100) and order of the PCE coefficients (2 – 5) on the accuracy of the temperature evaluations. The results were then benchmarked with the standard Monte Carlo (MC) method. The average temperature of the 4th-order global PCE showed good agreement with the MC results (which were positively skewed). Orders greater than 4 gave an underestimation of the temperatures while predictions for the peak temperature improved as the number of sampling points increased.

Fagbemi, Samuel [ORNL] (ORCID:0000000236995025)↗

DOE OSTI2022

Parallel-in-Time Integration for Nonlinear Hyperbolic Problems (Final Report)

The work for the subcontract was situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. Over the past two decades, much progress has been made in this area; however, this progress has largely been limited to diffusion-dominated PDEs, with some recent success in scalar linear hyperbolic PDEs. Given the ubiquity of numerical simulations of hyperbolic PDEs throughout the sciences, in particular, nonlinear hyperbolic systems, there is a strong need to develop efficient parallel-in-time techniques for hyperbolic problems beyond simple scalar and linear cases, which is the main topic of this subcontract. The main focus of the work was to further develop and perfect coarse-grid operators for the Multigrid Reduction-inTime (MGRIT) method applied to hyperbolic PDEs that were recently proposed in PhD thesis, based on a modified semi-Lagrangian approach.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Control simulations of many-body quantum systems by a synergism of discrete real-time learning and optimal control theory

We present a self-consistent algorithm for optimal control simulations of many-body quantum systems. The algorithm features a two-step synergism that combines discrete real-time machine learning (DRTL) with Quantum Optimal Control Theory (QOCT) using the time-dependent Schrödinger equation. Specifically, in step (1), DRTL is employed to identify a compact working space (i.e., the important portion of the Hilbert space) for the time evolution of the many-body quantum system in the presence of a control field (i.e., the initial or previously updated field), and in step (2), QOCT utilizes the DRTL-determined working space to find a newly updated control field for a chosen objective. Steps 1 and 2 are iterated until a self-consistent control objective value is reached such that the resulting optimal control field yields the same targeted objective value when the corresponding working space is systematically enlarged. Furthermore, to demonstrate this two-step self-consistent DRTL-QOCT synergistic algorithm, we perform optimal control simulations of strongly interacting 1D as well as 2D Heisenberg spin systems. In both scenarios, only a single spin (at the left end site for 1D and the upper left corner site for 2D) is driven by the time-dependent control fields to create an excitation at the opposite site as the target. It is found that, starting from all spin-down zero excitation states, the synergistic method is able to identify working spaces and convergence of the desired controlled dynamics with just a few iterations of the overall algorithm. In the cases studied, the dimensionality of the working space scales only quasi-linearly with the number of spins.

Artificial neural networks↗

DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗

DOE OSTI2026

The Q 10 of in situ microbial soil respiration varies with mean annual temperature, precipitation, pH, and plant cover: a meta-analysis and spatial prediction of Q 10

The temperature sensitivity of soil microbial respiration, commonly quantified using the Q 10 coefficient, is a key parameter in carbon cycle models. Uncovering how environmental factors affect in situ Q 10 values can therefore provide critical insight into potential shifts in global carbon stocks under climate change. We collected data from previously published field experiments that measured soil microbial respiration across a range of temperatures. We hypothesized that the Q 10 coefficient of in situ soil microbial respiration would vary based on environmental factors including mean annual temperature (MAT), mean annual precipitation (MAP), plant cover type, pH, soil C:N, and latitude. Linear regression revealed that Q 10 correlates negatively with MAT and MAP and positively with pH and absolute latitude. Additionally, average Q 10 varied significantly across different plant cover types; it was highest in mountain grasslands and lowest in tropical moist forests. Variation in microbial Q 10 across environmental factors may arise from underlying mechanisms such as enzyme kinetics, substrate availability and complexity, and microbial adaptation. To capture patterns in Q 10 more comprehensively, we developed a multiple linear regression model of Q 10 based on the most individually significant environmental drivers and applied it to public datasets to generate a global map of predicted Q 10 . Q 10 was higher in high-latitude and high-altitude regions, where large permafrost carbon stores are vulnerable to thawing and decomposition. We also compared fits between the Q 10 equation and a model produced from macromolecular rate theory (MMRT). We found that the MMRT model had the superior fit and may be better suited to model temperature sensitivity of complex biological reactions. Overall, our results emphasize that relationships between microbial Q 10 and environmental variables should be accounted for in climate models. Incorporating these variations in the Q 10 parameter, rather than using a fixed value, will help predict whether CO 2 emissions will be buffered or exacerbated by soil microbial respiration under climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

An Asymptotic Preserving Discontinuous Galerkin Method for a Linear Boltzmann Semiconductor Model

A key property of the linear Boltzmann semiconductor model is that as the collision frequency tends to infinity, the phase space density $f$ = $f$ ($x, v, t$) converges to an isotropic function $M (v)$$ρ$$(x, t)$, called the drift-diffusion limit, where $M$ is a Maxwellian and the physical density $ρ$ satisfies a second-order parabolic PDE known as the drift-diffusion equation. Numerical approximations that mirror this property are said to be asymptotic preserving. In this paper we build a discontinuous Galerkin method to the semiconductor model, and we show this scheme is both uniformly stable in $ε$, where 1/$ε$ is the scale of the collision frequency, and asymptotic preserving. Here in particular, we discuss what properties the discrete Maxwellian must satisfy in order for the schemes to converge in $ε$ to an accurate $h$-approximation of the drift-diffusion limit. Discrete versions of the drift-diffusion equation and error estimates in several norms with respect to $ε$ and the spacial resolution are also included.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Hugoniot, sound speed, and phase transitions of single-crystal sapphire for pressures 0.2–2.1 TPa

Sapphire (Al 2 O 3 ) is a major constituent of the Earth's mantle and has significant contributions to the field of high-pressure physics. Constraining its Hugoniot over a wide pressure range and identifying the location of shock-driven phase transitions allows for development of a multiphase equation of state and enables its use as an impedance-matching standard in shock physics experiments. In this paper we present measurements of the principal Hugoniot and sound velocity from direct impact experiments using magnetically launched flyers on the Z machine at Sandia National Laboratories. The Hugoniot was constrained for pressures from 0.2–2.1 TPa and a four-segment piecewise linear shock-velocity–particle-velocity fit was determined. First-principles molecular dynamics simulations were conducted and agree well with the experimental Hugoniot. Sound-speed measurements identified the onset of melt between 450 and 530 GPa, and the Hugoniot fit refined the onset to 525 ± 13 GPa. A phase diagram which incorporates literature diamond-anvil cell data and melting measurements is presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗