SEARCH · Engineering Papers
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.
Code for my PhD thesis that creates meshes using gmsh and code I wrote to create high-order hybrid meshes to use for testing a hybrid method between spectral element method and generalized finite differences. It also uses a tetrahedral mesh improvement method called "stellar" to improve mesh quality near a curved boundary.
PIAFS is a finite-difference code to solve the compressible Navier-Stokes equations with chemical heating on Cartesian grids. It models chemical reactions of air (oxygen and carbon dioxide) with ozone subject to radiation. It uses a high-order WENO spatial discretization and explicit Runge-Kutta time integration. It is capable of parallel simulations using MPI. The code is written in C/C++.
The mammalian ocular lens is an avascular multicellular organ that grows continuously throughout life. Traditionally, its cellular organization is investigated using dissected lenses, which eliminates in vivo environmental and structural support. Therefore, in vivo optical imaging methods for studying lenses in their native context in live animals are urgently needed. Here, we demonstrated that two-photon fluorescence microscopy can visualize lens cells in vivo. To maintain subcellular resolution at depth, we used adaptive optics to correct aberrations owing to ocular and lens tissues, which led to substantial signal and resolution improvements. Imaging lens cells up to 980 µm deep, we observed novel cellular organizations including suture-associated voids, enlarged vacuoles, and large cavities, contrary to the conventional view of a highly ordered organization. We tracked these features longitudinally over weeks and observed the incorporation of new cells during growth. Taken together, noninvasive longitudinal in vivo imaging of lens morphology using adaptive optics two-photon fluorescence microscopy will allow us to observe the development or alterations of lens cellular organization in living animals directly.
The objective of this paper is to develop a fully integrated model for the wheel loader, including the subsystem dynamics of the engine, drivetrain, working circuit, steering circuit, and vehicle. It leads to a high-order strongly nonlinear system, and all state variables are coupled together to form a Multi-Input and Multi-Output (MIMO) system. A control architecture is proposed to decouple the MIMO system into several Single-Input and Single-Output (SISO) systems. Here, a tracking problem has been formulated to validate this fully integrated model with the field test data. The accuracy of the model is verified by the 2.3% difference between the measured and simulated fuel consumption. Meanwhile, an energy distribution analysis is conducted to reveal the energy consumption and energy loss of each portion of the wheel loader. Such a model can be used to plan the working pattern, guide the driving habits of human operators, or refine the underlying architecture, leading to the ultimate goal of reducing total fuel consumption and improving productivity.
We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.
The advent of chirped-pulse amplification in the 1980s and femtosecond Ti:sapphire lasers in the 1990s enabled transformative advances in intense laser–matter interaction physics. Whereas most of experiments have been conducted in the limited near-infrared range of 0.8–1 μm, theories predict that many physical phenomena such as high harmonic generation in gases favor long laser wavelengths in terms of extending the high-energy cutoff. Significant progress has been made in developing few-cycle, carrier-envelope phase-stabilized, high-peak-power lasers in the 1.6–2 μm range that has laid the foundation for attosecond X ray sources in the water window. Even longer wavelength lasers are becoming available that are suitable to study light filamentation, high harmonic generation, and laser–plasma interaction in the relativistic regime. Long-wavelength lasers are suitable for sub-bandgap strong-field excitation of a wide range of solid materials, including semiconductors. In the strong-field limit, bulk crystals also produce high-order harmonics. In this review, we first introduce several important wavelength scaling laws in strong-field physics, then describe recent breakthroughs in short- (1.4–3 μm), mid- (3–8 μm), and long-wave (8–15 μm) infrared laser technology, and finally provide examples of strong-field applications of these novel lasers. Furthermore, some of the broadband ultrafast infrared lasers will have profound effects on medicine, environmental protection, and national defense, because their wavelengths cover the water absorption band, the molecular fingerprint region, as well as the atmospheric infrared transparent window.
Short wavelength high-harmonic sources are undergoing intense development for applications in spectroscopy and microscopy. Despite recent progress in peak and average power, spatial control over coherent extreme ultraviolet (XUV) beams remains a formidable challenge due to the lack of suitable optical elements for beam shaping and control. In this work, we demonstrate a robust and precise approach that structures XUV high-order harmonics in space as they are emitted from a nanostructured MgO crystal. Our demonstration paves the way for bridging the numerous applications of shaped light beams from the visible to the short wavelengths, with potential uses for applications in microscopy and nanoscale machining.
Non-linear processes such as four-wave-mixing have become instrumental in attosecond EUV spectroscopy. Using EUV high harmonics in conjunction with collinear near-infrared and mid-infrared fields, we extended this approach to high-order-mixing between three colors. Specifically, we find that atomic resonances in neon exhibit a significant cross section for six-wave mixing. The MIR probe frequency tunability in our multicolor scheme is employed to access several optically dark resonances and probe the quantum beat of a coherent electronic wavepacket using background-free EUV emission as a diagnostic. This technique can be easily extended to other atomic and molecular systems, and opens the door to multi-dimensional non-linear spectroscopy.
In this work, we study the point-by-point inscription of sapphire parallel fiber Bragg gratings (sapphire pFBGs) in a fully multimode system. A parallel FBG is shown to be critical in enabling detectable and reliable high-order grating signals. The impacts of modal volume, spatial coherence, and grating location on reflectivity are examined. Three cascaded seventh-order pFBGs are fabricated in one sapphire fiber for wavelength multiplexed temperature sensing. Using a low-cost, fully multimode 850-nm interrogator, reliable measurement up to 1500°C is demonstrated.
Resonant high-order harmonics, which result in quasimonochromatic extreme ultraviolet light with coherent intensity enhancement involving autoionizing resonances, have been demonstrated from laser-ablated plumes in the tunnel-ionization regime. Here, we demonstrate resonant harmonics in the previously unexplored multiphoton-ionization regime. We demonstrate an intense resonant harmonic from gallium with an intensity enhancement ratio of 714 relative to the neighboring harmonics, achieved without the need for extreme ultraviolet filtering methods, thus preventing a typical photon flux loss of more than 70%. Three-dimensional time-dependent Schrödinger equation calculations reveal that this increase in the enhancement ratio is due to the lowelectron wave packet spreading in the multiphoton-ionization regime. These results reveal a method for increasing the intensity and monochromaticity of intense multimicrojoule femtosecond extreme ultraviolet light and will also facilitate understanding of the involvement of autoionizing resonances in generating resonant harmonics in the multiphoton-ionization regime.
Fabry–Pérot (F–P) cavity and metal hole array are classic photonic devices. Integrating F–P cavity with holey metal typically enhances interfacial reflection and dampens wave transmission. In this work, a hybrid bound surface state is found within rectangular metal holes on a silicon substrate by merging an extraordinary optical transmission (EOT) mode and a high-order F–P cavity mode both spatially and spectrally. Transmission, Q-factor, and bandwidth can be enhanced significantly with respect to the classical EOT and F–P interference by simply sweeping the cavity length. This state can provide EOT properties and ten times broader EOT bandwidth well below the effective plasma frequency of the periodic metal holes, where the metal holes typically show evanescent properties and do not support EOT in theory. Furthermore, a large modulation range of 25 % and 39 % is demonstrated with various graphene patterns for the transmittance of this hybrid state at 500 and 582 GHz, respectively.
Low emittance electron beams are of high importance at facilities like the LCLS-II at SLAC. Emittance dilution effects due to off-axis beam transport for a TESLA-type cryomodule (CM) have been shown at the Fermilab Accelerator Science and Technology facility. The results showed the correlation between the electron beam-induced cavity high-order modes (HOMs) and submacropulse centroid slewing and oscillation downstream of the CM. Mitigation of emittance dilution can be achieved by reducing the HOM signals and the variances in the submacropulse beam positions downstream of the CM. Here we present a Machine Learning based optimization and model construction for HOM signal level reduction using Neural Networks and Gaussian Processes. To gather training data we performed experiments using single bunch and 50 bunch electron beams with charges up to 125 pC/b. We measured HOM signals of all cavities and beam position with a set of BPMs downstream of the CM. The beam trajectory was changed using V/H125 corrector set located upstream of the CM. The results presented here will inform the LCLS-II injector commissioning and will serve as a prototype for HOM reduction and emittance preservation.
Experiments at the Fermilab Accelerator Science and Technology* (FAST) facility detected electron beam-induced high order mode (HOM) signals from Tesla superconducting cavities. This paper describes some of the signal detection hardware used in this experiment, as well as measurements of the HOM signal magnitude versus beam trajectory. These measurements were made both with a single bunch and with a train of 50 bunches at bunch charges from 400 pC/b down to 10 pC/b. The detection hardware is designed for use with the Tesla superconducting cavities of LCLS-II at SLAC** and is based on a prototype already in use at Fermilab. The HOM signal passes through a bandpass filter that is centered on several cavity dipole modes and a zero bias Schottky diode detects its magnitude. Direct comparisons were made between the FNAL chassis and the SLAC prototype for identical beam steering conditions. To support measurements with bunch charges as low as 10 pC, the SLAC detector has RF amplification between the bandpass filter and the diode detector. With this hardware, usable HOM signal measurements are obtained with a single bunch of 10 pC in cryomodule cavities as will be needed for LCLS-II.
Jefferson Lab has an ongoing R&D program in plasma processing. The experimental program investi-gated processing using argon/oxygen and heli-um/oxygen gas mixtures. Plasma processing is a com-mon technique where the free oxygen produced by the plasma breaks down and removes hydrocarbons from surfaces. This increases the work function and reduces the secondary emission coefficient. The initial focus of the effort was processing C100 cavities by injecting RF power into the high order mode (HOM) coupler ports. We also developed the methods for establishing a plasma in C75 cryomodules where the RF power is injected via the fundamental power-coupler. Four C100 cryomodules were in situ processed in the CE-BAF accelerator in May 2023 with the cryomodules returning to an operational status in Sept. 2023. The overall operational energy gain for the four cryomod-ules was 49 MeV. Methods, systems and results from processing cryomodules in the CEBAF accelerator and vertical test results are presented. Current status and future plans are discussed.
The Electron Ion Collider (EIC) is a new nuclear phys-ics research facility to be constructed at Brookhaven Lab (BNL). It consists a few accelerator complexes, each re-quires some different SRF and RF systems. Among them, the electron storage ring?s (ESR) 591 MHz fundamental RF system is one of the most challenging. This 17 cavi-ties system will handle a variable beam current of up to 2.5 A, and replenish up to 10 MW of beam power losses from both synchrotron radiation (SR) and high-order modes (HOM). In this paper, we will report the progress and findings of the ongoing design and prototyping of this cavity.
Truncated Power Series Algebra (TPSA), or Differential Algebra (DA), is a well-established tool in accelerator physics, commonly used for generating high-order maps of dynamic systems, as well as in symplectic tracking, normal form analysis, verified integration, optimization, and fast multipole methods. This package is the first to perform symbolic DA computations, enabling traceability of initial condition contributions and runtime reduction for repeated DA calculations, potentially expanding DA’s applications.
One of the primary methods used by the United States to estimate the yields of atmospheric nuclear tests performed in the 1940’s, 50’s, and 60’s was to film each detonation with a high-speed camera, measure the blast radius as a function of time, and then use Taylor’s spherical blast wave equation to calculate the source energy (i.e., the yield). For surface detonations, the yield in Taylor’s equation was assumed to be twice as large as the actual yield to account for the energy reflected from the surface, assuming the ground was a perfect reflector. Recent simulations of surface detonations over water (referred to as barge shots) have shown that Taylor’s equation must be modified further to account for material entrainment and heat losses. We have computed this lost energy using a high-order radiation-hydrodynamics code to determine the reduction in blast radius vs. time. By analyzing simulation data in the same manner as film data, we have determined that a yield correction factor of 1.27± 0.04 should be applied to the barge shots.
The SPARC team has completed much of the workflow development needed for the coupled and uncoupled analyses needed to support our flight test validation work. The SPARC team has continued to make progress on development activities to support unsteady, full reentry vehicle aero analysis, which have focused on turbulence modeling, uniform mesh refinement, in-situ visualization (with the Paraview/Catalyst team), and performance analysis. The SPARC team is making solid progress towards our Q4 goal of demonstrating an improved Wall Modeled Large Eddy Simulation (WMLES) capability for standard geometries using high-order finite difference, discontinuous Galerkin, and low-dissipation finite volume methods. The SPARC team has met the minimum completion criteria for parts of the L1, including documenting performance and scaling of SPARC on Trinity, Sierra and Astra, and performing runs for the flight test validation case.