Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high-order”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2023

Towards critical and supercritical electromagnetic fields

Abstract The availability of ever stronger, laser-generated electromagnetic fields underpins continuing progress in the study and application of nonlinear phenomena in basic physical systems, ranging from molecules and atoms to relativistic plasmas and quantum electrodynamics. This raises the question: how far will we be able to go with future lasers? One exciting prospect is the attainment of field strengths approaching the Schwinger critical field ${E}_{\mathrm{cr}}$ in the laboratory frame, such that the field invariant ${E}^2-{c}^2{B}^2>{E}_{\mathrm{cr}}^2$ is reached. The feasibility of doing so has been questioned, on the basis that cascade generation of dense electron–positron plasma would inevitably lead to absorption or screening of the incident light. Here we discuss the potential for future lasers to overcome such obstacles, by combining the concept of multiple colliding laser pulses with that of frequency upshifting via a tailored laser–plasma interaction. This compresses the electromagnetic field energy into a region of nanometre size and attosecond duration, which increases the field magnitude at fixed power but also suppresses pair cascades. Our results indicate that laser facilities with peak power of tens of PW could be capable of reaching ${E}_{\mathrm{cr}}$ . Such a scenario opens up prospects for the experimental investigation of phenomena previously considered to occur only in the most extreme environments in the universe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Orbital-resolved calculations of two-center interferences in linear triatomic molecules

In this report we perform ab initio calculations of high-harmonic spectroscopy (HHS) of two-center interference phenomena in oriented carbon-dichalcogen molecules, using time-dependent density functional theory (TDDFT). We show that by resolving the harmonic response into contributions from individual Kohn-Sham orbitals, we can extract target-specific characteristics for both the spectral amplitude and phase. We also discuss that this extraction is predicated on a careful analysis and normalization of the harmonic spectrum. Finally, we present field-free scattering calculations that mimic the recollision step in high-order-harmonic generation and allow the calculation of recombination dipole matrix elements without explicitly calculating the scattering states of a molecule. We show that the orbital-resolved TDDFT HHS results and results based on our field-free scattering calculations are in very good agreement with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Wavelength scaling and multicolor operation of a plasma-driven attosecond x-ray source via harmonic generation

The generation of high-power coherent soft x-ray pulses of sub-100 as duration and 10 nm wavelength using beams from a GeV energy plasma wakefield accelerator has been recently investigated in Emma []. As a future upgrade to this concept, this contribution investigates scaling to shorter x-ray wavelengths by cascading undulators tuned to higher harmonics of the fundamental. We present two simulation studies for plasma-driven attosecond harmonic generation schemes with final photon wavelengths of 2 and 0.40 nm. We demonstrate in these schemes that using undulators with retuned fundamental frequencies can produce GW-scale pulses of sub-nm radiation with tens of attosecond-scale pulse lengths, an order of magnitude shorter than current state-of-the-art attosecond x-ray free electron lasers (XFELs). This multipulse multicolor operation will be broadly applicable to attosecond pump-probe experiments. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Smooth periodic gauge satisfying crystal symmetry and periodicity to study high-harmonic generation in solids

Intense lasers can easily drive nonadiabatic transitions of excited electron wave packets across the Brillouin zones, thus transition dipole moments (TDM) between energy bands of solids should be continuous, satisfying crystal symmetry, and periodic at zone boundaries. While current ab initio algorithms are powerful in calculating band structures of solids, they all introduced random phases into the eigenfunctions at each crystal momentum k. In this work, we show how to choose a “smooth-periodic” gauge where TDMs can be smooth versus k, preserving crystal symmetry, as well as maintaining periodic at boundaries. The symmetry properties of TDMs with respect to k ensure the absence of even-order harmonics from MgO with inversion symmetry, while the TDM in the “smooth-periodic” gauge for broken-symmetry ZnO is responsible for even harmonics that were underestimated in previous simulations. These results reveal the importance of correctly treating the complex TDMs that satisfy crystal symmetry and continuous across zone boundaries in nonlinear laser-solid interactions, which has been elusive in most theories so far.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analytical WKB theory for high-harmonic generation and its application to massive Dirac electrons

Here we propose an analytical approach to high-harmonic generation (HHG) for nonperturbative low-frequency and high-intensity fields based on the (Jeffreys-)Wentzel-Kramers-Brillouin (WKB) approximation. By properly taking into account Stokes phenomena of WKB solutions, we obtain wave functions that systematically include the repetitive dynamics of production and acceleration of electron-hole pairs and quantum interference due to phase accumulation between different pair production times (Stückelberg phase). Using the obtained wave functions without relying on any phenomenological assumptions, we explicitly compute electric current (including intra- and interband contributions) as the source of HHG for a massive Dirac system in (1+1) dimensions under an ac electric field. We demonstrate that the WKB approximation agrees well with numerical results obtained by solving the time-dependent Schrödinger equation and point out that the quantum interference is important in HHG. We also predict in the deep nonperturbative regime that (1) harmonic intensities oscillate with respect to electric-field amplitude E 0 and frequency Ω with a period determined by the Stückelberg phase, (2) the cutoff order of HHG is determined by 2eE 0 ℏΩ 2 , with e being the electron charge, and that (3) noninteger harmonics, controlled by the Stückelberg phase, appear as a transient effect. Our WKB theory is particularly suited for a parameter regime, where the Keldysh parameter γ = (Δ/2)Ω/eE 0 ,with Δ being the gap size, is small. This parameter regime corresponds to intense lasers in the terahertz regime for realistic massive Dirac materials. Our analysis implies that the so-called HHG plateau can be observed at the terahertz frequency within the current technology.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Light-Matter Interaction near the Schwinger Limit Using Tightly Focused Doppler-Boosted Lasers

Strong-field quantum electrodynamics (SF QED) is a burgeoning research topic dealing with electromagnetic fields comparable to the Schwinger field (≈1.32×10 18 V/m). While most past and proposed experiments rely on reaching this field in the rest frame of relativistic particles, the Schwinger limit could also be approached in the laboratory frame by focusing to its diffraction limit the light reflected by a plasma mirror irradiated by a multipetawatt laser. Here, we explore the interaction between such intense light and matter with particle-in-cell simulations. We find that the collision with a relativistic electron beam would enable the study of the nonperturbative regime of SF QED, while the interaction with a solid target leads to a profusion of SF QED effects that retroact on the interaction. In both cases, relativistic attosecond pair jets with high densities are formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Redundant roles of four ZIP family members in zinc homeostasis and seed development in Arabidopsis thaliana

Zinc (Zn) is essential for normal plant growth and development. The Zn-regulated transporter, iron-regulated transporter (IRT)-like protein (ZIP) family members are involved in Zn transport and cellular Zn homeostasis throughout the domains of life. In this study, we have characterized four ZIP transporters from Arabidopsis thaliana (IRT3, ZIP4, ZIP6, and ZIP9) to better understand their functional roles. The four ZIP proteins can restore the growth defect of a yeast Zn uptake mutant and are upregulated under Zn deficiency. Single and double mutants show no phenotypes under Zn-sufficient or Zn-limited growth conditions. In contrast, triple and quadruple mutants show impaired growth irrespective of external Zn supply due to reduced Zn translocation from root to shoot. All four ZIP genes are highly expressed during seed development, and siliques from all single and higher-order mutants exhibited an increased number of abnormal seeds and decreased Zn levels in mature seeds relative to wild type. The seed phenotypes could be reversed by supplementing the soil with Zn. Finally, our data demonstrate that IRT3, ZIP4, ZIP6, and ZIP9 function redundantly in maintaining Zn homeostasis and seed development in A. thaliana.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation

In this work, we propose and develop an arbitrary-order adaptive discontinuous Petrov--Galerkin (DPG) method for the nonlinear Grad--Shafranov equation. An ultraweak formulation of the DPG scheme for the equation is given based on a minimal residual method. The DPG scheme has the advantage of providing more accurate gradients compared to conventional finite element methods, which is desired for numerical solutions to the Grad--Shafranov equation. The numerical scheme is augmented with an adaptive mesh refinement approach, and a criterion based on the residual norm in the minimal residual method is developed to achieve dynamic refinement. Nonlinear solvers for the resulting system are explored and a Picard iteration with Anderson acceleration is found to be efficient to solve the system. Finally, the proposed algorithm is implemented in parallel on MFEM using a domain-decomposition approach, and our implementation is general, supporting arbitrary order of accuracy and general meshes. Furthermore, numerical results are presented to demonstrate the efficiency and accuracy of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Fourth-Order Embedded Boundary Finite Volume Method for the Unsteady Stokes Equations with Complex Geometries

A fourth-order finite volume embedded boundary (EB) method is presented for the unsteady Stokes equations. The algorithm represents complex geometries on a Cartesian grid using EB, employing a technique to mitigate the ``small cut-cell"" problem without mesh modifications, cell merging, or state redistribution. Spatial discretizations are based on a weighted least-squares technique that has been extended to fourth-order operators and boundary conditions, including an approximate projection to enforce the divergence-free constraint. Solutions are advanced in time using a fourth-order additive implicit-explicit Runge-Kutta method, with the viscous and source terms treated implicitly and explicitly, respectively. Formal accuracy of the method is demonstrated with several grid convergence studies, and results are shown for an application with a complex bio-inspired material. In conclusion, the developed method achieves fourth-order accuracy and is stable despite the pervasive small cells arising from complex geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The virtual element method for linear elastodynamics models: Design, analysis, and implementation

We design the conforming virtual element method for the numerical simulation of two dimensional time-dependent elastodynamics problems. We investigate the performance of the method both theoretically and numerically. We prove the stability and the convergence of the semi-discrete approximation in the energy norm and derive optimal error estimates. We also show the convergence in the L 2 norm. The performance of the virtual element method is assessed on a set of different computational meshes, including non-convex cells up to order four in the h-refinement setting. Exponential convergence is also experimentally seen in the p-refinement setting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High-Performance Electron Sources: Numerical Methods and Beam Dynamics at the Precision Frontier

Electron sources have a wide range of applications and there are many stakeholders that express continuing need for improvements and performance enhancements. Whether we consider ultra-cold, high-brightness, high-charge or high-average current source needs, there are some common themes from the point of view of the beam dynamics involved. These are the following ones: ability to model accurately the emission processes, including the presence of often complicated cathode and other boundary surfaces with a wide range of spatial scales; ability to deal accurately and efficiently with a large number of particles interacting pair-wise, including the stochastic part of these interactions with a wide range of spatial scales; and ability to propagate the particle distributions in time, including collisions with a wide range of temporal scales. These tasks require high precision and accuracy since the goal is usually generation, transport and preservation of very high-quality beams. This grant addressed one of the rem

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

A guide to the design of the virtual element methods for second- and fourth-order partial differential equations

Here we discuss the design and implementation details of two conforming virtual element methods for the numerical approximation of two partial differential equations that emerge in phase-field modeling of fracture propagation in elastic material. The two partial differential equations are: (i) a linear hyperbolic equation describing the momentum balance and (ii) a fourth-order elliptic equation modeling the damage of the material. Inspired by, we develop a new conforming VEM for the discretization of the two equations, which is implementation-friendly, i.e., different terms can be implemented by exploiting a single projection operator. We use C 0 and C 1 virtual elements for the second-and fourth-order partial differential equation, respectively. For both equations, we review the formulation of the virtual element approximation and discuss the details pertaining the implementation.

42 ENGINEERING↗