SEARCH · Engineering Papers

Results for “Nonlinear models”

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 1,297 records · Page 72

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Parallel iterative methods for sparse linear and nonlinear equations

As three-dimensional models are gaining importance, iterative methods will become almost mandatory. Among these, preconditioned Krylov subspace methods have been viewed as the most efficient and reliable, when solving linear as well as nonlinear systems of equations. There has been several different approaches taken to adapt iterative methods for supercomputers. Some of these approaches are discussed and the methods that deal more specifically with general unstructured sparse matrices, such as those arising from finite element methods, are emphasized.

Saad, Youcef↗

NASA NTRS1995

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗

NASA NTRS1983

Analysis and test for space shuttle propellant dynamics

This report presents the results of a study to develop an analytical model capable of predicting the dynamic interaction forces on the Shuttle External Tank, due to large amplitude propellant slosh during RTLS separation. The report details low-g drop tower and KC-135 test programs that were conducted to investigate propellant reorientation during RTLS. In addition, the development of a nonlinear finite element slosh model (LAMPS2, two dimensional, and one LAMPS3, three dimensional) is presented. Correlation between the model and test data is presented as a verification of the modeling approach.

Berry, R. L.↗

NASA NTRS2024

Finite Element Simulation of the Direct Energy Deposition using Comsol Multiphysics

Direct Energy Deposition (DED) is an emerging technology extensively employed in metal Additive Manufacturing (AM). Despite its widespread industrial application, mathematical modeling in this domain remains challenging. This complexity arises from the intricate nature of the modeling approach and the nonlinear behavior of material parameters across a broad temperature range. Consequently, experimentalists often resort to a trial-and-error method to achieve structures with desired properties. However, this approach can be time-consuming and may not always yield parts with the requisite characteristics, highlighting the necessity for mathematical modeling to enhance manufacturing success. This study focuses on thin-walled manufacturing with a single bead thickness, adjusting laser parameters to produce such structures. The laser cladding speed, typical for DED manufacturing, is considered to be on the order of centimeters per second. The various wall thicknesses and heights are explored to discern the general characteristics of walls manufactured via this DED approach. Our modeling approach diverges from the conventional DED modeling based on activation-deactivation of predefined mesh domains, commonly implemented in many finite element codes. Instead, a method is proposed that effectively models layer cladding and melt pool dynamics, enabling predictions of the microstructure in the resulting structures. The formation mechanisms of cellular, dendritic columnar, and stray (equiaxed) grains, which arise from the interplay between nucleation and growth from the surface is analyzed. The study also examines the feasibility of microstructure formation to demonstrate the various thermal characteristics inherent in wall manufacturing. Our results illustrate how different process parameters influence the temperature gradient and cooling rate of the molten pool, subsequently affecting the primary dendrite arm spacing (PDAS). This modeling technique allows for the investigation of diverse thermal conditions, facilitating the prediction of microstructure and residual stresses in the manufactured parts.

modeling↗

NASA NTRS1992

Paneling methods with vorticity effects and corrections for nonlinear compressibility

Supersonic panel methods and axisymmetric body-modeling singularity methods are presently combined with corrections for nonlinear flow phenomena to a complete missile, its airbreathing inlets, and wing-body combinations. The computer code LRCDM2 is used as an illustrative example of the methods in question. Attention is given to a preliminary method which employs panels to estimate additive drag and lift acting on supersonic rectangular inlets, as well as to the method used to correct off-body flowfields for the presence of a shock. Examples of missile applications of these methods with the appropriate nonlinear corrections are presented.

Dillenius, Marnix F. E.↗

NASA NTRS1995

Observational and numerical studies of extreme frontal scale contraction

The general objective of this effort is to increase understanding of how frontal scale contraction processes may create and sustain intense mesoscale precipitation along intensifying cold fronts. The five-part project (an expansion of the originally proposed two-part project) employed conventional meteorological data, special mesoscale data, remote sensing measurements, and various numerical models. First an idealized hydrostatic modeling study of the scale contraction effects of differential cloud cover on low-level frontal structure and dynamics was completed and published in a peer-reviewed journal. The second objective was to complete and publish the results from a three dimensional numerical model simulation of a cold front in which differential sensible heating related to cloud coverage patterns was apparently crucial in the formation of a severe frontal squall line. The third objective was to use a nonhydrostatic model to examine the nonlinear interactions between the transverse circulation arising from inhomogeneous cloud cover, the adiabatic frontal circulation related to semi-geostrophic forcing, and diabatic effects related to precipitation processes, in the development of a density current-like microstructure at the leading edge of cold fronts. Although the development of a frontal model that could be used to initialize such a primitive equation model was begun, we decided to focus our efforts instead on a project that could be successfully completed in this short time, due to the lack of prospects for continued NASA funding beyond this first year (our proposal was not accepted for future funding). Thus, a fourth task was added, which was to use the nonhydrostatic model to test tentative hypotheses developed from the most detailed observations ever obtained on a density current (primarily sodar and wind profiler data). These simulations were successfully completed, the findings were reported at a scientific conference, and the results have recently been submitted to a peer-reviewed journal. The fifth objective was to complete the analysis of data collected during the Cooperative Oklahoma Profiler Studies (COPS-91) field project, which was supported by NASA. The analysis of the mesoscale surface and sounding data, Doppler radar imagery, and other remote sensing data from multi frequency wind profiler, microwave radiometer, and the Radio Acoustic Sounding System has been completed. This study is a unique investigation of processes that caused the contraction of a cold front to a microscale zone exhibiting an undular bore-like structure. Results were reported at a scientific conference and are being prepared for publication. In summary, considerable progress has been achieved under NASA funding in furthering our understanding of frontal scale contraction and density current - gravity wave interaction processes, and in utilizing models and remotely sensed data in such studies.

Koch, Steven E.↗

NASA NTRS1993

Delay analysis of combinations of pass transistors and classical logic gates

Combinations of pass transistors and logic gates driving nonlinear capacitive loads are analyzed for the presence of characteristics that will permit easier and more accurate digital logic simulation. Accurate time delay models are developed by studying the nature of the response of simplified circuit models to variations of input waveform rise and fall times and output loading. The nonlinear effects of the CMOS logic devices are minimized to permit easier interpretation of the influence of nonlinear capacitive loads. The performance of a CMOS inverter with a complex nonlinear load consisting of a pass transistor that separates a range of capacitances is compared to the same inverter circuit with a linear capacitive load to develop an understanding of the unique requirements of modeling a nonlinear system. Several methods of modeling the delay of CMOS circuits are reviewed, and a multi-parameter linear model is described. General guidelines for designing CMOS circuits with complex load circuits are developed, emphasizing that the circuit output rise delays and fall delays must be separately analyzed.

Gibson, Jonathan S.↗

NASA NTRS1979

A statistical model for the initial stellar mass function

The nonlinear stages of the fragmentation of a collapsing molecular cloud are modeled by coagulation theory. Several distinct physical processes are discussed, including protostellar fragment coalescence, gas accretion, and binary formation. This work confirms and extends the earlier result of Nakano (1966) that an approximately self-similar limiting form of the mass spectrum develops after several mean initial collision times. An approximate solution to the velocity-averaged coagulation equation is given for an arbitrary power-law dependence of the coalescence rate on mass, with dimensional dependence proportional to m to the power lambda; i.e., the asymptotic mass spectrum varies as m to the -3 lambda/2 power at small masses and cuts off exponentially at large masses, the characteristic mass depending on the number of collision times elapsed. Simple physical arguments suggest that lambda may increase with increasing mass, but is restricted to the range from 2/3 to 4/3. A large fraction of collisions could result in binary formation.

Silk, J.↗

NASA NTRS1984

An elastic potential for the nonlinear response of unidirectional graphite composites

An elastic potential is proposed that is capable of modeling the reversible portion of the observed nonlinear response of unidirectional graphite fiber composites. The model includes both the stiffening stress-strain behavior as well as the softening Poisson's response for loading in the fiber direction. The model is compared with experimental results for Celion 6000/PMR-15 graphite-polyimide.

Pindera, M.-J.↗

NASA NTRS2013

Characterization of Swirl-Venturi Lean Direct Injection Designs for Aviation Gas-Turbine Combustion

Injector geometry, physical mixing, chemical processes, and engine cycle conditions together govern performance, operability and emission characteristics of aviation gas-turbine combustion systems. The present investigation explores swirl-venturi lean direct injection combustor fundamentals, characterizing the influence of key geometric injector parameters on reacting flow physics and emission production trends. In this computational study, a design space exploration was performed using a parameterized swirl-venturi lean direct injector model. From the parametric geometry, 20 three-element lean direct injection combustor sectors were produced and simulated using steady-state, Reynolds-averaged Navier-Stokes reacting computations. Species concentrations were solved directly using a reduced 18-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear κ-ε model. Results demonstrate sensitivities of the geometric perturbations on axially averaged flow field responses. Output variables include axial velocity, turbulent kinetic energy, static temperature, fuel patternation and minor species mass fractions. Significant trends have been reduced to surrogate model approximations, intended to guide future injector design trade studies and advance aviation gas-turbine combustion research.

Combustion Chamber↗

NASA NTRS1995

Parameter identification and modeling of longitudinal aerodynamics

Using a comprehensive flight test database and a parameter identification software program produced at NASA Ames Research Center, a math model of the longitudinal aerodynamics of the Harrier aircraft was formulated. The identification program employed the equation error method using multiple linear regression to estimate the nonlinear parameters. The formulated math model structure adhered closely to aerodynamic and stability/control theory, particularly with regard to compressibility and dynamic manoeuvring. Validation was accomplished by using a three degree-of-freedom nonlinear flight simulator with pilot inputs from flight test data. The simulation models agreed quite well with the measured states. It is important to note that the flight test data used for the validation of the model was not used in the model identification.

Aksteter, J. W.↗

NASA NTRS1982

Design of a helicopter autopilot by means of linearizing transformations

An automatic flight control systems design methods for aircraft that have complex characteristics and operational requirements, such as the powered lift STOL and V/STOL configurations are discussed. The method is effective for a large class of dynamic systems that require multiaxis control and that have highly coupled nonlinearities, redundant controls, and complex multidimensional operational envelopes. The method exploits the possibility of linearizing the system over its operational envelope by transforming the state and control. The linear canonical forms used in the design are described, and necessary and sufficient conditions for linearizability are stated. The control logic has the structure of an exact model follower with linear decoupled model dynamics and possibly nonlinear plant dynamics. The design method is illustrated with an application to a helicopter autopilot design.

Meyer, G.↗

NASA NTRS1978

Nonlinear transient analysis via energy minimization

The formulation basis for nonlinear transient analysis of finite element models of structures using energy minimization is provided. Geometric and material nonlinearities are included. The development is restricted to simple one and two dimensional finite elements which are regarded as being the basic elements for modeling full aircraft-like structures under crash conditions. The results indicate the effectiveness of the technique as a viable tool for this purpose.

Kamat, M. P.↗

NASA NTRS1994

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗

NASA NTRS1972

Digital modeling of human operator dynamics via class of Liapunov functions

A technique is described which is utilized for modeling human operator dynamics. The technique is based upon a model-reference system configuration and a class of Liapunov functions formulated for this purpose which possess variable characteristics. It is shown how such a class of Liapunov functions and their time derivatives is formulated. The crux of the formulation lies in three variable positive definite matrices used for the construction of the class of Liapunov functions. The form and order of these matrices depend upon the form and order of the error differential equation of the human operator model-reference system. These matrices can be modified to include nonlinear functions of the human operator model.

Mekel, R.↗

NASA NTRS1974

Helicopter derivative identification from analytic models and flight test data

Recent results of stability derivative identification from helicopter analytic models and flight test data are presented. Six and nine degree-of-freedom (DOF) linear models are identified from an analytic nonlinear helicopter simulation using a least square technique. The identified models are compared with the convectional partial differentiation method for obtaining derivatives to form the basis for interpretation of derivatives identified from flight data. Six degree-of-freedom models are identified from CH-53A and CH-54B flight data, using an extended Kalman filter modified to process several maneuvers simultaneously. The a priori derivative estimate is obtained by optimal filtering of the data and then using a least square method. The results demonstrate that a six DOF identified model is sufficient to determine the low frequency modes of motion, but a nine DOF rotor/body model is necessary for proper representation of short-term response.

Molusis, J. H.↗

NASA NTRS1982

On the interaction of surface heating anomalies with zonally symmetric and asymmetric atmospheric flows

Models of intermediate complexity have been used to study some aspects of the climatic effects of anomalous heating, but many aspects of the problem have yet to be explored thoroughly. The present study represents a preliminary investigated of the gaps in scientific understanding of the interaction of heating and atmospheric dynamics. The principle research tool is a model of intermediate complexity, including a time-dependent, nonlinear-two-layer quasi-geostrophic model of relatively high horizontal resolution which incorporates simple heating parameterizations. The model is used to examine systematically the interaction of heating arising from anomalies in surface temperature with zonally symmetric and zonally asymmetric flows characterized by different values of static stability and mean vertical wind shear.

Phillips, T. J.↗

NASA NTRS1985

Constitutive modeling for isotropic materials

Constitutive theories were evaluated against a large uniaxial and multiaxial data base that was generated as part of this work. The experimental approach was to determine the constitutive behavior of Rene 80 under a multitude of conditions that are important in the design of gas turbine blades and vanes. The experimental and analytical goals of this program were successfully accomplished. A new multiaxial constitutive model which can represent the complex nonlinear high temperature behavior of Rene 80 was developed. The model was extensively verified on data at several temperatures.

Ramaswamy, V. G.↗