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

Dynamics of electromechanical flow structures.

Free-surface gravity flows and capillary wicking provide examples of flow structures with fluid partially ducted at free surfaces by external forces. Wall-less electromechanical flow structures are developed which have a similar nature, but with polarization forces providing the orientation at free surfaces. Like their mechanical counterparts, these have the ability to ingest liquid or expel vapor through their walls. The structures consist of electrodes running in the flow direction z with slowly varying cross sections in a plane transverse to the flow. A formulation is given of the long-wave nonlinear (principal mode) dynamics, with use made of energy functions to represent a broad class of possible mechanical and electrical structure geometries.

Jones, T. B., Jr.↗
NASA NTRS1981

Dynamic behavior of continuous cantilevered pipes conveying fluid near critical velocities

The plane motion of a cantilevered pipe conveying fluid is examined when the flow velocity is in the neighborhood of that generating flutter. In contrast to previous studies, the flow velocity is not prescribed as a constant, but is determined from the laws of motion. A system of two nonlinear partial differential equations which are coupled through the nonlinear terms is thereby obtained. The solution is found by the use of the Krylov-Bogoliubov averaging method and the results are discussed indicating the effect of nonlinearities.

Rousselet, J.↗
NASA NTRS2015

Higher-Order Methods for Compressible Turbulent Flows Using Entropy Variables

Turbulent flows have a large range of spatial and temporal scales which need to be resolved in order to obtain accurate predictions. Higher-order methods can provide greater efficiency for simulations requiring high spatial and temporal resolution, allowing for solutions with fewer degrees of freedom and lower computational cost than traditional second-order computational fluid dynamics (CFD) methods.1 Higher-order methods have been widely used for turbulent flows. However, the reduced numerical stabilization present in higher-order schemes implies that special care needs to be taken in the development of numerical methods to suppress nonlinear instabilities.2–6 In this work we present the development of a higher-order space-time discontinuous Galerkin method with a focus on the aspects of our numerical scheme required for ensuring nonlinear stability for turbulent simulations at high Reynolds numbers.

Diosady, Laslo T.↗
NASA NTRS1994

Spectral Methods for Computational Fluid Dynamics

As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.

Zang, T. A.↗
NASA NTRS1978

Filtering of non-linear instabilities

For Courant numbers larger than one and cell Reynolds numbers larger than two, oscillations and in some cases instabilities are typically found with implicit numerical solutions of the fluid dynamics equations. This behavior has sometimes been associated with the loss of diagonal dominance of the coefficient matrix. It is shown that these problems can be related to the choice of the spatial differences, with the resulting instability related to aliasing or nonlinear interaction. Appropriate filtering can reduce the intensity of these oscillations and possibly eliminate the instability. These filtering procedures are equivalent to a weighted average of conservation and nonconservation differencing. The entire spectrum of filtered equations retains a three point character as well as second order spatial accuracy. Burgers equation was considered as a model.

Khosla, P. K.↗
NASA NTRS1989

Towards a rigorous approach to artificial dissipation

Artificial dissipation in modern computational fluid dynamics is something of an art. In spite of the dramatic strides made over the last decade, dissipation continues to lack a quantitative theory with predictive power. Such a theory, based on entropy production, is presented. An important side-effect is nonlinear stability for schemes which satisfy certain conditions.

Merriam, Marshal L.↗
NASA NTRS2001

Modeling Wave Driven Non-linear Flow Oscillations: The Terrestrial QBO and a Solar Analog

The Quasi Biennial Oscillation (QBO) of the zonal circulation observed in the terrestrial atmosphere at low latitudes is driven by wave mean flow interaction as was demonstrated first by Lindzen and Holton (1968), shown in a laboratory experiment by Plumb and McEwan (1978), and modeled by others (e.g., Plumb, Dunkerton). Although influenced by the seasonal cycle of solar forcing, the QBO, in principle, represents a nonlinear flow oscillation that can be maintained by a steady source of upward propagating waves. The wave driven non-linearity is of third or odd order in the flow velocity, which regenerates the fundamental harmonic itself to keep the oscillation going - the fluid dynamical analog of the displacement mechanism in the mechanical clock. Applying Hines' Doppler Spread Parameterization (DSP) for gravity waves (GW), we discuss with a global-scale spectral model numerical experiments that elucidate some properties of the QBO and its possible effects on the climatology of the atmosphere. Depending on the period of the QBO, wave filtering can cause interaction with the seasonal variations to produce pronounced oscillations with beat periods around 10 years. Since the seasonal cycle and its variability influence the period of the QBO, it may also be a potent conduit of solar activity variations to lower altitudes. Analogous to the terrestrial QBO, we propose that a flow oscillation may account for the 22-year periodicity of the solar magnetic cycle, potentially answering Dicke (1978) who asked, "Is there a chronometer hidden deep inside the Sun?" The oscillation would occur below the convection region, where gravity waves can propagate. Employing a simplified, analytic model, Hines' DSP is applied to estimate the flow oscillation. Depending on the adopted horizontal wavelengths of GW's, wave amplitudes less than 10 m/s can be made to produce oscillating zonal flows of about 20 m/s that should be large enough to generate a significant oscillation in the magnetic field. For the large length scales of the Sun, the flow cycle period tends to be very long. The period, however, can be made to be 22 years, provided the buoyancy frequency (stability) is sufficiently small, thus placing the proposed flow near the base of the convection zone where a dynamo is now believed to operate.

Mayr, Hans G.↗
NASA NTRS2013

Simulating Effects of High Angle of Attack on Turbofan Engine Performance

A method of investigating the effects of high angle of attack (AOA) flight on turbofan engine performance is presented. The methodology involves combining a suite of diverse simulation tools. Three-dimensional, steady-state computational fluid dynamics (CFD) software is used to model the change in performance of a commercial aircraft-type inlet and fan geometry due to various levels of AOA. Parallel compressor theory is then applied to assimilate the CFD data with a zero-dimensional, nonlinear, dynamic turbofan engine model. The combined model shows that high AOA operation degrades fan performance and, thus, negatively impacts compressor stability margins and engine thrust. In addition, the engine response to high AOA conditions is shown to be highly dependent upon the type of control system employed.

Liu, Yuan↗
NASA NTRS2001

Time-dependent oral absorption models

The plasma concentration-time profiles following oral administration of drugs are often irregular and cannot be interpreted easily with conventional models based on first- or zero-order absorption kinetics and lag time. Six new models were developed using a time-dependent absorption rate coefficient, ka(t), wherein the time dependency was varied to account for the dynamic processes such as changes in fluid absorption or secretion, in absorption surface area, and in motility with time, in the gastrointestinal tract. In the present study, the plasma concentration profiles of propranolol obtained in human subjects following oral dosing were analyzed using the newly derived models based on mass balance and compared with the conventional models. Nonlinear regression analysis indicated that the conventional compartment model including lag time (CLAG model) could not predict the rapid initial increase in plasma concentration after dosing and the predicted Cmax values were much lower than that observed. On the other hand, all models with the time-dependent absorption rate coefficient, ka(t), were superior to the CLAG model in predicting plasma concentration profiles. Based on Akaike's Information Criterion (AIC), the fluid absorption model without lag time (FA model) exhibited the best overall fit to the data. The two-phase model including lag time, TPLAG model was also found to be a good model judging from the values of sum of squares. This model also described the irregular profiles of plasma concentration with time and frequently predicted Cmax values satisfactorily. A comparison of the absorption rate profiles also suggested that the TPLAG model is better at prediction of irregular absorption kinetics than the FA model. In conclusion, the incorporation of a time-dependent absorption rate coefficient ka(t) allows the prediction of nonlinear absorption characteristics in a more reliable manner.

Non-NASA Center↗
NASA NTRS1989

A new rotational relaxation model for use in hypersonic computational fluid dynamics

The theoretical basis for the Landau-Teller equation commonly used to model rotational nonequilibrium is reviewed. Several assumptions underlying this model are indicated to be unrealistic for rarefied hypersonic flow. A new rotational nonequilibrium model based on recent measurements up to 2000K of quantum state transition rates is presented. The new model is applied to the continuum study of normal shock wave structure using both the Burnett equations and a simplified nonlinear constitutive relation. Comparisons are made between shock wave temperature profiles generated from the new model and those generated using the Landau-Teller model. Comparisons of shock reciprocal thicknesses between experimental data for nitrogen and continuum solutions using both rotational models are made. The new rotational model agrees well with experiment up to Mach 6, and under predicts shock thickness at higher Mach numbers. The Landau-Teller model agrees well with experimental shock thickness up to Mach 3, and over predicts the thickness at higher Mach numbers. A modification to the rotational collision number in the Landau-Teller model is found to give results which agree with experimental shock thicknesses at all Mach numbers up to 11.

Chapman, Dean E.↗
NASA NTRS2005

Mathematical Metaphors: Problem Reformulation and Analysis Strategies

This paper addresses the critical need for the development of intelligent or assisting software tools for the scientist who is working in the initial problem formulation and mathematical model representation stage of research. In particular, examples of that representation in fluid dynamics and instability theory are discussed. The creation of a mathematical model that is ready for application of certain solution strategies requires extensive symbolic manipulation of the original mathematical model. These manipulations can be as simple as term reordering or as complicated as discovery of various symmetry groups embodied in the equations, whereby Backlund-type transformations create new determining equations and integrability conditions or create differential Grobner bases that are then solved in place of the original nonlinear PDEs. Several examples are presented of the kinds of problem formulations and transforms that can be frequently encountered in model representation for fluids problems. The capability of intelligently automating these types of transforms, available prior to actual mathematical solution, is advocated. Physical meaning and assumption-understanding can then be propagated through the mathematical transformations, allowing for explicit strategy development.

Thompson, David E.↗
NASA NTRS1974

U.S. National Congress of Applied Mechanics, 7th, University of Colorado, Boulder, Colo., June 3-7, 1974, Proceedings

A series of experimental, theoretical, and review papers dealing with various aspects of applied mechanics, including structural analysis, mathematical methods, and fluid dynamics. Some of the topics covered are: the direct formulation of two-dimensional mechanics problems, the mechanical aspects of athletics, Laplace's tidal equations, cochlear mechanics, dynamic constitutive equations, optical techniques for the determination of displacement and strain fields, chemical and biochemical problems involving nonlinear diffusion processes, thermal stresses, the rotating rod viscometer, finite-element methods in plasticity, ships of minimum resistance, and wide-band random vibration of structures. Individual items are announced in this issue.

Source record↗
NASA NTRS1979

Filtering of non-linear instabilities

For Courant numbers larger than one and cell Reynolds numbers larger than two, oscillations and in some cases instabilities are typically found with implicit numerical solutions of the fluid dynamics equations. This behavior has sometimes been associated with the loss of diagonal dominance of the coefficient matrix. It is shown here that these problems can in fact be related to the choice of the spatial differences, with the resulting instability related to aliasing or nonlinear interaction. Appropriate 'filtering' can reduce the intensity of these oscillations and in some cases possibly eliminate the instability. These filtering procedures are equivalent to a weighted average of conservation and non-conservation differencing. The entire spectrum of filtered equations retains a three-point character as well as second-order spatial accuracy. Burgers equation has been considered as a model. Several filters are examined in detail, and smooth solutions have been obtained for extremely large cell Reynolds numbers.

Khosla, P. K.↗
NASA NTRS1994

Geometry of thin liquid sheet flows

Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Chubb, Donald L.↗
NASA NTRS2016

Parametric Modeling Investigation of a Radially-Staged Low-Emission Aviation Combustor

Aviation gas-turbine combustion demands high efficiency, wide operability and minimal trace gas emissions. Performance critical design parameters include injector geometry, combustor layout, fuel-air mixing and engine cycle conditions. The present investigation explores these factors and their impact on a radially staged low-emission aviation combustor sized for a next-generation 24,000-lbf-thrust engine. By coupling multi-fidelity computational tools, a design exploration was performed using a parameterized annular combustor sector at projected 100% takeoff power conditions. Design objectives included nitrogen oxide emission indices and overall combustor pressure loss. From the design space, an optimal configuration was selected and simulated at 7.1, 30 and 85% part-power operation, corresponding to landing-takeoff cycle idle, approach and climb segments. All results were obtained by solution of the steady-state Reynolds-averaged Navier-Stokes equations. Species concentrations were solved directly using a reduced 19-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear K-epsilon model. This research demonstrates revolutionary combustor design exploration enabled by multi-fidelity physics-based simulation.

design analysis↗
NASA NTRS1992

Role of shock dynamics in transonic flutter

A computational study of the influence of shock motion on flutter and divergence in transonic flow is presented. The numerical scheme models the entire fluid-structure system as a single continuum dynamics problem, by using a mixed Eulerian-Lagrangian formulation. No assumptions of small displacements are made, but the effect of viscosity is neglected. The results from this study indicate that the shock dynamics gives rise to limit cycles and highly nonlinear aeroelastic phenomena, such as weak divergence and flutter-divergence interactions. Although the shocks typically are destabilizing at the linear flutter boundary, they often have a strongly stabilizing effect for moderate-amplitude motions. The shocks are thus capable of quenching an emerging bending-torsion flutter motion and turning it into limit cycle flutter. The usefulness of classical flutter and divergence boundary diagrams is severely limited in transonic flow, because much of the global dynamic stability information is lost in such a presentation.

Bendiksen, Oddvar O.↗
NASA NTRS1994

Sensitivity derivatives for three dimensional supersonic Euler code using incremental iterative strategy

In a recent work, an incremental strategy was proposed to iteratively solve the very large systems of linear equations that are required to obtain quasianalytical sensitivity derivatives from advanced computational fluid dynamics (CFD) codes. The technique was sucessfully demonstrated for two large two-dimensional problems: a subsonic and a transonic airfoil. The principal feature of this incremental iterative stategy is that it allows the use of the identical approximate coefficient matrix operator and algorithm to solve the nonlinear flow and the linear sensitivity equations; at convergence, the accuracy of the sensitivity derivatives is not compromised. This feature allows a comparatively straightforward extension of the methodology to three-dimensional problems; this extension is successfully demonstrated in the present study for a space-marching solution of the three-dimensional Euler equations over a Mach 2.4 blended wing-body configuration.

Korivi, Vamshi Mohan↗
NASA NTRS2001

Simulation of Supersonic Jet Noise with the Adaptation of Overflow CFD Code and Kirchhoff Surface Integral

An acoustic prediction capability for supersonic axisymmetric jets was developed on the basis of OVERFLOW Navier-Stokes CFD (Computational Fluid Dynamics) code of NASA Langley Research Center. Reynolds-averaged turbulent stresses in the flow field are modeled with the aid of Spalart-Allmaras one-equation turbulence model. Appropriate acoustic and outflow boundary conditions were implemented to compute time-dependent acoustic pressure in the nonlinear source-field. Based on the specification of acoustic pressure, its temporal and normal derivatives on the Kirchhoff surface, the near-field and the far-field sound pressure levels are computed via Kirchhoff surface integral, with the Kirchhoff surface chosen to enclose the nonlinear sound source region described by the CFD code. The methods are validated by a comparison of the predictions of sound pressure levels with the available data for an axisymmetric turbulent supersonic (Mach 2) perfectly expanded jet.

Kandula, Max↗