Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential algebra”

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

Partial spectral flow in the D1D5 CFT

The two-dimensional 𝒩 = 4 superconformal algebra has a free field realization with four bosons and four fermions. There is an automorphism of the algebra called spectral flow. Under spectral flow, the four fermions are transformed together. In this paper, we study partial spectral flow where only two of the four fermions are transformed. Partial spectral flow is applied to the D1D5 CFT where a marginal deformation moves the CFT away from the free point. The partial spectral flow is broken by the deformation. We show that this effect can be studied due to a transformation of the deformation which is well-defined under partial spectral flow. As a result in the spectrum, we demonstrate how to compute the second-order energy lift of a D1D5P state through its partial spectral flowed state. We find that D1D5P states related by partial spectral flow do not have the same lift in general.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Finite-analytic numerical solution of heat transfer in two-dimensional cavity flow

Heat transfer in cavity flow is numerically analyzed by a new numerical method called the finite-analytic method. The basic idea of the finite-analytic method is the incorporation of local analytic solutions in the numerical solutions of linear or nonlinear partial differential equations. In the present investigation, the local analytic solutions for temperature, stream function, and vorticity distributions are derived. When the local analytic solution is evaluated at a given nodal point, it gives an algebraic relationship between a nodal value in a subregion and its neighboring nodal points. A system of algebraic equations is solved to provide the numerical solution of the problem. The finite-analytic method is used to solve heat transfer in the cavity flow at high Reynolds number (1000) for Prandtl numbers of 0.1, 1, and 10.

Chen, C.-J.↗

Bifurcations in unsteady aerodynamics

Nonlinear algebraic functional expansions are used to create a form for the unsteady aerodynamic response that is consistent with solutions of the time dependent Navier-Stokes equations. An enumeration of means of invalidating Frechet differentiability of the aerodynamic response, one of which is aerodynamic bifurcation, is proposed as a way of classifying steady and unsteady aerodynamic phenomena that are important in flight dynamics applications. Accomodating bifurcation phenomena involving time dependent equilibrium states within a mathematical model of the aerodynamic response raises an issue of memory effects that becomes more important with each successive bifurcation.

Tobak, M.↗

Bifurcations in unsteady aerodynamics

Nonlinear algebraic functional expansions are used to create a form for the unsteady aerodynamic response that is consistent with solutions of the time dependent Navier-Stokes equations. An enumeration of means of invalidating Frechet differentiability of the aerodynamic response, one of which is aerodynamic bifurcation, is proposed as a way of classifying steady and unsteady aerodynamic phenomena that are important in flight dynamics applications. Accommodating bifurcation phenomena involving time dependent equilibrium states within a mathematical model of the aerodynamic response raises an issue of memory effects that becomes more important with each successive bifurcation.

Tobak, M.↗

Algebraic Generation Of Two-Dimensional Grids

Two-Boundary Grid Generation (TBGG) computer program applies interactive algebraic grid-generation technique in two dimensions. Incorporates algebraic equations that relate uniform rectangular computational domain to arbitrary physical domain. Hermite cubic interpolation performed between fixed, nonintersecting boundaries. Applied to variety of problems in which discrete techniques used to solve such partial differential equations as governing equations of fluid flow. Written in FORTRAN 77.

Smith, Robert E.↗

Validation of Mars-GRAM and Planned New Features

For altitudes below 80 km, Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is based on output climatology from NASA Ames Mars General Circulation Model (MGCM). At COSPAR 2002, results were presented of validation tests of Mars-GRAM versus data from Mars Global Surveyor Thermal Emission Spectrometer (TES) and Radio Science (RS) experiment. Further validation tests are presented comparing Mars- GRAM densities with those from the European Mars Climate Database (MCD), and comparing densities from both Mars-GRAM and MCD against TES observations. Throughout most of the height and latitude range of TES data (040 km and 70s to 70N), good agreement is found between atmospheric densities from Mars-GRAM and MCD. However, at the season and latitude zone for Mars Phoenix arrival and landing (Ls = 65 to 80 degrees and latitude 65 to 75N), Mars-GRAM densities are about 30 to 45 percent higher than MCD densities near 40 km altitude. Further evaluation is warranted concerning potential impact of these model differences on planning for Phoenix entry and descent. Three planned features for Mars-GRAM update are also discussed: (1) new MGCM and Thermospheric General Circulation Model data sets to be used as a revised basis for Mars-GRAM mean atmosphere, (2) a new feature to represent planetary-scale traveling waves for upper altitude density variations (such as found during Mars Odyssey aerobraking), and (3) a new model for effects of high resolution topographic slope on winds near the surface (0 to 4.5 km above MOLA topography level). Mars-GRAM slope winds will be computed from a diagnostic (algebraic) relationship based on Ye, Segal, and Pielke (1990). This approach differs from mesoscale models (such as MRAMS and Mars MM5), which use prognostic, full-physics solutions of the time- and space-dependent differential equations of motion. As such, slope winds in Mars-GRAM will be consistent with its "engineering-level" approach, and will be extremely fast and easy to evaluate, compared with mesoscale model solutions. Mars-GRAM slope winds are not being suggested as a replacement for sophisticated, full-physics Mars mesoscale models, but may have value, particularly for preliminary screening of large numbers of candidate landing sites for future Mars missions, such as Phoenix and Mars Science Laboratory. Test output is presented from Mars-GRAM slope winds in the area of Gusev Crater and Valles Marineris.

Justus, C. G.↗

Sensitivity Analysis and Optimization of Aerodynamic Configurations with Blend Surfaces

A novel (geometrical) parametrization procedure using solutions to a suitably chosen fourth order partial differential equation is used to define a class of airplane configurations. Inclusive in this definition are surface grids, volume grids, and grid sensitivity. The general airplane configuration has wing, fuselage, vertical tail and horizontal tail. The design variables are incorporated into the boundary conditions, and the solution is expressed as a Fourier series. The fuselage has circular cross section, and the radius is an algebraic function of four design parameters and an independent computational variable. Volume grids are obtained through an application of the Control Point Form method. A graphic interface software is developed which dynamically changes the surface of the airplane configuration with the change in input design variable. The software is made user friendly and is targeted towards the initial conceptual development of any aerodynamic configurations. Grid sensitivity with respect to surface design parameters and aerodynamic sensitivity coefficients based on potential flow is obtained using an Automatic Differentiation precompiler software tool ADIFOR. Aerodynamic shape optimization of the complete aircraft with twenty four design variables is performed. Unstructured and structured volume grids and Euler solutions are obtained with standard software to demonstrate the feasibility of the new surface definition.

Thomas, A. M.↗

Shear free, twisting Einstein-Maxwell metrics in the Newman-Penrose formalism

The problem of finding algebraically special solutions to the vacuum Einstein-Maxwell equations was investigated using a spin coefficient formalism. The general case in which the degenerate null vectors are not hypersurface orthogonal is reduced to a problem of solving five coupled differential equations that are no longer dependent on the affine parameter along the degenerate null directions. It is shown that the most general regular, shear-free, nonradiating solution to these equations is the Kerr-Newman metric.

Lind, R. W.↗

Grid and design variables sensitivity analyses for NACA four-digit wing-sections

Two distinct parameterization procedures are developed for investigating the grid sensitivity with respect to design parameters of a wing-section example. The first procedure is based on traditional (physical) relations defining NACA four-digit wing-sections. The second is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B-Splines) for defining the wing-section geometry. An interactive algebraic grid generation technique, known as Hermite Cubic Interpolation, is employed to generate C-type grids around wing-sections. The grid sensitivity of the domain with respect to design and grid parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations. A comparison of the sensitivity coefficients with those obtained using a finite-difference approach has been made to verify the feasibility of the approach. The aerodynamic sensitivity coefficients are obtained using the compressible two-dimensional thin-layer Navier-Stokes equations.

Sadrehaghighi, Ideen↗

Periodic solutions to systems of reaction-diffusion equations

Necessary and sufficient conditions are derived for the existence of temporally periodic 'dissipative structure' solutions in weak diffusion with the reaction rate terms dominant in a generic system of reaction-diffusion differential equations. The enumerator index i of the equations denotes the density or concentration of the ith participating molecular or biological species, and D sub i is the diffusivity constant for the ith species while Q sub i (c), an algebraic function of the n-tuple c, expresses the local rate of production of the ith species due to chemical reactions or biological interactions.

Rosen, G.↗

Convergence of infinite dimensional sampled LQR problems - Theory and numerical results

A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.

Rosen, I. G.↗

Analysis and experimental evaluation of a Stewart platform-based force/torque sensor

The kinematic analysis and experimentation of a force/torque sensor whose design is based on the mechanism of the Stewart Platform are discussed. Besides being used for measurement of forces/torques, the sensor also serves as a compliant platform which provides passive compliance during a robotic assembly task. It consists of two platforms, the upper compliant platform (UCP) and the lower compliant platform (LCP), coupled together through six spring-loaded pistons whose length variations are measured by six linear voltage differential transformers (LVDT) mounted along the pistons. Solutions to the forward and inverse kinematics of the force sensor are derived. Based on the known spring constant and the piston length changes, forces/torques applied to the LCP gripper are computed using vector algebra. Results of experiments conducted to evaluate the sensing capability of the force sensor are reported and discussed.

Nguyen, Charles C.↗

Bleed-boundary conditions for numerically simulated mixed-compression supersonic inlet flow

A numerical study was performed to investigate boundary conditions for bleed openings that are needed in the numerical simulation of critical flow through turbojet inlets with bleed. Nine different boundary conditions (BCs) were devised and examined. Of the BCs examined, three were able to stabilize the terminal shock near the inlet's throat. This numerical study was based on the ensemble-averaged, unsteady, three-dimensional, compressible Navier-Stokes equations closed by the Baldwin-Lomax algebraic turbulence model. Solutions to the Navier-Stokes equations were obtained by using a hybrid implicit-explicit method (Beam and Warming/MacCormack). The multiblock grid system used in the solution procedure was generated by solving a set of elliptic partial differential equations. Solutions obtained were compared with available experimental data.

Chyu, W. J.↗

Grid and aerodynamic sensitivity analyses of airplane components

An algorithm is developed to obtain the grid sensitivity with respect to design parameters for aerodynamic optimization. The procedure is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B-Splines) for defining the wing-section geometry. An interactive algebraic grid generation technique, known as Two-Boundary Grid Generation (TBGG) is employed to generate C-type grids around wing-sections. The grid sensitivity of the domain with respect to geometric design parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations such as a wing or fuselage. The aerodynamic sensitivity coefficients are obtained by direct differentiation of the compressible two-dimensional thin-layer Navier-Stokes equations. An optimization package has been introduced into the algorithm in order to optimize the wing-section surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the wing-section.

Sadrehaghighi, Ideen↗

High accuracy solutions of incompressible Navier-Stokes equations

In recent years, high accuracy finite difference approximations were developed for partial differential equations of elliptic type, with particular emphasis on the convection-diffusion equation. These approximations are of compact type, have a local truncation error of fourth order, and allow the use of standard iterative schemes to solve the resulting systems of algebraic equations. These high accuracy approximations are extended to the solution of Navier-Stokes equations. Solutions are obtained for the model problem of driven cavity and are compared with solutions obtained using other approximations and those obtained by other authors. It is discovered that the high order approximations do indeed produce high accuracy solutions and have a potential for use in solving important problems of viscous fluid flows.

Gupta, Murli M.↗

Influence of cross-flow on nonlinear Tollmien-Schlichting/vortex interaction

The transition of an incompressible three-dimensional boundary layer with strong cross-flow is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two lower-branch Tollmien-Schlichting waves which mutually interact nonlinearly to induce a longitudinal vortex flow. The vortex motion in turn gives rise to significant wave modulation via wall-shear forcing. The characteristic Reynolds number is large and, as a consequence, the waves' and the vortex motion are governed primarily by triple deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Following analysis and computation over a wide range of parameters, three distinct responses are found to emerge in the nonlinear behavior of the flow solution downstream: an algebraic finite-distance singularity, far-downstream saturation or far-downstream wave decay leaving pure vortex flow. These depend on the input conditions, the wave angles and the size of the cross flow.

Davis, D. A. R.↗

The modified equation approach to the stability and accuracy analysis of finite-difference methods

The stability and accuracy of finite-difference approximations to simple linear partial differential equations are analyzed by studying the modified partial differential equation. Aside from round-off error, the modified equation represents the actual partial differential equation solved when a numerical solution is computed using a finite-difference equation. The modified equation is derived by first expanding each term of a difference scheme in a Taylor series and then eliminating time derivatives higher than first order by certain algebraic manipulations. The connection between 'heuristic' stability theory based on the modified equation approach and the von Neumann (Fourier) method is established. In addition to the determination of necessary and sufficient conditions for computational stability, a truncated version of the modified equation can be used to gain insight into the nature of both dissipative and dispersive errors.

Warming, R. F.↗

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity

The determination of the elastodynamic fields of an ellipsoidal inhomogeneity is studied in detail via the eigenstrain approach. A complete formulation and a treatment of both types of eigenstrains for equivalence between the inhomogeneity problem and the inclusion problem are given. This approach is shown to be mathematically identical to other approaches such as the direct volume integral formulation. Expanding the eigenstrains and applied strains in the polynomial form in the position vector and satisfying the equivalence conditions at every point, the governing simultaneous algebraic equations for the unknown coefficients in the eigenstrain expansion are derived. The elastodynamic field outside an ellipsoidal inhomogeneity in a linear elastic isotropic medium is given as an example. The angular and frequency dependence of the induced displacement field, as well as the differential and total cross sections are formally given in series expansion form for the case of uniformly distributed eigenstrains.

Fu, L. S.↗