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 91 records · Page 5

Radiative transfer calculated from a Markov chain formalism

The theory of Markov chains is used to formulate the radiative transport problem in a general way by modeling the successive interactions of a photon as a stochastic process. Under the minimal requirement that the stochastic process is a Markov chain, the determination of the diffuse reflection or transmission from a scattering atmosphere is equivalent to the solution of a system of linear equations. This treatment is mathematically equivalent to, and thus has many of the advantages of, Monte Carlo methods, but can be considerably more rapid than Monte Carlo algorithms for numerical calculations in particular applications. We have verified the speed and accuracy of this formalism for the standard problem of finding the intensity of scattered light from a homogeneous plane-parallel atmosphere with an arbitrary phase function for scattering. Accurate results over a wide range of parameters were obtained with computation times comparable to those of a standard 'doubling' routine. The generality of this formalism thus allows fast, direct solutions to problems that were previously soluble only by Monte Carlo methods. Some comparisons are made with respect to integral equation methods.

Esposito, L. W.↗

The general solution to the classical problem of the finite Euler-Bernoulli beam

An analytical solution is obtained for the problem of free and forced vibrations of a finite Euler-Bernoulli beam with arbitrary (partially fixed) boundary conditions. The effects of linear viscous damping, Winkler foundation, constant axial tension, a concentrated mass, and an arbitrary forcing function are included in the analysis. No restriction is placed on the values of the parameters involved, and the solution presented here contains all cited previous solutions as special cases.

Hussaini, M. Y.↗

Transformation of two and three-dimensional regions by elliptic systems

Grid smoothing and orthogonalization procedures were developed and implemented in the construction of two and three dimensional grids. The procedures are based on the variational methods of grid generation. The two-dimensional examples were computed using the MSU IRIS Graphics Workstation. It was demonstrated that the elliptic grid generation equations, with arbitrary forcing functions, can be solved, in their variational formulation, using a gradient method. Since gradient methods have a global convergence property, the divergence problems often encountered when using SOR iterative methods can be avoided. It is not to be concluded, however, that SOR methods should be abandoned, since gradient methods tend to converge very slowly. In fact, slow convergence was the major problem encountered in the three-dimensional grids. Further progress was made on the continuing effort to develop conservative interpolation formulas for overlapping grids.

Mastin, C. Wayne↗

The X-ray background and the evolution of quasars

The contribution of QSOs and active galaxies to the 2-10-keV X-ray background is calculated theoretically, assuming that the X-ray luminosity of each GSO has a power-law time evolution (dL/dt = AL exp alpha) and applying the continuity equation to derive the X-ray luminosity function at arbitrary redshift. The observed X-ray background is shown to require alpha greater than 1.2, ruling out pure luminosity evolution (alpha = 1).

Tucker, W. H.↗

Free vibration of rectangular plates with a small initial curvature

The method of matched asymptotic expansions is used to solve the transverse free vibration of a slightly curved, thin rectangular plate. Analytical results for natural frequencies and mode shapes are presented in the limit when the dimensionless bending rigidity, epsilon, is small compared with in-plane forces. Results for different boundary conditions are obtained when the initial deflection is: (1) a polynomial in both directions, and (2) the product of a polynomial and a trigonometric function, and arbitrary. For the arbitrary initial deflection case, the Fourier series technique is used to define the initial deflection. The results obtained show that the natural frequencies of vibration of slightly curved plates are coincident with those of perfectly flat, prestressed rectangular plates. However, the eigenmodes are very different from those of initially flat prestressed rectangular plates. The total deflection is found to be the sum of the initial deflection, the deflection resulting from the solution of the flat plate problem, and the deflection resulting from the static problem.

Adeniji-Fashola, A. A.↗

Systematic methods for calculation of the dielectric properties of arbitrary plasmas

In the novel approach presented for calculating the dispersion integrals needed for determining plasma dielectric properties, the dispersion integrals for an arbitrary distortion function with a continuous derivative are systematically expanded in terms of a set of orthogonal functions whose corresponding dispersion functions are already known. This general approach is, on the one hand, implemented for unmagnetized plasmas, and on the other generalized to treat relativistic and magnetized plasmas. The method allows the systematic and efficient calculation of dispersion integrals, for the cases of either real or complex arguments.

Robinson, P. A.↗

A finite element formulation for scattering from electrically large 2-dimensional structures

A finite element formulation is given using the scattered field approach with a fictitious material absorber to truncate the mesh. The formulation includes the use of arbitrary approximation functions so that more accurate results can be achieved without any modification to the software. Additionally, non-polynomial approximation functions can be used, including complex approximation functions. The banded system that results is solved with an efficient sparse/banded iterative scheme and as a consequence, large structures can be analyzed. Results are given for simple cases to verify the formulation and also for large, complex geometries.

Ross, Daniel C.↗

Fabrication techniques for very fast diffractive lenses

Aspheric lenses with arbitrary phase functions can be fabricated on thin light weight substrates via the binary optics fabrication technique. However, it is difficult and costly to fabricate a fast lens (f/number less than 1) for use as the shorter wavelengths. The pitch of the masks and the alignment accuracy must be very fine. For a large lens, the space-bandwidth product of the element can also become impractically large. In this paper, two alternate approaches for the fabrication of fast aspheric diffractive lenses are described. The first approach fabricates the diffractive lens interferometrically, utilizing a spherical wavefront to provide the optical power of the lens and a computer generated hologram to create the aspheric components. The second approach fabricates the aspheric diffractive lens in the form if a higher order kinoform which trades groove profile fidelity for coarser feature size. The design and implementation issues for these two fabrication techniques are discussed.

Tai, Anthony M.↗

The Schwinger Variational Method

Variational methods have proven invaluable in theoretical physics and chemistry, both for bound state problems and for the study of collision phenomena. The application of the Schwinger variational (SV) method to e-molecule collisions and molecular photoionization has been reviewed previously. The present chapter discusses the implementation of the SV method as applied to e-molecule collisions. Since this is not a review of cross section data, cross sections are presented only to server as illustrative examples. In the SV method, the correct boundary condition is automatically incorporated through the use of Green's function. Thus SV calculations can employ basis functions with arbitrary boundary conditions. The iterative Schwinger method has been used extensively to study molecular photoionization. For e-molecule collisions, it is used at the static exchange level to study elastic scattering and coupled with the distorted wave approximation to study electronically inelastic scattering.

Huo, Winifred M.↗

General Purpose Fortran Program for Discrete-Ordinate-Method Radiative Transfer in Scattering and Emitting Layered Media: An Update of DISORT

This update reports a state-of-the-art discrete ordinate algorithm for monochromatic unpolarized radiative transfer in non-isothermal, vertically inhomogeneous, but horizontally homogeneous media. The physical processes included are Planckian thermal emission, scattering with arbitrary phase function, absorption, and surface bidirectional reflection. The system may be driven by parallel or isotropic diffuse radiation incident at the top boundary, as well as by internal thermal sources and thermal emission from the boundaries. Radiances, fluxes, and mean intensities are returned at user-specified angles and levels. DISORT has enjoyed considerable popularity in the atmospheric science and other communities since its introduction in 1988. Several new DISORT features are described in this update: intensity correction algorithms designed to compensate for the 8-M forward-peak scaling and obtain accurate intensities even in low orders of approximation; a more general surface bidirectional reflection option; and an exponential-linear approximation of the Planck function allowing more accurate solutions in the presence of large temperature gradients. DISORT has been designed to be an exemplar of good scientific software as well as a program of intrinsic utility. An extraordinary effort has been made to make it numerically well-conditioned, error-resistant, and user-friendly, and to take advantage of robust existing software tools. A thorough test suite is provided to verify the program both against published results, and for consistency where there are no published results. This careful attention to software design has been just as important in DISORT's popularity as its powerful algorithmic content.

Tsay, Si-Chee↗

Quantum-Inspired Maximizer

A report discusses an algorithm for a new kind of dynamics based on a quantum- classical hybrid-quantum-inspired maximizer. The model is represented by a modified Madelung equation in which the quantum potential is replaced by different, specially chosen 'computational' potential. As a result, the dynamics attains both quantum and classical properties: it preserves superposition and entanglement of random solutions, while allowing one to measure its state variables, using classical methods. Such optimal combination of characteristics is a perfect match for quantum-inspired computing. As an application, an algorithm for global maximum of an arbitrary integrable function is proposed. The idea of the proposed algorithm is very simple: based upon the Quantum-inspired Maximizer (QIM), introduce a positive function to be maximized as the probability density to which the solution is attracted. Then the larger value of this function will have the higher probability to appear. Special attention is paid to simulation of integer programming and NP-complete problems. It is demonstrated that the problem of global maximum of an integrable function can be found in polynomial time by using the proposed quantum- classical hybrid. The result is extended to a constrained maximum with applications to integer programming and TSP (Traveling Salesman Problem).

Zak, Michail↗

Three-Dimensional, Space-Dependent Mesoscale Diffusivity: Derivation and Implications

Recently, we presented a parameterization of an arbitrary tracer 3D mesoscale flux that describes both diabatic and adiabatic regimes without using arbitrary tapering functions. However, we did not parameterize the mesoscale diffusivity, which is the subject of this work. A key difference between the present and previous diffusivity parameterizations is that in the latter, the two main ingredients, mesoscale drift velocity and eddy kinetic energy, were not parameterized but determined using present data, which deprives the models of predictive power. Since winds, stratification, etc., are predicted to change in the future, use of these parameterizations to study future climate scenarios becomes questionable. In this work, we parameterize drift velocity and eddy kinetic energy (vertical–horizontal components), which we first assess with data [WOCE, TOPEX/Poseidon (T/P), and North Atlantic Tracer Release Experiment (NATRE)] and then use in a coarse-resolution stand-alone ocean code under Coordinated Ocean-Ice Reference Experiment I (CORE-I) forcing. We present results for the global ocean temperature and salinity, Atlantic overturning circulation, meridional heat transport, and Drake Passage transport, which we compare with several previous studies. The temperature drift is less than that of five of seven previous OGCMs, and the salinity drift is among the smallest in those studies. The predicted winter Antarctic Circumpolar Current mixed layer depths (MLDs) are in good agreement with the data. Predicting the correct MLD is important in climate studies since models that predict very deep mixed layers transfer more of the radiative perturbation to the deep ocean, reducing surface warming (and vice versa).

Eddies; Mesoscale processes; Mixing; Ocean models;↗

The heating of a thermally conducting stratified medium. II - A simple plane model of an atmosphere

Exact solutions of the following theoretical problem are presented: A plane atmosphere is in hydrostatic equilibrium with a uniform gravity. The ideal gas law is assumed. Heat is generated everywhere at a rate proportional to the local density. The atmosphere is maintained in a steady state through cooling by thermal conduction and radiation. This problem is reducible to quadratures for a thermal conductivity which is an arbitrary, but prescribed, function of the temperature, and for a radiative loss which is expressible as the product of the density and an arbitrary, but prescribed, function of the pressure. The analysis is carried out for the case of power law thermal conductivity, and a radiative loss proportional to the square of the density and to the first power of the temperature. The radiative cooling function adopted here has the basic mathematical form for an optically thin medium. The solutions reproduce the macroscopic ordering of a hot 'corona' separated from a 'photosphere' by a layer of temperature minimum. The analytic solutions allow direct illustration of the interplay between steady energy transport and the requirements of hydrostatic equilibrium.

Lerche, I.↗

Random search optimization based on genetic algorithm and discriminant function

The general problem of optimization with arbitrary merit and constraint functions, which could be convex, concave, monotonic, or non-monotonic, is treated using stochastic methods. To improve the efficiency of the random search methods, a genetic algorithm for the search phase and a discriminant function for the constraint-control phase were utilized. The validity of the technique is demonstrated by comparing the results to published test problem results. Numerical experimentation indicated that for cases where a quick near optimum solution is desired, a general, user-friendly optimization code can be developed without serious penalties in both total computer time and accuracy.

Kiciman, M. O.↗

A New Merit Function for Evaluating the Flaw Tolerance of Composite Laminates: Arbitrary Size Holes and Center Cracks - Part 2

In a previous paper, a new merit function for determining the strength performance of flawed composite laminates was presented. This previous analysis was restricted to circular hole flaws that were large enough that failure could be predicted using the laminate stress concentration factor. In this paper, the merit function is expanded to include the flaw cases of an arbitrary size circular hole or center crack. Failure prediction for these cases is determined using the point stress criterion. An example application of the merit function is included for a wide range of graphite/epoxy laminates.

Martin, Mikulas M., Jr.↗

A New Merit Function for Evaluating the Flaw Tolerance of Composite Laminates: Arbitrary Size Holes and Center Cracks - Part 2

In a previous paper, a new merit function for determining the strength performance of flawed composite laminates was presented. This previous analysis was restricted to circular hole flaws that were large enough that failure could be predicted using the laminate stress concentration factor. In this paper, the merit function is expanded to include the flaw cases of an arbitrary size circular hole or a center crack. Failure prediction for these cases is determined using the point stress criterion. An example application of the merit function is included for a wide range of graphite/epoxy laminates.

Mikulas, Martin M., Jr.↗