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

Closed-Form Solutions for the Equations of Motion of the Heavy Symmetrical Top with One Point Fixed

The Equations of Motion (EOM) for the Heavy Symmetrical Top with One Point Fixed are highly non-linear. The literature describes the numerical methods that are used to resolve this Classical system including modern tools i.e. the Runge-Kutta Fourth Order method. It is more difficult to derivate closed-form solutions for the EOM and as mentioned in the literature it is not always possible to find the close-form solution for all the EOM. Fortunately, there are a few examples available that will serve as a guide to move further on in this topic. It is the purpose of this paper to find a methodology that will produce the solutions for a given subset of EOM that fulfill certain requisites. The report is organized as follows: it starts with a very short summary of the literature available on this topic and quickly follows into the derivation of the EOM using the Euler-Lagrange method. The Routhian will be used to reduce the size of the expression. It continues with the formulation of the Classical cubic function (f(u)) through a novel process. The roots of f(u) are of the outmost importance to be able to find the EOM closed-form solution, and once the final roots are selected the general method that will produce the closed-form solutions is presented. Two sets of examples are included to show the validity of the process and comparisons of the results from the closed-form solutions vs. the numerical results for these examples are shown.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Numerical solution of large scale Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size N b . We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most $\mathcal{O}$(N b 2 ) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with N b up to 2.88 × 10 6 , and the wall clock time is less than 100 s using 17 280 CPU cores. Finally, this enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Central vortex steady states and dynamics of Bose–Einstein condensates interacting with a microwave field

Here we study central vortex steady states and dynamics of a two-dimensional (2D), two-component, Gross–Pitaevskii equation (CGPE) system for two pseudo-spinor Bose Einstein condensates (BECs) interacting with an electromagnetic field (microwave) analytically and numerically. For the central vortex steady state at any given winding number S, we prove its existence in a reduced, single component detuning limit when contact-interaction strength β < β b , and nonexistence when β > β b , respectively, where β b is a threshold value, whose value is given in Theorem 3.1 in the paper. We extend the existence and nonexistence result to the general two pseudo-spinor case and prove that a central vortex steady state exists for any given S if β ≤ β b while it does not exist when β > 2β b . We then derive dynamical equations for some observables (expectations of the matter wave function) such as the center-of-mass, position of dispersion, the linear and angular momentum of the two pseudo-spinor CGPEs. Finally, numerical computations are brought in to validate and extend the existence of vortex steady state results to β b ≤ β < 2β b for the two pseudo-spinor case and to explore transient dynamics of the observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The equation of state and shock-driven decomposition of polymethylmethacrylate (PMMA)

We present new equations of state (EOS) for polymethylmethacrylate and its shock-driven decomposition products, generated in both tabular form and as linear U S – u p fits to all previously available shock data. Different approaches are taken for materials shocked above and below the apparent threshold for chemical reactions at pressures P ≈ 30 GPa. Both EOSs are calibrated to and compared with a broad array of experimental results, as well as new ab initio molecular dynamics simulations presented here. Comparisons are good in most cases, and we focus our particular attention on those observables that test the distinction between reactants and products.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2022

Polyethylene in Dead-End Silica Nanopores: Forces and Mobility from Non-Equilibrium Statistical Mechanics and Exchange Spectroscopy Nuclear Magnetic Resonance

Billions of tons of plastic have been produced, and only a small fraction of this has been recycled. Tennakoon et al. [Nature Catalysis 3, 893 (2020)] developed a catalyst that repeatedly cleaves C 10 –C 30 hydrocarbons from the end of a polyethylene chain. The reaction occurs at a Pt nanoparticle at the base of a cylindrical silica mesopore with a diameter of 2 nm and a length of 110 nm. Portions of the polymer situated inside the pore can be differentiated from those outside using 13 C nuclear magnetic resonance (NMR), allowing the dynamics and extent of polymer threading to be monitored using two-dimensional (2D) exchange spectroscopy NMR. Here we construct a Fokker–Planck equation for the polymer dynamics by assuming a reptation diffusivity and a graduated adsorption free energy that depends linearly on the depth of polymer penetration in the pore. The solutions allow us to predict the intensities of the 2D NMR resonances as a function of time. We use the solutions to extract a polymer diffusivity at each temperature and estimate the per-segment desorption free energy, enthalpy, and entropy. Random and systematic errors are examined to test key assumptions in the theory and interpretation of the experiments.

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

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Space-time dependent thermal conductivity in nonlocal thermal transport

Nonlocal thermal transport is generally described by the Peierls-Boltzmann transport equation (PBE). However, solving the PBE for a general space-time dependent problem remains a challenging task due to the high dimensionality of the integro-differential equation. In this work, we present a direct solution to the space-time dependent PBE with a linearized collision matrix using an eigendecomposition method. We show that there exists a generalized Fourier-type relation that links heat flux to the local temperature, and this constitutive relation defines a thermal conductivity that depends on both time and space. Combining this approach with ab initio calculations of phonon properties, we demonstrate that the space-time dependent thermal conductivity gives rise to an oscillatory response in temperature in a transient grating geometry in high thermal conductivity materials. The present solution method allows us to extend the reach of our computational capability for heat conduction to space-time dependent nondiffusive transport regimes. Here, this capability will not only enable a more accurate interpretation of thermal measurements that observe nonlocal thermal transport, but also enhance our physical understanding of nonlocal thermal transport in high thermal conductivity materials that are promising candidates for nanoscale thermal management applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Critical and near-critical relaxation of holographic superfluids

We investigate the relaxation of holographic superfluids after quenches, when the end state is either tuned to be exactly at the critical point, or very close to it. By solving the bulk equations of motion numerically, we demonstrate that in the former case the system exhibits a power law falloff, as well as an emergent discrete scale invariance. The latter case is in the regime dominated by critical slowing down, and we show that there is an intermediate time range before the onset of late-time exponential falloff, where the system behaves similarly to the critical point with its power law falloff. We further postulate a phenomenological Gross-Pitaevskii-like equation (corresponding to model F of Hohenberg and Halperin) that is able to make quantitative predictions for the behavior of the holographic superfluid after near-critical quenches into the superfluid and normal phase. Intriguingly, all parameters of our phenomenological equation, which describes the nonlinear time evolution, may be fixed with information from the static equilibrium solutions and linear response theory. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗
DOE OSTI2024

Exploiting Power Flow Manifold to Solve AC Optimal Power Flow

AC optimal power flow has proven difficult to solve with interior point methods on GPUs. This is largely due to challenging linear algebra problems that current state of the art massively parallel linear solvers struggle with. However, the advent of Riemannian optimization techniques and the fact that the power flow equations form a smooth manifold present an alternative approach. In this talk, we present the basics of Riemannian optimization techniques in which optimization is done directly on a manifold. Then we present computational results showing that Riemannian techniques are capable of producing solutions of comparable quality as interior point methods.

AC optimal power flow↗
DOE OSTI2024

Quasinormal modes of naked singularities in presence of nonlinear scalar fields

We study linear perturbations against static spherically symmetric background configurations of general relativity with a real scalar field (SF), which is minimally coupled with gravity; it is nonlinear due to the presence of the self-action potential. The background solutions have a naked singularity at the center of the configuration. The focus is on the stability of the background and fundamental frequencies of the quasinormal modes (QNM) of the axial perturbations in the Regge-Wheeler gauge. The problem is reduced to one hyperbolic master equation with an effective potential 𝑊eff, which turns out to be positive for a general non-negative SF potential; this ensures the linear stability with respect to this kind of perturbations. For numerical simulations, the SF potential was chosen in the power-law form 𝑉⁡(𝜙) ∼ 𝜙 2⁢𝑛 with 2 < 𝑛 ≤40. We extracted the fundamental frequencies of QNM for different 𝑛 and various sets of the background configuration parameters. The results show that even for a small background SF, there is a significant difference between the fundamental frequencies and ones in case of the Schwarzschild background. Furthermore, the results are also compared with the case of the Fisher-Janis-Newman-Winicour background dealing with a massless linear scalar field.

General relativity↗
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

New insights into supradense matter from dissecting scaled stellar structure equations

The strong-field gravity in general relativity (GR) realized in neutron stars (NSs) renders the equation of state (EOS) P(ε) of supradense neutron star matter to be essentially nonlinear and refines the upper bound for Φ ≡ P/ε to be much smaller than the special relativity (SR) requirement with linear EOSs, where P and ε are respectively the pressure and energy density of the system considered. Specifically, a tight bound Φ ≲ 0.374 is obtained by perturbatively anatomizing the intrinsic structures of the scaled Tolman–Oppenheimer–Volkoff (TOV) equations without using any input nuclear EOS. New insights gained from this novel analysis provide EOS-model-independent constraints on the properties (e.g., density profiles of the sound speed squared s 2 = dP/dε and trace anomaly Δ = 1/3 – Φ) of cold supradense matter in NS cores. Using the gravity-matter duality in theories describing NSs, we investigate the impact of gravity on supradense matter EOS in NSs. In particular, we show that the NS mass M NS , radius R, and compactness ξ ≡ M NS /R scale with certain combinations of its central pressure and energy density (encapsulating its central EOS). Thus, observational data on these properties of NSs can straightforwardly constrain NS central EOSs without relying on any specific nuclear EOS model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗