Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Mitigation of stimulated Raman scattering in the kinetic regime by external magnetic fields

Here, we show via particle-in-cell simulations that small normalized magnetic fields (ω c /ω p << 1) can significantly modify the evolution of backward stimulated Raman scattering (SRS) in the kinetic regime due to the enhanced dissipation of nonlinear electron plasma waves propagating perpendicular to magnetic fields. A magnetic field applied perpendicularly to the electron plasma wave (and driving light wave) increases the SRS threshold for kinetic inflation and decreases the amount of reflectivity when SRS is driven significantly above threshold. Analysis indicates that this arises because trapped electrons are accelerated as they surf across the wave, leading to the continual dissipation of the electron plasma waves over a wider range of wave amplitudes. The reduction in SRS reflectivity is most significant for a purely perpendicular field, although reduction also occurs for other angles; a parallel field can slightly increase single-speckle SRS but decreases multispeckle SRS. These simulations demonstrate the significance of magnetic-field contributions to nonlinear electron plasma wave damping with respect to nonlinear parametric decay instabilities; the simulation parameters are directly relevant for SRS in inertial confinement fusion devices and indicate that approximately 30 tesla magnetic fields might significantly reduce SRS backscatter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
NASA NTRS1957

Flight Investigation at Low Angles of Attack to Determine the Longitudinal Stability and Control Characteristics of a Cruciform Canard Missile Configuration with a Low-Aspect-Ratio Wing and Blunt Nose at Mach Numbers from 1.2 to 2.1

A full-scale rocket-powered model of a cruciform canard missile configuration with a low-aspect-ratio wing and blunt nose has been flight tested by the Langley Pilotless Aircraft Research Division. Static and dynamic longitudinal stability and control derivatives of this interdigitated canard-wing missile configuration were determined by using the pulsed-control technique at low angles of attack and for a Mach number range of 1.2 to 2.1. The lift-curve slope showed only small nonlinearities with changes in control deflection or angle of attack but indicated a difference in lift-curve slope of approximately 7 percent for the two control deflections of delta = 3.0 deg and delta = -0.3 deg. The large tail length of the missile tested was effective in producing damping in pitch throughout the Mach number range tested. The aerodynamic-center location was nearly constant with Mach number for the two control deflections but was shown to be less stable with the larger control deflection. The increment of lift produced by the controls was small and positive throughout the Mach number range tested, whereas the pitching moment produced by the controls exhibited a normal trend of reduced effectiveness with increasing Mach number. The effectiveness of the controls in producing angle of attack, lift, and pitching moment was good at all Mach numbers tested.

Brown, C. A., Jr.↗
NASA NTRS1972

Analysis of an electrohydraulic aircraft control surface servo and comparison with test results

An analysis of an electrohydraulic aircraft control-surface system is made in which the system is modeled as a lumped, two-mass, spring-coupled system controlled by a servo valve. Both linear and nonlinear models are developed, and the effects of hinge-moment loading are included. Transfer functions of the system and approximate literal factors of the transfer functions for several cases are presented. The damping action of dynamic pressure feedback is analyzed. Comparisons of the model responses with results from tests made on a highly resonant rudder control-surface servo indicate the adequacy of the model. The effects of variations in hinge-moment loading are illustrated.

Edwards, J. W.↗
NASA NTRS1974

Estimation of structural parameters from dynamic test data

Parameters for a lumped linear model approximating a distributed elastic structure are determined from dynamic test data comprised of several mode shapes and frequencies. Measurement errors, nonlinear response, and nonmeasurable quantities such as mode slope components are accommodated. Some mass and stiffness parameters may be known accurately, whereas the remainder are to be estimated. The method entails minimizing a quadratic function of the difference between corresponding modes and frequencies of the theoretical model and the test specimen. This technique was applied to some actual vibration test data, and the special techniques that are required to overcome convergence problems are described.

Hall, B. M.↗
NASA NTRS1977

The relative diffusion of strong magnetic fields and tenuous gases

General equations are given for treating the diffusion of a magnetic field in a fluid with electrical conductivity and a resistive diffusion coefficient. Consequences of nonnegligible gas motion are investigated by considering the one-dimensional problem where a tenuous isothermal gas with pressure expressed as a function of x-coordinate and time is enmeshed in a magnetic field in the y-direction. Approximate equations are obtained for the case of weak fluctuations in gas or field pressure, and the properties of a nonlinear diffusion equation are illustrated with the example of an initial thin slab of gas confined by a uniform field on either side. The advance of a gas front into a nonvanishing gas density is illustrated by using the progressive wave solution. Application of the results to the evolution of streamers and filaments observed in the solar corona is briefly discussed

Parker, E. N.↗
NASA NTRS1982

Development of a reduced basis technique for transient thermal analysis

A technique to reduce the degrees of freedom in static and dynamic problems, the reduced basis method, is described. The method combines the classical Rayleigh-Ritz approximation with contemporary finite element methods to retain modeling versatility as the degrees of freedom are reduced. Applications to a nonlinear dynamic response problem are discussed efforts to apply the method to nonlinear transient thermal response problems are summarized. The selection of basis vectors for reducing the system of equations is addressed.

Shore, C. P.↗
NASA NTRS1983

Multiple buoyancy driven flows in a vertical cylinder heated from below

The structure of axisymmetric buoyancy-driven convection in a vertical cylinder heated from below is probed by finite element solution of the Boussinesq equations coupled with computed-implemented perturbation techniques for detecting and tracking multiple flows and for determining flow stability. Results are reported for fluids with Prandtl number of one and for cylinders with aspect ratio (Lambda) (defined as the height to radius of the cylinder) between 0.5 and 2.25. Extensive calculations of the neutral stability curve for the static solution and of the nonlinear motions along the bifurcating flow families show a continuous evolution of the primary cellular motion from a single toroidal cell to two and three cells nested radially in the cylinder, instead of the sharp transitions found for a cylinder with shear-free sidewalls. The smooth transitions in flow structure with Rayleigh number and lambda are explained by nonlinear connectivity between the first two bifurcating flow families formed either by a secondary bifurcation point for Lambda or = Lambda * approximately 0.80 or by a limit point for Lambda Lambda *. The transition between these two modes may be described by the theory of multiple limit point bifurcation.

Yamaguchi, Y.↗
NASA NTRS1986

Mass fractionation during transonic escape and implications for loss of water from Mars and Venus

Hydrodynamic escape of hydrogen from a planetary atmosphere can remove heavier gases as well as hydrogen, provided that the escape rate is sufficiently large. Analytic approximations for the degree of mass fractionation of a trace species during hydrodynamic escape are compared with accurate numerical solutions for the case of transonic outflow. The analytic approximations are most accurate when the ratio of molecular weights of the heavier and lighter constituents is large so that nonlinear terms in the momentum equation for the heavy constituent become small. The simplest analytic formula is readily generalized to the case where a heavy constituent is also a major species. Application of the generalized formula to hypothetical episodes of hydrodynamic escape from Venus and Mars suggests that both hydrogen and oxygen could have escaped; thus, substantial quantities of water may have been lost without the need to oxidize large amounts of the crust.

Zahnle, Kevin J.↗
NASA NTRS1987

Model-size reduction for the non-linear dynamic analysis of quasi-symmetric structures

A numerical technique is developed to reduce the size of models describing the nonlinear dynamic response of quasi-symmetric structures (i.e., structures with unsymmetric geometry). The response vectors of the structure are approximated by a linear combination of the symmetric and antisymmetric vectors at each time step. The mathematical formulation and numerical implementation of the method are described in detail, and results for a shallow laminated anisotropic panel of quadrilateral planform are presented in graphs and normalized contour plots.

Noor, Ahmed K.↗
NASA NTRS1989

A hybrid perturbation-Galerkin method for differential equations containing a parameter

A two-step hybrid perturbation-Galerkin method to solve a variety of differential equations which involve a parameter is presented and discussed. The method consists of: (1) the use of a perturbation method to determine the asymptotic expansion of the solution about one or more values of the parameter; and (2) the use of some of the perturbation coefficient functions as trial functions in the classical Bubnov-Galerkin method. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is illustrated first with a simple linear two-point boundary value problem and is then applied to a nonlinear two-point boundary value problem in lubrication theory. The results obtained from the hybrid method are compared with approximate solutions obtained by purely numerical methods. Some general features of the method, as well as some special tips for its implementation, are discussed. A survey of some current research application areas is presented and its degree of applicability to broader problem areas is discussed.

Geer, James F.↗
NASA NTRS1999

On the Alignment of Strain, Vorticity and Scalar Gradient in Turbulent, Buoyant, Nonpremixed Flames

The alignment of vorticity and scalar gradient with the eigendirections of the rate of strain tensor is investigated in turbulent buoyant nonpremixed horizontal and vertical flames. The uniqueness of a buoyant nonpremixed flame is that it contains regions with distinct alignment characteristics. The strain-enstrophy angle Psi is used to identify these regions. Examination of the vorticity field and the vorticity production in these different regions indicates that Psi and consequently the alignment properties near the flame surface identified by the mixture fraction band F approximately equals F(sub st) differ from those in the fuel region, F > F(sub st) and the oxidizer region, F < F(sub st). The F approximately equals F(sub st) band shows strain-dominance resulting in vorticity/alpha alignment while F > F(sub st) (and F < F(sub st) for the vertical flame) band(s) show(s) vorticity/beta alignment. The implication of this result is that the scalar dissipation, epsilon(sub F), attains its maximum value always near F approximately equals F(sub st). These results are also discussed within the framework of recent dynamical results [Galanti et al., Nonlinearity 10, 1675 (1997)] suggesting that the Navier-Stokes equations evolved towards an attracting solution. It is shown that the properties of such an attracting solution are also consistent with our results of buoyant turbulent nonpremixed flames.

Boratav, O. N.↗
NASA NTRS2014

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗
NASA NTRS2016

Optic Nerve Sheath Mechanics in VIIP Syndrome

Visual Impairment Intracranial Pressure (VIIP) syndrome is a major concern in current space medicine research. While the exact pathology of VIIP is not yet known, it is hypothesized that the microgravity-induced cephalad fluid shift increases intracranial pressure (ICP) and drives remodeling of the optic nerve sheath. To investigate this possibility, we are culturing optic nerve sheath dura mater samples under different pressures and investigating changes in tissue composition. To interpret results from this work, it is essential to first understand the biomechanical response of the optic nerve sheath dura mater to loading. Here, we investigated the effects of mechanical loading on the porcine optic nerve sheath.Porcine optic nerves (number: 6) were obtained immediately after death from a local abattoir. The optic nerve sheath (dura mater) was isolated from the optic nerve proper, leaving a hollow cylinder of connective tissue that was used for biomechanical characterization. We developed a custom mechanical testing system that allowed for unconfined lengthening, twisting, and circumferential distension of the dura mater during inflation and under fixed axial loading. To determine the effects of variations in ICP, the sample was inflated (0-60 millimeters Hg) and circumferential distension was simultaneously recorded. These tests were performed under variable axial loads (0.6 grams - 5.6 grams at increments of 1 gram) by attaching different weights to one end of the dura mater. Results and Conclusions: The samples demonstrated nonlinear behavior, similar to other soft connective tissue (Figure 1). Large increases in diameter were observed at lower transmural pressures (approximately 0 to 5 millimeters Hg), whereas only small diameter changes were observed at higher pressures. Particularly interesting was the existence of a cross-over point at a pressure of approximately 11 millimeters Hg. At this pressure, the same diameter is obtained for all axial loads applied to the tissue; i.e., as the axial load is varied, the diameter of the dura mater remains constant. This cross-over in the pressure-diameter curves occurred in all optic nerve sheaths that were tested, and may correspond with in vivo ICP levels for pigs. These data suggest that diameter of the dura mater of the optic nerve remains nearly constant in vivo despite being stretched axially. This may be a homeostatic mechanism aimed at maintaining target stresses/strains on the cells in the dura mater, and deviations from these stresses may play an important role in optic nerve sheath remodeling. Future studies will involve subjecting the dura mater to varying pressures and axial tensions for extended periods of time, while monitoring changes in the biomechanical properties. The data can then be used to study the effects of changes in ICP on the remodeling of the dura mater.

biomechanics↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
DOE OSTI2024

Development of Segregated Thermal-Hydraulics Solvers in MOOSE

The simulation of fluid flows is an essential part of the design and analysis of nuclear systems. Algorithms able to simulate flows at different fidelity levels are available in the Multiphysics Object-Oriented Simulation Environment (MOOSE) and MOOSE-based applications such as Pronghorn \cite{novak2018pronghorn}, Pronghorn-Subchannel, RELAP-7, and SAM. Currently, significant effort is being invested in the development of coarse-mesh Computational Fluid Dynamics (CFD) capabilities within MOOSE and Pronghorn for the simulation of Generation IV nuclear reactors. Traditionally, the solution algorithms in MOOSE have relied on Newton or quasi-Newton methods (such as the preconditioned Jacobian-free Newton-Krylov method) where residuals and Jacobians (or approximations thereof) are constructed. Both Newton and quasi-Newton methods require the solution of a linear system at each nonlinear Newton iteration with the Jacobian as the system matrix. The Jacobian contains blocks originating from all variables in the problem (i.e., for thermal-hydraulics at least pressure, velocities, and temperature). Due to the formulation of the problem in a general multiphysics setting on unstructured mesh, creating a good preconditioner for the linear system can be challenging, thus many fluid applications have utilized direct solver-based methods such as LU factorization. However, with increasing system size and complexity in multi-dimensional problems, the direct solution of linear systems becomes computationally expensive both in execution time and and memory. For this reason, recent effort has focused on adapting segregated solution algorithms for CFD problems in MOOSE. These algorithms use fixed-point iteration between segregated systems whose assembly and preconditioning are easier those of the monolithic system. Initial results show that the segregated solution algorithm outperforms the monolithic approach in terms of memory usage and for large 3D problems in terms of CPU time as well.

42 ENGINEERING↗
NASA NTRS2002

Analysis of Three-Dimensional, Nonlinear Development of Wave-Like Structure in a Compressible Round Jet

An analysis of the nonlinear development of the large-scale structures or instability waves in compressible round jets was conducted using the integral energy method. The equations of motion were decomposed into two sets of equations; one set governing the mean flow motion and the other set governing the large-scale structure motion. The equations in each set were then combined to derive kinetic energy equations that were integrated in the radial direction across the jet after the boundary-layer approximations were applied. Following the application of further assumptions regarding the radial shape of the mean flow and the large structures, equations were derived that govern the nonlinear, streamwise development of the large structures. Using numerically generated mean flows, calculations show the energy exchanges and the effects of the initial amplitude on the coherent structure development in the jet.

Dahl, Milo D.↗
NASA NTRS1993

Active control of acoustic pressure fields using smart material technologies

An overview describing the use of piezoceramic patches in reducing noise in a structural acoustics setting is presented. The passive and active contributions due to patches which are bonded to an Euler-Bernoulli beam or thin shell are briefly discussed and the results are incorporated into a 2-D structural acoustics model. In this model, an exterior noise source causes structural vibrations which in turn lead to interior noise as a result of nonlinear fluid/structure coupling mechanism. Interior sound pressure levels are reduced via patches bonded to the flexible boundary (a beam in this case) which generate pure bending moments when an out-of-phase voltage is applied. Well-posedness results for the infinite dimensional system are discussed and a Galerkin scheme for approximating the system dynamics is outlined. Control is implemented by using linear quadratic regulator (LQR) optimal control theory to calculate gains for the linearized system and then feeding these gains back into the nonlinear system of interest. The effectiveness of this strategy for this problem is illustrated in an example.

Banks, H. T.↗