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 415 records · Page 23

Universal Location of the Yang-Lee Edge Singularity in O( N ) Theories

We determine a previously unknown universal quantity, the location of the Yang-Lee edge singularity for the O(N) theories in a wide range of N and various dimensions. At large N, we reproduce the N→∞ analytical result on the location of the singularity and, additionally, we obtain the mean-field result for the location in d=4 dimensions. In order to capture the nonperturbative physics for arbitrary N, d and complex-valued external fields, we use the functional renormalization group approach.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cutting Corners: Curvilinear-Surface-Based Gravity Models for Asteroids and Comets

We report that asteroids and comets are often highly irregular in shape and can contain large density inhomogenities. For such bodies, the simple point-mass gravity model does not capture the dynamics of the environment and this has driven a need for higher fidelity alternatives. A number of approaches have been developed with one of the most popular being the analytic polyhedral model. The analytic polyhedral model approximates the irregular shaped body as a polyhedron, typically consisting of triangular facets, with some assumed internal density distribution. Several analytic formulations have been derived for the gravitational fields of a constant-density polyhedron. A thorough history of the model is provided in Ref. [6] as well as a comparison of the implementations of Refs. [2–4]. The three versions provide near identical accuracy with the implementation of Werner requiring the fewest transcendental function evaluations suggesting its superior efficiency. Implementations with linear density contrasts and arbitrary polynomial density contrasts have also been developed.

79 ASTRONOMY AND ASTROPHYSICS↗

Fully unsteady subsonic and supersonic potential aerodynamics of complex aircraft configurations for flutter applications

The paper presents a general theory of steady, oscillatory or fully unsteady, subsonic and supersonic aerodynamics for aircraft of arbitrary shape. The finite element method is used to obtain the time functional relationship between aerodynamic potential and its normal derivative (normal wash) in a form suitable for computational analysis. The matrix of the aerodynamic influence coefficients, as necessary for flutter calculations, is obtained as a result. The theory is implemented in a computer program, and some typical results obtained with this program are presented.

Tseng, K.↗

Gravitational spectra from the tracking of planetary spacecraft in eccentric orbits

Two dimensional gravitational spectra are derived from simple harmonic analysis of range rate tracking data on planetary orbiters. The eccentricity of the orbit is arbitrary and results are shown to vary substantially with the aspect angle of the tracking line of sight with the orbit plane. The development for arbitrary start with stop times (with respect to periapsis) uses modified eccentricity functions evaluated by quadrature. Simulations with a point-masses model of Venus using tracking data on the Pioneer Venus Orbiter show excellent predictions of the average orbiter spectrum over one Venus day. The Venus gravitational signal should be above the tracking noise level for arc lengths longer than 40 deg (in true anomaly) about periapsis and for terms as high as 55th degree. analysis has been made of tracking residuals from a short arc fit to Mariner Mars 9 data over the Hellas Basin (using a complete 6th degree field). Results are most consistent with higher residual gravitational power than predicted from Kaula's rule for Mars.

Wagner, C. A.↗

The axisymmetric elasticity problem for a laminated plate containing a circular hole

The elasticity problem for a laminated thick plate which consists of two bonded dissimilar layers and which contains a circular hole is considered. The problem is formulated for arbitrary axisymmetric tractions on the hole surface by using the Love strain function. Through the expansion of the boundary conditions into Fourier series the problem is reduced to an infinite system of algebraic equations which is solved by the method of reduction. Of particular interest in the problem are the stresses along the interface as they relate to the question of delamination failure of the composite plate. These stresses are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singular behavior of the stress state is presented and the stress intensity factors are calculated.

Delale, F.↗

Radiative transfer in multilayered random medium with laminar structure - Green's function approach

For a multilayered random medium with a laminar structure a Green's function approach is introduced to obtain the emitted intensity due to an arbitrary point source. It is then shown that the approach is applicable to both active and passive remote sensing. In active remote sensing, the computed radar backscattering cross section for the multilayered medium includes the effects of both volume multiple scattering and surface multiple scattering at the layer boundaries. In passive remote sensing, the brightness temperature is obtained for arbitrary temperature profiles in the layers. As an illustration the brightness temperature and reflectivity are calculated for a bounded layer and compared with results in the literature.

Karam, M. A.↗

On the generation of magnetohydrodynamic waves in a stratified and magnetized fluid. I - Vertical propagation

The generation of MHD waves by turbulent motions in a stratified medium with an embedded uniform magnetic field, a topic which is relevant to the study of the solar atmosphere, is considered. Both compressible and incompressible MHD waves are treated in a one-dimensional approach; however, the direction of the background magnetic field is permitted to vary in an arbitrary direction. Theoretical expressions for MHD energy fluxes are obtained as a function of wave frequency and multipole coefficients. It is shown that monopole, dipole, and quadrupole emissions are responsible for the generation of the compressible components of the fast and slow modes. However, the incompressible components and the Alfven modes can be generated by the dipole emission only. Specific results obtained for special magnetic field geometries are discussed for the fast and slow modes.

Musielak, Z. E.↗

Automated generation of influence functions for planar crack problems

A numerical procedure for the generation of influence functions for Mode I planar problems is described. The resulting influence functions are in a form for convenient evaluation of stress-intensity factors for complex stress distributions. Crack surface displacements are obtained by a least-squares solution of the Williams eigenfunction expansion for displacements in a cracked body. Discrete values of the influence function, evaluated using the crack surface displacements, are curve fit using an assumed functional form. The assumed functional form includes appropriate limit-behavior terms for very deep and very shallow cracks. Continuous representation of the influence function provides a convenient means for evaluating stress-intensity factors for arbitrary stress distributions by numerical integration. The procedure is demonstrated for an edge-cracked strip and a radially cracked disk. Comparisons with available published results demonstrate the accuracy of the procedure.

Sire, Robert A.↗

Stress Intensity Factors for Part-Through Surface Cracks in Hollow Cylinders

Flaws resulting from improper welding and forging are usually modeled as cracks in flat plates, hollow cylinders or spheres. The stress intensity factor solutions for these crack cases are of great practical interest. This report describes some recent efforts at improving the stress intensity factor solutions for cracks in such geometries with emphasis on hollow cylinders. Specifically, two crack configurations for cylinders are documented. One is that of a surface crack in an axial plane and the other is a part-through thumb-nail crack in a circumferential plane. The case of a part-through surface crack in flat plates is used as a limiting case for very thin cylinders. A combination of the two cases for cylinders is used to derive a relation for the case of a surface crack in a sphere. Solutions were sought which cover the entire range of the geometrical parameters such as cylinder thickness, crack aspect ratio and crack depth. Both the internal and external position of the cracks are considered for cylinders and spheres. The finite element method was employed to obtain the basic solutions. Power-law form of loading was applied in the case of flat plates and axial cracks in cylinders and uniform tension and bending loads were applied in the case of circumferential (thumb-nail) cracks in cylinders. In the case of axial cracks, the results for tensile and bending loads were used as reference solutions in a weight function scheme so that the stress intensity factors could be computed for arbitrary stress gradients in the thickness direction. For circumferential cracks, since the crack front is not straight, the above technique could not be used. Hence for this case, only the tension and bending solutions are available at this time. The stress intensity factors from the finite element method were tabulated so that results for various geometric parameters such as crack depth-to-thickness ratio (a/t), crack aspect ratio (a/c) and internal radius-to-thickness ratio (R/t) or the crack length-to-width ratio (2c/W) could be obtained by interpolation and extrapolation. Such complete tables were then incorporated into the NASA/FLAGRO computer program which is widely used by the aerospace community for fracture mechanics analysis.

Mettu, Sambi R.↗

User's Manual for FEMOM3DR

FEMoM3DR is a computer code written in FORTRAN 77 to compute radiation characteristics of antennas on 3D body using combined Finite Element Method (FEM)/Method of Moments (MoM) technique. The code is written to handle different feeding structures like coaxial line, rectangular waveguide, and circular waveguide. This code uses the tetrahedral elements, with vector edge basis functions for FEM and triangular elements with roof-top basis functions for MoM. By virtue of FEM, this code can handle any arbitrary shaped three dimensional bodies with inhomogeneous lossy materials; and due to MoM the computational domain can be terminated in any arbitrary shape. The User's Manual is written to make the user acquainted with the operation of the code. The user is assumed to be familiar with the FORTRAN 77 language and the operating environment of the computers on which the code is intended to run.

Reddy, C. J.↗

Finite Element Analysis of Reverberation Chambers

The primary motivating factor behind the initiation of this work was to provide a deterministic means of establishing the validity of the statistical methods that are recommended for the determination of fields that interact in -an avionics system. The application of finite element analysis to reverberation chambers is the initial step required to establish a reasonable course of inquiry in this particularly data-intensive study. The use of computational electromagnetics provides a high degree of control of the "experimental" parameters that can be utilized in a simulation of reverberating structures. As the work evolved there were four primary focus areas they are: 1. The eigenvalue problem for the source free problem. 2. The development of a complex efficient eigensolver. 3. The application of a source for the TE and TM fields for statistical characterization. 4. The examination of shielding effectiveness in a reverberating environment. One early purpose of this work was to establish the utility of finite element techniques in the development of an extended low frequency statistical model for reverberation phenomena. By employing finite element techniques, structures of arbitrary complexity can be analyzed due to the use of triangular shape functions in the spatial discretization. The effects of both frequency stirring and mechanical stirring are presented. It is suggested that for the low frequency operation the typical tuner size is inadequate to provide a sufficiently random field and that frequency stirring should be used. The results of the finite element analysis of the reverberation chamber illustrate io-W the potential utility of a 2D representation for enhancing the basic statistical characteristics of the chamber when operating in a low frequency regime. The basic field statistics are verified for frequency stirring over a wide range of frequencies. Mechanical stirring is shown to provide an effective frequency deviation.

Bunting, Charles F.↗

A Fast Butterfly-Compressed Hadamard–Babich Integrator for High-Frequency Helmholtz Equations in Inhomogeneous Media with Arbitrary Sources

Here we present a butterfly-compressed representation of the Hadamard-Babich (HB) ansatz for the Green's function of the high-frequency Helmholtz equation in smooth inhomogeneous media. For a computational domain discretized with Nv discretization cells, the proposed algorithm first solves and tabulates the phase and HB coefficients via eikonal and transport equations with observation points and point sources located at the Chebyshev nodes using a set of much coarser computation grids, and then butterfly compresses the resulting HB interactions from all Nv cell centers to each other. The overall CPU time and memory requirement scale as O(Nv log2 Nv) for any bounded two-dimensional (2D) domains with arbitrary excitation sources. A direct extension of this scheme to bounded 3D domains yields an O(Nv4/3) CPU complexity, which can be further reduced to quasi-linear complexities with proposed remedies. The scheme can also efficiently handle scattering problems involving inclusions in inhomogeneous media. Although the current construction of our HB integrator does not accommodate caustics, the resulting HB integrator itself can be applied to certain sources, such as concave-shaped sources, to produce caustic effects. Compared to finite-difference frequency domain methods, the proposed HB integrator is free of numerical dispersion and requires fewer discretization points per wavelength. As a result, it can solve wave propagation problems well beyond the capability of existing solvers. Remarkably, the proposed scheme can accurately model wave propagation in 2D domains with 640 wavelengths per direction and in 3D domains with 54 wavelengths per direction on a state-of-the-art supercomputer at Lawrence Berkeley National Laboratory.

Hadamard--Babich ansatz↗

Techniques for estimating Space Station aerodynamic characteristics

A method was devised and calculations were performed to determine the effects of reflected molecules on the aerodynamic force and moment coefficients for a body in free molecule flow. A procedure was developed for determining the velocity and temperature distributions of molecules reflected from a surface of arbitrary momentum and energy accommodation. A system of equations, based on momentum and energy balances for the surface, incident, and reflected molecules, was solved by a numerical optimization technique. The minimization of a 'cost' function, developed from the set of equations, resulted in the determination of the defining properties of the flow reflected from the arbitrary surface. The properties used to define both the incident and reflected flows were: average temperature of the molecules in the flow, angle of the flow with respect to a vector normal to the surface, and the molecular speed ratio. The properties of the reflected flow were used to calculate the contribution of multiply reflected molecules to the force and moments on a test body in the flow. The test configuration consisted of two flat plates joined along one edge at a right angle to each other. When force and moment coefficients of this 90 deg concave wedge were compared to results that did not include multiple reflections, it was found that multiple reflections could nearly double lift and drag coefficients, with nearly a 50 percent increase in pitching moment for cases with specular or nearly specular accommodation. The cases of diffuse or nearly diffuse accommodation often had minor reductions in axial and normal forces when multiple reflections were included. There were several cases of intermediate accommodation where the addition of multiple reflection effects more than tripled the lift coefficient over the convex technique.

Thomas, Richard E.↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Noniterative estimation of a nonlinear parameter

An algorithm is described which solves the parameters X = (x1,x2,...,xm) and p in an approximation problem Ax nearly equal to y(p), where the parameter p occurs nonlinearly in y. Instead of linearization methods, which require an approximate value of p to be supplied as a priori information, and which may lead to the finding of local minima, the proposed algorithm finds the global minimum by permitting the use of series expansions of arbitrary order, exploiting an a priori knowledge that the addition of a particular function, corresponding to a new column in A, will not improve the goodness of the approximation.

Bergstroem, A.↗

Reflection and interference of electromagnetic waves in inhomogeneous media

Solutions were obtained of the wave equation for a plane horizontally polarized electro-magnetic wave incident on a semi infinite two dimensional inhomogeneous medium. Two problems were considered: An inhomogeneous half space, and an inhomogeneous layer of arbitrary thickness. Solutions of the wave equation were obtained in terms of Hankel functions with complex arguments. Numerical calculations were made of the reflection coefficient R at the interface of the homogeneous medium. The startling results show that the reflection coefficient for a complex dielectric constant with gradient, can be less than that of the same medium with zero gradient.

Geiger, F. E.↗

Radiation and scattering by thin-wire structures in the complex frequency domain

Piecewise-sinusoidal expansion functions and Galerkin's method are employed to formulate a solution for an arbitrary thin-wire configuration in a homogeneous conducting medium. The analysis is performed in the real or complex frequency domain. In antenna problems, the solution determines the current distribution, impedance, radiation efficiency, gain and far-field patterns. In scattering problems, the solution determines the absorption cross section, scattering cross section and the polarization scattering matrix. The electromagnetic theory is presented for thin wires and the forward-scattering theorem is developed for an arbitrary target in a homogeneous conducting medium.

Richmond, J. H.↗

Approximate solutions to the Navier-Stokes initial value problem

A Galerkin-Ritz procedure for any arbitrary system of field equations is shown to follow generically from 'two-functional' variational conditions. The initial-value problem for boundary-free incompressible Navier-Stokes flow is solved analytically in the two-parameter-function approximation.

Rosen, G.↗