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.
In this work, the formation of a two-phase surface molecular overlayer that transitions from isolated propene molecules to a highly ordered 1D chain structure on Cu(111) is elucidated through combined high-resolution STM imaging and DFT-based calculations. These models reveal how disordered molecules present in-between the 1D chains stabilizes the system as a whole.
Magnetic kagome materials provide a fascinating playground for exploring the interplay of magnetism, correlation and topology. Many magnetic kagome systems have been reported including the binary Fe m X n (X = Sn, Ge; m:n = 3:1, 3:2, 1:1) family and the rare earth RMn 6 Sn 6 (R = rare earth) family, where their kagome flat bands are calculated to be near the Fermi level in the para magnetic phase. While partially filling a kagome flat band is predicted to give rise to a Stoner-type ferromagnetism, experimental visualization of the mag netic splitting across the ordering temperature has not been reported for any of these systems due to the high ordering temperatures, hence leaving the nature of magnetism in kagome magnets an open question. Here, we probe the electronic structure with angle-resolved photoemission spectroscopy in a kagome magnet thin film FeSn synthesized using molecular beam epitaxy. We identify the exchange-split kagome flat bands, whose splitting persists above the magnetic ordering temperature, indicative of a local moment picture. Such local moments in the presence of the topological flat band are consistent with the compact molecular orbitals predicted in theory. We further observe a large spin-orbital selective band renormalization in the Fe d xy + d x 2 -y 2 spin majority channel reminiscent of the orbital selective correlation effects in the iron based superconductors. Our discovery of the coexistence of local moments with topological flat bands in a kagome system echoes similar findings in magic-angle twisted bilayer graphene, and provides a basis for theoretical effort towards modeling correlation effects in magnetic flat band systems.
The aim of this work is to present a detailed code-to-code comparison of two Large-Eddy Simulation (LES) solvers. Corresponding experimental measurements are used as a reference to validate the quality of the CFD simulations. The comparison highlights the effects of solver order on the solutions, and it tries to answer the question of whether a high order solver is necessary to capture the main characteristics of a wind farm. Both solvers were used on different grids to study their convergence behavior. While both solvers show a good match with experimental measurements, it appears that the low order solver is more accurate and substantially cheaper in terms of computational cost.
We propose a novel design of a crab cavity (CC) for the ILC project with low parasitic High Order Modes (HOM) losses and preserving the beam emittance, which is critical for operation with high beam current intensity. Multiple electrodes immersed in the hollow waveguide form a trapped-mode resonator. The transverse components of the electromagnetic field of the trapped dipole mode induce a transverse kick and efficiently deflect charged particles passing through the cavity. We present a scalable design of a superconducting Quasi-waveguide Multicell Resonator (QMiR) seamlessly connected with a beam vacuum chamber. The cavity is completely open at both ends, which significantly reduces the maximum loaded quality factor of the higher order modes (HOM), avoids complex HOM couplers, and thus simplifies the mechanical design of the cavity. The same port is used to feed RF power to the operating mode and to extract the same order modes (SOM). Finally, we estimate the expected cryogenic losses, HOM impedance limits, RF input power required, and frequency tuning for a QMiR cavity designed to operate at 2.6 GHz.
Polarizable nanoparticles are of interest in materials science because of their rich and complex phase behavior that can be used to engineer nanostructured materials with long-range crystalline order. To understand and rationally navigate the design space of polarizable nanoparticles for self-assembling highly ordered superlattices, we developed a coarse-grained computational model to describe the nanoparticle-nanoparticle interactions in implicit solvent and employ the computationally efficient image method to model many-body polarization interactions. We conducted high-throughput virtual screening over a five-dimensional particle design space spanned by temperature, particle size, particle charge, particle dielectric, and solvent dielectric using enhanced sampling molecular dynamics calculations within an active learning framework to efficiently map out the regions of thermodynamic stability of the self-assembled aggregates. We validate our predictions in comparisons against small angle x-ray scattering measurements of gold nanoparticles surface functionalized with metal chalcogenide ligands. Lastly, we use our validated phase maps to computationally design switchable nanostructured materials capable of triggered assembly and disassembly as a function of temperature and solvent dielectric with potential applications as sensors, smart windows, optoelectronic devices, and in medical diagnostics.
Supramolecular polymers can emulate some of the physical properties of covalent polymers but offer new opportunities given the possibility of designing monomers that will form highly ordered assemblies with defined shapes. Internally ordered supramolecular polymers formed through nucleation–elongation self-assembly are well-known but highly crystalline examples which exhibit important properties such as light harvesting, charge transport, and ferroelectricity are not common. We report here on a detailed study of supramolecular polymers formed in water by carboxylated naphtho-p-quinodimethane amphiphiles. We found that supramolecular polymerization of these amphiphiles in aqueous media yields crystalline assemblies with morphologies that included ribbons, helically rolled ribbons, and twisted filaments. This polymorphism was found to be controlled exclusively by repulsive electrostatic interactions controlled by the degree of protonation of the carboxylic head groups which also dictates the nature of supramolecular packing. Substoichiometric amounts of base lead to highly crystalline ribbons due to a decreased surface charge density and less electrostatic repulsion. Increasing deprotonation results in helically rolled ribbons with a different polymorph crystal lattice, whereas excessive deprotonation leads to twisted filaments with maximum surface charge density. Ribbons, helical rolled ribbons, and twisted filaments revealed an increasing red shift in their visible absorption maxima. These crystalline assemblies could be potential candidates for solar energy materials and photocatalytic systems.
We report solids irradiated by energetic photons can be eroded in two modes, depending on the radiation intensity. High average, low-peak power sources, e.g., synchrotron radiation and high-order harmonics, induce desorption of the material at a low etch rate. In contrast, high-peak-power radiation from extreme ultraviolet and x-ray lasers usually causes a massive removal of the material even by a single shot. In this contribution, an effective material erosion is reported in PMMA exposed to multiple accumulated pulses generated by the free-electron x-ray-laser Linac Coherent Light Source (LCLS, tuned at a photon energy of 830 eV in this study, operated in Menlo Park at Stanford, CA, USA) at a fluence below the single-pulse ablation threshold. The effect is caused by polymer-chain scissions initiated by single photons carrying enough energy to break the C—C bounds. High efficiency of the erosion is supposed to occur due to a correlation of the single-photon effects. The subthreshold damage exhibits a nonlinear dose dependence resulting from a competition between chain scissions and cross-linking processes. The cross-linking is proven by Raman spectroscopy of the irradiated polymer. Two theoretical models of the x-ray free-electron-laser-induced erosion are suggested, which provide an excellent agreement with the experimental results.
We investigate identifying differential equations in the frequency domain. Fourier analysis is an important tool in theoretical analysis and numerical solvers of differential equations, yet there is limited work in exploring this connection in the identification of differential equations. This paper aims to identify the underlying differential equation in the frequency domain, from a given single realization of the differential equation perturbed by noise. Such setting imposes difficulties which are different from other identification methods where computation is carried out in the physical domain. We propose several ways to mitigate the challenges arising from noise in data and large differences in the magnitudes of frequency responses. The main takeaways are that identifying differential equations solely in the frequency domain is challenging, the method we propose is based on a form of domain partitions in the frequency domain, and this method shows benefits for complex data even with high level of noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain and to enhance the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms, and shows advantages on complex data with many frequency modes, even under high level of noise.
High-energy-density (HED) hydrodynamics studies such as those relevant to inertial confinement fusion and astrophysics require highly disparate densities, temperatures, viscosities, and other diffusion parameters over relatively short spatial scales. This presents a challenge for high-order accurate methods to effectively resolve the hydrodynamics at these scales, particularly in the presence of highly disparate diffusion. A significant volume of engineering and physics applications use an unstructured discontinuous Galerkin (DG) method developed based on the finite element mesh generation and algorithmic framework. This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection–diffusion equations. However, to the best of the authors’ knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. As a result, the algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.
Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.
Exciting solids with intense femtosecond laser pulses prompts electrons of the interrogated material to respond in a highly nonlinear manner, as is evident in the emission of high-order harmonic radiation and photoelectrons with kinetic energies well above that of the driving photons. Such high-field interactions can be resolved, for example, in above-threshold multiphoton photoemission (ATP) spectroscopy. In this work, we interrogate the nonlinear photoelectric responses of the pristine copper, silver, and gold noble metal surfaces in (111) and (100) crystal orientations in the perturbative regime. Using multiphoton photoemission spectroscopy (mPP) excited by finely tuned optical fields, we characterize enhancement of the mPP and ATP yields from (111) surfaces in selected k||-momentum ranges when the occupied Shockley surface (SS) states are (near-)resonantly coupled by multiphoton transitions to image potential (IP) intermediate states in the excitation process. The ATP signal from the IP states of (111) surfaces is largely defined by their formation through polarization of SS electrons; this observation is contrasted with ATP experiments from the Ag(100) surface, for which the SS becomes an unoccupied resonance and the IP states can only be excited from bands with significantly more bulk character. In addition, based on the optical power and nonlinear order-dependent mPP spectra, we provide evidence for ATP being a one-step, rather than a sequential process, as previously postulated.
Context. 3C 84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of very-long-baseline interferometry (VLBI) above the hitherto available maximum frequency of 86 GHz.Aims. Using ultrahigh resolution VLBI observations at the currently highest available frequency of 228 GHz, we aim to perform a direct detection of compact structures and understand the physical conditions in the compact region of 3C 84.Methods. We used Event Horizon Telescope (EHT) 228 GHz observations and, given the limited (u, v)-coverage, applied geometric model fitting to the data. Furthermore, we employed quasi-simultaneously observed, ancillary multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure.Results. We report the detection of a highly ordered, strong magnetic field around the central, supermassive black hole of 3C 84. The brightness temperature analysis suggests that the system is in equipartition. We also determined a turnover frequency of νm = (113 ± 4) GHz, a corresponding synchrotron self-absorbed magnetic field of BSSA = (2.9 ± 1.6) G, and an equipartition magnetic field of Beq = (5.2 ± 0.6) G. Three components are resolved with the highest fractional polarisation detected for this object (mnet = (17.0 ± 3.9)%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017–2018. We report a steeply negative slope of the spectrum at 228 GHz. We used these findings to test existing models of jet formation, propagation, and Faraday rotation in 3C 84.Conclusions. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C 84. However, systematic uncertainties due to the limited (u, v)-coverage, however, cannot be ignored. Our upcoming work using new EHT data, which offer full imaging capabilities, will shed more light on the compact region of 3C 84.Key words: techniques: high angular resolution / techniques: interferometric / galaxies: active / galaxies: individual: NGC 1275 / galaxies: jets
Explore the source record for details and available documents.
Copper and its oxides are key catalytic materials, on which reactions often occur at the metal/oxide interface. Here, in this work, we directly connect the induction period observed during methanol-driven reduction of thin-film copper oxides to their atomic-scale structural order. Using temperature-programmed desorption (TPD) methanol titrations combined with scanning tunneling microscopy, we show that highly ordered oxide phases – particularly the “29” structure with its low defect density – exhibit long induction periods and initially low reactivity. The induction period, defined as the number of methanol TPD cycles required to reach half of the maximum formaldehyde yield, scales with oxide order and oxygen coverage. Enhanced reactivity of well-ordered oxides emerges only after repeated methanol adsorption/desorption cycles generate oxygen vacancies and new Cu(111)/Cu x O interfacial sites. In contrast, disordered or sub-stoichiometric oxides, which contain more intrinsic defects and interfaces, are active from the first TPD cycle. We further examine how dilute Pt and Rh dopants influence oxide order and reactivity: 1% Pt increases defect density and catalytic activity, while 1% Rh promotes oxide ordering and longer induction periods. These findings demonstrate that dilute alloying provides a potential method for tuning the structure and reactivity of Cu(111)/Cu x O interfaces.
For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10 6 times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).
The structural rigidity of fused units in the polymer backbone, in addition to the resulting stabilizing effect of the quinoidal structure, and tunable electronic properties have played a key role in promoting highly-ordered π-stacking moieties, exhibiting promising charge carrier mobilities. The electron-deficient thiazole moiety shows high planarity and effective π–π stacking, which leads to the reduction in the energy levels of the highest occupied and lowest unoccupied molecular orbitals (HOMO/LUMO), and ideally enhances the electron charge mobility. Four heterocycle-based monomers BTzS, BTzSe, DTS, and DTG based on fused bithiazole and dithiophene units incorporated with sulfur, selenium, silicon, and germanium as the bridging atoms were synthesized and characterized. Here, the monomers were copolymerized with the electron-rich alkylated thienylenevinylene (TV) unit to afford copolymers P1–P4. The thermal, optical, and electrochemical properties and crystallinity of the copolymers were thoroughly investigated. Extensive OFET device optimization using different solvents and annealing temperatures resulted in the best charge mobility of 0.09 cm2 V –1 s –1 for the electron-deficient bithiazole BTzS copolymer P1 and 0.36 cm2 V –1 s –1 for the DTS copolymer P3.
In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.
Although still a very active area of research with many open questions, there exist highly-accurate and efficient algorithms for numerically estimating solutions to the equations of fluid motion. Similarly, great strides have been made in numerically modeling the motion of solids. However, there exist many practical applications where one must solve both fluid and structure in tandem. There is a gap in the literature for high order unconditionally stable partitioned algorithms for problems of this type. We investigate a new partitioned scheme which shows promise in filling this gap. We use a finite-element-based model with an added mass approach on problems coupling the Euler beam equation with the incompressible Navier-Stokes equations.