Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hyperbolic systems”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Superior discretizations and AMG solvers for extremely anisotropic diffusion via hyperbolic operators [Slides]

Diffusion in magnetic confinement fusion is extremely anisotropic in the direction of field lines. Rewrote diffusion system based on directional gradients, apply discretization and solver techniques developed for advection. Orders of magnitude decrease in error and solve wallclock time vs. traditional methods. The next steps include: (1) incorporate into larger MHD simulations, (2) better solvers for closed field lines or mixed regimes, and (3) possibly other extremely anisotropic equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Colossal Optical Anisotropy from Atomic-Scale Modulations

Materials with large birefringence (Δn, where n is the refractive index) are sought after for polarization control (e.g., in wave plates, polarizing beam splitters, etc.), nonlinear optics, micromanipulation, and as a platform for unconventional light–matter coupling, such as hyperbolic phonon polaritons. Layered 2D materials can feature some of the largest optical anisotropy; however, their use in most optical systems is limited because their optical axis is out of the plane of the layers and the layers are weakly attached. This work demonstrates that a bulk crystal with subtle periodic modulations in its structure—Sr 9/8 TiS 3 —is transparent and positive-uniaxial, with extraordinary index n e = 4.5 and ordinary index n o = 2.4 in the mid- to far-infrared. The excess Sr, compared to stoichiometric SrTiS 3 , results in the formation of TiS 6 trigonal-prismatic units that break the chains of face-sharing TiS 6 octahedra in SrTiS 3 into periodic blocks of five TiS 6 octahedral units. The additional electrons introduced by the excess Sr form highly oriented electron clouds, which selectively boost the extraordinary index n e and result in record birefringence (Δn > 2.1 with low loss). The connection between subtle structural modulations and large changes in refractive index suggests new categories of anisotropic materials and also tunable optical materials with large refractive-index modulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Environmental damping and vibrational coupling of confined fluids within isolated carbon nanotubes

Abstract Because of their large surface areas, nanotubes and nanowires demonstrate exquisite mechanical coupling to their surroundings, promising advanced sensors and nanomechanical devices. However, this environmental sensitivity has resulted in several ambiguous observations of vibrational coupling across various experiments. Herein, we demonstrate a temperature-dependent Radial Breathing Mode (RBM) frequency in free-standing, electron-diffraction-assigned Double-Walled Carbon Nanotubes (DWNTs) that shows an unexpected and thermally reversible frequency downshift of 10 to 15%, for systems isolated in vacuum. An analysis based on a harmonic oscillator model assigns the distinctive frequency cusp, produced over 93 scans of 3 distinct DWNTs, along with the hyperbolic trajectory, to a reversible increase in damping from graphitic ribbons on the exterior surface. Strain-dependent coupling from self-tensioned, suspended DWNTs maintains the ratio of spring-to-damping frequencies, producing a stable saturation of RBM in the low-tension limit. In contrast, when the interior of DWNTs is subjected to a water-filling process, the RBM thermal trajectory is altered to that of a Langmuir isobar and elliptical trajectories, allowing measurement of the enthalpy of confined fluid phase change. These mechanisms and quantitative theory provide new insights into the environmental coupling of nanomechanical systems and the implications for devices and nanofluidic conduits.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Core-shell metallic alloy nanopillars-in-dielectric hybrid metamaterials with magneto-plasmonic coupling

Combining plasmonic and magnetic properties, namely magneto-plasmonic coupling, inspires great research interest and the search for magneto-plasmonic nanostructure becomes considerably critical. In this work, we designed a nanopillar-in-matrix structure with core–shell alloyed nanopillars for both BaTiO 3 (BTO)-Au 0.5 Co 0.5 (AuCo) and BTO-Au 0.25 Cu 0.25 Co 0.25 Ni 0.25 (AuCuCoNi) hybrid systems, i.e., ferromagnetic alloy cores (e.g., Co or CoNi) with plasmonic shells (e.g., Au or Au/Cu). These core–shell alloy nanopillars are uniformly embedded into a dielectric BTO matrix to form a vertically aligned nanocomposite (VAN) structure. Both hybrid systems present excellent epitaxial quality and interesting multi-functionality, e.g., high magnetic anisotropy, magneto-optical coupling response, tailorable plasmonic resonance wavelength, tunable hyperbolic properties and strong optical anisotropy. These alloyed nanopillars-in-matrix designs provide enormous potential for complex hybrid material designs with multi-functionality and demonstrate strong interface enabled magneto-plasmonic coupling along with plasmonic and magnetic performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nemchinov-Dyson Solutions of the Two-Dimensional Axisymmetric Inviscid Compressible Flow Equations [Review]

We investigate the two-dimensional (2D) inviscid compressible flow equations in axisymmetric coordinates, constrained by an ideal gas equation of state (EOS). Beginning with the assumption that the 2D velocity field is space-time separable and linearly variable in each corresponding spatial coordinate, we proceed to derive an in nite family of elliptic or hyperbolic, uniformly expanding or contracting \gas cloud" solutions. Construction of specific example solutions belonging to this family is dependent on the solution of a system of nonlinear, coupled, second-order ordinary differential equations, and the prescription of an additional physical process of interest (e.g., uniform temperature or uniform entropy flow). The physical and computational implications of these solutions as pertaining to quantitative code verification or model qualification studies are discussed in some detail.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Geometry and entanglement in the scattering matrix

A formulation of nucleon–nucleon scattering is developed in which the S-matrix, rather than an effective-field theory (EFT) action, is the fundamental object. Spacetime plays no role in this description: the S-matrix is a trajectory that moves between RG fixed points in a compact theory space defined by unitarity. This theory space has a natural operator definition, and a geometric embedding of the unitarity constraints in four-dimensional Euclidean space yields a flat torus, which serves as the stage on which the S-matrix propagates. Trajectories with vanishing entanglement are special geodesics between RG fixed points on the flat torus, while entanglement is driven by an external potential. The system of equations describing S-matrix trajectories is in general complicated, however the very-low-energy S-matrix –that appears at leading-order in the EFT description– possesses a UV/IR conformal invariance which renders the system of equations integrable, and completely determines the potential. In this geometric viewpoint, inelasticity is in correspondence with the radius of a three-dimensional hyperbolic space whose two-dimensional boundary is the flat torus. This space has a singularity at vanishing radius, corresponding to maximal violation of unitarity. The trajectory on the flat torus boundary can be explicitly constructed from a bulk trajectory with a quantifiable error, providing a simple example of a holographic quantum error correcting code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Observation of Universal Quench Dynamics and Townes Soliton Formation from Modulational Instability in Two-Dimensional Bose Gases

We experimentally study universal nonequilibrium dynamics of two-dimensional atomic Bose gases quenched from repulsive to attractive interactions. Here, we observe the manifestation of modulational instability that, instead of causing collapse, fragments a large two-dimensional superfluid into multiple wave packets universally around a threshold atom number necessary for the formation of Townes solitons. We confirm that the density distributions of quench-induced solitary waves are in excellent agreement with the stationary Townes profiles. Furthermore, our density measurements in the space and time domain reveal detailed information about this dynamical process, from the hyperbolic growth of density waves, the formation of solitons, to the subsequent collision and collapse dynamics, demonstrating multiple universal behaviors in an attractive many-body system in association with the formation of a quasistationary state.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Non-intrusive nonlinear model reduction via machine learning approximations to low-dimensional operators

Abstract Although projection-based reduced-order models (ROMs) for parameterized nonlinear dynamical systems have demonstrated exciting results across a range of applications, their broad adoption has been limited by their intrusivity: implementing such a reduced-order model typically requires significant modifications to the underlying simulation code. To address this, we propose a method that enables traditionally intrusive reduced-order models to be accurately approximated in a non-intrusive manner. Specifically, the approach approximates the low-dimensional operators associated with projection-based reduced-order models (ROMs) using modern machine-learning regression techniques. The only requirement of the simulation code is the ability to export the velocity given the state and parameters; this functionality is used to train the approximated low-dimensional operators. In addition to enabling nonintrusivity, we demonstrate that the approach also leads to very low computational complexity, achieving up to $$10^3{\times }$$ 10 3 × in run time. We demonstrate the effectiveness of the proposed technique on two types of PDEs. The domain of applications include both parabolic and hyperbolic PDEs, regardless of the dimension of full-order models (FOMs).

42 ENGINEERING↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantum dynamics of cold atomic gas with SU(1,1) symmetry

Motivated by recent advances in quantum dynamics, we investigate the dynamics of the system with SU(1,1) symmetry. Instead of performing the time-ordered integral for the evolution operator of the time-dependent Hamiltonian, we show that the time evolution operator can be expressed as an SU(1,1) group element. Since the SU(1,1) group describes the “rotation” on a hyperbolic surface, the dynamics can be visualized on a Poincaré disk, a stereographic projection of the upper hyperboloid. As an example, we present the trajectory of the revival of Bose-Einstein condensation and that of the scale-invariant Fermi gas on the Poincaré disk. Further considering quantum gases in an oscillating lattice, we also study the dynamics of the system with time-dependent single-particle dispersion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows

Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.

42 ENGINEERING↗
DOE OSTI2020

Design of 3D Oxide–Metal Hybrid Metamaterial for Tailorable Light–Matter Interactions in Visible and Near-Infrared Region

Dielectric–metallic hybrid metamaterials exhibit extraordinary optical properties due to the light–matter interactions at the dielectric–metallic interfaces. The ability in precision control of the light–matter interactions in nanoscale is key to tailor the optical properties of hybrid metamaterials. In this work, a complex 3D framework of multilayered self-assembled BaTiO 3 (BTO)-Au hybrid thin films is demonstrated with such precision control of the light–matter interaction in nanoscale. Here, unique “bamboo-like” Au nanostructures are formed via the bilayer and trilayer stacking of BTO-Au hybrid layers with interlayers of SrTiO 3 , CeO 2 , or MgO. Different film strain states introduced by the three interlayers result in variable diameter and density of Au nanopillars. Both simulated and experimental optical data demonstrate the localized surface plasmon resonance change and hyperbolic dispersion wavelength shift in visible to near-infrared because of the effective tuning of the Au nanopillar aspect ratio and free electron density. The highly tunable optical properties along with the ferroelectric behavior and thermal robustness of the 3D hybrid film enable it to be a great candidate for multifunctional applications. This study demonstrates a unique 3D approach for precision optical property tuning and combined functionalities in oxide–metal metamaterial systems toward future integrated photonic and electronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Good plasmons in a bad metal

Correlated metals may exhibit unusually high resistivity that increases linearly in temperature, breaking through the Mott-Ioffe-Regel bound, above which coherent quasiparticles are destroyed. The fate of collective charge excitations, or plasmons, in these systems is a subject of debate. Several studies have suggested that plasmons are overdamped, whereas other studies have detected propagating plasmons. Here, in this work, we present direct nano-optical images of low-loss hyperbolic plasmon polaritons (HPPs) in the correlated van der Waals metal MoOCl 2 . HPPs are plasmon-photon modes that waveguide through extremely anisotropic media and are remarkably long-lived in MoOCl 2 . Photoemission data presented here reveal a highly anisotropic Fermi surface, reconstructed and made partly incoherent, likely through electronic interactions as explained by many-body theory. HPPs remain long-lived despite this, revealing previously unseen imprints of many-body effects on plasmonic collective modes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced magnetic and optical properties of Y 3 Fe 5 O 12 (YIG) films with Au nanoinclusions

Y 3 Fe 5 O 12 (YIG) thin films are well known for their ferrimagnetic insulating property and low Gilbert damping coefficient (α), allowing them to be used for various spintronic applications and as magneto-optical isolators for photonic devices. Instead of doping, incorporation of plasmonic metals as nanoinclusions could be a promising route for improved magneto-optical coupling properties. In this work, YIG–Au nanocomposites have been deposited with ferrimagnetic insulating YIG as the matrix and Au nanoinclusions which introduce plasmonic absorption, optical anisotropy, and hyperbolic properties. Films with varying Au nanoinclusion densities have been processed and annealed to compare with the as-deposited ones. The films that had low Au nanoinclusion density and were annealed presented a lower magnetic damping coefficient of 2.84 × 10 −4 than the pure YIG film (9.66 × 10 −4 ). The as-deposited film with the highest Au density shows the strongest hyperbolic properties among all samples. These results demonstrate that both magnetic damping and optical properties can be tuned through deposition conditions in YIG–Au nanocomposite thin films, allowing for a balance of both properties. This YIG–Au nanocomposite system presents promising potential in next-generation opto-spintronic devices.

Quigley, Lizabeth [Purdue Univ., West Lafayette, I↗
DOE OSTI2024

ZnO-Au Hybrid Metamaterial Thin Films with Tunable Optical Properties

ZnO-Au nanocomposite thin films have been previously reported as hybrid metamaterials with unique optical properties such as plasmonic resonance properties and hyperbolic behaviors. In this study, Au composition in the ZnO-Au nanocomposites has been effectively tuned by target composition variation and thus resulted in microstructure and optical property tuning. Specifically, all the ZnO-Au nanocomposite thin films grown through the pulsed laser deposition (PLD) method show obvious vertically aligned nanocomposite (VAN) structure with the Au nanopillars embedded in the ZnO matrix. Moreover, the average diameter of Au nanopillars increases as Au concentration increases, which also leads to the redshifts in the surface plasmon resonance (SPR) wavelength and changes in the hyperbolic behaviors of the films. As a whole, this work discusses how strain-driven tuning of optical properties and microstructure resulted through a novel Au concentration variation approach which has not been previously attempted in the ZnO-Au thin film system. These highly ordered films present great promise in the areas of sensing, waveguides, and nanophotonics to name a few.

Bhatt, Nirali A.↗