SEARCH · Engineering Papers
Results for “partial differential equation”
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.
The partial derivatives of the fundamental thermodynamic functions in terms of the variables T and P
Partial derivatives of fundamental thermodynamic functions in terms of temperature and pressure variables
Application of Krylov-Bogoliubov method to the solution of the stellar three body problem.
Higher order perturbations of elements in three- body problem determined using methods of Krylov- Bogoliubov and Poincare
Longitudinal shear loading of a unidirectional composite.
Fiber reinforced composite material behavior examined from shear loading tests, using finite difference analysis for representative PDEs
A study analysis of cable-body systems totally immersed in a fluid stream
A general stability analysis of a cable-body system immersed in a fluid stream is presented. The analytical portion of this analysis treats the system as being essentially a cable problem, with the body dynamics giving the end conditions. The mathematical form of the analysis consists of partial differential wave equations, with the end and auxiliary conditions being determined from the body equations of motion. The equations uncouple to give a lateral problem and a longitudinal problem as in first order airplane dynamics. A series of tests on a tethered wind tunnel model provide a comparison of the theory with experiment.
Calculation of sheath and wake structure about a pillbox-shaped spacecraft in a flowing plasma
A computer program was used for studies of the disturbed zones around bodies in flowing plasmas, particularly spacecraft and their associated sheaths and wakes. The program solved a coupled Poisson-Vlasov system of nonlinear partial differential integral equations to obtain distributions of electric potential and ion and electron density about a finite length cylinder in a plasma flow at arbitrary ion Mach numbers. The approach was applicable to a larger range of parameters than other available approaches. In sample calculations, bodies up to 100 Debye lengths in radius were treated, that is, larger than any previously treated realistically. Applications were made to in-situ satellite experiments.
Rotor blade aeroelastic stability and response in forward flight
The aeroelastic stability and response problem of the coupled flap-lag-torsional dynamics of a hingeless rotor blade in forward flight is treated in a comprehensive manner. The spatial dependence of the partial differential, nonlinear, equations of motion is discretized using a multimodal Galerkin method. The aeroelastic problem is coupled with the trim state of the helicopter obtained from improved, representative, trim procedures. The nonlinear time dependent equilibrium position, or response, about which the equations are linearized is obtained by solving a sequence of linear periodic response problems, using quasi-linearization. Numerous results illustrating blade behavior in forward flight are presented.
Coupled flap-lag-torsional dynamics of hingeless rotor blades in forward flight
The aeroelastic stability and response problem of the coupled flap-lag-torsional dynamics of an isolated hingeless rotor blade in forward flight is considered. Linear, quasi-steady aerodynamics below stall is included. The spatial dependence of the partial differential nonliner equations of motion is discretized using a multimodal Galerkin method. The nonlinear time dependent equilibrium position (steady state response) about which the equations are linearized is obtained by solving a sequence of linear periodic response problems (quasilinearization). Results illustrating blade behavior in forward flight for both soft-in-plane and stiff-in-plane designs are presented. Quasilinearization provides a clear indication of the cases when nonlinear terms due to moderate deflections are important. The results indicate that nonlinearities affect system stability much more than system response. This implies that for obtaining blade vibrations and loads, approximate analytical models based on linearized formulations could be employed, whereas for determining blade stability a more accurate nonlinear analysis would be required.
The fully nonlinear development of Goertler vortices in growing boundary layers
The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.
A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model
The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment.
Analyses of composite fuselage structure under various loading conditions
With the characteristics of high strength/weight and stiffness/weight ratios in advanced composite materials, composite plate structures have been successfully applied to secondary load-carrying structural components in the aerospace industry for the past two decades. Recently, filament wound composite shells are being considered for design of primary fuselage structures. To implement the design the structural response under various loading conditions needs to be predicted accurately. Therefore, it is essential to establish a rigorous analytical solution in the area of composite laminated shells. A closed form solution is presented that predicts the response of a composite shell subjected to internal pressure, axial tension, bending, and torsion. The material of the shell is assumed to be general cylindrically anisotropic. Based on the theory of cylindrical anisotropic elasticity coupled partial differential governing equations are developed using Lehknitskii's stress function approach. The general expressions for the stresses and displacements in the composite cylinders under these loading conditions is discussed. Three examples (45) off-axis unidirectional, (45/45) unsymmetric, and (45/45)(sub s) symmetric angle-ply fiber-reinforced laminated shells are shown to illustrate the effect of radius-to-thickness ratios, coupling, and stacking sequence.
Time as an Observable in Nonrelativistic Quantum Mechanics
The argument follows from the viewpoint that quantum mechanics is taken not in the usual form involving vectors and linear operators in Hilbert spaces, but as a boundary value problem for a special class of partial differential equations-in the present work, the nonrelativistic Schrodinger equation for motion of a structureless particle in four- dimensional space-time in the presence of a potential energy distribution that can be time-as well as space-dependent. The domain of interest is taken to be one of two semi-infinite boxes, one bounded by two t=constant planes and the other by two t=constant planes. Each gives rise to a characteristic boundary value problem: one in which the initial, input values on one t=constant wall are given, with zero asymptotic wavefunction values in all spatial directions, the output being the values on the second t=constant wall; the second with certain input values given on both z=constant walls, with zero asymptotic values in all directions involving time and the other spatial coordinates, the output being the complementary values on the z=constant walls. The first problem corresponds to ordinary quantum mechanics; the second, to a fully time-dependent version of a problem normally considered only for the steady state (time-independent Schrodinger equation). The second problem is formulated in detail. A conserved indefinite metric is associated with space-like propagation, where the sign of the norm of a unidirectional state corresponds to its spatial direction of travel.
Consistent Boundary Conditions for Multicomponent Real Gas Mixtures Based on Characteristic Waves
Previously developed characteristic-wave-based boundary conditions for multicomponent perfect gas mixtures are here extended to account for real gases.
Oxidation Behavior of Carbon Fiber-Reinforced Composites
OXIMAP is a numerical (FEA-based) solution tool capable of calculating the carbon fiber and fiber coating oxidation patterns within any arbitrarily shaped carbon silicon carbide composite structure as a function of time, temperature, and the environmental oxygen partial pressure. The mathematical formulation is derived from the mechanics of the flow of ideal gases through a chemically reacting, porous solid. The result of the formulation is a set of two coupled, non-linear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined at each time step using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The non-linear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual finite element method, allowing for the solution of the differential equations numerically.
Numerical solution of flow problems using body-fitted coordinate systems
The paper deals with numerically generated boundary-fitted coordinate systems. This procedure eliminates the shape of the boundaries as a complicating factor and allows the flow about arbitrary boundaries to be treated essentially as easily as that about simple boundaries. The technique of boundary-fitted coordinate systems is based on a method of automatic numerical generation of a general curvilinear coordinate system having a coordinate line coincident with each boundary of a general multiconnected region involving any number of arbitrarily shaped boundaries. Once the curvilinear coordinate system is generated, any partial differential system of interest may be solved on the coordinate system by transforming the equations and solving the resulting system in finite-difference approximation on the rectangular transformed plane. Attention is given to the types of boundary-fitted coordinate systems, coordinate system control, operation of the coordinate codes, solution of partial differential equations, application to free-surface flow, and other applications of interest.
Modeling and control of flexible space stations (slew maneuvers)
Large orbiting space structures are expected to experience mechanical vibrations arising from several disturbing forces such as those induced by shuttle takeoff or docking and crew movements. The problem is considered of modeling and control of large space structures subject to these and other disturbing forces. The system consists of a (rigid) massive body, which may play the role of experimental modules located at the center of the space station and flexible configurations, consisting of several beams, forming the space structure. A complete dynamic model of the system was developed using Hamilton's principle. This model consists of radial equations describing the translational motion of the central body, rotational equations describing the attitude motions of the body and several beam equations governing the vibration of the flexible members (platform) including appropriate boundary conditions. In summary, the dynamics of the space structure is governed by a complex system of interconnected partial and ordinary differential equations. Using Lyapunov's approach the asymptotic stability of the space structure is investigated. For asymptotic stability of the rest state (nominal trajectory), feedback controls are suggested. In the investigation, stability of the slewing maneuvers is also considered. Several numerical results are presented for illustration of the impact of coupling and the effectiveness of the stabilizing controls. Some insight is provided into the complexity of modeling, analysis and stabilization of actual space structures.
Non-linear heat transfer of solids in orthogonal coordinate systems.
Nonlinear partial differential heat conduction equation /Poisson type/ obtained by using Kirchhoff transformation
Galactic civilizations - Population dynamics and interstellar diffusion
A model is developed of the interstellar diffusion of galactic civilizations which takes into account the population dynamics of such civilizations. The problem is formulated in terms of potential theory, with a family of nonlinear partial differential and difference equations specifying population growth and diffusion for an organism with advantageous genes that undergoes random dispersal while increasing in population locally, and a population at zero population growth. In the case of nonlinear diffusion with growth and saturation, it is found that the colonization wavefront from the nearest independently arisen galactic civilization can have reached the earth only if its lifetime exceeds 2.6 million years, or 20 million years if discretization can be neglected. For zero population growth, the corresponding lifetime is 13 billion years. It is concluded that the earth is uncolonized not because interstellar spacefaring civilizations are rare, but because there are too many worlds to be colonized in the plausible colonization lifetime of nearby civilizations, and that there exist no very old galactic civilizations with a consistent policy of the conquest of inhabited worlds.