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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

On the peridynamic effective force state and multiphase constitutive correspondence principle

This article concerns modeling unsaturated deformable porous media as an equivalent single-phase and single-force state peridynamic material through the effective force state. The balance equations of linear momentum and mass of unsaturated porous media are presented by defining relevant peridynamic states. The energy balance of unsaturated porous media is utilized to derive the effective force state for the solid skeleton that is an energy conjugate to the nonlocal deformation state of the solid, and the suction force state. Through an energy equivalence, a multiphase constitutive correspondence principle is built between classical unsaturated poromechanics and peridynamic unsaturated poromechanics. The multiphase correspondence principle provides a means to incorporate advanced constitutive models in classical unsaturated porous theory directly into unsaturated peridynamic poromechanics. Finally, numerical simulations of localized failure in unsaturated porous media under different matric suctions are presented to demonstrate the feasibility of modeling the mechanical behavior of such three-phase materials as an equivalent single-phase peridynamic material through the effective force state concept.

42 ENGINEERING↗
DOE OSTI2021

A general solution for accelerating screw dislocations in arbitrary slip systems with reflection symmetry

Solutions to the differential equations of linear elasticity in the continuum limit in arbitrary crystal symmetry are known only for steady-state dislocations of arbitrary character, i.e. line defects moving at constant velocity. Troubled by singularities at certain ‘critical’ velocities (typically close to certain sound speeds), these dislocation fields are thought to be too idealized, and divergences are usually attributed to neglecting the finite size of the core and to the restriction to constant velocity. In the isotropic limit, accelerating pure screw and edge dislocations were studied some time ago. A generalization to anisotropic crystals has been attempted for pure screw and edge dislocations only for some special cases. This work aims to fill the gap of deriving a general anisotropic solution for pure screw dislocations applicable to slip systems featuring a reflection symmetry, a prerequisite to studying pure screw dislocations without mixing with edge dislocations. Finally, further generalizations to arbitrary mixed dislocations as well as regularizations of the dislocation core are beyond the scope of this paper and are left for future work.

42 ENGINEERING↗
DOE OSTI2020

Pore-scale simulation of gas flow in microscopic permeable media with complex geometries

Lattice Boltzmann method (LBM) is an efficient tool to perform direct numerical simulation of gas flow in micro-size pores. Despite great advances in LBM and its use in describing gas flow in micro-size channels in slip and transitional flow regimes, exploring formulations that are able to capture flow behavior in domains with complex boundary geometries remains a challenging task. As a result, the impact of complexities in pore structures on gas flow may not be fully explored using LBM. In this study, we propose to use the combined bounce-back and Maxwellian diffusive reflection scheme to capture the slip velocity with a new set of slip coefficients. Optimal values of slip coefficients are determined using a gradient based method, where the results from the linearized Boltzmann equation is used as the reference solution. The proposed formulation is further validated against predictions from molecular dynamics simulation in the presence of complex geometries by introducing obstacles of different shapes in straight channels. LBM with the proposed boundary treatment is then used to investigate gas flow in a synthetic microscopic permeable medium. Here, results indicate that gas exhibits different flow configurations as Knudsen number varies, and apparent gas permeability appears to have up to a second-order dependency on the reciprocal mean pressure in slip and transitional flow regimes.

03 NATURAL GAS↗
DOE OSTI2024

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An Open‐Source, Physics‐Based, Tropical Cyclone Downscaling Model With Intensity‐Dependent Steering

Abstract An open‐source, physics‐based tropical cyclone (TC) downscaling model is developed, in order to generate a large climatology of TCs. The model is composed of three primary components: (a) a random seeding process that determines genesis, (b) an intensity‐dependent beta‐advection model that determines the track, and (c) a non‐linear differential equation set that determines the intensification rate. The model is entirely forced by the large‐scale environment. Downscaling ERA5 reanalysis data shows that the model is generally able to reproduce observed TC climatology, such as the global seasonal cycle, genesis locations, track density, and lifetime maximum intensity distributions. Inter‐annual variability in TC count and power‐dissipation is also well captured, on both basin‐wide and global scales. Regional TC hazard estimated by this model is also analyzed using return period maps and curves. In particular, the model is able to reasonably capture the observed return period curves of landfall intensity in various sub‐basins around the globe. The incorporation of an intensity‐dependent steering flow is shown to lead to regionally dependent changes in power dissipation and return periods. Advantages and disadvantages of this model, compared to other downscaling models, are also discussed.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Pseudospin-triplet pairing in iron-chalcogenide superconductors

Abstract Understanding the pairing symmetry is a crucial theoretical aspect in the study of unconventional superconductivity for interpreting experimental results. Here we study superconductivity of electron systems with both spin and pseudospin-1/2 degrees of freedom. By solving linearized gap equations, we derive a weak coupling criterion for the even-parity spin-singlet pseudospin-triplet pairing. It can generally mix with the on-site s -wave pairing since both of them belong to the same symmetry representation ( A 1 g ) and their mixture could naturally give rise to anisotropic intra-band pairing gap functions with or without nodes. This may directly explain why some of the iron-chalcogenide superconductors are fully gapped (e.g. FeSe thin film) and some have nodes (e.g. LaFePO and LiFeP). We also find that the anisotropy of gap functions can be enhanced when the principal rotation symmetry is spontaneously broken in the normal state such as nematicity, and the energetic stabilization of pseudospin-triplet pairings indicates the coexistence of nematicity and superconductivity. This could be potentially applied to bulk FeSe, where gap anisotropy has been experimentally observed.

Physics↗
DOE OSTI2023

Molecular hydrodynamic theory of the velocity autocorrelation function

The velocity autocorrelation function (VACF) encapsulates extensive information about a fluid’s molecular-structural and hydrodynamic properties. We address the following fundamental question: How well can a purely hydrodynamic description recover the molecular features of a fluid as exhibited by the VACF? To this end, we formulate a bona fide hydrodynamic theory of the tagged-particle VACF for simple fluids. Our approach is distinguished from previous efforts in two key ways: collective hydrodynamic modes and tagged-particle self-motion are modeled by linear hydrodynamic equations; the fluid’s spatial velocity power spectrum is identified as a necessary initial condition for the momentum current correlation. This formulation leads to a natural physical interpretation of the VACF as a superposition of products of quasinormal hydrodynamic modes weighted commensurately with the spatial velocity power spectrum, the latter of which appears to physically bridge continuum hydrodynamical behavior and discrete-particle kinetics. The methodology yields VACF calculations quantitatively on par with existing approaches for liquid noble gases and alkali metals. Furthermore, we obtain a new, hydrodynamic form of the self-intermediate scattering function whose description has been extended to low densities where the Schmidt number is of order unity; various calculations are performed for gaseous and supercritical argon to support the general validity of the theory. Excellent quantitative agreement is obtained with recent MD calculations for a dense supercritical Lennard-Jones fluid.

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

Closure theory for high-collisionality multi-ion plasmas

A general formalism is developed to construct and solve a system of linearized moment equations for parallel and perpendicular closures in high-collisionality plasmas. It is applicable for multiple ion species with arbitrary masses, temperatures, charges, and densities. The convergence of closure coefficients is evaluated by increasing the number of moments from 2 to 32 for scalar, vector, and rank-2 tensor moments. As an example, the complete set of closure coefficients for a deuterium-carbon plasma over the entire Hall parameter range is presented. Furthermore, the closure coefficients at various temperature ratios show that the one-temperature closure coefficients can differ significantly from the two-temperature coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An energy-dependent electro-thermal response model of CUORE cryogenic calorimeter

The Cryogenic Underground Observatory for Rare Events (CUORE) is the most sensitive experiment searching for neutrinoless double-beta decay (0νββ) in 130 Te. CUORE uses a cryogenic array of 988 TeO 2 calorimeters operated at ~10 mK with a total mass of 741 kg. To further increase the sensitivity, the detector response must be well understood. Here, we present a non-linear thermal model for the CUORE experiment on a detector-by-detector basis. We have examined both equilibrium and dynamic electro-thermal models of detectors by numerically fitting non-linear differential equations to the detector data of a subset of CUORE channels which are well characterized and representative of all channels. Further, we demonstrate that the hot-electron effect and electric-field dependence of resistance in NTD-Ge thermistors alone are inadequate to describe our detectors' energy-dependent pulse shapes. We introduce an empirical second-order correction factor in the exponential temperature dependence of the thermistor, which produces excellent agreement with energy-dependent pulse shape data up to 6 MeV. We also present a noise analysis using the fitted thermal parameters and show that the intrinsic thermal noise is negligible compared to the observed noise for our detectors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Randomized Adiabatic Quantum Linear Solver Algorithm with Optimal Complexity Scaling and Detailed Running Costs

Solving linear systems of equations is a fundamental problem with a wide variety of applications across many fields of science, and there is increasing effort to develop quantum linear solver algorithms. Subaşı et al. [Phys. Rev. Lett. 122, 060504 (2019)] proposed a randomized algorithm inspired by adiabatic quantum computing, based on a sequence of random Hamiltonian simulation steps, with suboptimal scaling in the condition number 𝜅 of the linear system and the target error 𝜖. Here we go beyond these results in several ways. Firstly, using filtering [Lin and Tong, Quantum 4, 361 (2020)] and Poissonization techniques [Cunningham and Roland, ArXiv:2406.03972 (2024)], the algorithm complexity is improved to the optimal scaling 𝑂⁡(𝜅⁢log (1/𝜖))—an exponential improvement in 𝜖, and a shaving of a log 𝜅 scaling factor in 𝜅. Secondly, the algorithm is further modified to achieve constant factor improvements, which are vital as we progress towards hardware implementations on fault-tolerant devices. We introduce a cheaper randomized walk operator method replacing Hamiltonian simulation—which also removes the need for potentially challenging classical precomputations; randomized routines are sampled over optimized random variables; circuit constructions are improved. We obtain a closed formula rigorously upper bounding the expected number of times one needs to apply a block-encoding of the linear system matrix to output a quantum state encoding the solution to the linear system. The upper bound is 837⁢𝜅 at 𝜖 = 10 −10 for Hermitian matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

3D theory of microscopic instabilities driven by space-charge forces

Microscopic, or short-wavelength, instabilities are known for a drastic reduction of the beam quality and strong amplification of the noise in a beam. Space charge and coherent synchrotron radiation are known to be the leading causes of such instabilities. In this paper, we present a rigorous 3D theory of such instabilities driven by the space-charge forces. We define the condition when our theory is applicable to an arbitrary accelerator system with 3D coupling. Finally, we derive a linear integral equation describing such instability and identify conditions when it can be reduced to an ordinary second-order differential equation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Accurate relativistic observables from postprocessing light cone catalogs

We introduce and study a new scheme to construct relativistic observables from postprocessing light cone data. This construction is based on a novel approach, LC-metric, which takes general light cone or snapshot output generated by arbitrary N-body simulations or emulations and solves the linearized Einstein equations to determine the spacetime metric on the light cone. We find that this scheme is able to determine the metric to high precision, and subsequently generate accurate mock cosmological observations sensitive to effects such as post-Born lensing and nonlinear integrated Sachs-Wolfe contributions. By comparing to conventional methods in quantifying those general relativistic effects, we show that this scheme is able to accurately construct the lensing convergence signal. We also find that the accuracy of this method in quantifying the integrated Sachs-Wolfe effects in the highly nonlinear regime outperforms conventional methods by an order of magnitude. Furthermore, this scheme opens a new path for exploring and modeling higher-order and nonlinear general relativistic contributions to cosmological observables, including mock observations of gravitational lensing and the moving lens and Rees-Sciama effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Emulating 2D Materials with Magnons

Spin waves (magnons) in two-dimensional (2D) materials have received increasing interest due to their unique states and potential for tunability. However, many interesting features of these systems, including Dirac points and topological states, occur at high frequencies, where experimental probes are limited. Here, we study a crystal formed by patterning a hexagonal array of holes in a perpendicularly magnetized thin film. Through simulation, we find that the magnonic band structure imitates that of graphene, but additionally has some kagomelike character and includes a few flat bands. Surprisingly, its nature can be understood using a nine-band tight-binding Hamiltonian. This clear analogy to 2D materials enables band-gap engineering in 2D, topological magnons along 1D phase boundaries, and spectrally isolated modes at 0D point defects. Interestingly, the 1D phase boundaries allow access to the valley degree of freedom through a magnonic analog of the quantum valley Hall insulator. These approaches can be extended to other magnonic systems, but are potentially more general due to the simplicity of the model, which resembles existing results from electron, phonon, photon, and cold-atom systems. This finding brings the physics of spin waves in 2D materials to more experimentally accessible scales, augments it, and outlines a few principles for controlling magnonic states.

Ferrimagnets↗
DOE OSTI2022

Phase transformation and growth mechanism of RF sputtered ferroelectric lead scandium tantalate (PbSc 0.5 Ta 0.5 O 3 ) films

Lead scandium tantalate (PbSc 0.5 Ta 0.5 O 3 , PST), an order/disorder ferroelectric, is a potential candidate for electrocaloric cooling and pyroelectric infrared (IR) detector. In this work, we report the phase transformation kinetics from two series of samples containing pure amorphous and mixture of amorphous and pyrochlore to desired perovskite phase using postdeposition rapid thermal processing (RTP) as well as growth mechanism of RF sputtered PST thin films using excess lead target on platinized silicon (Pt/Ti/SiO 2 /Si) substrates. We find that small changes in the temperature ramp have a large effect on the degree of perovskite conversion (ferroelectric phase), orientation (crystallographic texture), and long-range order parameter (< S 111 >). Through isothermal annealing, we obtained optimal perovskite phase at ≥700°C temperature. The phase transformation is characterized by spontaneous formation of center-type in-plane radial rosette-like structures revealed by scanning electron microscopy. The PST perovskite crystallites were found to coexist with pyrochlore in RTP annealed films. The volume fractions for perovskite and pyrochlore phase were obtained from the analysis of “rosettes” and respective X-ray diffraction intensities which helped to determine various parameters associated with phase kinetics (n, k, and activation energy, Ea) and accompanying growth. The effective activation energies of perovskite transition and growth were found to be 332 ± 11 kJ/mol (345 ± 11 kJ/mol) and 114 ± 10 kJ/mol (122 ± 10 kJ/mol), respectively, for pure amorphous only (and mixed amorphous and pyrochlore) phase following nucleation-growth controlled Avrami's equation. A linear growth rate (n~1) for the perovskite phase indicates predominant interface-controlled process and diffusion-limited phenomena thus inhibiting rosette size owing to reactant depletion and soft impingement at the grain boundary. However, the growth behavior is isotropic in two-dimension parallel to the plane of the substrates for both sample series. Furthermore, lead loss was severe for in-situ growth and RTP combined with conventional furnace annealing than those of RTP only films, which were closer to stoichiometric albeit with excess lead and marginal oxygen vacancies (V o ).

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Analytic solutions of the DGLAP evolution and theoretical uncertainties

The energy dependence for the singlet sector of Parton Distributions Functions (PDFs) is described by an entangled pair of ordinary linear differential equations. Although there are no exact analytic solutions, it is possible to provide approximated results depending on the assumptions and the methodology adopted. These results differ in their sub-leading, neglected terms and ultimately they are associated with different treatments of the theoretical uncertainties. In this work, a novel analytic approach in Mellin space is presented and a new methodology for obtaining closed and exponentiated analytic solutions is devised. Different results for the DGLAP evolution at Next-Leading-Order are compared, discussing advantages and disadvantages for each solution. The generalizations to higher orders are addressed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A Simple, Scalable Large Deformation Solid Mechanics Implementation in the MOOSE Framework

This article describes a large deformation solid mechanics solver implemented as part of the freely available and open source MOOSE finite element simulation framework. The article documents the choices made in developing the solid mechanics framework and describes novel formulations for the gradient operator and constitutive modeling framework made to simplify implementations of different coordinate systems, stabilized gradient operators, and different constitutive model inputs and outputs. In the process, the article describes a new formulation that casts objective integration of the Cauchy stress as a linear transformation of the small stress rate. Finally, the article presents key implementation details and examines the parallel efficiency of the solid mechanics solver implemented in MOOSE. The implementation retains a good weak scaling efficiency beyond 1,000 parallel processes. The article includes a discussion of the factors limiting the parallel efficiency of implicit, large deformation solid mechanics codes on current high-performance computers, with the main current limitation being the scalability of the algebraic multigrid methods used to solve the linearized equilibrium equations.

Applied computing → Computer-aided design↗