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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in the NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency. In contrast to the legacy FUN3D iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear Courant-Friedrichs-Lewy number (CFL) ramping, HANIM is based upon a hierarchy of modules including preconditioner, generalized conjugate residual, realizability check, nonlinear control, and CFL adaption modules. HANIM performance is systematically compared with the performance of the legacy solver of FUN3D and a baseline solver based on a preconditioner alone. Iterative solutions are compared for three benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-fuselage juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmaras model and its nonlinear extension based on quadratic constitutive relations.

CFD↗

NASA NTRS1976

An efficient procedure for computing partial ionization of a gas in numerical fluid dynamics

A numerical procedure is described that allows efficient solution for the three nonlinear and implicit equations of state relating the equilibrium properties of a partially (and singly) ionized neutral gas. In the computation, the values for the intensive variables, that is pressure, temperature and ionization coefficient, are iteratively obtained with a given accuracy starting from assigned values for energy and volume per unit mass. Iterative computation is performed only when the value for the ionization coefficient lies within the lower and the upper bounds fixed by the required number of significant digits; otherwise the two explicit equations of state for nonionized or for fully ionized gases are employed.

Garribba, S.↗

DOE OSTI2022

Additional Degree of Freedom for WEC Model

'Additional Degree of Freedom for WEC' - WEC-Sim numerical model from RFTS 1 TEAMER project. An increase in wave energy converter (WEC) efficiency requires not only consideration of the nonlinear effects in the WEC dynamics and the power take-off (PTO) mechanisms, but also more integrated treatment of the whole system, i.e., the buoy dynamics, the PTO system, and the control strategy. It results in an optimization formulation that has a nonquadratic and nonstandard cost functional. This model presents the application of real-time nonlinear model predictive controller (NMPC) to two degrees of freedom point absorber type WEC with highly nonlinear PTO characteristics. The nonlinear effects, such as the fluid viscous drag, are also included in the plant dynamics. The controller may be implemented on a real-time target machine, and the WEC device is emulated in real-time using the WECSIM toolbox. Please see the ‘ReadMe’ within the root directory for an explanation of how to set-up and run the model. An article covering the development of this model is attached as "Real-Time Nonlinear Model Predictive Controller for Multiple Degrees of Freedom Wave Energy Converters with Non-Ideal Power Take-Off"

16 TIDAL AND WAVE POWER↗

NASA NTRS1998

A Method for Generating Reduced-Order Linear Models of Multidimensional Supersonic Inlets

Simulation of high speed propulsion systems may be divided into two categories, nonlinear and linear. The nonlinear simulations are usually based on multidimensional computational fluid dynamics (CFD) methodologies and tend to provide high resolution results that show the fine detail of the flow. Consequently, these simulations are large, numerically intensive, and run much slower than real-time. ne linear simulations are usually based on large lumping techniques that are linearized about a steady-state operating condition. These simplistic models often run at or near real-time but do not always capture the detailed dynamics of the plant. Under a grant sponsored by the NASA Lewis Research Center, Cleveland, Ohio, a new method has been developed that can be used to generate improved linear models for control design from multidimensional steady-state CFD results. This CFD-based linear modeling technique provides a small perturbation model that can be used for control applications and real-time simulations. It is important to note the utility of the modeling procedure; all that is needed to obtain a linear model of the propulsion system is the geometry and steady-state operating conditions from a multidimensional CFD simulation or experiment. This research represents a beginning step in establishing a bridge between the controls discipline and the CFD discipline so that the control engineer is able to effectively use multidimensional CFD results in control system design and analysis.

Chicatelli, Amy↗

NASA NTRS1998

Parallel Implicit Algorithms for CFD

The main goal of this project was efficient distributed parallel and workstation cluster implementations of Newton-Krylov-Schwarz (NKS) solvers for implicit Computational Fluid Dynamics (CFD.) "Newton" refers to a quadratically convergent nonlinear iteration using gradient information based on the true residual, "Krylov" to an inner linear iteration that accesses the Jacobian matrix only through highly parallelizable sparse matrix-vector products, and "Schwarz" to a domain decomposition form of preconditioning the inner Krylov iterations with primarily neighbor-only exchange of data between the processors. Prior experience has established that Newton-Krylov methods are competitive solvers in the CFD context and that Krylov-Schwarz methods port well to distributed memory computers. The combination of the techniques into Newton-Krylov-Schwarz was implemented on 2D and 3D unstructured Euler codes on the parallel testbeds that used to be at LaRC and on several other parallel computers operated by other agencies or made available by the vendors. Early implementations were made directly in Massively Parallel Integration (MPI) with parallel solvers we adapted from legacy NASA codes and enhanced for full NKS functionality. Later implementations were made in the framework of the PETSC library from Argonne National Laboratory, which now includes pseudo-transient continuation Newton-Krylov-Schwarz solver capability (as a result of demands we made upon PETSC during our early porting experiences). A secondary project pursued with funding from this contract was parallel implicit solvers in acoustics, specifically in the Helmholtz formulation. A 2D acoustic inverse problem has been solved in parallel within the PETSC framework.

Keyes, David E.↗

NASA NTRS2004

Computation of Sensitivity Derivatives of Navier-Stokes Equations using Complex Variables

Accurate computation of sensitivity derivatives is becoming an important item in Computational Fluid Dynamics (CFD) because of recent emphasis on using nonlinear CFD methods in aerodynamic design, optimization, stability and control related problems. Several techniques are available to compute gradients or sensitivity derivatives of desired flow quantities or cost functions with respect to selected independent (design) variables. Perhaps the most common and oldest method is to use straightforward finite-differences for the evaluation of sensitivity derivatives. Although very simple, this method is prone to errors associated with choice of step sizes and can be cumbersome for geometric variables. The cost per design variable for computing sensitivity derivatives with central differencing is at least equal to the cost of three full analyses, but is usually much larger in practice due to difficulty in choosing step sizes. Another approach gaining popularity is the use of Automatic Differentiation software (such as ADIFOR) to process the source code, which in turn can be used to evaluate the sensitivity derivatives of preselected functions with respect to chosen design variables. In principle, this approach is also very straightforward and quite promising. The main drawback is the large memory requirement because memory use increases linearly with the number of design variables. ADIFOR software can also be cumber-some for large CFD codes and has not yet reached a full maturity level for production codes, especially in parallel computing environments.

Vatsa, Veer N.↗

DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗

DOE OSTI2023

Simulating dielectric spectra: A demonstration of the direct electric field method and a new model for the nonlinear dielectric response

Here we demonstrate a method to compute the dielectric spectra of fluids in molecular dynamics (MD) by directly applying electric fields to the simulation. We obtain spectra from MD simulations with low magnitude electric fields (≈0.01 V/Å) in agreement with spectra from the fluctuation–dissipation method for water and acetonitrile. We examine this method’s trade-off between noise at low field magnitudes and the nonlinearity of the response at higher field magnitudes. We then apply the Booth equation to describe the nonlinear response of both fluids at low frequency (0.1 GHz) and high field magnitude (up to 0.5 V/Å). We develop a model of the frequency-dependent nonlinear response by combining the Booth description of the static nonlinear dielectric response of fluids with the frequency-dependent linear dielectric response of the Debye model. We find good agreement between our model and the MD simulations of the nonlinear dielectric response for both acetonitrile and water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

OC6 Phase Ia Definition Document: Validation of Nonlinear Hydrodynamic Loading on the DeepCwind Semisubmersible

The objective of OC6 Phase I is to investigate the persistent under-prediction (about 20% on average) of the structural loads in the OC5-DeepCwind floating semisubmersible wind system, as observed within the previous OC5 Phase II project. The OC5 results indicate that much of this under-prediction originates from the low-frequency response of the system at the surge/pitch natural frequencies, resulting from nonlinear hydrodynamic loading and/or overprediction of hydrodynamic damping. This work package therefore focuses on better understanding this low-frequency response behavior, the applicability of computational-fluid-dynamics (CFD) models and engineering-level hydrodynamic models for predicting the nonlinear hydrodynamic loading, and the level of uncertainty in the measured response characteristics. To address these objectives, two experimental campaigns were developed and performed in 2017 and 2018 by a sub-group of the OC5 project: Floating test campaign: A simplified OC5-DeepCwind semisubmersible model (no turbine, rigid tower) was moored to three taut-spring lines and subjected to motion decay and wave tests (regular and irregular waves), during which motion amplitudes of the model were measured (see Figure 1a). The purpose of this test campaign was to examine the repeatability of the low-frequency response of the system to wave excitation and determine the level of uncertainty in the response behavior.Constrained test campaign: The same model, but without the tower present, was fixed to a carriage and subjected to towing, forced oscillation and wave tests (regular and irregular waves), during which forces exerted upon the model were measured (see Figure 1b). The purpose of this test campaign was to examine the hydrodynamic loading contributions by separating the wave excitation and radiation forces. Both model test campaigns were performed in the concept basin of the Maritime Research Institute Netherlands (MARIN) in Wageningen, the Netherlands and the first test was within the framework of the MARINET2 project. Both experimental campaigns were examined within Phase Ia of the OC6 project. The purpose of this report is to supply the needed information to participants of the OC6 project to build and simulate the experimental campaigns.

17 WIND ENERGY↗

NASA NTRS1976

Status and prospects of computational fluid dynamics

The use and limitations on using computational aerodynamics in approximating inviscid linear, inviscid nonlinear, vicous time averaged, and viscous time dependent flow past airfoils, wings, and aircraft is reviewed. The current status of two- and three-dimensional time averaged Navier-Stokes equation is discussed and possible applications for the 1980 and 1985 to 1990 period is projected for three-dimensional applications.

Chapman, D. R.↗

NASA NTRS1992

A new method for transonic static aeroelasticity problems

A new method has been developed to calculate the steady flow and structural deformations for fluid/structure interaction problems. The discretized fluid dynamic and structural equations are regarded as a single set of coupled, nonlinear, algebraic equations. The equilibrium solution is directly obtained using Newton's method. The governing equations used for the fluid flow are the two-dimensional Navier-Stokes equations, and a finite-element model is used to represent the structure. This paper describes the analytical method and presents sample calculations demonstrating the technique. The results show rapid convergence and good agreement with experimental data.

Felker, Fort F.↗

DOE OSTI2022

Phenomenological Model of Nonlinear Dynamics and Deterministic Chaotic Gas Migration in Bentonite: Experimental Evidence and Diagnostic Parameters

Understanding gas migration in compacted clay materials, e.g., bentonite and claystone, is important for the design and performance assessment of an engineered barrier system of a radioactive waste repository system, as well as many practical applications. Existing field and laboratory data on gas migration processes in low-permeability clay materials demonstrate the complexity of flow and transport processes, including various types of instabilities, caused by nonlinear dynamics of coupled processes of liquid–gas exchange, dilation, fracturing, fracture healing, etc., which cannot be described by classical models of fluid dynamics in porous media. We here show that the complexity of gas migration processes can be explained using a phenomenological concept of nonlinear dynamics and deterministic chaos theory. To do so, we analyzed gas pressure and gas influx (i.e., input) and outflux (i.e., output), recorded during the gas injection experiment in the compact Mx80-D bentonite sample, and calculated a set of the diagnostic parameters of nonlinear dynamics and chaos, such a global embedding dimension, a correlation dimension, an information dimension, and a spectrum of Lyapunov exponents, as well as plotted 2D and 3D pseudo-phase-space strange attractors, based on the univariate influx and outflux time series data. These results indicate the presence of phenomena of low-dimensional deterministic chaotic behavior of gas migration in bentonite. In particular, during the onset of gas influx in the bentonite core, before the breakthrough, the development of gas flow pathways is characterized by the process of chaotic gas diffusion. After the breakthrough, with inlet-to-outlet movement of gas, the prevailing process is chaotic advection. During the final phase of the experiment, with no influx to the sample, the relaxation pattern of gas outflux is resumed back to a process of chaotic diffusion. The types of data analysis and a proposed phenomenological model can be used to establish the basic principles of experimental data-gathering, modeling predictions, and a research design.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Real-Time Nonlinear Model Predictive Controller for Multiple Degrees of Freedom Wave Energy Converters with Non-Ideal Power Take-Off

An increase in wave energy converter (WEC) efficiency requires not only consideration of the nonlinear effects in the WEC dynamics and the power take-off (PTO) mechanisms, but also more integrated treatment of the whole system, i.e., the buoy dynamics, the PTO system, and the control strategy. It results in an optimization formulation that has a nonquadratic and nonstandard cost functional. This article presents the application of real-time nonlinear model predictive controller (NMPC) to two degrees of freedom point absorber type WEC with highly nonlinear PTO characteristics. The nonlinear effects, such as the fluid viscous drag, are also included in the plant dynamics. The controller is implemented on a real-time target machine, and the WEC device is emulated in real-time using the WECSIM toolbox. The results for the successful performance of the design are presented for irregular waves under linear and nonlinear hydrodynamic conditions.

Haider, Ali S. (ORCID:0000000246513157)↗

DOE OSTI2022

Midfidelity model verification for a point-absorbing wave energy converter with linear power take-off

In the preliminary design stage of a wave energy converter (WEC), researchers need fast and reliable simulation tools. High-fidelity numerical models are usually employed to study the wave-structure interaction, but the computational cost is demanding. As an alternative, midfidelity models can provide simulations in the order of real time. In this study, we operate Uppsala University’s WEC in a relatively mild sea state and model it using WEC-Sim. The model is verified based on OpenFOAM simulations. To analyze the ability of the midfidelity model to capture WEC dynamics, we investigate the system separately with 1, 2, and 3 degrees of freedom. We examine the contribution of viscous phenomena, and study both linear and weakly nonlinear solutions provided by WEC-Sim. Our results indicate that the viscous effects can be neglected in heave and surge motion, but not for pitch. We also find that the weakly nonlinear WEC-Sim solution successfully agrees with the computational fluid dynamics, whereas the linear solution could suggest misleading results.

16 TIDAL AND WAVE POWER↗

NASA NTRS2016

A Nonlinear Modal Aeroelastic Solver for FUN3D

A nonlinear structural solver has been implemented internally within the NASA FUN3D computational fluid dynamics code, allowing for some new aeroelastic capabilities. Using a modal representation of the structure, a set of differential or differential-algebraic equations are derived for general thin structures with geometric nonlinearities. ODEPACK and LAPACK routines are linked with FUN3D, and the nonlinear equations are solved at each CFD time step. The existing predictor-corrector method is retained, whereby the structural solution is updated after mesh deformation. The nonlinear solver is validated using a test case for a flexible aeroshell at transonic, supersonic, and hypersonic flow conditions. Agreement with linear theory is seen for the static aeroelastic solutions at relatively low dynamic pressures, but structural nonlinearities limit deformation amplitudes at high dynamic pressures. No flutter was found at any of the tested trajectory points, though LCO may be possible in the transonic regime.

Goldman, Benjamin D.↗

NASA NTRS2000

A Low Order and a High Order Solution for a Converging-Diverging Nozzle Problem

The Computational Aeroacoustics Workshops on Benchmark problems are conducted in order to generate solutions with a variety of methods for problems that are chosen both to be representative of significant and relevant acoustic phenomena. It is generally recognized that CAA applications require the resolution in time and in space of solution details over a significantly broader range of scales than typical computational fluid dynamics applications. These requirements are forcing the development of new methodologies for CAA applications. This paper presents numerical results from a third and a seventh order algorithm for the propagation of an acoustic signal through a converging-diverging nozzle. These two algorithms are members of a new family of high accuracy methods that have the same order accuracy in both space and time, and are an extension of previous work for linearized Euler equations to fully nonlinear time dependent problems. The simulations are all with the fully nonlinear quasi-1D Euler equations for the total solution, which includes both the fluid dynamics and the acoustics. The acoustic solution is obtained from the time dependent nonlinear solution by subtracting the steady solution.

Goodrich, John W.↗

NASA NTRS1983

Spectral methods for partial differential equations

Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid dynamical applications are emphasized.

Hussaini, M. Y.↗

DOE OSTI2020

Bichromatic Wave Selection for Validation of the Difference-Frequency Transfer Function for the OC6 Validation Campaign: Preprint

The focus of the Offshore Code Comparison Collaboration, Continuation, with Correlation and unCertainity (OC6) project, which operates under the International Energy Agency (IEA) Wind Task 30, is to refine the accuracy of engineering tools used to design offshore wind turbines. In support of this work, a new validation campaign is being developed that seeks to better understand the nonlinear wave loading that excites floating wind systems at their low-frequency, rigid-body modes in surge and pitch. The validation data will be used in a three-way validation between simplified engineering tools and higher-fidelity tools, such as computational fluid dynamics (CFD). Irregular wave spectrums, which are traditionally used to examine the nonlinear wave interaction with offshore structures, are too computationally expensive to be simulated in CFD tools, and so bichromatic wave cases will be used instead. This paper reviews the process used to choose the bichromatic wave pairs to be used in the campaign to validate the second-order difference-frequency quadratic and potential loads at the surge and pitch natural frequencies of a floating semisubmersible.

17 WIND ENERGY↗