Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear equations”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Material characterization of structural adhesives in the lap shear mode

A general method for characterizing structual adhesives in the bonded lap shear mode is proposed. Two approaches in the form of semiempirical and theoretical approaches are used. The semiempirical approach includes Ludwik's and Zhurkov's equations to describe respectively, the failure stresses in the constant strain rate and constant stress loading modes with the inclusion of the temperature effects. The theoretical approach is used to describe adhesive shear stress-strain behavior with the use of viscoelastic or nonlinear elastic constitutive equations. Two different model adhesives are used in the single lap shear mode with titanium adherends. These adhesives (one of which was developed at NASA Langley Research Center) are currently considered by NASA for possible aerospace applications. Use of different model adhesives helps in assessment of the generality of the method.

Sancaktar, E.↗
NASA NTRS1983

Material characterization of structural adhesives in the lap shear mode

A general method for characterizing structural adhesives in the bonded lap shear mode is proposed. Two approaches in the form of semi-empirical and theoretical approaches are used. The semi-empirical approach includes Ludwik's and Zhurkov's equations to describe respectively, the failure stresses in the constant strain rate and constant stress loading modes with the inclusion of the temperature effects. The theoretical approach is used to describe adhesive shear stress-strain behavior with the use of viscoelastic or nonlinear elastic constitutive equations. Three different model adhesives are used in the simple lap shear mode with titanium adherends. These adhesives (one of which was developed at NASA Langley Research Center) are currently considered by NASA for possible aerospace applications. Use of different model adhesives helps in assessment of the generality of the method.

Schenck, S. C.↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2022

Nonlinear spectral analysis of ion acoustic solitons arising from a streaming charged object using the numerical inverse scattering transform

Observations of nonlinear coherent plasma phenomena in spacecraft and terrestrial experiments often rely on visual identification of solitary modes because nonlinear coherent modes are broadband in a Fourier transform based analysis. We implement an alternative spectral decomposition known as the inverse scattering transform and demonstrate its ability to successfully isolate nonlinear modes on the homogeneous and forced Korteweg–De Vries equation which models these nonlinear modes. We also demonstrate for the first time that this decomposition is useful when forcing is applied. This is because the stable modes generated by localized forcing are similar to the homogeneous solutions where the inverse scattering transform is a rigorous decomposition. This spectral technique is then applied to simulations of ion acoustic waves generated by a 10 mm spherical debris object interacting with the ionospheric plasma orbiting at 2000 km altitude. The algorithm is found to successfully detect the resultant solitons. This demonstrates the feasibility of using this spectral technique as a real time analysis tool for screening spacecraft data for nonlinear solitary modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Pinching solutions of slender cylindrical jets

Simplified equations for slender jets are derived for a circular jet of one fluid flowing into an ambient second fluid, the flow being confined in a circular tank. Inviscid flows are studied which include both surface tension effects and Kelvin-Helmholtz instability. For slender jets a coupled nonlinear system of equations is found for the jet shape and the axial velocity jump across it. The equations can break down after a finite time and similarity solutions are constructed, and studied analytically and numerically. The break-ups found pertain to the jet pinching after a finite time, without violation of the slender jet ansatz. The system is conservative and admissible singular solutions are those which conserve the total energy, mass, and momentum. Such solutions are constructed analytically and numerically, and in the case of vortex sheets with no surface tension certain solutions are given in closed form.

Papageorgiou, Demetrios T.↗
NASA NTRS1994

Progress in modeling atmospheric propagation of sonic booms

The improved simulation of sonic boom propagation through the real atmosphere requires greater understanding of how the transient acoustic pulses popularly termed sonic booms are affected by humidity and turbulence. A realistic atmosphere is invariably somewhat turbulent, and may be characterized by an ambient fluid velocity v and sound speed c that vary from point to point. The absolute humidity will also vary from point to point, although possibly not as irregularly. What is ideally desired is a relatively simple scheme for predicting the probable spreads in key sonic boom signature parameters. Such parameters could be peak amplitudes, rise times, or gross quantities obtainable by signal processing that correlate well with annoyance or damage potential. The practical desire for the prediction scheme is that it require a relatively small amount of knowledge, possibly of a statistical nature, concerning the atmosphere along, the propagation path from the aircraft to the ground. The impact of such a scheme, if developed, implemented, and verified, would be that it would give the persons who make planning decisions a tool for assessing the magnitude of environmental problems that might result from any given overflight or sequence of overflights. The technical approach that has been followed by the author and some of his colleagues is to formulate a hierarchy of simple approximate models based on fundamental physical principles and then to test these models against existing data. For propagation of sonic booms and of other types of acoustic pulses in nonturbulent model atmospheres, there exists a basic overall theoretical model that has evolved as an outgrowth of geometrical acoustics. This theoretical model depicts the sound as propagating within ray tubes in a manner analogous to sound in a waveguide of slowly varying cross-section. Propagation along the ray tube is quasi-one-dimensional, and a wave equation for unidirectional wave propagation is used. A nonlinear term is added to this equation to account for nonlinear steepening, and the formulation has been carried through to allow for spatially varying sound speed, ambient density, and ambient wind velocities. The model intrinsically neglects diffraction, so it cannot take into account what has previously been mentioned in the literature as possibly important mechanisms for turbulence-related distortion. The model as originally developed could predict an idealized N-waveform which often agrees with data in terms of peak amplitude and overall positive phase duration. It is possible, moreover, to develop simple methods based on the physics of relaxation processes for incorporating molecular relaxation into the quasi-one-dimensional model of nonlinear propagation along ray tubes.

Pierce, Allan D.↗
NASA NTRS1999

Photochemistry of Pluto's Atmosphere

This work include studies of two problems: (1) Modeling thermal balance, structure. and escape processes in Pluto's upper atmosphere. This study has been completed in full. A new method, of analytic solution for the equation of hydrodynamic flow from in atmosphere been developed. It was found that the ultraviolet absorption by methane which was previously ignored is even more important in Pluto's thermal balance than the extreme ultraviolet absorption by nitrogen. Two basic models of the lower atmosphere have been suggested, with a tropopause and a planetary surface at the bottom of the stellar occultation lightcurve, respectively, Vertical profiles, of temperature, density, gas velocity, and the CH4 mixing ratio have been calculated for these two models at low, mean, and high solar activity (six models). We prove that Pluto' " s atmosphere is restricted to 3060-4500 km, which makes possible a close flyby of future spacecraft. Implication for Pluto's evolution have also been discussed. and (2) Modeling of Pluto's photochemistry. Based on the results of (1), we have made some changes in the basic continuity equation and in the boundary conditions which reflect a unique can of hydrodynamic escape and therefore have not been used in modeling of other planetary atmospheres. We model photochemistry of 44 neutral and 23 ion species. This work required solution of a set of 67 second-order nonlinear ordinary differential equations. Two models have been developed. Each model consists of the vertical profiles for 67 species, their escape and precipitation rates. These models predict the chemical structure and basic chemical processes in the current atmosphere and possible implication of these processes for evolution. This study has also been completed in full.

Krasnopolsky, Vladimir A.↗
NASA NTRS1986

Self-adaptive difference method for the effective solution of computationally complex problems of boundary layer theory

An implicit difference procedure for the solution of equations for a chemically reacting hypersonic boundary layer is described. Difference forms of arbitrary error order in the x and y coordinate plane were used to derive estimates for discretization error. Computational complexity and time were minimized by the use of this difference method and the iteration of the nonlinear boundary layer equations was regulated by discretization error. Velocity and temperature profiles are presented for Mach 20.14 and Mach 18.5; variables are velocity profiles, temperature profiles, mass flow factor, Stanton number, and friction drag coefficient; three figures include numeric data.

Schoenauer, W.↗
DOE OSTI2021

Fully Implicit Monte Carlo (FIMC) (September 2019-September 2020 Progress Report)

This research is focused on applying the NLem method to the full thermal radiative transfer (TRT) equations in the context of a Monte Carlo solution strategy. A framework has been developed for obtaining a solution to the TRT equations using the nonlinear elimination method (NLem) method. Current efforts are focused of the development and verification of the required Monte Carlo algorithm, diffusion algorithm, and material energy solver. Once all three pieces have been verified, the algorithms will then be coupled into a fully functioning method for applying NLem to the TRT equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Aerodynamic Influence Coefficient Computations Using Euler/Navier-Stokes Equations on Parallel Computers

Modem design requirements for an aircraft push current technologies used in the design process to their limit or sometimes require more advanced technologies to meet the requirement. New design requirements always demand to improve the operational performance. Accurate prediction of aerodynamic coefficients is essential to improve the performance. For example, in the design of an advanced subsonic civil transport, since the fluid flow at transonic regime shows strong nonlinearities, high fidelity equations, such as the Euler or Navier-Stokes equations predict flow characteristics more accurately than the linear aerodynamics, which are widely used in the current design process However, high fidelity flow equations are computationally expensive and require an order of magnitude longer time to obtain aerodynamic coefficients required in the design. Parallel computing is one possibility to cut down the computational turn-around time in using high fidelity equations so that high fidelity equations would be incorporated into the design process. By doing so, high fidelity equations would be used in the routine design process. This work will demonstrate the feasibility of using high fidelity flow equations in a design process by computing aerodynamic influence coefficients of a wing-body-empennage configuration on a multiple-instruction, multiple-data parallel computer.

Byun, Chansup↗
NASA NTRS1983

Transient responses of phosphoric acid fuel cell power plant system

An analytical and computerized study of the steady state and transient response of a phosphoric acid fuel cell (PAFC) system was completed. Parametric studies and sensitivity analyses of the PAFC system's operation were accomplished. Four non-linear dynamic models of the fuel cell stack, reformer, shift converters, and heat exchangers were developed based on nonhomogeneous non-linear partial differential equations, which include the material, component, energy balance, and electrochemical kinetic features. Due to a lack of experimental data for the dynamic response of the components only the steady state results were compared with data from other sources, indicating reasonably good agreement. A steady state simulation of the entire system was developed using, nonlinear ordinary differential equations. The finite difference method and trial-and-error procedures were used to obtain a solution. Using the model, a PAFC system, that was developed under NASA Grant, NCC3-17, was improved through the optimization of the heat exchanger network. Three types of cooling configurations for cell plates were evaluated to obtain the best current density and temperature distributions. The steady state solutions were used as the initial conditions in the dynamic model. The transient response of a simplified PAFC system, which included all of the major components, subjected to a load change was obtained. Due to the length of the computation time for the transient response calculations, analysis on a real-time computer was not possible. A simulation of the real-time calculations was developed on a batch type computer. The transient response characteristics are needed for the optimization of the design and control of the whole PAFC system. All of the models, procedures and simulations were programmed in Fortran and run on IBM 370 computers at Cleveland State University and the NASA Lewis Research Center.

Lu, Cheng-Yi↗
NASA NTRS1995

High-order essentially non-oscillatory methods for computational aeroacoustics

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. In addition, it is desirable to highly resolve discontinuities that can develop in the solutions of the Euler or Navier-Stokes equations. The class of essentially non-oscillatory (ENO) shock-capturing schemes has been designed to have both of these properties. The dual capacity of ENO schemes for high-order accuracy and non-oscillatory shock-capturing is achieved through the use of adaptive stenciling, which makes these schemes highly nonlinear. These schemes are briefly described and referenced herein. A fourth-order algorithm is then applied to the solution of an acoustic wave in a quasi-one-dimensional converging-diverging nozzle.

Casper, Jay↗
NASA NTRS1983

New discretization and solution techniques for incompressible viscous flow problems

Several topics arising in the finite element solution of the incompressible Navier-Stokes equations are considered. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. The role of artificial viscosity in viscous flow calculations is studied, emphasizing work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some modifications are mentioned.

Gunzburger, M. D.↗
NASA NTRS1986

Modification of a variational objective analysis model for new equations for pressure gradient and vertical velocity in the lower troposphere and for spatial resolution and accuracy of satellite data

Since late 1982 NASA has supported research to develop a numerical variational model for the diagnostic assimilation of conventional and space-based meteorological data. In order to analyze the model components, four variational models are defined dividing the problem naturally according to increasing complexity. The first of these variational models (MODEL I), the subject of this report, contains the two nonlinear horizontal momentum equations, the integrated continuity equation, and the hydrostatic equation. This report summarizes the results of research (1) to improve the way the large nonmeteorological parts of the pressure gradient force are partitioned between the two terms of the pressure gradient force terms of the horizontal momentum equations, (2) to generalize the integrated continuity equation to account for variable pressure thickness over elevated terrain, and (3) to introduce horizontal variation in the precision modulus weights for the observations.

Achtemeier, G. L.↗
DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

New discretization and solution techniques for incompressible viscous flow problems

This paper considers several topics arising in the finite element solution of the incompressible Navier-Stokes equations. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. Following this, the role of artificial viscosity in viscous flow calculations is studied, emphasizing recent work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some recent modifications are mentioned.

Gunzburger, M. D.↗
NASA NTRS1985

Elastohydrodynamic lubrication calculations used as a tool to study scuffing

A new Reynolds equation was developed that takes into account the nonlinear viscous behavior of the fluid. The new Reynolds equation considers the nonlinear viscous fluid model of Eyring, the equilibrium equation, the constant mass flow, and the kinematic boundary condition. The new Reynolds equation and the elasticity equation are solved simultaneously by using a system approach and a Newton-Raphson technique. Comparisons are made with results obtained from the classical Reynolds equation. The effects of sliding speed and introducing a bump or a groove within the conjunction are studied. Results are shown for both moderate and heavy loads.

Houpert, L. G.↗
NASA NTRS2004

A Constitutive Equation for Stratospheric Balloon Materials

The selection of a suitable material for use as a reliable stratospheric balloon gas barrier and structural component is based on a variety of properties. Due to a more desirable combination of properties, the low density polyethylene that has been used for the last half century has been replaced during the last decade by linear low density polyethylene (LLDPE). This paper describes the effort to characterize the time dependent properties of a 38 micron coextrusion of LLDPE. The nonlinear viscoelastic constitutive equation presented may be used to accurately describe the creep and/or relaxation of this film when subjected to a biaxial state of stress, such as might be required for an extended balloon flight. Recent laboratory data have been used to mod@ an existing model of LLDPE to account for differences caused by the coextrusion process. The new model will facilitate structural design optimization and reliability assessment, and may be further utilized as a predictive tool to benefit in-flight operations. Current structural analysis tech&ques based on linear elastic properties have predicted stresses in excess of those which would actually exist.

Rand, J. L.↗