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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Fluid hydrogen at high density - Pressure dissociation

A model for the Helmholtz free energy of fluid hydrogen at high density and high temperature is developed. This model aims at describing both pressure and temperature dissociation and ionization and bears directly on equations of state of partially ionized plasmas, as encountered in astrophysical situations and high-pressure experiments. This paper focuses on a mixture of hydrogen atoms and molecules and is devoted to the study of the phenomenon of pressure dissociation at finite temperatures. In the present model, the strong interactions are described with realistic potentials and are computed with a modified Weeks-Chandler-Andersen fluid perturbation theory that reproduces Monte Carlo simulations to better than 3 percent. Theoretical Hugoniot curves derived from the model are in excellent agreement with experimental data.

Saumon, Didier↗
DOE OSTI2024

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.

42 ENGINEERING↗
NASA NTRS1993

Nonlinear electrostatic drift Kelvin-Helmholtz instability

Nonlinear analysis of electrostatic drift Kelvin-Helmholtz instability is performed. It is shown that the analysis leads to the propagation of the weakly nonlinear dispersive waves, and the nonlinear behavior is governed by the nonlinear Burger's equation.

Sharma, Avadhesh C.↗
NASA NTRS1996

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗
DOE OSTI2020

Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon

The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Stacking Independence and Resonant Interlayer Excitation of Monolayer WSe 2 /MoSe 2 Heterostructures for Photocatalytic Energy Conversion

Here, we report a comparison of the photocatalytic performance of WSe 2 -on-MoSe 2 and MoSe 2 -on-WSe 2 heterostructures. While built-in electric fields exist in these heterostructures on the order of 100 kV/cm due to band offsets between these two materials, the photocatalytic performance (i.e., photocurrent) is independent of the stacking order of the two materials. Solving Poisson’s equation under these conditions, we find that the built-in electric field produced in the heterostructure is at least 1 order of magnitude smaller than that produced in the electrochemical double layer (i.e., Helmholtz layer). Mott–Schottky measurements indicate that transition metal dichalcogenides (TMDCs) on ITO electrodes have similar capacitance to that of bare ITO, providing further evidence that the interfacial electric fields produced in the solid state heterostructure are negligible compared to the fields generated by the ions in solution. The photocatalytic performance of these heterostructures provided the largest relative enhancement in the heterojunction region under 920 and 785 nm irradiation compared with 532 and 633 nm wavelength excitation. Here, the 920 nm (1.35 eV) photons lie below the band gaps and produce very little photocurrent in the constituent monolayer materials but resonantly excite the interlayer optical transition in the heterostructure, producing a 5-fold enhancement in the measured photocurrent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Finite Element Prediction of Acoustic Scattering and Radiation from Submerged Elastic Structures

A finite element formulation is derived for the scattering and radiation of acoustic waves from submerged elastic structures. The formulation uses as fundamental unknowns the displacement in the structure and a velocity potential in the field. Symmetric coefficient matrices result. The outer boundary of the fluid region is terminated with an approximate local wave-absorbing boundary condition which assumes that outgoing waves are locally planar. The finite element model is capable of predicting only the near-field acoustic pressures. Far-field sound pressure levels may be determined by integrating the surface pressures and velocities over the wet boundary of the structure using the Helmholtz integral. Comparison of finite element results with analytic results show excellent agreement. The coupled fluid-structure problem may be solved with general purpose finite element codes by using an analogy between the equations of elasticity and the wave equation of linear acoustics.

Everstine, G. C.↗
DOE OSTI2025

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

A transverse Kelvin-Helmholtz instability in a magnetized plasma

An analysis is conducted of the transverse Kelvin-Helmholtz instability in a magnetized plasma for unstable flute modes. The analysis makes use of a two-fluid model. Details regarding the instability calculation are discussed, taking into account the ion continuity and momentum equations, the solution of a zero-order and a first-order component, and the properties of the solution. It is expected that the linear calculation conducted will apply to situations in which the plasma has experienced no more than a few growth periods.

Kintner, P.↗
NASA NTRS1989

The nonlinear breakup of the sun's toroidal field

There are good reasons for believing that the sun has a strong toroidal magnetic field in the stably stratified region of convective overshoot sandwiched between the radiative zone and convective zone proper. The magnetic field in this region is modeled by studying the behavior of a layer of uniform field embedded in a subadiabatic atmosphere. Since the field can support extra mass, such a configuration is top-heavy, and instabilities of the Rayleigh-Taylor type can occur. Numerical integration of the two-dimensional compressible MHD equations makes it possible to follow the evolution of this instability into the nonlinear regime. The initial buoyancy-driven instability of the magnetic field gives rise to strong shearing motions, thereby exciting secondary Kelvin-Helmholtz instabilities which wrap the gas into regions of intense vorticity. The somewhat surprising subsequent motions are determined primarily by the strong interactions between vortices.

Hughes, D. W.↗
NASA NTRS2007

3D Relativistic MHD Simulations of Magnetized Spine-Sheath Relativistic Jets

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using the RAISHIN code. In the numerical simulations a jet with Lorentz factor gamma=2.5 is precessed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode-analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized moving sheath. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

Mizuno, Yosuke↗
NASA NTRS2008

3D RMHD Simulations of Magnetized Spine-Sheath Relativistic Jets

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath flow using the RAISHIN code. In the numerical simulations a jet with Lorentz factor gamma=2.5 is processed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized sheath flow. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

Mizuno, Yosuke↗
NASA NTRS2008

3D RMHD Simulations of Magnetized Spine-sheath Relativistic Jets

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic'(0.5c) sheath flow using the RAISHIN code. In the numerical simulations a jet with Lorentz factor gamma=2.5 is processed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized sheath flow. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

Mizuno, Yosuke↗
NASA NTRS1994

Origin of flux ropes in Venus' ionosphere

The joule dissipation inside flux ropes in Venus' ionosphere is so great that they must be formed near, and maintained at, the place where they are observed. Thus ropes are not formed by a Kelvin-Helmholtz instability of the ionopause. The hypothesis that ropes may be formed by the dynamo action of internal gravity waves in Venus' thermosphere (Luhmann and Elphic, 1985; Cole, 1993) is strengthened by discussion of a magnetic evolution equation which includes neutral air motion. However, the dynamo process would work only at altitudes at which v(sub in) is greater than or equal to omega(sub i). At altitudes or parts of a rope where v(sub in) is much less than omega(sub i), the process does not work. A solar wind dynamo is therefore examined to account for the ropes. Thereby a major new heat source for ions of the Venus ionosphere associated with the ropes is uncovered.

Cole, Keith D.↗
NASA NTRS1984

On the Kelvin-Helmholtz instabilities of supersonic shear layers

Analytical consideration is given to a Kelvin-Helmholtz instability in a finite width supersonic shear layer with linear and sinusoidal velocity profiles in the transition region. An analogy is explored between the pressure perturbation and the bound and scattering wave states of the one-dimensional Schroedinger equation. A numerical model is obtained for generating contour maps of the real and imaginary parts of the perturbation frequency in the (A,B)-plane. A comparison is made between the numerical results and analytical expressions derived for very small Mach numbers. Two unstable regimes are characterized in the (A,B)-plane. The method is concluded useful for both laboratory and astrophysics problems.

Choudhury, S. R.↗
DOE OSTI2020

A discontinuous Galerkin method for general relativistic hydrodynamics in thornado

Discontinuous Galerkin (DG) methods provide a means to obtain high-order accurate solutions in regions of smooth fluid flow while still resolving strong shocks. These and other properties make DG methods attractive for solving problems involving hydrodynamics; e.g., the core-collapse supernova problem. With that in mind we are developing a DG solver for the general relativistic, ideal hydrodynamics equations under a 3+1 decomposition of spacetime, assuming a conformally-flat approximation to general relativity. With the aid of limiters we verify the accuracy and robustness of our code with several difficult test-problems: a special relativistic Kelvin-Helmholtz instability problem, a two-dimensional special relativistic Riemann problem, and a one- and two-dimensional general relativistic standing accretion shock (SAS) problem. We find good agreement with published results, where available. We also establish sufficient resolution for the 1D SAS problem and find encouraging results regarding the standing accretion shock instability (SASI) in 2D.

Dunham, Samuel↗
NASA NTRS1978

Impedance measurement using a two-microphone, random-excitation method

The feasibility of using a two-microphone, random-excitation technique for the measurement of acoustic impedance was studied. Equations were developed, including the effect of mean flow, which show that acoustic impedance is related to the pressure ratio and phase difference between two points in a duct carrying plane waves only. The impedances of a honeycomb ceramic specimen and a Helmholtz resonator were measured and compared with impedances obtained using the conventional standing-wave method. Agreement between the two methods was generally good. A sensitivity analysis was performed to pinpoint possible error sources and recommendations were made for future study. The two-microphone approach evaluated in this study appears to have some advantages over other impedance measuring techniques.

Seybert, A. F.↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗