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.

SourceResearch connectionFrom yearThrough year

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

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Thermosolutal convection during dendritic solidification of a binary alloy

The FEM model of Heinrich (1988) is presently extended to the case of dendritic vertical solidification under steady-state conditions. An analysis is conducted of the nonlinear governing equations for linear stability, and FEM calculations are performed in cases where instability has been predicted for Pb-Sn melt systems; the results obtained show that the model is more stable than the plane front liquid-solid model, and shows similar, finger-like covection when unstable. It is found that convection in the 'mush' zone is confined to the top of that region; convective motion is driven by convection in the bulk fluid, in the absence of natural convection.

Heinrich, Juan C.↗

NASA NTRS1993

Dynamics of aging magnetic clouds

The dynamics of radially expanding magnetic clouds is rigorously analyzed within the framework of ideal MHD. The cloud is modelled as a cylindrically symmetric magnetic flux rope. In the force balance we include the gas pressure gradient and the Lorentz force. Interaction with the ambient solar wind due to expansion of the magnetic cloud is represented by a drag force proportional to the bulk velocity. We consider the self-similar expansion of a polytrope, and reduce the problem to an ordinary nonlinear differential equation for the evolution function. Analyzing the asymptotic behavior of the evolution function, we formulate theoretical expectations for the long-term behavior of cloud parameters. We focus on the temporal evolution of (1) the magnetic field strength; (2) the twist of the field lines; (3) the asymmetry of the total field profile; and (4) the bulk flow speed. We present data from two magnetic clouds observed at 1 AU and 2 AU, respectively, and find good agreement with theoretical expectations. For a peak magnetic field strength at 1 AU of 25 nT and a polytropic index of 0.5, we find that a magnetic cloud can be distinguished from the background interplanetary field up to a distance of about 5 AU. Taking larger magnetic fields and bigger polytropic indices this distance can double.

Osherovich, V. A.↗

NASA NTRS1949

Two-dimensional compressible flow in turbomachines with conic flow surfaces

A general method of analysis is developed for two-dimensional, steady, compressible flow in stators or rotors of radial and mixed flow turbomachines with conic flow surfaces (surfaces of right circular cones generated by center line of flow passage in the axial-radial plane). The variables taken into account are: (1) tip speed of the rotor, (2) flow rate, (3) blade shape, (4) variation in passage height with radius, (5) number of blades, and (6) cone angle of the flow surface. Relaxation methods are used to solve the nonlinear differential equation for the stream function. Two numerical examples are presented; one for compressible and the other for incompressible flow in a centrifugal compressor with thin, straight blades. The results of these examples are given by plots of the streamlines, constant velocity-ratio lines, and constant pressure-ratio line.

Stanitz, John D↗

NASA NTRS1953

On transonic flow past a wave-shaped wall

This report is an extension of a previous investigation (described in NACA rep. 1069) concerned with the solution of the nonlinear differential equation for transonic flow past a wavy wall. In the present work several new notions are introduced which permit the solution of the recursion formulas arising from the method of integration in series. In addition, a novel numerical tests of convergence, applied to the power series (in transonic similarity parameter) representing the local Mach number distribution at the boundary, indicates that smooth symmetrical potential flow past the wavy wall is no longer possible once the critical value of the stream Mach number has been exceeded.

Kaplan, Carl↗

NASA NTRS1993

Marangoni instability in a liquid layer with two free surfaces

The onset of the Marangoni instability in a liquid layer with two free nearly insulating surfaces heated from below is studied. Linear stability analysis yields a condition for the emergence of a longwave or a finite wavelength instability from the quiescent equilibrium state. Using the method of asymptotic expansions, a weakly nonlinear evolution equation describing the spatiotemporal behavior of the velocity and temperature fields at the onset of the longwave instability is derived. The latter is given by delta(M) = 24, delta(M) being the difference between the upper and the lower Marangoni numbers. It is shown that in some parametric range one convective cell forms across the layer, while in other parametric domains two convective cells emerge between the two free surfaces.

Deissler, Robert J.↗

NASA NTRS1991

MHD shocks in coronal mass ejections

The primary objective of this research program is the study of the magnetohydrodynamic (MHD) shocks and nonlinear simple waves produced as a result of the interaction of ejected lower coronal plasma with the ambient corona. The types of shocks and nonlinear simple waves produced for representative coronal conditions and disturbance velocities were determined. The wave system and the interactions between the ejecta and ambient corona were studied using both analytic theory and numerical solutions of the time-dependent, nonlinear MHD equations. Observations from the SMM coronagraph/polarimeter provided both guidance and motivation and are used extensively in evaluating the results. As a natural consequence of the comparisons with the data, the simulations assisted in better understanding the physical interactions in coronal mass ejections (CME's).

Steinolfson, R. S.↗

NASA NTRS1993

Boundary conditions and the simulation of low Mach number flows

The problem of accurately computing low Mach number flows, with the specific intent of studying the interaction of sound waves with incompressible flow structures, such as concentrations of vorticity is considered. This is a multiple time (and/or space) scales problem, leading to various difficulties in the design of numerical methods. Concentration is on one of these difficulties - the development of boundary conditions at artificial boundaries which allow sound waves and vortices to radiate to the far field. Nonlinear model equations are derived based on assumptions about the scaling of the variables. Then these are linearized about a uniform flow and exact boundary conditions are systematically derived using transform methods. Finally, useful approximations to the exact conditions which are valid for small Mach number and small viscosity are computed.

Hagstrom, Thomas↗

NASA NTRS1993

On the estimation of the correlation dimension and its application to radar reflector discrimination

Recently, system theorists have recognized that low order systems of nonlinear differential equations can give rise to solutions which are neither periodic, constant, nor predictable in steady state, but which are nonetheless bounded and deterministic. This behavior, which was first described in the study of weather systems, has been termed 'chaotic.' Much study of chaotic systems has concentrated on analysis of the systems' phase space attractors. It has been recognized that invariant measures of the attractor possess inherent information about the system. One such measure is the dimension of the attractors. The dimension of a chaotic attractor has been shown to be noninteger, leading to the term 'strange attractor;' the attractor is said to have a fractal structure. The correlation dimension has become one of the most popular measures of dimension. However, many problems have been identified in correlation dimension estimation from time sequences. The most common methods for obtaining the correlation dimension have been least squares curves fitting to find the slope of the correlation integral and the Takens Estimator. However, these estimates show unacceptable sensitivity to the upper limit on the distance chosen. Here, a new method is proposed which is shown to be rather insensitive to the upper limit and to perform in a very stable manner, at least in the absence of noise. The correlation dimension is also shown to be an effective discriminant in distinguishing between radar returns resulting from weather and those from the ground. The weather returns are shown to have a correlation dimension generally between 2.0 and 3.0, while ground returns have a correlation dimension exceeding 3.0.

Barnett, Kevin D.↗

NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗

NASA NTRS1994

Microwave heating and joining of ceramic cylinders: A mathematical model

A thin cylindrical ceramic sample is placed in a single mode microwave applicator in such a way that the electric field strength is allowed to vary along its axis. The sample can either be a single rod or two rods butted together. We present a simple mathematical model which describes the microwave heating process. It is built on the assumption that the Biot number of the material is small, and that the electric field is known and uniform throughout the cylinder's cross-section. The model takes the form of a nonlinear parabolic equation of reaction-diffusion type, with a spatially varying reaction term that corresponds to the spatial variation of the electromagnetic field strength in the waveguide. The equation is analyzed and a solution is found which develops a hot spot near the center of the cylindrical sample and which then propagates outwards until it stabilizes. The propagation and stabilization phenomenon concentrates the microwave energy in a localized region about the center where elevated temperatures may be desirable.

Booty, Michael R.↗

NASA NTRS1995

Preliminary control system design and analysis for the Space Station Furnace Facility thermal control system

This report presents the Space Station Furnace Facility (SSFF) thermal control system (TCS) preliminary control system design and analysis. The SSFF provides the necessary core systems to operate various materials processing furnaces. The TCS is defined as one of the core systems, and its function is to collect excess heat from furnaces and to provide precise cold temperature control of components and of certain furnace zones. Physical interconnection of parallel thermal control subsystems through a common pump implies the description of the TCS by coupled nonlinear differential equations in pressure and flow. This report formulates the system equations and develops the controllers that cause the interconnected subsystems to satisfy flow rate tracking requirements. Extensive digital simulation results are presented to show the flow rate tracking performance.

Jackson, M. E.↗

NASA NTRS1995

Sensitivity analysis of frictional contact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the static frictional contact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of nodal displacements, and Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented for the frictional contact of a composite spherical cap pressed against a rigid plate.

Karaoglan, I.↗

NASA NTRS1995

Dynamic sensitivity analysis of frictional contact/impact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Karaoglan, Levent↗

NASA NTRS1993

Kinetic theory model predictions compared with low-thrust axisymmetric nozzle plume data

A system of nonlinear integral equations equivalent to the steady-state Krook kinetic equation was used to model the flow from a low-thrust axisymmetric nozzle. The mathematical model was used to numerically calculate the number density, temperature, and velocity of a simple gas as it expands into a near vacuum. With these quantities the gas pressure and flow directions of the gas near the exit plane were calculated and compared with experimental values for a low-thrust nozzle of the same geometry and mass flow rate.

Riley, B. R.↗

NASA NTRS1994

Cyclogenesis in a saturated environment

The dynamics of baroclinic wave growth in a saturated environment is examined using linear and nonlinear models employing a parameterization of latent heat release that assumes all rising air is saturated, and saturation equivalent potential temperature is conserved on ascent. Piecewise potential vorticity (PV) diagnostics are used to interpret the results. When the stability to vertical displacements in saturated air is allowed to increase with height, as it must in an atmosphere with a constant, positive lapse rate of potential temperature, the growth rates of the most unstable modes of the Eady problem grow only marginally faster than the modes of the dry problem. The vertical variation of moist static stability produces a gradient of moist potential vorticity in the rising air, eliminating the short-wave cutoff present in the dry Eady problem. The destabilization of the short waves is shown to be associated with the interaction between surface potential temperature anomalies and diabatically generated lower-tropospheric potential vorticity anomalies. Nonlinear primitive equation simulations, starting from normal-mode initial conditions, show that while the dry wave grows at nearly the linear growth rate until maximum amplitude is reached, the moist wave grows significantly faster than the linear growth rate at finite amplitude. This enhanced growth is associated with the rapid amplification of a mesoscale PV anomaly generated by latent heat release at the warm front. The rapid amplification of the surface cyclone results from the superposition of the circulation associated with this misoscale PV anomaly upon the circulation associated with the surface and upper boundary potential temperature anomalies. Additional integrations with finite-amplitude initial conditions more typical of atmospheric conditions exhibit similar behavior. It is suggested that many of the rapid cyclogenesis events that occur as upper-troposheric PV anomalies cross the east coasts of continents may arise from the rapid generation of PV anomalies by condensational heating in the moist maritime lower troposphere.

Whitaker, Jeffrey S.↗

NASA NTRS1994

Parasitic instabilities in magnetized, differentially rotating disks

Velikhov, Chandrasekhar, Balbus & Hawley have discovered a MHD instability of differentially rotating fluids that may explain the effective viscosity of accretion disks. If the unperturbed magnetic energy density is much less than the gas pressure, this `magnetorotational' instability (MRI) arises on small scales and is approximately incompressible. If the unperturbed field is suffiently, strong, the MRI is suppressed. We therefore ask whether the MRI mechanism can amplify even very weak initial fields until they reach equipartition with the gas pressure. We show that in the total incompressible limit, the MRI modes are exact solutions of the nonlinear fluid equations, even if the perturbed magnetic field is much larger than the unperturbed field. Also, we present a new exact solution in the opposite limit that the magnetic pressure is much larger than that of the gas. On the other hand, we show that the incompressible MRI modes are themselves subject to parasitic instabilities with instantaneous growth rates proportional to the MRI amplitude. Some of the parasitic instabilities are of the Kelvin-Helmholtz type, while others are less familiar. The eigenfunctions of the latter group suggest that they may promote rapid reconnection of the MRI field. Thus, parasitic modes may stop MRI growth at subequipartition amplitudes if it developes from a sufficiently weak initital field.

Goodman, Jeremy↗

NASA NTRS1995

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗

NASA NTRS1995