Engineering Papers⌕ Search

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.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

A theory for predicting composite laminate warpage resulting from fabrication

Linear laminate theory is used in conjunction with the moment-curvature relationship to derive equations for predicting end deflections due to warpage without solving the coupled fourth-order partial differential equations of the plate. Using these equations, it is found that a 1 deg error in the orientation angle of one ply is sufficient to produce warpage end deflection equal to two laminate thicknesses in a 10 inch by 10 inch laminate made from 8-ply Mod-I/epoxy. From a sensitivity analysis on the governing parameters, it is found that a 3 deg fiber migration or a void volume ratio of three percent in some plies is sufficient to produce laminate warpage corner deflection equal to several laminate thicknesses. Tabular and graphical data are presented which can be used to identify possible errors contributing to laminate warpage and/or to obtain an a priori assessment when unavoidable errors during fabrication are anticipated.

Chamis, C. C.↗
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given; the equations are used to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the intial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given, in order to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1976

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitational and rotational terms in the equations are of first order in the space variables, the pressure-gradient terms are of second order, and the turbulent-viscosity term is of third order. The presence of turbulent viscosity ensures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial flow is always inward and allows collapse to occur (axially) even when the rotation is large. An approximate solution of the governing partial differential equations is also given in order to study the spatial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1986

Optimal control of aeroassisted plane change maneuver using feedback expansions

A guidance law for an aeroassisted plane change maneuver is developed by an asymptotic expansion technique using a small parameter which essentially represents the ratio of the inertial forces to the atmospheric forces. This guidance law minimizes the energy loss while meeting terminal constraints on the altitude, flight path angle, and heading angle. By neglecting the inertial forces, the resulting optimization problem is integrable and can be determined in closed form. This zeroth-order solution is the first term in an asymptotic series solution of the Hamilton-Jacobi-Bellman equation. The remaining terms are determined from the solution of a first-order, linear partial differential equation whose solution requires only quadrature integration. Our initial results in using this guidance scheme are encouraging.

Mishne, D.↗
NASA NTRS1989

Numerical computation of shock wave-turbulent boundary layer interaction in transonic flow over an axisymmetric curved hill

A control-volume based finite difference computation of a turbulent transonic flow over an axisymmetric curved hill is presented. The numerical method is based on the SIMPLE algorithm, and hence the conservation of mass equation is replaced by a pressure correction equation for compressible flows. The turbulence is described by a k-epsilon turbulence model supplemented by a near-wall turbulence model. In the method, the dissipation rate in the region very close to the wall is obtained from an algebraic equation and that for the rest of the flow domain is obtained by solving a partial differential equation for the dissipation rate. The other flow equations are integrated up to the wall. It is shown that the present turbulence model yields the correct location of the compression shock. The other computational results are also in good agreement with experimental data.

Kim, S.-W.↗
NASA NTRS1990

Quadrature imposition of compatibility conditions in Chebyshev methods

This paper handles a problem encountered in the solution by way of Chebyshev spectral collocation discretization of the Poisson equation with homogeneous Neumann boundary conditions. When pseudospectral Chebyshev methodology is utilized to discretize the partial differential equation these integrals have to be approximated by an appropriate quadrature formula. The correct quadrature formula is developed for these problems, taking into account the degree of the polynomials involved. It is shown that this formula leads to a well-conditioned Chebyshev approximation to the differential equations and that the compatibility condition is automatically satisfied.

Gottlieb, D.↗
NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗
NASA NTRS2000

Viscoelastic flow modeling in the extrusion of a dough-like fluid

This work attempts to investigate the effect of viscoelasticity and three-dimensional geometry in screw channels. The Phan-Thien Tanner (PTT) constitutive equation with simplified model parameters was solved in conjunction with the flow equations. Polyflow, a commercially available finite element code was used to solve the resulting nonlinear partial differential equations. The PTT model predicted one log scale lower pressure buildup compared to the equivalent Newtonian results. However, the velocity profile did not show significant changes for the chosen PTT model parameters. Past Researchers neglected viscoelastic effects and also the three dimensional nature of the flow in extruder channels. The results of this paper provide a starting point for further simulations using more realistic model parameters, which may enable the food engineer to more accurately scale-up and design extrusion processes.

NASA Discipline Life Support Systems↗
NASA NTRS1973

Nonlinear Equations for Bending of Rotating Beams with Application to Linear Flap-Lag Stability of Hingeless Rotors

The nonlinear partial differential equations for the flapping and lead-lag degrees of freedom of a torisonally rigid, rotating cantilevered beam are derived. These equations are linearized about an equilibrium condition to study the flap-lag stability characteristics of hingeless helicopter rotor blades with zero twist and uniform mass and stiffness in the hovering flight condition. The results indicate that these configurations are stable because the effect of elastic coupling more than compensates for the destabilizing flap-lag Coriolis and aerodynamic coupling. The effect of higher bending modes on the lead-lag damping was found to be small and the common, centrally hinged, spring restrained, rigid blade approximation for elastic rotor blades was shown to be resonably satisfactory for determining flap-lag stability. The effect of pre-cone was generally stabilizing and the effects of rotary inertia were negligible.

Hodges, D. H.↗
NASA NTRS1973

Unsteady turbulent boundary layer analysis

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of turbulent energy equations which include upstream history efforts, are solved numerically. An explicit finite difference analog to the partial differential equations is formulated and developed into a computer program. Calculations were made for a variety of unsteady flows in both two and three dimensions but primarily for two dimensional flow fields in order to first understand some of the fundamental physical aspects of unsteady turbulent boundary layers. Oscillating free stream flows without pressure gradient, oscillating retarded free stream flows and monotonically time-varying flows have all been studied for a wide frequency range. It was found that to the lowest frequency considered, the lower frequency bound being determined by economic considerations (machine time), there were significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗
NASA NTRS1981

A Quasi-Parabolic technique for computation of three-dimensional viscous flows

A computational technique is presented for obtaining flowfield solutions to a parabolic form of the Navier-Stokes equations. The point of departure is the General Interpolant Method (GIM) which provides a discretization for partial differential equations on arbitrary three-dimensional geometries. The new scheme, termed Quasi-Parabolic, treats the parabolized equations but with 'time-like' terms appended. Addition of these extra terms, which are relaxed by iteration, avoid many of the singularities inherent in classical parabolic Navier-Stokes methods. Streamwise derivatives are approximated by three-point backward differences and the cross plane operators use an alternating forward-backward sweep. A two-step sequence is used to implement the difference scheme in the spatial dimensions and a time-like relaxation converges the quasi-parabolic procedure at each plane. Solutions are presented for flows in two and three dimensions. Inviscid flows are solved for internal and external applications and viscous flows in boundary layers and free shear layers are also computed with the GIM/Quasi-Parabolic scheme.

Spradley, L. W.↗
NASA NTRS1983

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem.

Phillips, T. N.↗
NASA NTRS1983

Model error estimation for distributed systems described by elliptic equations

A function space approach is used to develop a theory for estimation of the errors inherent in an elliptic partial differential equation model for a distributed parameter system. By establishing knowledge of the inevitable deficiencies in the model, the error estimates provide a foundation for updating the model. The function space solution leads to a specification of a method for computation of the model error estimates and development of model error analysis techniques for comparison between actual and estimated errors. The paper summarizes the model error estimation approach as well as an application arising in the area of modeling for static shape determination of large flexible systems.

Rodriguez, G.↗
NASA NTRS1984

Numerical solution of a coupled pair of elliptic equations from solid state electronics

Iterative methods are considered for the solution of a coupled pair of second order elliptic partial differential equations which arise in the field of solid state electronics. A finite difference scheme is used which retains the conservative form of the differential equations. Numerical solutions are obtained in two ways, by multigrid and dynamic alternating direction implicit methods. Numerical results are presented which show the multigrid method to be an efficient way of solving this problem. Previously announced in STAR as N83-30109

Phillips, T. N.↗
NASA NTRS1986

Reduction of the effects of the communication delays in scientific algorithms on message passing MIMD architectures

The efficient implementation of algorithms on multiprocessor machines requires that the effects of communication delays be minimized. The effects of these delays on the performance of a model problem on a hypercube multiprocessor architecture is investigated and methods are developed for increasing algorithm efficiency. The model problem under investigation is the solution by red-black Successive Over Relaxation YOUN71 of the heat equation; most of the techniques described here also apply equally well to the solution of elliptic partial differential equations by red-black or multicolor SOR methods. Methods for reducing communication traffic and overhead on a multiprocessor are identified and results of testing these methods on the Intel iPSC Hypercube reported. Methods for partitioning a problem's domain across processors, for reducing communication traffic during a global convergence check, for reducing the number of global convergence checks employed during an iteration, and for concurrently iterating on multiple time-steps in a time-dependent problem. Empirical results show that use of these models can markedly reduce a numewrical problem's execution time.

Saltz, J. H.↗