Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Helmholtz equation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The development of thermodynamically consistent and physics-informed equation-of-state model through machine learning

Ab initio molecular dynamics (AIMD) simulations have become an important tool used in the construction of equations of state (EOS) tables for warm dense matter. Due to computational costs, only a limited number of system state conditions can be simulated, and the remaining EOS surface must be interpolated for use in radiation-hydrodynamic simulations of experiments. In this work, we develop a thermodynamically consistent EOS model that utilizes a physics-informed machine learning approach to implicitly learn the underlying Helmholtz free-energy from AIMD generated energies and pressures. The model, referred to as PIML-EOS, was trained and tested on warm dense polystyrene producing a fit within a 1% relative error for both energy and pressure and is shown to satisfy both the Maxwell and Gibbs–Duhem relations. In addition, we provide a path toward obtaining thermodynamic quantities, such as the total entropy and chemical potential (containing both ionic and electronic contributions), which are not available from current AIMD simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Scalar mixing in a Kelvin-Helmholtz shear layer and implications for Reynolds-averaged Navier-Stokes modeling of mixing layers

Large-eddy simulation of a temporally evolving Kelvin-Helmholtz (KH) mixing layer is performed with the tenth-order compact difference code miranda to examine the steady-state behavior of a passive scalar in a shear-driven mixing layer. It is shown that the integral behavior of scalar variance in a KH mixing layer behaves similarly to the integral behavior of scalar variance in a Rayleigh-Taylor (RT) mixing layer, and mixedness of the simulated KH shear layer tends towards a value of about 0.8. It is further shown that if the k-L-a-V Reynolds-averaged Navier-Stokes (RANS) model [B. E. Morgan et al., Phys. Rev. E 98, 033111 (2018)], calibrated to reproduce steady-state mixing in an RT layer, is applied to simulate a KH mixing layer, the RANS model will significantly overpredict the magnitude of scalar variance in the KH layer. A straightforward addition to the k-L-a-V model is then suggested, and self-similarity analysis is applied to determine constraints on model coefficients. Furthermore, it is shown that with the addition of a buoyancy production term in the model equation for scalar variance, it becomes possible to eliminate the model deficiency and match steady-state mixedness in simulations of both RT and KH mixing layers with a single model calibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Incompressible magnetohydrodynamic surface waves - Nonlinear aspects

The nonlinear properties of MHD surface waves in the solar atmosphere are investigated analytically, assuming that the fluid is incompressible and that the waves are confined to a single surface, with semiinfinite regions on both sides. The governing equations are derived in detail, and qualitative results are presented in a graph. For propagating waves, second-order terms in the wave amplitude are found to lead to wave steepening at leading or trailing edges, the steepening rate becoming very large as the threshold for the linear Kelvin-Helmholtz instability is approached. Second-order effects on standing waves include crest and trough sharpening (increasing with time), a current independent of distance on the surface but decreasing exponentially with distance from the surface, and pressure-field fluctuations of infinite extent. It is suggested that these effects could account for a large fraction of solar-atmosphere heating.

Hollweg, Joseph V.↗
NASA NTRSDate not supplied

TPSAS-NF1676L-33433-DND

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction, and are characterized by a complex valued frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic scattering of geometric bodies treated with acoustic liners. This work considers using either an impedance boundary condition or admittance boundary condition to allow acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces; admittance is defined to be the inverse of impedance. Three different acoustic liner models will be discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and a Broadband Model. In the Extended Helmholtz Resonator Model and Three-Parameter Impedance Model, the liner impedance is specified at a single frequency; the Broadband Model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for each model will be derived and coupled with a time-domain boundary integral equation. The scattering solution will be obtained iteratively using spatial and temporal basis functions. Stability will be demonstrated through eigenvalue analysis. Moreover, fast algorithms and high performance computing are considered to both assess and reduce the computational cost of the simulation.

Michelle E Rodio↗
NASA NTRS1984

On the stability of the ionopause of Venus

The stability of the Venus ionopause is examined in light of the importance of gravitation and curvature. Using a one-fluid approximation for the equation of motion of the plasma, and ignoring the effects of neutrals, a dispersion relation is obtained that includes the effects of the magnetic field, sheared plasma flow, buoyancy, centrifugal force and magnetic tension due to boundary curvature. It is found that buoyancy acts to neutralize the flute instability. As expected, the Kelvin-Helmholtz mode is the dominant instability over most of the dayside ionopause. The expected growth times of this mode are short in comparison with the wave-convection time over the boundary; the waves can grow and saturate quickly, producing a turbulent boundary that may affect electrodynamic coupling between the solar wind and ionospheric plasmas.

Elphic, R. C.↗
NASA NTRS1994

The effect of parallel currents on auroral micropulsations

Field aligned currents play an important role in the global coupling between the magnetosphere and the ionosphere and in their relationship to the auroral phenomena. Moreover, there exists evidence that resonant oscillations are related to large-scale Birkeland currents. The spatial confinement of the field-aligned currents forms an inhomogeneous system susceptible to low-frequency oscillations, which can be excited due to periodic variations in the solar wind pressure or to the Kelvin-Helmholtz (KH) instability. In this paper we present a study of ultralow-frequency (ULF) oscillations in an inhomogeneous magnetic field formed by a large-scale current. We investigate the effects of the field-aligned currents on the generation of localized Alfven waves. The field oscillations are described by an eigenvalue wave equation which includes the effects of the field aligned currents and which produces a discrete spectrum of Alfven waves. These waves are observed mainly in three regions of the magnetosphere: in the magnetosheath, in the polar cusp, and in the plasmasphere. In the present study we limit our investigation to the auroral region.

Tavares, M.↗
NASA NTRS1977

Fluid mechanical model of the Helmholtz resonator

A semi-empirical fluid mechanical model of the acoustic behavior of Helmholtz resonators is presented which predicts impedance as a function of the amplitude and frequency of the incident sound pressure field and resonator geometry. The model assumes that the particle velocity approaches the orifice in a spherical manner. The incident and cavity sound fields are connected by solving the governing oscillating mass and momentum conservation equations. The model is in agreement with the Rayleigh slug-mass model at low values of incident sound pressure level. At high values, resistance is predicted to be independent of frequency, proportional to the square root of the amplitude of the incident sound pressure field, and virtually independent of resonator geometry. Reactance is predicted to depend in a very complicated way upon resonator geometry, incident sound pressure level, and frequency. Nondimensional parameters are defined that divide resonator impedance into three categories corresponding to low, moderately low, and intense incident sound pressure amplitudes. The two-microphone method was used to measure the impedance of a variety of resonators. The data were used to refine and verify the model.

Hersh, A. S.↗
NASA NTRS1992

Computing interface motion in compressible gas dynamics

An analysis is conducted of the coupling of Osher and Sethian's (1988) 'Hamilton-Jacobi' level set formulation of the equations of motion for propagating interfaces to a system of conservation laws for compressible gas dynamics, giving attention to both the conservative and nonconservative differencing of the level set function. The capabilities of the method are illustrated in view of the results of numerical convergence studies of the compressible Rayleigh-Taylor and Kelvin-Helmholtz instabilities for air-air and air-helium boundaries.

Mulder, W.↗
NASA NTRSDate not supplied

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstructs the acoustic signal at each microphone in the array. The mode structure of the lowest-order mode is shown to be significantly impacted by shear flow refraction effects, and higher-order mode structures are also affected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest-order mode, but discrepancies arise when higher-order modes are the dominant component of the acoustic field.

Duct acoustics↗
NASA NTRSDate not supplied

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstruct the acoustic signals at each microphone in the array. The mode structure of the lowest order mode is shown to be significantly impacted by shear flow refraction effects, and higher order mode structures are also effected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest order mode, but discrepancies arise when higher order modes are the dominant component of the acoustic field.

Duct acoustics↗
NASA NTRSDate not supplied

Helmholtz Motor: A Novel Electric Machine That Enables High Temperature Superconducting Armatures

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Second generation high temperature superconductors (HTS) are the ideal conductor for the field winding (typically the rotor) of superconducting machines, and they have attractive features for the armature winding (typically the stator) of these machines. However, their use in armature windings has been severely limited due to excessive AC losses in conventional electric machines. This presentation will present the working principle, design, and a performance trade study of a novel electric machine – termed a Helmholtz machine – that enables HTS armatures by significantly reducing AC losses. The working principle of the patent pending Helmholtz motor is to optimize the rotor-produced magnetic field (i.e., the field winding) so that it produces a magnetic field almost purely in the plane of the HTS. This is accomplished using 2 rotors that contain a set of matched Helmholtz coil pairs (or matched permanent magnet Halbach arrays). The armature (stator) HTS coils are positioned between the rotors and oriented to minimize the out of plane magnetic field. The armature HTS coils are also designed to minimize the out of plane component of their self field. The presentation will include the lessons learned in the sizing and design of this motor, including equations for how to match the magnet arrays. A performance trade study will be presented for a partially superconducting, radial flux version of the motor. The trade study involved an optimization of a motor design code based on analytical calculations (electromagnetic, AC loss, thermofluid, mechanical, and motor sizing) that determine the motor’s total efficiency and total specific power. The design code will be described in the presentation, with an emphasis on the calculation of magnetic fields, AC loss, and temperature distribution.

Superconducting electric machines↗
NASA NTRS2000

Nonlinear Stability and Structure of Compressible Reacting Mixing Layers

The parabolized stability equations (PSE) are used to investigate issues of nonlinear flow development and mixing in compressible reacting shear layers. Particular interest is placed on investigating the change in flow structure that occurs when compressibility and heat release are added to the flow. These conditions allow the 'outer' instability modes- one associated with each of the fast and slow streams-to dominate over the 'central', Kelvin-Helmholtz mode that unaccompanied in incompressible nonreacting mixing layers. Analysis of scalar probability density functions in flows with dominant outer modes demonstrates the ineffective, one-sided nature of mixing that accompany these flow structures. Colayer conditions, where two modes have equal growth rate and the mixing layer is formed by two sets of vortices, offer some opportunity for mixing enhancement. Their extent, however, is found to be limited in the mixing layer's parameter space. Extensive validation of the PSE technique also provides a unique perspective on central- mode vortex pairing, further supporting the view that pairing is primarily governed perspective sheds insight on how linear stability theory is able to provide such an accurate prediction of experimentally-observed, fully nonlinear flow phenomenon.

Day, M. J.↗
NASA NTRS2005

Numerical Simulation of Liquid Jet Atomization Including Turbulence Effects

This paper describes numerical implementation of a newly developed hybrid model, T-blob/T-TAB, into an existing computational fluid dynamics (CFD) program for primary and secondary breakup simulation of liquid jet atomization. This model extend two widely used models, the Kelvin-Helmholtz (KH) instability of Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O'Rourke and Amsden to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. Several assessment studies are presented and the results indicate that the existing KH and TAB models tend to under-predict the product drop size and spray angle, while the current model provides superior results when compared with the measured data.

Trinh, Huu P.↗
DOE OSTI2025

Femtosecond temperature measurements of laser-shocked copper deduced from the intensity of the x-ray thermal diffuse scattering

We present 50-fs, single-shot measurements of the x-ray thermal diffuse scattering (TDS) from copper foils that have been shocked via nanosecond laser ablation up to pressures above ∼135 GPa. We hence deduce the x-ray Debye–Waller factor, providing a temperature measurement. The targets were laser-shocked with the DiPOLE 100-X laser at the High Energy Density endstation of the European X-ray Free-Electron Laser. Single x-ray pulses, with a photon energy of 18 keV, were scattered from the samples and recorded on Varex detectors. Despite the targets being highly textured (as evinced by large variations in the elastic scattering) and with such texture changing upon compression, the absolute intensity of the azimuthally averaged inelastic TDS between the Bragg peaks is largely insensitive to these changes, and allowing for both Compton scattering and the low-level scattering from a sacrificial ablator layer provides a reliable measurement of $T/Θ^2_D$, where Θ D is the Debye temperature. We compare our results with the predictions of the SESAME 3336 and LEOS 290 equations of state for copper and find good agreement within experimental errors. We, thus, demonstrate that single-shot temperature measurements of dynamically compressed materials can be made via thermal diffuse scattering of XFEL radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

A self-consistent linear-mode model of stellar convection

A normal-mode expansion of the linearized fluid equations in terms of small subset of spherical harmonics can provide a foundation for a physically motivated, self-consistent description of a solar-type convection zone. In the absence of dissipation, a second-order differential equation governs the radial dependence of the modes, so that interpretation of the effects on convection quantities of the normal-form 'potential well' is straightforward. The philosophy is quite different from the more recent work of Narasimha and Antia (1982): all envelopes presented here differ substantially from MLT envelopes, and therefore, from theirs, which are constructed to be consistent with MLT. The amplitude of all modes is set by a Kelvin-Helmholtz-('shear'-) instability argument unrelated to solar observations, with the result that the convection description may be considered to arise from 'first-hueristic-principles'. The thermodynamics modelled vaguely resemble the sun's, and more vigorously convective envelopes show some phenomena qualitatively like solar observations (e.g., atmospheric velocity spectra).

Macauslan, J.↗
NASA NTRS1993

Mach wave emission from a high temperature supersonic jet

The paper considers the compressible Rayleigh equation as a model for the Mach wave emission mechanism associated with high temperature supersonic jets. Solutions to the compressible Rayleigh equation reveal the existence of several families of supersonically convecting instability waves. These waves directly radiate noise to the jet far field. The predicted noise characteristics are compared to previously acquired experimental data for an axisymmetric Mach 2 fully pressure balanced jet operating over a range of jet operating total temperatures from ambient to 1370 K. The results of this comparison show that the first order supersonic instability wave and the Kelvin-Helmholtz first, second, and third order modes have directional radiation characteristics that are in agreement with observed data. The assumption of equal initial amplitudes for all of the waves leads to the conclusion that the flapping mode of instability dominates the noise radiation process of supersonic jets. At a jet temperature of 1370 K, supersonic instability waves are predicted to dominate the noise radiated at high frequency at narrow angles to the jet axis.

Seiner, J. M.↗
NASA NTRS1994

Parasitic instabilities in magnetized, differentially rotating disks

Velikhov, Chandrasekhar, Balbus & Hawley have discovered a MHD instability of differentially rotating fluids that may explain the effective viscosity of accretion disks. If the unperturbed magnetic energy density is much less than the gas pressure, this `magnetorotational' instability (MRI) arises on small scales and is approximately incompressible. If the unperturbed field is suffiently, strong, the MRI is suppressed. We therefore ask whether the MRI mechanism can amplify even very weak initial fields until they reach equipartition with the gas pressure. We show that in the total incompressible limit, the MRI modes are exact solutions of the nonlinear fluid equations, even if the perturbed magnetic field is much larger than the unperturbed field. Also, we present a new exact solution in the opposite limit that the magnetic pressure is much larger than that of the gas. On the other hand, we show that the incompressible MRI modes are themselves subject to parasitic instabilities with instantaneous growth rates proportional to the MRI amplitude. Some of the parasitic instabilities are of the Kelvin-Helmholtz type, while others are less familiar. The eigenfunctions of the latter group suggest that they may promote rapid reconnection of the MRI field. Thus, parasitic modes may stop MRI growth at subequipartition amplitudes if it developes from a sufficiently weak initital field.

Goodman, Jeremy↗
NASA NTRSDate not supplied

Flow field Reconstruction for Inhomogeneous Turbulence using Data and Physics Driven Models

A methodology combining Large Eddy Simulation (LES) trained data and a physics driven wave packet model to obtain a reduced order reconstruction for broadband, three-dimensional, temporally stationary but spatially inhomogeneous, incompressible turbulence. Wake turbulence generated by an axisymmetric dragging disk with a turbulent co-flow serves as the benchmark test case. We begin by studying the proper-orthogonal decomposition of the turbulent fluctuations taken from a high-resolution LES to first identify whether the fields demonstrate a low-rank character. It is argued that the presence of the turbulent co-flow results in a largely broadband character lacking any tonal properties. This is especially true for Strouhal numbers greater than 1 and only a small fraction of energy is contained in the leading order Kelvin-Helmholtz modes. As such reconstructions and reduced order modeling purely relying on data from LES does not appear to be a lucrative solution - contrary to problems with strongly tonal character. To supplement the missing energy from a low order truncated mode expansion, we utilize a physics based super-resolution (enrichment) algorithm that relies on spatio-temporally localized Gabor wave packets whose time evolution is described using a set of ordinary differential equations. The reconstructed flow has single- and two-point correlations that are consistent with the reference high resolution simulation data.

SLS↗