Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Continuous derivative”

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

Implementation of structural response sensitivity calculations in a large-scale finite-element analysis system

The methodology used to implement structural sensitivity calculations into a major, general-purpose finite-element analysis system (SPAR) is described. This implementation includes a generalized method for specifying element cross-sectional dimensions as design variables that can be used in analytically calculating derivatives of output quantities from static stress, vibration, and buckling analyses for both membrane and bending elements. Limited sample results for static displacements and stresses are presented to indicate the advantages of analytically calculating response derivatives compared to finite difference methods. Continuing developments to implement these procedures into an enhanced version of SPAR are also discussed.

Giles, G. L.↗

Effect of Liquid Droplets on Turbulence Structure in a Round Gaseous Jet

A second-order model which predicts the modulation of turbulence in jets laden with uniform size solid particles or liquid droplets is discussed. The approach followed is to start from the separate momentum and continuity equations of each phase and derive two new conservation equations. The first is for the carrier fluid's kinetic energy of turbulence and the second for the dissipation rate of that energy. Closure of the set of transport equations is achieved by modeling the turbulence correlations up to a third order. The coefficients (or constants) appearing in the modeled equations are then evaluated by comparing the predictions with LDA-measurements obtained recently in a turbulent jet laden with 200 microns solid particles. This set of constants is then used to predict the same jet flow but laden with 50 microns solid particles. The agreement with the measurement in this case is very good.

Elghobashi, S. E.↗

Incompressible Navier-Stokes solutions for a sharp-edged double-delta wing

An implicit finite-difference scheme is used to compute the incompressible laminar vortical flow around a double-delta wing with an aspect ratio of 2.06. By adding a time derivative of the pressure to the continuity equation, the unsteady incompressible Navier-Stokes equations can be integrated like a conventional parabolic time-dependent system of equations. The flux-difference split scheme combines approximate factorization in crossflow planes with a symmetric planar Gauss-Seidel relaxation in the remaining spatial direction. The trajectory of the vortical core is well predicted in comparison with experimental data.

Hsu, Chung-Hao↗

Fully elliptic incompressible flow calculations on regular grid by a new pressure substitution method

A new method is presented for the solution of incompressible flow in generalized coordinates. This method is based on the substitution of the pressure weighted form of the momentum equations into the continuity equation. The algorithm is rigorously derived and a Fourier analysis is used to assess its suitability to act as an error smoother. Linear stability analysis results indicate that the performance of the new pressure substitution method (PSM) and the pressure correction method (PCM) is about the same at low Reynolds numbers, with no significant pressure gradient. At high Reynolds numbers the PSM shows much faster convergence. Likewise prediction of various flows indicate that the PSM has better accuracy for high Reynolds number flows with significant pressure gradients. Since most practical aerodynamic flows have significant pressure gradients, the PSM seems to be attractive for such flows. Solutions for both laminar and turbulent flow are compared with the experimental data. A two-equation low Reynolds number turbulence model is used to resolve the turbulent flowfield.

Hobson, G. V.↗

Prediction of vortical flows on wings using incompressible Navier-Stokes equations

Numerical solutions to the incompressible Navier-Stokes equations are in greater demand than ever before as the field of computational fluid dynamics (CFD) increases its impact as an engineering tool. An implicit upwind-relation finite-difference scheme is developed for solving unsteady incompressible Navier-Stokes equations incorporating an artificial time derivative of the pressure to the continuity equation. The essential features observed in experiments for complicated vortical flows past three round-edged low-aspect-ratio wings at high angles of attack are reproduced both qualitatively and quantitatively.

Hsu, C.-H.↗

Acoustic radiation damping of flat rectangular plates subjected to subsonic flows

The acoustic radiation damping for various isotropic and laminated composite plates and semi-infinite strips subjected to a uniform, subsonic and steady flow has been predicted. The predictions are based on the linear vibration of a flat plate. The fluid loading is characterized as the perturbation pressure derived from the linearized Bernoulli and continuity equations. Parameters varied in the analysis include Mach number, mode number and plate size, aspect ratio and mass. The predictions are compared with existing theoretical results and experimental data. The analytical results show that the fluid loading can significantly affect realistic plate responses. Generally, graphite/epoxy and carbon/carbon plates have higher acoustic radiation damping values than similar aluminum plates, except near plate divergence conditions resulting from aeroelastic instability. Universal curves are presented where the acoustic radiation damping normalized by the mass ratio is a linear function of the reduced frequency. A separate curve is required for each Mach number and plate aspect ratio. In addition, acoustic radiation damping values can be greater than or equal to the structural component of the modal critical damping ratio (assumed as 0.01) for the higher subsonic Mach numbers. New experimental data were acquired for comparison with the analytical results.

Lyle, Karen Heitman↗

Pseudospectral collocation methods for fourth order differential equations

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Malek, Alaeddin↗

Ground Based, Millimeter Wave Measurement of Ozone in the Middle Atmosphere

There is a need for highly reliable measurements of stratospheric ozone. Policy makers worldwide concerned with public health rely oil a clear consensus from the scientific community as a basis for ozone-related environmental policy that has a significant impact oil national economies. Tile latest Such consensus was presented in WMO, and used in a 1999 meeting of the parties considering amendments to the Montreal Protocol oil Substances that Deplete the Ozone Layer. Tile scientific community, in turn, needs highly precise and accurate measurements of ozone levels, and small time derivatives of these levels, both in continued - development of its understanding of the physical and chemical processes involved and as clear evidence that these processes are occurring as stated. Over most of the world, changes in ozone levels are small. For example, over the heavily populated northern midlatitudes, the linearized rate of ozone decline is between 0.2% per year and 0.7% per year, depending on altitude. These values are small enough to make measurement requirements technically challenging. Data quality may suffer from imperfections in individual instruments. In one instance, early results from a satellite-borne ozone sensor were later found to be invalid because of calibration drift. Even in the absence of drift, tile absolute calibration of a new sensor may differ slightly from that of its predecessor in service. Most ozone remote sensing instruments operate at ultraviolet or infrared wavelengths where scattering from dust and aerosols must be taken into account; results from these systems may be or are affected following a major volcanic eruption, such as tile one at Mt. Pinatubo in 1991. Given these difficulties, a consensus of measurements from several independent systems is required to insure a reliable understanding of stratospheric ozone levels. Because of the above-described need for highly precise and accurate ozone measurements using several independent techniques, there was interest in developing several techniques which were known but not highly developed in the 1980's into systems capable of being used in ail operational manner to make measurements with the level of quality needed to-detect small trends in ozone levels. A ground-based microwave instrument capable of remotely sensing stratospheric ozone had been designed by tile Principal Investigator of the present project. This instrument was built at tile Millitech Corporation in South Deerfield, Massachusetts before tile present work began. (Funding for design and construction of the instrument came from sources other than the present grant.) Tile instrument measured the spectrum of one of the many emission lines produced by purely rotational transitions of ozone, one at a frequency of 110.8 GHz. The altitude distribution of ozone can. in principle, be retrieved from the details of the pressure-broadened spectrum of the ozone transition. However, the level of contamination of the spectral measurement by instrumentally induced artifacts must be very low in order to retrieve a ozone profile of useful quality from it. The Millitech instrument demonstrated spectral purity at ail adequate level, and there were promising ideas for instrumental improvements and for further development of the technique. The initial objectives of the present project, then, were to develop techniques for calibrating the Millitech instrument, to minimize artifacts in tile spectra it produces, to optimally retrieve ozone profiles from tile spectra, to test tile quality of the microwave profiles by comparing them with profiles obtained using several other, independent techniques over both short and periods of time, and to perform research using the ozone data gathered with the instrument.

Parrish, Alan↗

An Analysis of the NOAA Satellite-Derived Snow-Cover Record, 1972 - Present

The large-scale distribution of snow cover over northern hemisphere lands has been a topic of increasing attention in recent years. This interest has been spurred, at least in part, by concerns associated with potential changes in the global climate system associated with anthropogenic and natural causes. Satellite observations using visible satellite imagery permit a hemispheric analysis of snow extent. For almost three decades the National Oceanic and Atmospheric Administration (NOAA) has been using visible imagery to produce weekly charts depicting the extent of snow cover over northern hemisphere lands. These charts constitute the longest satellite-derived environmental dataset available on a continuous basis and produced in a consistent manner. We will briefly describe the NOAA charts and then provide an update on the variability of snow extent over the hemisphere from January 1972 through August 1995. Concentration will be on snow kinematics.

Robinson, David A.↗

Modeling Turbulent Combustion for Variable Prandtl and Schmidt Number

This report consists of two abstracts submitted for possible presentation at the AIAA Aerospace Science Meeting to be held in January 2005. Since the submittal of these abstracts we are continuing refinement of the model coefficients derived for the case of a variable Turbulent Prandtl number. The test cases being investigated are a Mach 9.2 flow over a degree ramp and a Mach 8.2 3-D calculation of crossing shocks. We have developed an axisymmetric code for treating axisymmetric flows. In addition the variable Schmidt number formulation was incorporated in the code and we are in the process of determining the model constants.

Hassan, H. A.↗

Ejecta- and Size-Scaling Considerations from Impacts of Glass Projectiles into Sand

One of the most promising means of learning how initial impact conditions are related to the processes leading to the formation of a planetary-scale crater is through scaling relationships.1,2,3 The first phase of deriving such relationships has led to great insight into the cratering process and has yielded predictive capabilities that are mathematically rigorous and internally consistent. Such derivations typically have treated targets as continuous media; in many, cases, however, planetary materials represent irregular and discontinuous targets, the effects of which on the scaling relationships are still poorly understood.4,5 We continue to examine the effects of varying impact conditions on the excavation and final dimensions of craters formed in sand. Along with the more commonly treated variables such as impact speed, projectile size and material, and impact angle,6 such experiments also permit the study of changing granularity and friction angle of the target materials. This contribution presents some of the data collected during and after the impact of glass spheres into a medium-grained sand.

Anderson J. L. B.↗

Full Multigrid Flow Solver

FMG3D (full multigrid 3 dimensions) is a pilot computer program that solves equations of fluid flow using a finite difference representation on a structured grid. Infrastructure exists for three dimensions but the current implementation treats only two dimensions. Written in Fortran 90, FMG3D takes advantage of the recursive subroutine feature, dynamic memory allocation, and structured-programming constructs of that language. FMG3D supports multi-block grids with three types of block-to-block interfaces: periodic, C-zero, and C-infinity. For all three types, grid points must match at interfaces. For periodic and C-infinity types, derivatives of grid metrics must be continuous at interfaces. The available equation sets are as follows: scalar elliptic equations, scalar convection equations, and the pressure-Poisson formulation of the Navier-Stokes equations for an incompressible fluid. All the equation sets are implemented with nonzero forcing functions to enable the use of user-specified solutions to assist in verification and validation. The equations are solved with a full multigrid scheme using a full approximation scheme to converge the solution on each succeeding grid level. Restriction to the next coarser mesh uses direct injection for variables and full weighting for residual quantities; prolongation of the coarse grid correction from the coarse mesh to the fine mesh uses bilinear interpolation; and prolongation of the coarse grid solution uses bicubic interpolation.

Mineck, Raymond E.↗

A Theoretical and Computational Revisit of Conversions Between Whitham’s F-Function and Equivalent Area

This paper proves mathematically that the integral transforms between Whitham’s F-function and equivalent area are the inverse transforms of each other if and only if the slope of the equivalent area at the origin is zero. This mathematical fact contradicts the accepted unconditional inverse relation between Whitham’s F-function and equivalent area in the sonic boom research literature. Piecewise linear approximations of an F-function and of the second derivative of an equivalent area are used to derive numerical formulas for conversions between Whitham’s F-function and equivalent area. Numerical results are included to show convergence of the numerical conversions as the maximum length of the segments for piecewise linear approximations goes to zero. These numerical conversions are approximately the inverse transforms of each other when the second derivative of an equivalent area is continuous and the slope of the equivalent area at the origin is zero.

Whitham's F-function↗

Explicit asymmetric bounds for robust stability of continuous and discrete-time systems

The problem of robust stability in linear systems with parametric uncertainties is considered. Explicit stability bounds on uncertain parameters are derived and expressed in terms of linear inequalities for continuous systems, and inequalities with quadratic terms for discrete-times systems. Cases where system parameters are nonlinear functions of an uncertainty are also examined.

Gao, Zhiqiang↗

Light-ray wave functions and integrability

Using integrability, we construct (to leading order in perturbation theory) the explicit form of twist-three light-ray operators in planar $\mathcal{N}$ = 4 SYM. This construction allows us to directly compute analytically continued CFT data at complex spin. We derive analytically the “magic” decoupling zeroes previously observed numerically. Using the Baxter equation, we also show that certain Regge trajectories merge together into a single unifying Riemann surface. Perhaps more surprisingly, we find that this unification of Regge trajectories is not unique. If we organize twist-three operators differently into what we call “cousin trajectories” we find infinitely more possible continuations. We speculate about which of these remarkable features of twist-three operators might generalize to other operators, other regimes and other theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Interface crack in a nonhomogeneous elastic medium

The linear elasticity problem for an interface crack between two bonded half planes is reconsidered. It is assumed that one of the half planes is homogeneous and the second is nonhomogeneous in such a way that the elastic properties are continuous throughout the plane and have discontinuous derivatives along the interface. The problem is formulated in terms of a system of integral equations and the asymptotic behavior of the stress state near the crack tip is determined. The results lead to the conclusion that the singular behavior of stresses in the nonhomogeneous medium is identical to that in a homogeneous material provided the spacial distribution of material properties is continuous near and at the crack tip. The problem is solved for various values of the nonhomogeneity parameter and for four different sets of crack surface tractions, and the corresponding stress intensity factors are tabulated.

Delale, F.↗

The Effect of Functional Groups in Bio-Derived Fuel Candidates

Interest in developing renewable fuels is continuing to grow and biomass represents a viable source of renewable carbon with which to replace fossil-based components in transportation fuels. During our own work, we noticed that chemists think in terms of functional groups whereas fuel engineers think in terms of physical fuel properties. In this Concept article, we discuss the effect of carbon and oxygen functional groups on potential fuel properties. This serves as a way of informing our own thinking and provides us with a basis with which to design and synthesize molecules from biomass that could provide useful transportation fuels.

09 BIOMASS FUELS↗

Further information on structure of the atmosphere of Venus derived from the Vega Venus balloon and lander mission

Continued analysis of the pressure and temperature data returned by the two Vega mission balloons has revealed an apparently significant difference in mean atmospheric static stability between the two data sets. Furthermore, the stability is time dependent within each data set. The 6.5 K temperature contrast between the two balloons remains, and appears to have a counterpart in the contrast between two of the Pioneer Venus probes at these levels, which has been attributed to planetary scale waves. Comparisons of the Vega 2 lander data with those of the Pioneer Venus Large Probe shows relatively close agreement in the state properties and in the atmospheric static stability profiles as well.

Seiff, A.↗