Search NASA⌕ Search

SEARCH · Search NASA

Results for “WAVE 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An Investigation into High Voltage Spiral Generators Utilizing Thyristor Input Switches

High voltage nanosecond pulses are widely used in scientific research, but their wider adoption in industry requires compact, cost effective and easy to use generators to be developed. This paper presents modelling and experimental investigations into one method of producing such pulses – a spiral generator with a solid-state, thyristor-based input switch. It includes how the pulses are formed within the spiral, why a high speed input switch is required, how the geometry of the spiral dictates its output characteristics, and the effects of different loads. Using thyristors, often connected in series to increase the operating voltage of the spiral, enables the spiral generators to have low jitter, high repetition rate, and long lifetime. Modelling of the circuit used a combination of telegraph equations to account for the wave propagation along the spiral and a lumped circuit exchanging charge between the spiral, the input switch, and the load. The model is verified by detailed experimental results with the relative error being < 10% in most cases. Here, the output voltage pulse was often observed to have an initial peak of much lower magnitude than the subsequent peak(s) – which can only be fully explained by considering wave propagation effects. Lower input switch inductance, shorter switching time, larger mean diameter of the spiral, and increasing the width of the copper tape that makes up the spiral can all increase the voltage multiplication efficiency. Though increasing the number of turns that make up the spiral can increase the output voltage, it can also lower the multiplication efficiency. By understanding the effects of different geometries the spiral can be optimized to drive different loads - three applications of such spiral generators are then presented: pulses with 10 kV amplitude and 10 kHz repetition rate for driving DBD plasma, pulses with amplitude of 10 kV and 10 kV/ns rising rate for triggering of advanced solid state switches, and pulses with -50 kV amplitude and 50 ns rising time for triggering high current gas switches through field distortion.1 Index Terms—spiral generator, nanosecond pulses, telegraph equations, series-connected thyristor module, pulse forming process

42 ENGINEERING↗
DOE OSTI2025

Fast meta-solvers for 3D complex-shape scatterers using neural operators trained on a non-scattering problem

Three-dimensional target identification using scattering techniques requires high accuracy solutions and very fast computations for real-time predictions in some critical applications. We first train a deep neural operator (DeepONet) to solve wave propagation problems described by the Helmholtz equation in a domain without scatterers but at different wavenumbers and with a complex absorbing boundary condition. We then design two classes of fast meta-solvers by combining DeepONet with either relaxation methods, such as Jacobi and Gauss-Seidel, or with Krylov methods, such as GMRES and BiCGStab, using the trunk basis of DeepONet as a coarse-scale preconditioner. We leverage the spectral bias of neural networks to account for the lower part of the spectrum in the error distribution while the upper part is handled inexpensively using relaxation methods or fine-scale preconditioners. The meta-solvers are then applied to solve scattering problems with different shape of scatterers, at no extra training cost. We first demonstrate that the resulting meta-solvers are shape-agnostic, fast, and robust, whereas the standard standalone solvers may even fail to converge without the DeepONet. We then apply both classes of meta-solvers to scattering from a submarine, a complex three-dimensional problem. We achieve very fast solutions, especially with the DeepONet-Krylov methods, which require orders of magnitude fewer iterations than any of the standalone solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A second-order distributed memory parallel fast sweeping method for the Eikonal equation

The Eikonal equation is used to calculate wave propagation and distance fields, and due to its complexity requires numerical treatment for its solution. In this work, we present a second-order distributed memory parallel fast sweeping method. The second-order solution switches on a two-point stencil when two upwind points are available, and reverts to first-order otherwise. In all examples, the second-order method improves the solution over the first-order, allowing for significant savings in memory while achieving the same accuracy. Parallelization over distributed memory saw good weak scaling with optimal convergence. The computational time for second-order was approximately 2.5 times slower than first-order, where the largest amount of mesh points ran on 144 cores (512 GB) was ≈20 billion. The savings in memory from the second-order method combined with the distributed memory algorithm result in the ability to solve problems much larger than are possible with the serial first-order method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Role of Phytoplankton Biomacromolecules in Controlling Ocean Surface Roughness

Satellite altimetric data routinely map sea surface topography by measuring the ocean return signal. One source of altimeter measurement contamination occurs when the radar ocean backscatter becomes unusually large, a situation termed a Sigma-0 bloom. Past research suggests Sigma-0 blooms are associated with weak wind and natural surface slick conditions where capillary waves at the air–sea interface are suppressed. To date, no explicit connection between these conditions and Sigma-0 bloom presence has been provided. Using a series of simplified equations, our reduced model determines capillary wave heights from estimates of planktonic carbon concentrations and regional wind speed. Our results suggest that the radar signal reflection increases as capillary wave height decreases. This relationship depends on surfactant concentration, surfactant composition, and wind speed. Model sensitivity analysis indicates that the interface reflectivity depends on biological activity and wind speed. Our proposed simplified model provides a method to identify potential Sigma-0 bloom regions. We conclude that because of the demonstrated impact of biological surfactants on ocean roughness, it is necessary to consider the biological activity, i.e., phytoplankton bloom events, when interpreting signals from radar altimetry and when developing ocean hydrology models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On exact-WKB analysis, resurgent structure, and quantization conditions

There are two well-known approaches to studying nonperturbative aspects of quantum mechanical systems: saddle point analysis of the partition functions in Euclidean path integral formulation and the exact-WKB analysis based on the wave functions in the Schrödinger equation. In this work, based on the quantization conditions obtained from the exact-WKB method, we determine the relations between the two formalism and in particular show how the two Stokes phenomena are connected to each other: the Stokes phenomenon leading to the ambiguous contribution of different sectors of the path integral formulation corresponds to the change of the “topology” of the Stoke curves in the exact-WKB analysis. We also clarify the equivalence of different quantization conditions including Bohr-Sommerfeld, path integral and Gutzwiller’s ones. In particular, by reorganizing the exact quantization condition, we improve Gutzwiller’s analysis in a crucial way by bion contributions (incorporating complex periodic paths) and turn it into an exact result. Furthermore, we argue the novel meaning of quasi-moduli integral and provide a relation between the Maslov index and the intersection number of Lefschetz thimbles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Modeling and optimization of steady flow of natural gas and hydrogen mixtures in pipeline networks

Here, we extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity.

08 HYDROGEN↗
DOE OSTI2021

Photoinduced anomalous Hall effect in two-dimensional transition metal dichalcogenides

A circularly polarized ac pump field illuminated near resonance on two-dimensional transition metal dichalcogenides (TMDs) produces an anomalous Hall effect in response to a dc bias field. In this work, we develop a theory for this photoinduced anomalous Hall effect in undoped TMDs irradiated by a strong coherent laser field. The strong field renormalizes the equilibrium bands and opens up a dynamical energy gap where single-photon resonance occurs. The resulting photon dressed states, or Floquet states, are treated within the rotating-wave approximation. A quantum kinetic equation approach is developed to study the nonequilibrium density matrix and time-averaged transport currents under the simultaneous influence of the strong ac pump field and the weak dc probe field. Dissipative effects are taken into account in the kinetic equation that captures relaxation and dephasing. The photoinduced longitudinal and Hall conductivities display notable resonant signatures when the pump field frequency reaches the spin-split interband transition energies. Rather than valley polarization, we find that the anomalous Hall current is mainly driven by the intraband response of photon-dressed electron populations near the dynamical gap at both valleys, accompanied by a smaller contribution due to the interband response. In conclusion, these findings highlight the importance of photon-dressed bands and nonequilibrium distribution functions in achieving a proper understanding of the photoinduced anomalous Hall effect in a strong pump field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Total energy in supernova neutrinos and the tidal deformability and binding energy of neutron stars

The energy radiated in supernova neutrinos is a fundamental quantity that is closely related to the gravitational binding energy of a neutron star. Recently the tidal deformability of neutron stars was constrained by gravitational wave observations. By considering several equations of state, we find a strong correlation between the tidal deformability and neutron star binding energy. We use this correlation to sharpen predictions of the binding energy of neutron stars and the total neutrino energy in supernovae. We find a minimum binding energy for a neutron star formed in a supernova of ~ 1.5 × 10 53 ergs . Should the neutrino energy in a supernova be significantly below this value, it would strongly suggest new unobserved particles are carrying away some of the supernova energy. Alternatively, if the neutrino energy is observed above ~ 6 × 10 53 ergs , it would strongly imply the formation of a (perhaps surprisingly) massive neutron star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Low Mach number fluctuating hydrodynamics model for ionic liquids

We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

42 ENGINEERING↗
DOE OSTI2023

Exploring regular and turbulent flow states in active nematic channel flow via Exact Coherent Structures and their invariant manifolds

This work is a unified study of stable and unstable steady states of 2D active nematic channel flow using the framework of Exact Coherent Structures (ECS). ECS are stationary, periodic, quasiperiodic, or traveling wave solutions of the governing equations that, together with their invariant manifolds, organize the dynamics of nonlinear continuum systems. We extend our earlier work on ECS in the preturbulent regime by performing a comprehensive study of stable and unstable ECS for a wide range of activity values spanning the preturbulent and turbulent regimes. In the weakly turbulent regime, we compute more than 200 unstable ECS that co-exist at a single set of parameters, and uncover the role of symmetries in organizing the phase space geometry. We provide conclusive numerical evidence that in the preturbulent regime, generic trajectories shadow a series of unstable ECS before settling onto an attractor. Lastly, our studies hint at shadowing of quasiperiodic type ECS in the turbulent regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic “active turbulence.” Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Overall, our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Array directivity enhancement by leveraging angle-dependent scattering

The quality of a sonar array's localization capabilities, often expressed as directivity, is limited by the sonar's aperture, that is, the length of the sonar array. Previous attempts to improve directivity, without increasing array size, have been moderately successful. Wave scattering within a nontraditional array, such as an array fabricated from a non-homogeneous material, could provide additional information to the localization calculations and improve array directivity without increasing the size of the array. An investigation of array directivity improvement through wave scattering is performed. This paper modifies existing localization and directivity calculations to consider the scattered waves and uses the derived equations to explain why previous proposed scattering was incapable of increasing directivity. Finally, a scattering relationship capable of enhancing array localization without increasing array size is proposed, and the directivity improvement claims are verified with beamform plot comparisons and directivity index calculations.

42 ENGINEERING↗
DOE OSTI2020

Two-component axionic dark matter halos

We consider a two-component dark matter halo (DMH) of a galaxy containing ultra-light axions (ULA) of different mass. The DMH is described as a Bose–Einstein condensate (BEC) in its ground state. In the mean-field (MF) limit, we have derived the integro-differential equations for the spherically symmetrical wave functions of the two DMH components. We studied, numerically, the radial distribution of the mass density of ULA and constructed the parameters which could be used to distinguish between the two- and one-component DMH. Here, we also discuss an interesting connection between the BEC ground state of a one-component DMH and Black Hole temperature and entropy, and Unruh temperature.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A comparison of explosively driven shock wave radius versus time scaling approaches

Abstract Explosively driven shock wave radius versus time profiles are frequently used to document energy release and relative explosive performance. Recently, two universal shock wave radius versus time profiles have been presented in the literature, which demonstrate the ability to represent explosively driven shock wave profiles for all explosive sources in any fluid environment. These two universal shock wave profiles are examined here relative to each other and relative to a commonly used nonlinear shock wave profile, which is fit to experimental data for individual explosive materials. The nonlinear profile, originally developed by Dewey, is examined here, and a universal non-dimensional form of the equation is proposed. The universal shock wave profiles are all found to be relatively similar, but with slight variations in a transition region of non-dimensional radii $$0.15\lesssim R^*\lesssim 2$$ 0.15 ≲ R ∗ ≲ 2 . The variations in this region result in different estimations of energy release or blast strength between the curve fits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Neutron Star Equation of State in Light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. Here, we construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗