Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical differentiation”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Analytical results for postbuckling behavior of plates in compression and shear

The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.

Stein, M.↗
NASA NTRS1984

Singular asymptotic expansions in nonlinear rotordynamics

During hot firing ground testing of the Space Shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.

Day, W. B.↗
NASA NTRS1984

Chemoviscosity modeling for thermosetting resins - I

A new analytical model for chemoviscosity variation during cure of thermosetting resins was developed. This model is derived by modifying the widely used WLF (Williams-Landel-Ferry) Theory in polymer rheology. Major assumptions involved are that the rate of reaction is diffusion controlled and is linearly inversely proportional to the viscosity of the medium over the entire cure cycle. The resultant first order nonlinear differential equation is solved numerically, and the model predictions compare favorably with experimental data of EPON 828/Agent U obtained on a Rheometrics System 4 Rheometer. The model describes chemoviscosity up to a range of six orders of magnitude under isothermal curing conditions. The extremely non-linear chemoviscosity profile for a dynamic heating cure cycle is predicted as well. The model is also shown to predict changes of glass transition temperature for the thermosetting resin during cure. The physical significance of this prediction is unclear at the present time, however, and further research is required. From the chemoviscosity simulation point of view, the technique of establishing an analytical model as described here is easily applied to any thermosetting resin. The model thus obtained is used in real-time process controls for fabricating composite materials.

Hou, T. H.↗
NASA NTRS1985

Singular asymptotic expansions in nonlinear rotordynamics

During hot firing ground testing of the Space shuttle's Main Engine, vibrations of the liquid oxygen pump occur at frequencies which cannot be explained by the linear Jeffcott model of the rotor. The model becomes nonlinear after accounting for deadband, side forces, and rubbing. Two phenomena present in the numerical solutions of the differential equations are unexpected periodic orbits of the rotor and tracking of the nonlinear frequency. A multiple scale asymptotic expansion of the differential equations is used to give an analytic explanation of these characteristics.

Day, W. B.↗
NASA NTRS1985

Analytical results for post-buckling behaviour of plates in compression and in shear

The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.

Stein, M.↗
NASA NTRS1986

Modelling of jet- and swirl-stabilized reacting flows in axisymmetric combustors

Turbulent nonreactive and reactive flows with and without swirl are analyzed, with particular attention given to the flow fields of a gas-fueled nonpremixed swirl-stabilized combustor and a premixed opposed-jet combustor. Local mean flow properties, including velocity, temperature, and major species concentrations, are calculated by solving numerically the governing partial differential equations with associated submodels for turbulence and combustion. The results of the study indicate that the constant-density k-epsilon turbulence model provides a satisfactory representation of the aerodynamics in most practical combustor flows. The exception is the case of jet-stabilized combustor flow, due to the fact that the k-epsilon model cannot replicate the highly dissipative phenomenon found in such flows.

Nikjooy, M.↗
NASA NTRS1986

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis

A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.

Warming, Robert F.↗
NASA NTRS1988

Partitioning of regular computation on multiprocessor systems

Problem partitioning of regular computation over two dimensional meshes on multiprocessor systems is examined. The regular computation model considered involves repetitive evaluation of values at each mesh point with local communication. The computational workload and the communication pattern are the same at each mesh point. The regular computation model arises in numerical solutions of partial differential equations and simulations of cellular automata. Given a communication pattern, a systematic way to generate a family of partitions is presented. The influence of various partitioning schemes on performance is compared on the basis of computation to communication ratio.

Lee, Fung Fung↗
NASA NTRS1990

Two-dimensional convection and radiation with scattering from a Poiseuille flow

Two-dimensional combined convection and radiation heat transfer from a gray scattering fluid in a reflecting channel is considered. The model, represented by a set of simultaneous nonlinear integro-partial differential equations, is solved numerically. The effects of aspect ratio, conduction-radiation parameter, scattering albedo, and wall emissivity, are systematically investigated. It is found that these parameters have a significant influence on the temperature field and alter the radiative and convective fluxes at the hot and cold walls. In particular, when radiation effects are considerable, the heat-transfer characteristics of the fluid at the hot and cold walls are very different.

Kassemi, M.↗
NASA NTRS1990

Partitioning of regular computation on multiprocessor systems

Problem partitioning of regular computation over two dimensional meshes on multiprocessor systems is examined. The regular computation model considered involves repetitive evaluation of values at each mesh point with local communication. The computational workload and the communication pattern are the same at each mesh point. The regular computation model arises in numerical solutions of partial differential equations and simulations of cellular automata. Given a communication pattern, a systematic way to generate a family of partitions is presented. The influence of various partitioning schemes on performance is compared on the basis of computation to communication ratio.

Lee, Fung F.↗
NASA NTRS1991

Experiments and modeling

The usual approach in establishing the correctness and accuracy of turbulence models is to numerically solve the modeled differential equations and then compare the results with the experiment. However, in the case of a discrepancy, this procedure does not pinpoint where in the model the drawback lies. It is also possible that the model overcompensates one physical phenomenon and undercompensates the other so that the net result is a good agreement between the two. Therefore, a more desirable approach is to directly compare the individual terms in the equations with their models. To achieve this objective, primary physical experiments were used to carry out the second moment budgets. These can then be used to analyze and assess various models and closure assumptions and seek improvements/modifications where models prove deficient.

Shabbir, Aamir↗
NASA NTRS1992

On constructing three-dimensional overlapping grids with CMPGRD

Techniques for the construction of three-dimensional composite overlapping grids, using the grid construction program CMPGRD, are described. The overlapping approach can be used to generate grids for regions of complicated geometry. The grids can be constructed to be smooth and free from coordinate singularities. The ability to create smooth grids for complicated regions is an important first step towards the accurate numerical solution of partial differential equations. The creation of grids for surfaces defined by cross-sections such as an airplane wing is described. A method for integrating the patched surfaces generated by a computer aided design (CAD) package with the CMPGRD program is described. Additionally, the creation of grids in regions where surfaces intersect is described.

Henshaw, William D.↗
NASA NTRS1993

An efficient model for coupling structural vibrations with acoustic radiation

The scattering of an incident wave by a flexible panel is studied. The panel vibration is governed by the nonlinear plate equations while the loading on the panel, which is the pressure difference across the panel, depends on the reflected and transmitted waves. Two models are used to calculate this structural-acoustic interaction problem. One solves the three dimensional nonlinear Euler equations for the flow-field coupled with the plate equations (the fully coupled model). The second uses the linear wave equation for the acoustic field and expresses the load as a double integral involving the panel oscillation (the decoupled model). The panel oscillation governed by a system of integro-differential equations is solved numerically and the acoustic field is then defined by an explicit formula. Numerical results are obtained using the two models for linear and nonlinear panel vibrations. The predictions given by these two models are in good agreement but the computational time needed for the 'fully coupled model' is 60 times longer than that for 'the decoupled model'.

Frendi, Abdelkader↗
NASA NTRS1993

Nonlinear evolution of magnetic flux ropes. I - Low-beta limit

We study the nonlinear self-similar evolution of a cylindrical magnetic flux tube with two components of the magnetic field, axial and azimuthal. We restrict ourselves to the case of a plasma of low beta. Introducing a special class of configurations we call 'separable fields', we reduce the problem to an ordinary differential equation. Two cases are to be distinguished: (1) when the total field minimizes on the symmetry axis, the magnetic configuration inexorably collapses, and (2) when, on the other hand, the total field maximizes on the symmetry axis, the magnetic configuration behaves analogously to a nonlinear oscillator. Here we focus on the latter case. The effective potential of the motion contains two terms: a strong repulsive term and a weak restoring term associated with the pinch. We solve the nonlinear differential equation of motion numerically and find that the period of oscillations grows exponentially with the energy of the oscillator. Our treatment emphasizes the role of the force-free configuration as the lowest potential energy state about which the system oscillates.

Osherovich, V. A.↗
NASA NTRS1950

Approximate Turbulent Boundary-Layer Development in Plane Compressible Flow Along Thermally Insulated Surfaces with Application to Supersonic-Tunnel Contour Correction

Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

BOUNDARY LAYER - WIND TUNNEL WALLS↗
NASA NTRS1993

The Chebyshev-Legendre method: Implementing Legendre methods on Chebyshev points

We present a new collocation method for the numerical solution of partial differential equations. This method uses the Chebyshev collocation points, but because of the way the boundary conditions are implemented, it has all the advantages of the Legendre methods. In particular, L2 estimates can be obtained easily for hyperbolic and parabolic problems.

Don, Wai Sun↗
NASA NTRS1994

Ionospheric corrections to precise time transfer using GPS

The free electrons in the earth's ionosphere can retard the time of reception of GPS signals received at a ground station, compared to their time in free space, by many tens of nanoseconds, thus limiting the accuracy of time transfer by GPS. The amount of the ionospheric time delay is proportional to the total number of electrons encountered by the wave on its path from each GPS satellite to a receiver. This integrated number of electrons is called Total Electron Content, or TEC. Dual frequency GPS receivers designed by Allen Osborne Associates, Inc. (AOA) directly measure both the ionospheric differential group delay and the differential carrier phase advance for the two GPS frequencies and derive from this the TEC between the receiver and each GPS satellite in track. The group delay information is mainly used to provide an absolute calibration to the relative differential carrier phase, which is an extremely precise measure of relative TEC. The AOA Mini-Rogue ICS-4Z and the AOA TurboRogue ICS-4000Z receivers normally operate using the GPS P code, when available, and switch to cross-correlation signal processing when the GPS satellites are in the Anti-Spoofing (A-S) mode and the P code is encrypted. An AOA ICS-Z receiver has been operated continuously for over a year at Hanscom AFB, MA to determine the statistics of the variability of the TEC parameter using signals from up to four different directions simultaneously. The 4-channel ICS-4Z and the 8-channel ICS-4000Z, have proven capabilities to make precise, well calibrated, measurements of the ionosphere in several directions simultaneously. In addition to providing ionospheric corrections for precise time transfer via satellite, this dual frequency design allows full code and automatic codeless operation of both the differential group delay and differential carrier phase for numerous ionospheric experiments being conducted. Statistical results of the data collected from the ICS-4Z during the initial year of ionospheric time delay in the northeastern U.S., and initial results with the ICS-4000Z, will be presented.

Snow, Robert W.↗
NASA NTRS1994

Analysis of wavelet technology for NASA applications

The purpose of this grant was to introduce a broad group of NASA researchers and administrators to wavelet technology and to determine its future role in research and development at NASA JSC. The activities of several briefings held between NASA JSC scientists and Rice University researchers are discussed. An attached paper, 'Recent Advances in Wavelet Technology', summarizes some aspects of these briefings. Two proposals submitted to NASA reflect the primary areas of common interest. They are image analysis and numerical solutions of partial differential equations arising in computational fluid dynamics and structural mechanics.

Wells, R. O., Jr.↗