SEARCH · Engineering Papers

Results for “nonlinear hydrodynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Nonlinearity Effects of Lateral Density Diffusion Coefficient on Gain-Guided VCSEL Performance

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. In this paper, we apply a hydrodynamic model developed from the semiconductor Bloch equations to numerically study the effects of nonlinearity in the diffusion coefficient on single mode operation and direct modulation of a gain-guided InGaAs/GaAs multiple quantum well laser, operating not too far from threshold. We found that a small diffusion coefficient is advantageous for lowering the threshold current and increasing the modulation bandwidth. Most importantly, the effects of nonlinearity in the coefficient can be approximately reproduced by replacing the coefficient with an effective constant diffusion coefficient, which corresponds roughly to the half height density of the density distribution.

Li, Jian-Zhong↗

NASA NTRS1986

Unstable transition properties of the driven magnetohydrodynamic sheet pinch

The unstable transition behavior of a bounded current-carrying two-dimensional magnetofluid is explored, using the hydrodynamic theory developed for parallel shear flows as a guide. Nonlinear excitation of the higher wavenumbers results in the development of electric current sheets of finite extent, as well as the formation of 'attraction currents' centered at the magnetic O-points. A secondary instability mechanism, the dynamic tearing of the electric current sheet, is also observed. This dynamic tearing leads to sawtoothlike temporal oscillations in certain global quantities. The long-time state of the system resembles a nonlinearly saturated state with significant excitation of many wavenumbers. Some features of this state can be understood by means of a Landau nonlinear stability theory based on certain assumptions about the perturbation energy balance.

Dahlburg, R. B.↗

NASA NTRS2008

Tide Corrections for Coastal Altimetry: Status and Prospects

Knowledge of global oceanic tides has markedly advanced over the last two decades, in no small part because of the near-global measurements provided by satellite altimeters, and especially the long and precise Topex/Poseidon time series e.g. [2]. Satellite altimetry in turn places very severe demands on the accuracy of tidal models. The reason is clear: tides are by far the largest contributor to the variance of sea-surface elevation, so any study of non-tidal ocean signals requires removal of this dominant tidal component. Efforts toward improving models for altimetric tide corrections have understandably focused on deep-water, open-ocean regions. These efforts have produced models thought to be generally accurate to about 2 cm rms. Corresponding tide predictions in shelf and near-coastal regions, however, are far less accurate. This paper discusses the status of our current abilities to provide near-global tidal predictions in shelf and near-coastal waters, highlights some of the difficulties that must be overcome, and attempts to divine a path toward some degree of progress. There are, of course, many groups worldwide who model tides over fairly localized shallow-water regions, and such work is extremely valuable for any altimeter study limited to those regions, but this paper considers the more global models necessary for the general user. There have indeed been efforts to patch local and global models together, but such work is difficult to maintain over many updates and can often encounter problems of proprietary or political nature. Such a path, however, might yet prove the most fruitful, and there are now new plans afoot to try again. As is well known, tides in shallow waters tend to be large, possibly nonlinear, and high wavenumber. The short spatial scales mean that current mapping capabilities with (multiple) nadir-oriented altimeters often yield inadequate coverage. This necessitates added reliance on numerical hydrodynamic models and data assimilation, which in turn necessitates very accurate bathymetry with high spatial resolution. Nonlinearity means that many additional compound tides and overtides must be accounted for in our predictions, which increases the degree of modeling effort and increases the amounts of data required to disentangle closely aliased tides.

Ray, Richard D.↗

DOE OSTI2026

Lyman-$α$ forest holography: 3D predictions from 1D measurements

Cosmological analyses of Lyman-$α$ forest clustering rely on either one-dimensional correlations along individual sightlines or three-dimensional correlations between different sightlines. Because these observables probe the matter distribution on very different scales, they have traditionally been analyzed independently. In this work, we bridge this gap using ForestFlow, an emulator trained on a suite of cosmological hydrodynamical simulations that provides a unified description of Lyman-$α$ forest clustering from linear to nonlinear scales. This framework enables us to determine the range of three-dimensional clustering models compatible with the DESI one-dimensional flux power spectrum ($P_{\rm 1D}$). The resulting predictions successfully reproduce the large-scale clustering measured by the DESI BAO analysis and provide physically motivated priors on nonlinear clustering that are used in a companion paper presenting the full-shape analysis of the DESI DR2 Lyman-$α$ forest. We validate our methodology using the large-volume, high-resolution hydrodynamical simulation ACCEL-2, demonstrating excellent agreement across the full range of scales considered. Finally, we combine constraints from the $P_{\rm 1D}$ and BAO analyses on the parameter combinations $b_δσ_8$ and $b_ηf σ_8$, finding that the two probes provide comparable constraining power while exhibiting complementary parameter degeneracies. Our results establish a direct connection between one- and three-dimensional Lyman-$α$ forest measurements through ForestFlow, an approach we term Lyman-$α$ holography by analogy with the reconstruction of higher-dimensional structure from lower-dimensional information.

Chaves-Montero, J. [Barcelona, IFAE] (ORCID:000000↗

DOE OSTI2020

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗

DOE OSTI2024

Interplay between gain and loss in arrays of nonlinear plasmonic nanoparticles: toward parametric downconversion and amplification

With the help of a theoretical model and finite-difference time-domain (FDTD) simulations based on the hydrodynamic–Maxwell model, we examine the effect of difference-frequency generation (DFG) in an array of L-shaped metal nanoparticles (MNPs) characterized by intrinsic plasmonic nonlinearity. The outcomes of the calculations reveal the spectral interplay between gain and loss in the vicinity of the fundamental frequency of the localized surface plasmon resonances. Subsequently, we identify different array thicknesses and pumping regimes facilitating parametric amplification and spontaneous parametric downconversion. Our results suggest that the parametric amplification regime becomes feasible on a scale of hundreds of nanometers and spontaneous parametric downconversion on the scale of tens of nanometers, opening up new exciting opportunities for developing building blocks of photonic metasurfaces.

Shah, Syed A.↗

DOE OSTI2020

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Nonlinear evolution of tidally distorted accretion disks: Two-dimensional simulations

According to a previously published linear analysis, the tidal distortion of accretion disks in binary star systems produces a local hydrodynamic instability to m = 1 internal waves, which may have arbitrarily small wavelengths in the absence of viscosity. The instability is three-dimensional and approximately incrompressible. To explore the nonlinear outcome of this instability, we develop a shearing-sheet approximation on scales comparable to the disk thickness. The large-scale azimuthal variation of the disk is represented by varying the local metric with time (local orbital phase). The hydrodynamic equations can then be posed two-dimensionally on local meridional planes. We solve these equations with a second-order gasdynamical code based on the Total-Variation-Diminishing scheme. Our simulations confirm the predicted linear growth rate. The modes saturate chaotically at velocities scaling as the product of the linear growth rate and the wavelength. If the wavelength is small compared with the disk thickness, the modes remain nearly incompressible even when nonlinear. The two-dimensional power spectrum of velocities after saturation is roughly isotropic and extends over a broad range of scales in an approximately power-law fashion. We measure the heating rate associated with the nonlinear dissipation of the modes. The dissipation implies a secular torque on the disk and a return of angular momentum to the secondary star via the tidal potential. The estimated torque is somewhat larger than the tidal torque produced by maximal disk viscosity (alpha approximately 1). At least in these two-dimensional simulations, however, there is no significant angular momentum flux within the disk.

Ryu, Dongsu↗

DOE OSTI2020

Thermomagnetic instability of plasma composition gradients

In this paper, we show that, under Braginskii magneto-hydrodynamics, anti-parallel gradients in an average ion charge state and electron temperature can be unstable to the growth of self-generated magnetic fields. The instability is analogous to the field-generating thermomagnetic instability, although it is driven by the collisional thermal force magnetic source term rather than the Biermann battery term. The gradient in ion charge state causes a gradient in collisionality, which couples with temperature perturbations to create a self-generated magnetic field. This magnetic field deflects the electron heat flux in a way that reinforces the temperature perturbation. The derived linearized growth rate, typically on hydrodynamic timescales, includes resistive and thermal smoothing. It increases with large ion composition gradients and electron heat flux, conditions typical of the hohlraum walls or contaminant mix jets in inertial confinement fusion implosions. However, extended magneto-hydrodynamic simulations indicate that the instability is usually dominated and stabilized by nonlinear Nernst advection, in a similar manner to the standard thermomagnetic instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Electromagnetic pinned solitons for space debris detection

Electromagnetic “pinned” solitons in the form of stationary nonlinear waves are studied within the framework of an inertial magneto-hydrodynamic model. These structures, that can arise when a charged source moves in a magnetized plasma, have a velocity that is equal to the source velocity and, hence, appear as “pinned” structures that envelope the source. We investigate the excitation of such solitons in the Low Earth Orbit region due to the passage of charged orbital debris objects. The spatial size of these electromagnetic solitons, typically of the order of a few ion skin depths, can be very large in this region. Such solitons can be detected using a variety of ground- or orbit-based radio sounding techniques and may provide a convenient additional means of tracking small sized orbital debris objects that are difficult to spot optically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

The stability of cosmic-ray-dominated shocks - A secondary instability

Previous studies in one spatial dimension showed that the precursors of cosmic-ray-dominated shocks are unstable against traveling acoustic disturbances. Here we report that a secondary, Rayleigh-Taylor type instability can exist in cosmic-ray-dominated media influenced by the acoustic instability. Using the local WKB approximation, the growth rate of the secondary instability is shown to be comparable to that of the one-dimensional acoustic instability itself in the cases we have considered. The nonlinear development has been followed numerically with a two-dimensional PPM hydrodynamics code that also incorporates the two-fluid cosmic-ray energy equation. We show that the secondary instability may cause the precursor and postshock flows to become highly turbulent.

Ryu, Dongsu↗

NASA NTRS2007

Hydrodynamic Force on a Cylinder Oscillating at Low Frequency

The hydrodynamic force on a cylinder oscillating transversely to its axis is a nonlinear function of the displacement amplitude x0. We report measurements and numerical calculations of the force at frequencies low enough that delta > R, where delta is the viscous penetration length and R is the cylinder radius. For small amplitudes, the numerically calculated Fourier transform of the force per unit length, F(sub small), agrees with Stokes' analytical calculation. For larger amplitudes, the force per unit length found by both calculation and measurement is F = F(sub small)C (x(sub 0)/delta,R/delta). The complex function C depends only weakly on R/delta, indicating that x0/delta is more appropriate as a scaling variable than the Keulegan-Carpenter number KC = pi*x(sub 0)/R. The measurements used a torsion oscillator driven at frequencies from 1 to 12 Hz while immersed in dense xenon. The oscillator comprised cylinders with an effective radius of R = 13.4 micron and oscillation amplitudes as large as x(sub 0)/delta = 4 (corresponding to KC as large as 71). The calculations used similar conditions except that the amplitudes were as large as x0/delta = 28.

Berg, Robert F.↗

DOE OSTI2024

High‐Order Nonlinear Frequency Conversion in Transparent Conducting Oxide Thin Films

The study of conductive oxides has gained momentum within the photonics community due to their unique linear and nonlinear optical properties. Despite recent experiments reporting on high harmonic generation from thin films, the optical/electronic behavior of these compounds at the nanoscale is still not fully understood due to the lack of a suitable theoretical model. In the present work, aluminum zinc oxide is excited near its epsilon-near-zero crossing point using incident femtosecond pulses having peak power densities in the 1 TW cm −2 range. A relatively efficient frequency up-conversion including even and odd harmonics up to the seventh order is observed. A hydrodynamic-Maxwell theoretical approach is adopted, capable of simultaneously taking into account linear and nonlinear dispersions, nonlocal effects, surface, magnetic, and bulk nonlinearities in a spectral region that spans over two and a half octaves from the UV to the NIR region. The study enables a deeper understanding of the fundamental material parameters regulating optical nonlinearities, providing important insights to engineer this class of materials for applications in sensing, ultra-fast physics, and spectroscopy.

material modeling↗

DOE OSTI2024

Stability optimization of energetic particle driven modes in nuclear fusion devices: the FAR3d gyro-fluid code

The development of reduced models provide efficient methods that can be used to perform short term experimental data analysis or narrow down the parametric range of more sophisticated numerical approaches. Reduced models are derived by simplifying the physics description with the goal of retaining only the essential ingredients required to reproduce the phenomena under study. This is the role of the gyro-fluid code FAR3d, dedicated to analyze the linear and nonlinear stability of Alfvén Eigenmodes (AE), Energetic Particle Modes (EPM) and magnetic-hydrodynamic modes as pressure gradient driven mode (PGDM) and current driven modes (CDM) in nuclear fusion devices. Such analysis is valuable for improving the plasma heating efficiency and confinement; this can enhance the overall device performance. The present review is dedicated to a description of the most important contributions of the FAR3d code in the field of energetic particles (EP) and AE/EPM stability. FAR3d is used to model and characterize the AE/EPM activity measured in fusion devices as LHD, JET, DIII-D, EAST, TJ-II and Heliotron J. In addition, the computational efficiency of FAR3d facilitates performing massive parametric studies leading to the identification of optimization trends with respect to the AE/EPM stability. This can aid in identifying operational regimes where AE/EPM activity is avoided or minimized. This technique is applied to the analysis of optimized configurations with respect to the thermal plasma parameters, magnetic field configuration, external actuators and the effect of multiple EP populations. In addition, the AE/EPM saturation phase is analyzed, taking into account both steady-state phases and bursting activity observed in LHD and DIII-D devices. The nonlinear calculations provide: the induced EP transport, the generation of zonal structures as well as the energy transfer towards the thermal plasma and between different toroidal/helical families. Finally, FAR3d is used to forecast the AE/EPM stability in operational scenarios of future devices as ITER, CFETR, JT60SA and CFQS as well as possible approaches to optimization with respect to variations in the most important plasma parameters.

Alfv én Eigenmodes↗

NASA NTRS1983

Flare loop radiative hydrodynamics. IV - Dynamic evolution of unstable semiempirical loop models

The evolution of the unstable solar atmosphere into the nonlinear phase, in response to various perturbations, is followed. The initial dynamic evolution of the atmosphere follows the predictions of linear stability analysis. In the nonlinear phase, rapid changes are confined to the transition region; these changes are manifested as a propagation of the transition region through the plasma, i.e., chromospheric evaporation or condensation. Global evolution therefore proceeds on the coronal conductive time scale. The rate of propagation of the transition region is determined by the imbalance between the energy supplied by thermal conduction from the corona and radiative cooling within the transition region itself. Flow velocities in the lower corona during evaporation or condensation are, in the cases studied, of order 3 km/s. The observed dynamic evolution is consistent with the existence of relatively long-lived coronal loops whose brightnesses vary on the evaporative time scale.

An, C.-H.↗

DOE OSTI2024

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1981

Galaxy formation in an intergalactic medium dominated by explosions

The evolution of galaxies in an intergalactic medium dominated by explosions of star systems is considered analogously to star formation by nonlinearly interacting processes in the interstellar medium. Conditions for the existence of a hydrodynamic instability by which galaxy formation leads to more galaxy formation due to the propagation of the energy released at the death of massive stars are examined, and it is shown that such an explosive amplification is possible at redshifts less than about 5 and stellar system masses between 10 to the 8th and 10 to the 12th solar masses. Explosions before a redshift of about 5 are found to lead primarily to the formation of massive stars rather than galaxies, while those at a redshift close to 5 will result in objects of normal galactic scale. The model also predicts a dusty interstellar medium preventing the detection of objects of redshift greater than 3, numbers and luminosities of protogalaxies comparable to present observations, unvirialized groups of galaxies lying on two-dimensional surfaces, and a significant number of black holes in the mass range 1000-10,000 solar masses.

Ostriker, J. P.↗

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↗