Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear dynamics in fluids”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mean Flow from Phase Averages in the 2D Boussinesq Equations

The atmosphere and ocean are described by highly oscillatory PDEs that challenge both our understanding of their dynamics and their numerical approximation. This paper presents a preliminary numerical study of one type of phase averaging applied to mean flows in the 2D Boussinesq equations that also has application to numerical methods. The phase averaging technique, well-known in dynamical systems theory, relies on a mapping using the exponential operator, and then an averaging over the phase. The exponential operator has connections to the Craya–Herring basis pioneered by Jack Herring to study the fluid dynamics of oscillatory, nonlinear fluid dynamics. In this paper, we perform numerical experiments to study the effect of this averaging technique on the time evolution of the solution. We explore its potential as a definition for mean flows. We also show that, as expected from theory, the phase-averaging method can reduce the magnitude of the time rate of change in the PDEs, making them potentially suitable for time stepping methods.

Wingate, Beth A.↗
DOE OSTI2022

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, nonlinear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these nonlinear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modeled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.

Young, David P.↗
NASA NTRS1977

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗
NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗
DOE OSTI2022

Physics-constrained deep learning of nonlinear normal modes of spatiotemporal fluid flow dynamics

In this study, we present a physics-constrained deep learning method to discover and visualize from data the invariant nonlinear normal modes (NNMs) which contain the spatiotemporal dynamics of the fluid flow potentially containing strong nonlinearity. Specifically, we develop a NNM-physics-constrained convolutional autoencoder (NNM-CNN-AE) integrated with a multi-temporal-step dynamics prediction block to learn the nonlinear modal transformation, the NNMs containing the spatiotemporal dynamics of the flow, and reduced-order reconstruction and long-time future-state prediction of the flow fields, simultaneously. In test cases, we apply the developed method to analyze different flow regimes past a cylinder, including laminar flows with low Reynolds number in transient and steady states (RD = 100) and high Reynolds number flow (RD = 1000), respectively. The results indicate that the identified NNMs are able to reveal the nonlinear spatiotemporal dynamics of these flows, and the NNMs-based reduced-order modeling consistently achieves better accuracy with orders of magnitudes smaller errors in construction and prediction of the nonlinear velocity and vorticity fields, compared to the linear proper orthogonal decomposition (POD) method and the Koopman-constrained-CNN-AE using the same number or dimension of modes. We perform an analysis of the modal energy distribution of NNMs and find that compared to POD modes, the few fundamental NNMs capture a very high level of total energy of the flow, which is advantageous for reduced-order modeling and representation of the complex flows. Finally, we discuss the potentials and limitations of the presented method.

Mechanics↗
NASA NTRS1996

Turbulent Fluid Motion 6: Turbulence, Nonlinear Dynamics, and Deterministic Chaos

Several turbulent and nonturbulent solutions of the Navier-Stokes equations are obtained. The unaveraged equations are used numerically in conjunction with tools and concepts from nonlinear dynamics, including time series, phase portraits, Poincare sections, Liapunov exponents, power spectra, and strange attractors. Initially neighboring solutions for a low-Reynolds-number fully developed turbulence are compared. The turbulence is sustained by a nonrandom time-independent external force. The solutions, on the average, separate exponentially with time, having a positive Liapunov exponent. Thus, the turbulence is characterized as chaotic. In a search for solutions which contrast with the turbulent ones, the Reynolds number (or strength of the forcing) is reduced. Several qualitatively different flows are noted. These are, respectively, fully chaotic, complex periodic, weakly chaotic, simple periodic, and fixed-point. Of these, we classify only the fully chaotic flows as turbulent. Those flows have both a positive Liapunov exponent and Poincare sections without pattern. By contrast, the weakly chaotic flows, although having positive Liapunov exponents, have some pattern in their Poincare sections. The fixed-point and periodic flows are nonturbulent, since turbulence, as generally understood, is both time-dependent and aperiodic.

Deissler, Robert G.↗
NASA NTRS2003

The Effect of the Density Ratio on the Nonlinear Dynamics of the Unstable Fluid Interface

Here we report multiple harmonic theoretical solutions for a complete system of conservation laws, which describe the large-scale coherent dynamics in RTI and RMI for fluids with a finite density ratio in the general three-dimensional case. The analysis yields new properties of the bubble front dynamics. In either RTI or RMI, the obtained dependencies of the bubble velocity and curvature on the density ratio differ qualitatively and quantitatively from those suggested by the models of Sharp (1984), Oron et al. (2001), and Goncharov (2002). We show explicitly that these models violate the conservation laws. For the first time, our theory reveals an important qualitative distinction between the dynamics of the RT and RM bubbles.

Abarzhi, S. I.↗
NASA NTRS1988

The nonlinear dynamics of a spacecraft coupled to the vibration of a contained fluid

The dynamics of a linear spacecraft mode coupled to a nonlinear low gravity slosh of a fluid in a cylindrical tank is investigated. Coupled, nonlinear equations of motion for the fluid-spacecraft dynamics are derived through an assumed mode Lagrangian method. Unlike linear fluid slosh models, this nonlinear slosh model retains two fundamental slosh modes and three secondary modes. An approximate perturbation solution of the equations of motion indicates that the nonlinear coupled system response involves fluid-spacecraft modal resonances not predicted by either a linear, or a nonlinear, uncoupled slosh analysis. Experimental results substantiate the analytical predictions.

Peterson, Lee D.↗
NASA NTRS1992

Sonic boom predictions using a modified Euler code

The environmental impact of a next generation fleet of high-speed civil transports (HSCT) is of great concern in the evaluation of the commercial development of such a transport. One of the potential environmental impacts of a high speed civilian transport is the sonic boom generated by the aircraft and its effects on the population, wildlife, and structures in the vicinity of its flight path. If an HSCT aircraft is restricted from flying overland routes due to excessive booms, the commercial feasibility of such a venture may be questionable. NASA has taken the lead in evaluating and resolving the issues surrounding the development of a high speed civilian transport through its High-Speed Research Program (HSRP). The present paper discusses the usage of a Computational Fluid Dynamics (CFD) nonlinear code in predicting the pressure signature and ultimately the sonic boom generated by a high speed civilian transport. NASA had designed, built, and wind tunnel tested two low boom configurations for flight at Mach 2 and Mach 3. Experimental data was taken at several distances from these models up to a body length from the axis of the aircraft. The near field experimental data serves as a test bed for computational fluid dynamic codes in evaluating their accuracy and reliability for predicting the behavior of future HSCT designs. Sonic boom prediction methodology exists which is based on modified linear theory. These methods can be used reliably if near field signatures are available at distances from the aircraft where nonlinear and three dimensional effects have diminished in importance. Up to the present time, the only reliable method to obtain this data was via the wind tunnel with costly model construction and testing. It is the intent of the present paper to apply a modified three dimensional Euler code to predict the near field signatures of the two low boom configurations recently tested by NASA.

Siclari, Michael J.↗
DOE OSTI2021

Beyond Ohm's law: Bernoulli effect and streaming in electron hydrodynamics

Recent observations of nonlocal transport in ultraclean two-dimensional materials raised the tantalizing possibility of accessing hydrodynamic correlated transport of a many-electron state. However, it has been pointed out that nonlocal transport can also arise from impurity scattering rather than interaction. At the crux of the ambiguity is the focus on linear effects, i.e., Ohm's law, which cannot easily differentiate among different modes of transport. Furthermore we propose experiments that can reveal rich hydrodynamic features in the system by tapping into the nonlinearity of the Navier-Stokes equation. Three experiments we propose will each manifest a unique phenomenon that is well known in classical fluids: the Bernoulli effect, Eckart streaming, and Rayleigh streaming. Analysis of known parameters confirms that the proposed experiments are feasible and the hydrodynamic signatures are within reach of graphene-based devices. Experimental realization of any one of the three phenomena will provide a stepping stone to formulating and exploring the notions of nonlinear electron fluid dynamics with an eye to celebrated examples from classical nonlaminar flows, e.g., pattern formation and turbulence.

2-dimensional systems↗
DOE OSTI2020

Experimental determination of shock speed versus exciter speed in a two-dimensional dusty plasma

A shock that is continuously driven by a moving exciter will propagate at a speed that depends on the exciter speed. We obtained this dependence experimentally, in a strongly coupled dusty plasma that was prepared as a single two-dimensional layer of charged microparticles. Attaining this result required an experimental advance, developing a method of driving a shock continuously, which we did using an exciter moving at a constant supersonic speed, analogous to a piston in a cylinder. The resulting compressional pulse was a shock that propagated steadily without weakening, ahead of the moving exciter. In this work, we compare our experimental results to an empirical form $M_{\text{shock}} = 1 + sM_{\text{exciter}}$, and to the prediction of a recent simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Inclusive curvaturelike framework for describing dissipation: Metriplectic 4-bracket dynamics

An inclusive framework for joined Hamiltonian and dissipative dynamical systems that are thermodynamically consistent, i.e., preserve energy and produce entropy, is given. The dissipative dynamics of the framework is based on the metriplectic 4-bracket , a quantity like the Poisson bracket defined on phase space functions, but unlike the Poisson bracket has four slots with symmetries and properties motivated by Riemannian curvature. Metriplectic 4-bracket dynamics is generated using two generators, the Hamiltonian and the entropy, with the entropy being a Casimir of the Hamiltonian part of the system. The formalism includes known previous binary bracket theories for dissipation or relaxation as special cases. Rich geometrical significance of the formalism and methods for constructing metriplectic 4-brackets are explored. Many examples of both finite and infinite dimensions are given.

Physics↗
DOE OSTI2022

Effective drift velocity from turbulent transport by vorticity

Here, we highlight the differing roles of vorticity and strain in the transport of coarse-grained scalars at length scales larger than ℓ by smaller scale (subscale) turbulence. We use the first term in a multiscale gradient expansion due to which exhibits excellent correlation with the exact subscale physics when the partitioning length ℓ is any scale smaller than that of the spectral peak. We show that unlike subscale strain, which acts as an anisotropic diffusion/antidiffusion tensor, subscale vorticity's contribution is solely a conservative advection of coarse-grained quantities by an eddy-induced nondivergent velocity, v * , that is proportional to the curl of vorticity. Therefore, material (Lagrangian) advection of coarse-grained quantities is accomplished not by the coarse-grained flow velocity, u¯ ℓ , but by the effective velocity, u¯ ℓ + v * , the physics of which may improve commonly used LES models.

42 ENGINEERING↗
DOE OSTI2026

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗