Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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

Radiative transfer in a sphere illuminated by a parallel beam - An integral equation approach

The problem of multiple scattering of nonpolarized light in a planetary body of arbitrary shape illuminated by a parallel beam is formulated using the integral equation approach. There exists a simple functional whose stationarity condition is equivalent to solving the equation of radiative transfer and whose value at the stationary point is proportional to the differential cross section. The analysis reveals a direct relation between the microscopic symmetry of the phase function for each scattering event and the macroscopic symmetry of the differential cross section for the entire planetary body, and the interconnection of these symmetry relations and the variational principle. The case of a homogeneous sphere containing isotropic scatterers is investigated in detail. It is shown that the solution can be expanded in a multipole series such that the general spherical problem is reduced to solving a set of decoupled integral equations in one dimension. Computations have been performed for a range of parameters of interest, and illustrative examples of applications to planetary problems as provided.

Shia, R.-L.↗

Time domain unsteady incompressible cascade airfoil theory for helicopter rotors in hover

A detailed derivation of a finite-time arbitrary-motion incompressible cascade theory is presented for both Laplace and frequency domains. The generalized cascade lift deficiency function is shown to be consistent with the generalized Theodorsen's lift deficiency function when the wake spacing approaches infinity or when the reduced frequency tends to infinity; it also yields a correct value for the zero reduced frequency limit. Efficient numerical procedures are presented for the evaluation of the cascade lift deficiency function. Numerical examples comparing the cascade lift deficiency function with Loewy's (1957) rotary-wing lift deficiency function in frequency domain are presented. Pade approximants of the cascade lift deficiency function are constructed using a Bode plot approach that allows for complex poles. A general-purpose optimization program is used to determine the coefficients of the approximant.

Dinyavari, M. A. H.↗

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗

Some Experiences with High Order Discretization Schemes for Compressible Fluid Flow

It is widely recognized that the structure and complexity of compressible fluid flow places high demands on numerical discretization techniques for the fluid flow equations. Fluid flows arising in external aerodynamics often contain both flow field discontinuities and fluid boundary-layers. Both must be accurately resolved to provide useful information to aerodynamic design and analysis engineers. These accuracy requirements motivated the present author to examine a class of finite-volume techniques on arbitrary triangulated domains based on linear or quadratic reconstruction of integral-averaged data followed by upwind flux function evaluation and small time evolution. More recently, we have considered some new upwind techniques which yield compact discretizations while maintaining higher order accuracy. In the mini-symposium talk we will discuss both of these techniques as well as demonstrate the relative merits of each method by computing a number of aerodynamic flows containing shock waves and boundary-layers.

Barth, Timothy J.↗

A Systematic Methodology for Constructing High-Order Energy-Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter (AIAA 2008-2876, 2008) was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables \energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗

High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables 'energy stable' modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗

A Systematic Methodology for Constructing High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non{Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter [1] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables "energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗

Chordwise and compressibility corrections for arbitrary planform slender wings

The Lomax and Sluder method for adapting slender-wing theory to delta or rectangular wings by making chordwise and compressibility corrections is extended to cover wings of any arbitrary planform in subsonic and supersonic flows. The numerical accuracy of the present work is better than that of the Lomax-Sluder results. Comparison of the results of this work with those of the vortex-lattice method and Kernel function method for a family of Gothic and arrowhead wings shows good agreement. A universal curve is proposed for the evaluation of the lift coefficient of a low aspect ratio wing of an arbitrary planform in subsonic flow. The location of the center of pressure can also be estimated.

Levin, D.↗

Sum rules and static local-field corrections of electron liquids in two and three dimensions

The dielectric functions of electron liquids which take into account short-range electron-electron correlations via the static local-field corrections are examined in the light of the frequency-moment sum rules. The formation is given for degenerate as well as classical electron liquids in arbitrary (d) spatial dimensions, which is suitable for comparison between the two- and three-dimensional cases. By using the virial equations of state it is shown that such dielectric functions cannot satisfy the compressibility sum rule and the third-frequency-moment sum rule simultaneously. In the degenerate case, the plasmon, single-pair, and multipair contributions to the sum rules are analyzed, and the reason for this incompatibility is discussed.

Iwamoto, N.↗

On hierarchical solutions to the BBGKY hierarchy

It is thought that the gravitational clustering of galaxies in the universe may approach a scale-invariant, hierarchical form in the small separation, large-clustering regime. Past attempts to solve the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy in this regime have assumed a certain separable hierarchical form for the higher order correlation functions of galaxies in phase space. It is shown here that such separable solutions to the BBGKY equations must satisfy the condition that the clustered component of the solution has cluster-cluster correlations equal to galaxy-galaxy correlations to all orders. The solutions also admit the presence of an arbitrary unclustered component, which plays no dyamical role in the large-clustering regime. These results are a particular property of the specific separable model assumed for the correlation functions in phase space, not an intrinsic property of spatially hierarchical solutions to the BBGKY hierarchy. The observed distribution of galaxies does not satisfy the required conditions. The disagreement between theory and observation may be traced, at least in part, to initial conditions which, if Gaussian, already have cluster correlations greater than galaxy correlations.

Hamilton, A. J. S.↗

A half plane and a strip with an arbitrarily located crack

A technique is presented for dealing with the problem of an elastic domain containing an arbitrarily oriented internal crack. The problem is formulated as a system of integral equations for a fictitious layer of body forces imbedded in the plane along a closed smooth curve encircling the original domain. The problems of a half plane with a crack in the neighborhood of its free boundary and of an infinite strip containing a symmetrically located internal crack with an arbitrary orientation are considered as examples. In each case the stress intensity factors are computed and are given as functions of the crack angle.

Erdogan, F.↗

Proton dose approximation in arbitrary convex geometry

An expansion is derived for the solution to the transport equation in two dimensions subject to boundary conditions given for an arbitrary convex region. Questions of high-energy transport are considered along with the properties of the dose response function. The expansion of the solution of the transport equation is presented in terms of a parameter which measures the lateral dispersion of an unidirectional beam. This parameter is usually small and the expansion is expected to converge rapidly. The dominant term in the expansion is related to fluence-to-dose conversion factors in a semiinfinite slab for normal incidence. A convenient parameterization of the conversion factors is provided along with numerical examples.

Wilson, J. W.↗

A half plane and a strip with an arbitrarily located crack

The paper introduces a technique to deal with the problem of an elastic domain containing an arbitrarily oriented internal crack. The problem is formulated as a system of integral equations for a fictitious layer of body forces imbedded in the plane along a closed smooth curve encircling the original domain. The problems of a half plane with a crack in the neighborhood of its free boundary and of an infinite strip containing a symmetrically located internal crack with an arbitrary orientation are considered as examples. In each case the stress intensity factors are computed and are given as functions of the crack angle.

Erdogan, F.↗

A new determination of the halo luminosity density of the galaxy

The distribution of tangential velocities for arbitrary halo velocity dispersions and rotation rates are explicitly computed, assuming a particular form of the halo velocity distribution function near the sun. It is found that for a specific model of the halo velocity dispersion tensor that Schmidt's (1975) estimate of the mass carried by halo stars of 0.25-0.70 solar masses should be revised by a factor of between 1.0 and 2.0, depending on the rotation rate of the halo. Although obtaining a mass density in all stars from the stellar density in that mass range is extremely uncertain, it is found, in agreement with earlier work, that even under rather extreme assumptions about the remnants the observed stars and expected remnants fail to provide the theoretically expected mass. Comparison of the luminosity to the theoretical mass yields a mass to light ratio consistent with that of rich clusters of galaxies.

Richstone, D. O.↗

An accurate method for evaluating the kernel of the integral equation relating lift to downwash in unsteady potential flow

The method is capable of generating approximations of arbitrary accuracy. It is based on approximating the algebraic part of the nonelementary integrals in the kernel by exponential functions and then integrating termwise. The exponent spacing in the approximation is a geometric sequence. The coefficients and exponent multiplier of the exponential approximation are computed by least squares so the method is completely automated. Exponential approximates generated in this manner are two orders of magnitude more accurate than the exponential approximation that is currently most often used for this purpose. The method can be used to generate approximations to attain any desired trade-off between accuracy and computing cost.

Desmarais, R. N.↗

NLTE model stellar atmospheres with line blanketing near the series limits

In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.

Hubeny, I.↗

An Electrothermodynamic Model with Distributed Properties for Effective Permittivities of Sea Ice

This paper presents a model to calculate the temperature dependence of effective permittivities for sea ice, a heterogeneous medium containing multiphase scatterers. With the strong permittivity fluctuation approach the model accounts for the electrodynamic scattering effect together with the quasi-static characteristics of multiple species and subspecies of inhomogeneities with distributed orientations, sizes, and shapes. Because of a preferential direction in the orientation distribution, the medium is effectively anisotropic. The size distribution is described with a probability density function in terms of normalized volumetric sizes. Scatterer shapes are nonuniform and have a general ellipsoidal form characterized by arbitrary axial ratios of correlation lengths which are related to physical geometries of the scatterers. In this formulation, sea ice consisting of solid ice, liquid brine, and gaseous inclusions is modeled to derive effective permittivities with thermodynamic phase redistribution and structural metamorphism. Theoretical results are in good agreement with experimental data at the C band frequency of 4.8 GHz for saline ice undergoing warming and cooling cycles. A competitive effect between the increase of liquid brine and the shape rounding of ellipsoidal scatterers at increasing temperatures explains the trend observed in measured data. Sensitivities of effective permittivities to structural and physical parameters characterizing sea ice are also studied.

Nghiem, S. V.↗