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.
Some of the significant interference fields that may affect stability of aircraft at supersonic speeds are briefly summarized. Illustrations and calculations are presented to indicate the importance of interference fields created by wings, bodies, wing-body combinations, jets, and nacelles.
We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.
We develop for the first time a microscopic global nucleon-nucleus optical potential with quantified uncertainties suitable for analyzing nuclear reaction experiments at next-generation rare-isotope beam facilities. Within the improved local density approximation and without any adjustable parameters, we begin by computing proton-nucleus and neutron-nucleus optical potentials from a set of five nuclear forces from chiral effective field theory for 1800 target nuclei in the mass range 12 ≤ A ≤ 242 for energies between 0 MeV < E ≲ 150 MeV. We then parameterize a global optical potential for each chiral force that depends smoothly on the projectile energy as well as the target nucleus mass number and isospin asymmetry. Uncertainty bands for elastic scattering observables are generated from a full covariance analysis of the parameters entering in the description of our global optical potential and benchmarked against existing experimental data for stable target nuclei. Since our approach is purely microscopic, we anticipate a similar quality of the model for nucleon scattering on unstable isotopes.
This article describes the design, assembly, and implementation of a hand-held, magnetic-field-based sensor system that can be adapted for a variety of sensing applications. The miniaturized system is based on Chemical Identification by Magneto-Elastic Sensing (ChIMES) technology, which uses three concentric solenoid coils to wirelessly interrogate a sensor body comprised of a response material coupled to a magnetoelastic wire. The response material expands when it encounters a target, imposing mechanical stress on the wire and altering its magnetic permeability. The sensor bodies are passive, requiring no external power source, and they are small, measuring about 15 mm in length and 3.0 mm in diameter. Up to four sensor bodies can be configured as an evenly-spaced linear array. The sensor system operates by applying a low-frequency, current-stabilized, filtered triangle wave to a uniform-density excitation coil to switch the magnetic domains within the wire. Further, the responses from the sensors are picked up by a detection coil as stress-induced changes in the Faraday voltage, and the strong magnetic field induced by the excitation coil in the detection coil is nullified by a cancellation coil reverse-wound in series with the detection coil. The responses of the sensors in an array are separated in time by a linear gradient dc biasing coil. The sensors can be interrogated through metallic and nonmetallic barriers. The signals from the detection coil and the excitation coil are digitized by a pair of bipolar analog-to-digital converters (ADCs). A Raspberry Pi single-board computer (SBC) and associated software perform data acquisition and control all aspects of the sensor system hardware. The program allows the user to select the number of sensors in the array, the type of signal that is being collected, and the number of samples to take. The program also allows for signal processing of the sensor data, such as baseline correction. The program can differentiate sensor peaks from each other and calculate the magnitude of each sensor response with less than 1% error. The data are then displayed along with a graph of the signal.
This research program investigates boundary-layer receptivity in the leading-edge region for bodies with blunt leading edges. Receptivity theory provides the link between the unsteady distrubance environment in the free stream and the initial amplitudes of the instability waves in the boundary layer. This is a critical problem which must be addressed in order to develop more accurate prediction methods for boundary-layer transition. The first phase of this project examines the effects of leading-edge bluntness and aerodynamic loading for low Mach number flows. In the second phase of the project, the investigation is extended to supersonic Mach numbers. Singular perturbation techniques are utilized to develop an asymptotic theory for high Reynolds numbers. In the first year, the asymptotic theory was developed for leading-edge receptivity in low Mach number flows. The case of a parabolic nose is considered. Substantial progress was made on the Navier-Sotkes computations. Analytical solutions for the steady and unsteady potential flow fields were incorporated into the code, greatly expanding the types of free-stream disturbances that can be considered while also significantly reducing the the computational requirements. The time-stepping algorithm was modified so that the potential flow perturbations induced by the unsteady pressure field are directly introduced throughout the computational domain, avoiding an artificial 'numerical diffusion' of these from the outer boundary. In addition, the start-up process was modified by introducing the transient Stokes wave solution into the downstream boundary conditions.
Results of extensive studies on CFD algorithms for 2D inviscid flows in Cartesian and body fitted coordinates geometries are reviewed. These studies represent part of an ongoing investigation of combustion instabilities involving the interactions of high-speed nonlinear acoustic waves. Four numerical methods for the treatment of high speed flows are compared, namely, Roe-Sweby TVD, Yee symmetric TVD; Osher-Chakravarthy TVD; and the Colella's multi-dimensional Godunov method.
In the present study, the far wake was examined numerically using an implicit, upwind, finite-volume, compressible Navier-Stokes code. The numerical grid started at 500 equivalent circular cylinder diameters in the wave, and extended to 4000 equivalent diameters. By concentrating only on the far wake, the numerical difficulties and fine mesh requirements near the wake-generating body were eliminated. At the time of this writing, results for the K-epsilon and K-omega turbulence models at low Mach number have been completed and show excellent agreement with previous incompressible results and far-wake similarity solutions. The code is presently being used to compare the performance of various other turbulence models, including Reynolds stress models and the new anisotropic two-equation turbulence models being developed at NASA Langley. By increasing our physical understanding of the deficiencies and limits of these models, it is hoped that improvements to the universality of the models can be made. Future plans include examination of two-dimensional momentumless wakes as well.
The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.
XMbodyinfo was designed to evaluate potential reference trajectories, providing a proficient way to assess the quality of all satellite body flybys for a Cassini type mission tour. It is autonomous and will generate a variety of ORS (optical remote sensing) and FPW (fields, particles, and waves) plots that aid in the evaluation, qualification, selection, and improvement of a potential tour.
There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory (NREL) has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that can be neglected for wind turbines. These include buoyant loads, added mass and inertia loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbine systems at a speed comparable to real time. While efficient for long simulations, large sets of load cases, and design studies, mid-fidelity codes cannot capture all of the potentially important physical phenomenon impacting marine turbine systems. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system designed by NREL was modeled both with OpenFAST and with the commercial CFD code, STARCCM+. The CFD model used a 3-D unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modelled with a custom coupling of the CFD solver with the open-source code, MoorDyn. This improves upon the existing quasi-static catenary solver in STARCCM+, which lacks seabed contact or line-to-line connections. Spatial and temporal convergence studies were conducted. The simulation results for a combined current and wave condition are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.
Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here, in this paper, we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.
While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024
There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.
The neutron total cross sections σ tot of 16,18 O, 58 64 Ni, 103 Rh, and 112,124 Sn have been measured at the Los Alamos Neutron Science Center from low to intermediate energies (3 ≤ E lab ≤ 450 MeV) by leveraging wave-form-digitizer technology. The σ tot relative differences between isotopes are presented, revealing additional information about the isovector components needed for an accurate optical-model description away from stability. Digitizer-enabled σ tot -measurement techniques are discussed and a series of uncertainty-quantified dispersive optical model (DOM) analyses using these new data is presented, validating the use of the DOM for modeling light systems ( 16,18 O) and systems with open neutron shells ( 58,64 Ni and 112,124 Sn). Here, the valence-nucleon spectroscopic factors extracted for each isotope reaffirm the usefulness of high-energy proton reaction cross sections for characterizing depletion from the mean-field expectation.
Abstract We employ scattering amplitudes in curved space to model the dynamics of a light probe particle with massmorbiting in the background spacetime induced by a heavy gravitational source with massM. Observables are organized as an expansion inm/Mto all orders inG— the gravitational self-force expansion. An essential component of our analysis is the backreaction of the heavy source which we capture by including the associated light degrees of freedom. As illustration we consider a Schwarzschild background and verify geodesic motion as well as the first-order self-force correction to two-body scattering through$$\mathcal{O}$$(G 3 ). Amplitudes in curved space offer several advantages, and further developments along these lines may advance the computation of gravitational-wave signals for extreme-mass-ratio inspirals.
An optical lattice is a periodic light crystal constructed from the standing-wave interference patterns of laser beams. It can be used to store and manipulate quantum degenerate atoms and is an ideal platform for the quantum simulation of many-body physics. A principal feature is that optical lattices are flexible and possess a variety of multidimensional geometries with modifiable band structure. An even richer landscape emerges when control functions can be applied to the lattice by modulation of the position or amplitude with Floquet driving. However, the desire of realizing high-modulation bandwidths while preserving extreme lattice stability has been difficult to achieve. In this paper, we demonstrate an effective solution that consists of overlapping two counterpropagating lattice beams and controlling the phase and intensity of each with independent acousto-optic modulators. Our phase controller mixes sampled light from both lattice beams with a common optical reference. This dual heterodyne locking method allows exquisite determination of the lattice position, while also removing parasitic phase noise accrued as the beams travel along separate paths. Here, we report up to 50 dB suppression in lattice phase noise in the 0.1–1 Hz band, along with significant suppression spanning more than 4 decades of frequency. The absolute phase diffusion of the lattice position is only 10 Å when integrated over 10 s. This method permits precise, high-bandwidth modulation (greater than 50 kHz) of the optical lattice intensity and phase. We demonstrate the efficacy of this approach by executing intricate time-varying phase profiles for atom interferometry.
The full configuration interaction (FCI) wave function provides the exact solution to the Schrödinger equation in a given basis set. While FCI is intractable due to exponential computational costs, the many-body expansion (or the method of increments) can reduce scaling to a low-order polynomial with system size. This project entailed advances in the incremental FCI (iFCI) approach, designed to allow iFCI to reach larger system sizes and maintain its intrinsic high accuracy. This document describes advances in solvers for iFCI, strategies to treat multiple charge and spin states, and virtual state management methods to reduce memory requirements. Overall, this project allows iFCI to correlate (for the first time) 142 valence electrons in 444 orbitals in a realistic model of a transition metal complex.
Two different versions of the Green's function for the scalar wave equation in weakly curved spacetime (one due to DeWitt and DeWitt, the other to Thorne and Kovacs) are compared and contrasted; and their mathematical equivalence is demonstrated. Then the DeWitt-DeWitt Green's function is used to construct several alternative versions of the Thorne-Kovacs postlinear formalism for gravitational-wave generation. Finally it is shown that, in calculations of gravitational bremsstrahlung radiation, some of our versions of the postlinear formalism allow one to treat the interacting bodies as point masses, while others do not.