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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Geostrophic adjustment in a shallow-water numerical model as it relates to thermospheric dynamics

The theory of geostrophic adjustment and its application to the dynamics of the high latitude thermosphere have been discussed in previous papers based on a linearized treatment of the fluid dynamical equations. However, a linearized treatment is only valid for small Rossby numbers given by Ro = V/fL, where V is the wind speed, f is the local value of the Coriolis parameter, and L is a characteristic horizontal scale for the flow. For typical values in the auroral zone, the approximation is not reasonable for wind speeds greater than 25 m/s or so. A shallow-water (one layer) model was developed that includes the spherical geometry and full nonlinear dynamics in the momentum equations in order to isolate the effects of the nonlinearities on the adjustment process. A belt of accelerated winds between 60 deg and 70 deg latitude was used as the initial condition. The adjustment process was found to proceed as expected from the linear formulation, but that an asymmetry between the response for an eastward and westward flow results from the nonlineawr curvature (centrifugal) terms. In general, the amplitude of an eastward flowing wind will be less after adjustment than a westward wind. For instance, if the initial wind velocity is 300 m/s, the linearized theory predicts a final wind speed of 240 m/s, regardless of the flow direction. However, the nonlinear curvature terms modify the response and produce a final wind speed of only 200 m/s for an initial eastward wind and a final wind speed of almost 300 m/s for an initial westward flow direction. Also, less gravity wave energy is produced by the adjustment of the westward flow than by the adjustment of the eastward flow. The implications are that the response of the thermosphere should be significantly different on the dawn and dusk sides of the auroral oval. Larger flow velocities would be expected on the dusk side since the plasma will accelerate the flow in a westward direction in that sector.

Larsen, M. F.↗
NASA NTRS1990

Computational technology for flight vehicles; Proceedings of the Symposium, Washington, DC, Nov. 5-7, 1990

Recent advances in computational fluid mechanics are discussed in reviews and reports. Sections are devoted to (1) the modeling of local phenomena and edge effects in solids, (2) stochastic modeling and simulation of fracture toughness, and (3) partitioning strategy and new finite elements. Particular attention is given to global and local finite-element/spectral-boundary-element techniques for failure analysis, simulations of microfracture in metal-matrix composites, fatigue analysis of cracked anisotropic plates under stochastic loading, mathematical modeling for the analysis of nonlinear aircraft dynamics, physical and mathematical modeling of wave propagation in the Ariane 5 VEB structure, partitioning based on symmetry transformations, an FEM approach to adaptive reliability assurance, and time-domain FEMs for the large rotational dynamics of multibody systems.

Noor, Ahmed K.↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Validation of CFD simulations of the moored DeepCwind offshore wind semisubmersible in irregular waves

This article examines the use of computational fluid dynamics (CFD) simulations to predict the response of a floating wind platform to irregular-wave excitation. This work was conducted as part of the Reproducible CFD Joint Industry Project for Floating Offshore Wind Applications and involved verification across several participants and modeling tools and validation against experiments. The authors pay special attention to the uncertainties in both CFD results and experimental measurements. In this work, we perform detailed comparisons of the incident waves and the motion of a moored structure. The nonlinear, low-frequency resonance motion is of particular interest because it potentially drives the mooring and tower-base loads. The verification and validation study is partially successful in that the CFD simulations capture the low-frequency slow-drift motion well but underpredict the low-frequency pitch resonance. This underprediction can be attributed in part to the CFD incident waves, which showed some discrepancies with the experimental waves, especially around extreme events. The effects of the wave discrepancies are also estimated and investigated using a mid-fidelity OpenFAST model. Overall, the present study increases our confidence in using CFD simulations to predict the global performance of offshore wind platforms in irregular waves and to produce data for the calibration of lower-fidelity models.

17 WIND ENERGY↗
NASA NTRS1991

Numerical models of extragalactic radio sources

When supercomputer-implemented numerical simulations analyzing the nonlinear physics inherent in the hydrodynamic and MHD equations are applied to extragalactic radio sources, many of the complex structures observed on telescopic images are reproduced. Attention is presently given to recently obtained results from 2D and 3D numerical simulations of the formation and evolution of extended radio morphologies; these numerical models allow the exploration of such physical phenomena as the role of magnetic fields in the dynamics and emissivity of extended radio galaxies, intermittent outflow from the cores of active galaxies, fluid-jet instabilities, and the bending of collimated outflows by motion through the intergalactic medium.

Burns, Jack O.↗
DOE OSTI2020

CU-BENs: A structural modeling finite element library

The present work discusses capabilities within the finite element library CU-BENs. CU-BENs focuses on applying the finite element method to structural mechanics problems encountered within the context of inverse problems and partitioned fluid–structure interaction; thus element formulations are primarily of a structural type — truss, frame, and triangular discrete Kirchhoff theory shells. CU-BENs defaults to the skyline sparse storage scheme for the system matrix, but also supports other storage schemes when using external libraries such as LAPACK and UMFPACK. CU-BENs includes built-in nonlinear solution strategies, such as the Newton Raphson method and modified spherical arc length method, that are available within static analyses as well within the context of transient dynamic analyses involving a generalized-α implementation of the Newmark implicit time integration scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
NASA NTRS1983

Nonlinear aerodynamic modeling of flap oscillations in transonic flow - A numerical validation

The regime of validity of a nonlinear aerodynamic force and moment formulation, based on concepts from nonlinear functional analysis and applicable to a transonic airfoil with a deflecting flap, is investigated. A time-dependent finite difference technique is used to evaluate the aerodynamic data of the formulation in terms of specified, characteristic motions. Flap-motion histories are generated from the flap inertial equations of motion, with aerodynamic reactions specified by the moment formulation. The motion histories depicting the cases of decaying and growing flap oscillations are compared with histories generated through simultaneous, coupled solution of the fluid-dynamic equations and flap inertial equations of motion. The range of applicability of the formulation is discussed.

Chyu, W. J.↗
NASA NTRS1990

A finite element solver for 3-D compressible viscous flows

Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.

Reddy, K. C.↗
NASA NTRS1992

Propagating confined states in phase dynamics

Theoretical treatment is given to the possibility of the existence of propagating confined states in the nonlinear phase equation by generalizing stationary confined states. The nonlinear phase equation is set forth for the case of propagating patterns with long wavelengths and low-frequency modulation. A large range of parameter values is shown to exist for propagating confined states which have spatially localized regions which travel on a background with unique wavelengths. The theoretical phenomena are shown to correspond to such physical systems as spirals in Taylor instabilities, traveling waves in convective systems, and slot-convection phenomena for binary fluid mixtures.

Brand, Helmut R.↗
NASA NTRS1980

Possible fluid dynamical interpretation of some reported features in the Jovian atmosphere

A fluid dynamical interpretation is presented of the two major types of disturbance found in the southern hemisphere of Jupiter by the Voyager 1 imaging data. The observed features always occur together, and consist of a compact elliptically shaped formation having an anticyclonic flow which is poleward of a pair of more elongated cyclonic structures, as in the Great Red Spot and the white ovals. It is noted that the anticyclonic features at 41 deg S may be described by the cnoidal wave solutions to the appropriate nonlinear evolution equation, and that flow patterns derived in the vicinity of the Great Red Spot and white ovals are strikingly similar to those obtained for the flow around a solitary wave of the type than can exist in a zonal flow such as that found in the Jupiter atmosphere. Results of computations in terms of solitary wave theory of flow fields in the atmospheric structure and zonal velocity profiles determined from Voyager infrared spectroscopy and radiometry data are then presented which show that the pattern must be a singular solitary wave mode, the east-west structure of which is best described by the Korteweg-de-Vries equation

Maxworthy, T.↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

CMLM (Co-Optimized Machine-Learned Manifolds) [SWR-23-41]

Co-optimized Machine-Learned Manifolds (CMLM) is a data-driven approach for developing reduced-order manifold models for high-dimensional chemically reacting systems. It involves a specially designed neural network, the training of which simultaneously optimizes linear combinations of species that define the manifold, nonlinear mapping to outputs of interest such as reaction rates, and (optionally) subfilter closure for large eddy simulation. This software package provides an implementation of the CMLM approach in Python using the PyTorch machine learning library. A few example cases are included, showing how the tool can be applied to different types of data from 0D and 1D reacting simulations performed using Cantera. The neural networks can be saved in a format that is readable by the Pele suite of combustion solvers for use in reacting computational fluid dynamics simulations. This software repository contains several python scripts to perform various tasks associated with the Co-optimized Machine Learned Manifolds (CMLM) model, which is described in Perry, Henry de Frahan, and Yellapantula, CNF, 2022 (https://doi.org/10.1016/j.combustflame.2022.112286). This includes not only the code that defines the CMLM model, but also scripts to generate suitable training data, scripts to pre-process the data, scripts to train the CMLM model, and scripts to plot the output, as well as various other helper files. The scripts depend on several commonly used python libraries for data analysis and chemical reaction computations. The trained models that result from this tool are designed to work with the an interface being implemented in the Pele suite of reacting flow solvers (https://github.com/AMReX-Combustion).

Perry, Bruce↗
NASA NTRS1989

Influence of impeller shroud forces on turbopump rotor dynamics

The shrouded-impeller leakage path forces calculated by Childs (1987) have been analyzed to answer two questions. First, because of certain characteristics of the results of Childs, the forces could not be modeled with traditional approaches. Therefore, an approach has been devised to include the forces in conventional rotordynamic analyses. The forces were approximated by traditional stiffness, damping and inertia coefficients with the addition of whirl-frequency-dependent direct and cross-coupled stiffness terms. The forces were found to be well-modeled with this approach. Finally, the effect these forces had on a simple rotor-bearing system was analyzed, and, therefore, they, in addition to seal forces, were applied to a Jeffcott rotor. The traditional methods of dynamic system analysis were modified to incorporate the impeller forces and yielded results for the eigenproblem, frequency response, critical speed, transient response and an iterative technique for finding the frequency of free vibration as well as system stability. All results lead to the conclusion that the forces have little influence on natural frequency but can have appreciable effects on system stability. Specifically, at higher values of fluid swirl at the leakage path entrance, relative stability is reduced. The only unexpected response characteristics that occurred are attributed to the nonlinearity of the model.

Williams, Jim P.↗
NASA NTRS1993

Influence of impeller shroud forces on turbopump rotor dynamics

The shrouded-impeller leakage path forces calculated by Childs have been analyzed to answer two questions. First, because of certain characteristics or the results of Childs, the forces could not be modeled with traditional approaches. Therefore, an approach has been devised to include the forces in conventional rotordynamic analyses. The forces were found to be well-modeled with this approach. Finally, the effect these forces had on a simple rotor-bearing system was analyzed, and, therefore, they, in addition to seal forces, were applied to a Jeffcott rotor. The traditional methods of dynamic system analysis were modified to incorporate the impeller forces and yielded results for the eigenproblem, frequency response, critical speed, transient response, and an iterative technique for finding the frequency of free vibration as well as system stability. All results lead to the conclusion that the forces have little influence on natural frequency but can have appreciable effects on system stability. Specifically, at higher values of fluid swirl at the leakage path entrance, relative stability is reduced. The only unexpected response characteristics that occurred are attributed to the nonlinearity of the model.

Williams, J. P.↗
NASA NTRS1983

Influence of unbalance on the nonlinear dynamical response and stability of flexible rotor-bearing systems

In this paper, some of the effects of unbalance on the nonlinear response and stability of flexible rotor-bearing systems is presented from both a theoretical and experimental standpoint. In a linear system, operating above its stability threshold, the amplitude of motion grows exponentially with time and the orbits become unbounded. In an actual system, this is not necessarily the case. The actual amplitudes of motion may be bounded due to various nonlinear effects in the system. These nonlinear effects cause limit cycles of motion. Nonlinear effects are inherent in fluid film bearings and seals. Other contributors to nonlinear effects are shafts, couplings and foundations. In addition to affecting the threshold of stability, the nonlinear effects can cause jump phenomena to occur at not only the critical speeds, but also at stability onset or restabilization speeds.

Gunter, E. J.↗
NASA NTRS1993

Nonlinear waves and solitons in molecular clouds

We begin a study of nonlinear wave phenomena in molecular clouds. These clouds exhibit highly nonlinear structure that is often described in terms of 'clumps' and 'filaments' which are bouncing around, twisting, and colliding within the cloud. These clouds are important because they ultimately produce the initial conditions for the star formation process. Our motivation is to explore the possibility that solitons (i.e., spatially localized, single-hump wave entities which often exhibit remarkable stability) can live in these molecular clouds and produce their observed structure. In this paper we focus on the case of one spatial dimension, and we show that a rich variety of nonlinear waves can exist in molecular cloud fluid systems (where self-gravity is included). We show that in the absence of magnetic fields no true soliton solutions are allowed, although highly nonlinear waves (whose crests become widely spaced and thus soliton-like) do exist. For clouds with embedded magnetic fields, we derive a model equation which describes the behavior of wave phenomena; this model equation allows solutions which correspond to nonlinear waves, solitons, and topological solitons. We briefly consider the stability of these wave entities and discuss the possible role they play in molecular cloud dynamics.

Adams, Fred C.↗
NASA NTRSDate not supplied

Linear Parameter-Varying Control Compatible with Excitation for Real-Time Modeling

The Modeling and Control for Agile Aircraft Development project (MCAAD) seeks to reduce the time spent on wind tunnel testing and computational fluid dynamic modeling when developing control laws for new commercial and military aircraft. MCAAD uses advanced dynamic modeling procedures coupled with real-time control design which can reduce the time and cost spent on traditional aircraft development methods. This is done in-flight by using a real-time modeling procedure where programmed test inputs (PTI) deflect the aircraft control surfaces, and the responses are measured in order to create a dynamic model. This paper presents an 𝐻∞-based linear parameter-varying (LPV) control law designed on a nonlinear rigid-body dynamic model which commands control surface deflections at frequencies above the PTI frequency band. Thus the controller will not suppress the response due to the PTI. The LPV controller is compared with a previously designed Nonlinear Dynamic Inversion (NDI) controller and the coherence between the commanded control surface deflections and the PTI is plotted for both controllers. The LPV controller is shown to distort the PTI less than the NDI controller.

Luke J. Miller↗