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.
Static torque and no load test data for a 1hp, 300rpm axial-flux magnetically geared generator prototype developed by Texas A&M EMPE Lab.
We theoretically investigate modulation instability (MI) in a two-core nonlinear oppositely directed coupler with a negative-index material channel in the presence of a non-kerr-like nonlinear refractive-index change. Using linear stability analysis, we obtain an expression for the instability gain. It is found that the MI in the nonlinear two-core oppositely directed coupler is significantly influenced by the ratio of the forward- to backward-propagating wave power, nonlinearity, power and saturation. In the normal dispersion regime, the septic nonlinearity is found to suppress the MI by reducing both the gain and width of the instability nonlinearity. But in anomalous dispersion regime, the situation is different because septic nonlinearity enhances the MI by increasing both the gain and width of MI. The increase of septic nonlinearity both in normal and anomalous dispersion regimes increases the gain of the instability nonlinearity. Thus, we report new ways to generate and manipulate the MI and solitons in two-core oppositely directed couplers with a particular emphasis on a negative-index material (NIM) channel in presence of septic nonlinearity.
It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.
In this study, we explore the role of heave plate topology on fluid reaction forces using three, representative shapes: a hexagonal flat plate, a hexagonal conic with an open top, and the same with a closed top that encloses a fluid mass. We force each test article sinusoidally in a quiescent tank and decompose the reaction force using forms of the Morison equation for phase-invariant and phase-dependent parameterizations. We find that a flat plate generates 5.3 % more fluid reaction force than the open conic topology, and 21.4% more than the enclosed conic. Asymmetric topologies generate asymmetric reaction forces, but the magnitude of asymmetry is limited by nearly symmetric fluid inertia forces, which dominate over drag for these test articles. Additionally, we observe asymmetric vortex dynamics for the flat plate when the Keulegan-Carpenter number is between 1 and 2, accompanied by a shift in the phase of the peak force by of the oscillation period. As a consequence of this shift, the hydrodynamic coefficients estimated from the phase-dependent Morison equation decomposition are asymmetric, suggesting that phase-dependent representations may not provide physical insight in some hydrodynamic regimes.
In order to improve their long-term viability, wave energy converters (WECs) need to be able to shed loads when a threshold wave condition is exceeded. As shown by Tom et al. (2016) [1], provision of adjustable flaps within the body of an oscillating surge wave energy converter (OSWEC) allows wave energy to pass through the device. A control system may then be able to open and close the flaps when waves approaching the device exceed preset thresholds. The variable-geometry OSWEC (VG-OSWEC) concept studied in this paper is a bottom-hinged, rectangular wave paddle with five flaps of elliptical cross-section embedded into the face of the paddle. System ID tests were conducted on this VG-OSWEC device at a 1:14 scale in a wave basin. Free decay tests showed that the damping was distinctly nonlinear when the flaps were fully open, and the natural frequency increased by 40% when compared with the flaps in a fully closed configuration. Tests with regular wave conditions were used to develop the response amplitude operator for the rotational motion about the hinge. These response amplitude operator results when compared with numerical simulations run using WEC-Sim/WAMIT and ANSYS AQWA, show strong agreement with the flap open and closed conditions. The regular-wave condition measurements also show that the wave excitation moment about the hinge was reduced by up to 60% when the flaps were fully open. The experiments serve to demonstrate the potential of the variable geometry design to shed loads and survive harsh ocean environments.
Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.
During the 9 March 2018 event with two consecutive interplanetary shocks compressing the dayside magnetosphere, the azimuthal mode structure and frequency spectrum of ultra low frequency magnetic pulsations are resolved using a cross-spectral analysis based on high-fidelity multi-probe Magnetospheric Multiscale Mission (MMS) magnetometer data. The results based on the MMS 4 and MMS 3 pair of measurements show that shock arrival leads to low mode (|m| ≤ 3) magnetic fluctuations in the Pc4-5 regimes, and smaller spatial scale fluctuations implied by the dominant high mode numbers are observed after both shock signatures hit and passed the magnetosphere. Detailed evolution of the mode structure is also shown for the first shock to reveal the development of high mode structure from a bump-on-tail distribution at m≈20 to a dominant peak at m≈50 in about 10 min. In addition, an interesting change of sign in m from negative to positive is observed as MMS crosses ~11 MLT pre-noon, which is consistent with the picture of wave generation by dayside magnetopause compression and then anti-sunward propagation. Finally, for both shocks, the contribution of higher frequency waves (Pc-4 range compared with Pc-5) to the total wave power is found to be negligible before and after the shock impact, but it becomes more significant during the shock impact.
Abstract Baroclinic or extratropical cyclones (ETCs) transport heat and moisture to higher latitudes, making it fundamentally important to understand how their influence changes as Earth's climate evolves. A 2–8‐day Lanzcos bandpass filter is applied to European Center for Medium Range Weather Forecasting 5th Generation Reanalysis latent energy (LE) and kinetic energy (KE) data to assess how ETCs have changed from 1940 to 2020 relative to full‐scale changes in LE and KE. Full‐scale KE trends are more positive at high latitudes relative to mid‐latitudes, confirming several previous studies that ETCs have shifted poleward. LE increases have occurred globally, and trends in both full‐scale LE and KE are statistically significant in the southern high latitudes. The high relative fractional contribution of 2–8‐day LE wave power and trend clearly suggest that ETCs have an increasingly important role in poleward moisture transport but are not solely responsible for the observed statistically significant increases.
Here, we compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.
Linear mode conversion of Langmuir/z waves to electromagnetic radiation near the plasma and upper hybrid frequency in the presence of density gradients is potentially relevant to type II and III solar radio bursts, ionospheric radar experiments, pulsars, and continuum radiation for planetary magnetospheres. Here, we study mode conversion in warm, magnetized plasmas using a numerical electron fluid simulation code when the density gradient has a wide range of angle, δ, to the ambient magnetic field, B{sub 0}, for a range of incident Langmuir/z wavevectors. Our results include: (1) Left-handed polarized ordinary (oL) and right-handed polarized extraordinary (xR) mode waves are produced in various ranges of δ for Ω{sub 0} = (ωL/c){sup 1/3}(ω{sub ce}/ω) < 1.5, where ω{sub ce} is the (angular) electron cyclotron frequency, ω is the angular wave frequency, L is the length scale of the (linear) density gradient, and c is the speed of light; (2) the xR mode is produced most strongly in the range, 40° < δ < 60°, for intermediately magnetized plasmas with Ω{sub 0} = 1.0 and 1.5, while it is produced over a wider range, 0° ≤ δ ≤ 90°, for weakly magnetized plasmas with Ω{sub 0} = 0.1 and 0.7; (3) the maximum total conversion efficiencies for wave power from the Langmuir/z mode to radiation are of order 50%–99% and the corresponding energy conversion efficiencies are 5%–14% (depending on the adiabatic index γ and β = T{sub e}/m{sub e}c{sup 2}, where T{sub e} is the electron temperature and m{sub e} is the electron) for various Ω{sub 0}; (4) the mode conversion window becomes wider as Ω{sub 0} and δ increase. Hence, the results in this paper confirm that linear mode conversion under these conditions can explain the weak total circular polarization of interplanetary type II and III solar radio bursts because a strong xR mode can be generated via linear mode conversion near δ ∼ 45°.
The accurate and efficient description of strong electronic correlations remains an important objective in electronic structure theory. Projected Hartree-Fock theory, where symmetries of the Hamiltonian are deliberately broken and projectively restored, all with a mean-field computational scaling, shows considerable promise in this regard. However, the method is neither size extensive nor size consistent; in other words, the correlation energy per particle beyond broken-symmetry mean field vanishes in the thermodynamic limit, and the dissociation limit of a molecule is not the sum of the fragment energies. These two problems are closely related. Recently, Neuscamman [Phys. Rev. Lett. 109, 203001 (2012)] has proposed a method to cure the lack of size consistency in the context of the antisymmetrized geminal power wave function (equivalent to number-projected Hartree-Fock-Bogoliubov) by using a Jastrow-type correlator in Hilbert space. Here, we apply the basic idea in the context of projected Hartree-Fock theory, linearizing the correlator for computational simplicity but extending it to include spin fluctuations. Results are presented for the Hubbard Hamiltonian and for some simple molecular systems.
In this work we show the possibility to extract Kohn-Sham orbitals, orbital energies, and exchange correlation potentials from accurate Quantum Monte Carlo (QMC) densities for atoms (He, Be, Ne) and molecules (H{sub 2}, Be{sub 2}, H{sub 2}O, and C{sub 2}H{sub 4}). The Variational Monte Carlo (VMC) densities based on accurate Jastrow Antisymmetrised Geminal Power wave functions are calculated through different estimators. Using these reference densities, we extract the Kohn-Sham quantities with the method developed by Zhao, Morrison, and Parr (ZMP) [Phys. Rev. A 50, 2138 (1994)]. We compare these extracted quantities with those obtained form CISD densities and with other data reported in the literature, finding a good agreement between VMC and other high-level quantum chemistry methods. Our results demonstrate the applicability of the ZMP procedure to QMC molecular densities, that can be used for the testing and development of improved functionals and for the implementation of embedding schemes based on QMC and Density Functional Theory.
Here, a theory is proposed that describes mutual conversion of two electromagnetic modes in cold low-density plasma, specifically, in the high-frequency limit where the ion response is negligible. In contrast to the classic (Landau–Zener-type) theory of mode conversion, the region of resonant coupling in low-density plasma is not necessarily narrow, so the coupling matrix cannot be approximated with its first-order Taylor expansion; also, the initial conditions are set up differently. For the case of strong magnetic shear, a simple method is identified for preparing a two-mode wave such that it transforms into a single-mode wave upon entering high-density plasma. The theory can be used for reduced modeling of wave-power input in fusion plasmas. In particular, applications are envisioned in stellarator research, where the mutual conversion of two electromagnetic modes near the plasma edge is a known issue.
We show that for both single-Slater-Jastrow and Jastrow geminal power wave functions the formal cost scaling of Hilbert space variational Monte Carlo can be reduced from fifth to fourth order in the system size, thus bringing it in line with the long-standing scaling of its real space counterpart. While traditional quantum chemistry methods can reduce costs related to the two-electron integral tensor through various tensor decomposition methods, we show that such approaches are ineffective in the presence of Hilbert space Jastrow factors. Instead, we develop a simple semi-stochastic approach that can take similar advantage of the near-sparsity of this four-index tensor. Through demonstrations on alkanes of increasing length, we show that accuracy and overall statistical uncertainty are not meaningfully affected and that a total cost crossover is reached as early as 50 electrons when using a minimal basis. Further study will be needed to assess where the crossover occurs in more compact molecular geometries and larger basis sets and to explore how in that context the crossover can be accelerated.
A carefully designed waveguide-based millimeter-wave notch filter, operating at 140 GHz, safeguards plasma diagnostic instruments from gyrotron leakage. Here, utilizing cylindrical cavity resonators with aperture coupling, the filter efficiently resonates 140 GHz wave-power into the TE 11p mode, optimizing various geometrical parameters for practical fabrication and high-yield production. Thorough thermal analysis ensures its ability to handle power. The filter achieves outstanding performance with over 90 dB rejection at 140 GHz while providing low insertion loss over the passband (110–138 GHz), which is ideally suited for system-on-chip approach F-band diagnostic system applications.
Laser diode absorption spectroscopy is used to experimentally measure Ar(1s 5 ) metastable density and translational gas temperature within a 94 GHz microplasma. A square two-dimensional photonic crystal (PhC) at this resonance frequency serves to ignite and sustain the plasma from 20 to 200 Torr (2.7 × 10 3 –2.7 × 10 4 Pa) by using millimeter wave power from 300 to 1000 mW. Metastable density within the plasma is estimated from the absorption line shape of the laser traversing the PhC. The metastable density reaches an order of 10 19 m –3 at lower pressure and decreases as pressure increases. From the Lorentzian line shape of the absorption profile at 811.53 nm, the gas temperature is extracted and found to increase from 500 K at 20 Torr to 1300 K at 200 Torr. Furthermore, these data are analyzed and compared with a zero-dimensional plasma model and with previous experimental plasma results at 43 GHz.
Here, noninductive plasma current start-up using 2nd harmonic electron cyclotron resonance heating (ECRH) with oblique radio frequency (RF) injection is demonstrated in a Q-shu University experiment with steady-state spherical tokamak. A strong transition was observed in the heating and plasma current ramp-up. The initial bulk electron heating regime exhibits T ebulk ~ 140 eV and no hard x-ray (HXR) emission with a low I p of ~15 kA; it abruptly transitions to a regime that exhibits a low T ebulk of ~10 eV and a strong HXR emission with a high I p of ~50 kA. This behavior is distinctly different from that observed in previous fundamental ECRH experiments. The mechanism of the heating and current drive transition are investigated considering wave power absorption and plasma power balance. The results indicate that the transition is caused by the favorable heating of tail electrons where the RF power absorption at the 2nd harmonic increases nearly linearly with T etail , while the power transfer from the tail electrons to the bulk electrons decreases with 1/T etail 0.5 . This causes a rapid transition to a state with high T etail while reducing T ebulk towards colder ion temperature. The understanding of the transition mechanism helps to consider plasma current start-up using 2nd harmonic ECRH for tokamak reactors such as JT-60 SA and ITER.
Despite high demonstrated efficiency, lower-hybrid current drive (LHCD) has not been considered localized enough for neoclassical tearing mode (NTM) stabilization in tokamaks. This assessment must be reconsidered in view of the radiofrequency current condensation effect. Here we show that an island with a central hot spot induces significant localization of LHCD. Furthermore, in steady state tokamaks where a significant amount of current is provided by LHCD, passive stabilization of NTMs may occur automatically, particularly as islands become large, without requiring precise aiming of the wave power.