Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAPLACE EQUATION”

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 181 records · Page 10

Heat transfer characteristics of an emergent strand

A mathematical model was developed to describe the heat transfer characteristics of a hot strand emerging into a surrounding coolant. A stable strand of constant efflux velocity is analyzed, with a constant (average) heat transfer coefficient on the sides and leading surface of the strand. After developing a suitable governing equation to provide an adequate description of the physical system, the dimensionless governing equation is solved with Laplace transform methods. The solution yields the temperature within the strand as a function of axial distance and time. Generalized results for a wide range of parameters are presented, and the relationship of the results and experimental observations is discussed.

Simon, W. E.↗

Coupling finite and boundary element methods for two-dimensional potential problems

A finite element-boundary element (FE-BE) coupling method based on a weighted residual variational method is presented for potential problems, governed by either the Laplace or the Poisson equations. In this method, a portion of the domain of interest is modeled by finite elements (FE) and the remainder of the region by boundary elements (BE). Because the BE fundamental solutions are valid for infinite domains, a procedure that limits the effects of the BE fundamental solution to a small region adjacent to the FE region, called the transition region (TR), is developed. This procedure involves a judicious choice of functions called the transition (T) functions that have unit values on the BE-TR interface and zero values on the FE-TR interface. The present FE-BE coupling algorithm is shown to be independent of the extent of the transition region and the choice of the transition functions. Therefore, transition regions that extend to only one layer of elements between FE and BE regions and the use of simple linear transition functions work well.

Krishnamurthy, T.↗

Feeling Gravity's Pull: Gravity Modeling. The Gravity Field of Mars

Most people take the constant presence of gravitys pull for granted. However, the Earth's gravitational strength actually varies from location to location. This variation occurs because mass, which influences an object's gravitational pull, is not evenly distributed within the planet. Changes in topography, such as glacial movement, an earthquake, or a rise in the ocean level, can subtly affect the gravity field. An accurate measurement of the Earth's gravity field helps us understand the distribution of mass beneath the surface. This insight can assist us in locating petroleum, mineral deposits, ground water, and other valuable substances. Gravity mapping can also help notice or verify changes in sea surface height and other ocean characteristics. Such changes may indicate climate change from polar ice melting and other phenomena. In addition, gravity mapping can indicate how land moves under the surface after earthquakes and other plate tectonic processes. Finally, changes in the Earth's gravity field might indicate a shift in water distribution that could affect agriculture, water supplies for population centers, and long-term weather prediction. Scientists can map out the Earth's gravity field by watching satellite orbits. When a satellite shifts in vertical position, it might be passing over an area where gravity changes in strength. Gravity is only one factor that may shape a satellite's orbital path. To derive a gravity measurement from satellite movement, scientists must remove other factors that might affect a satellite's position: 1. Drag from atmospheric friction. 2. Pressure from solar radiation as it heads toward Earth and. as it is reflected off the surface of the Earth 3. Gravitational pull from the Sun, the Moon, and other planets in the Solar System. 4. The effect of tides. 5. Relativistic effects. Scientists must also correct for the satellite tracking process. For example, the tracking signal must be corrected for refraction through the atmosphere of the Earth. Supercomputers can calculate the effect of gravity for specific locations in space following a mathematical process known as spherical harmonics, which quantifies the gravity field of a planetary body. The process is based on Laplace's fundamental differential equation of gravity. The accuracy of a spherical harmonic solution is rated by its degree and order. Minute variations in gravity are measured against the geoid, a surface of constant gravity acceleration at mean sea level. The geoid reference gravity model strength includes the central body gravitational attraction (9.8 m/sq s) and a geopotential variation in latitude partially caused by the rotation of the Earth. The rotational effect modifies the shape of the geoid to be more like an ellipsoid, rather than a perfect, circle. Variations of gravity strength from the ellipsoidal reference model are measured in units called milli-Galileos (mGals). One mGal equals 10(exp -5) m/sq s. Research projects have also measured the gravity fields of other planetary bodies, as noted in the user profile that follows. From this information, we may make inferences about our own planet's internal structure and evolution. Moreover, mapping the gravity fields of other planets can help scientists plot the most fuel-efficient course for spacecraft expeditions to those planets.

Lemoine, Frank↗

Tidal Friction in the Earth-Moon System and Laplace Planes: Darwin Redux

The dynamical evolution of the Earth-Moon system due to tidal friction is treated here. George H. Darwin used Laplace planes (also called proper planes) in his study of tidal evolution. The Laplace plane approach is adapted here to the formalisms of W.M. Kaula and P. Goldreich. Like Darwin, the approach assumes a three-body problem: Earth, Moon, and Sun, where the Moon and Sun are point-masses. The tidal potential is written in terms of the Laplace plane angles. The resulting secular equations of motion can be easily integrated numerically assuming the Moon is in a circular orbit about the Earth and the Earth is in a circular orbit about the Sun. For Earth-Moon distances greater than ∼10 Earth radii, the Earth's approximate tidal response can be characterized with a single parameter, which is a ratio: a Love number times the sine of a lag angle divided by another such product. For low parameter values it can be shown that Darwin's low-viscosity molten Earth, M. Ross's and G. Schubert's model of an Earth near melting, and Goldreich's equal tidal lag angles must all give similar histories. For higher parameter values, as perhaps has been the case at times with the ocean tides, the Earth's obliquity may have decreased slightly instead of increased once the Moon's orbit evolved further than 50 Earth radii from the Earth, with possible implications for climate. This is contrast to the other tidal friction models mentioned, which have the obliquity always increasing with time. As for the Moon, its orbit is presently tilted to its Laplace plane by 5.2deg. The equations do not allow the Moon to evolve out of its Laplace plane by tidal friction alone, so that if it was originally in its Laplace plane, the tilt arose with the addition of other mechanisms, such as resonance passages.

Darwin Redux↗

ISAC - A tool for aeroservoelastic modeling and analysis

This paper discusses the capabilities of the Interaction of Structures, Aerodynamics, and Controls (ISAC) system of program modules. The major modeling, analysis, and data management components of ISAC are identified. Equations of motion are displayed for a Laplace-domain representation of the unsteady aerodynamic forces. Options for approximating a frequency-domain representation of unsteady aerodynamic forces with rational functions of the Laplace variable are shown. Linear time invariant state-space equations of motion that result are discussed. Model generation and analyses of stability and dynamic response characteristics are shown for an aeroelastic vehicle which illustrate some of the capabilities of ISAC as a modeling and analysis tool for aeroelastic applications.

Adams, William M., Jr.↗

ISAC: A tool for aeroservoelastic modeling and analysis

The capabilities of the Interaction of Structures, Aerodynamics, and Controls (ISAC) system of program modules is discussed. The major modeling, analysis, and data management components of ISAC are identified. Equations of motion are displayed for a Laplace-domain representation of the unsteady aerodynamic forces. Options for approximating a frequency-domain representation of unsteady aerodynamic forces with rational functions of the Laplace variable are shown. Linear time invariant state-space equations of motion that result are discussed. Model generation and analyses of stability and dynamic response characteristics are shown for an aeroelastic vehicle which illustrates some of the capabilities of ISAC as a modeling and analysis tool for aeroelastic applications.

Adams, William M., Jr.↗

Transfer Functions Via Laplace- And Fourier-Borel Transforms

Approach to solution of nonlinear ordinary differential equations involves transfer functions based on recently-introduced Laplace-Borel and Fourier-Borel transforms. Main theorem gives transform of response of nonlinear system as Cauchy product of transfer function and transform of input function of system, together with memory effects. Used to determine responses of electrical circuits containing variable inductances or resistances. Also possibility of doing all noncommutative algebra on computers in such symbolic programming languages as Macsyma, Reduce, PL1, or Lisp. Process of solution organized and possibly simplified by algebraic manipulations reducing integrals in solutions to known or tabulated forms.

Can, Sumer↗

Renormalization and scale evolution of the soft-quark soft function

Soft functions defined in terms of matrix elements of soft fields dressed by Wilson lines are central components of factorization theorems for cross sections and decay rates in collider and heavy-quark physics. While in many cases the relevant soft functions are defined in terms of gluon operators, at subleading order in power counting soft functions containing quark fields appear. We present a detailed discussion of the properties of the soft-quark soft function consisting of a quark propagator dressed by two finite-length Wilson lines connecting at one point. This function enters in the factorization theorem for the Higgs-boson decay amplitude of the h → γγ process mediated by light-quark loops. We perform the renormalization of this soft function at one-loop order, present a conjecture for its two-loop anomalous dimension and discuss solutions to its renormalization-group evolution equation in momentum space, in Laplace space and in the “diagonal space”, where the evolution is strictly local in the momentum variable.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electromagnetic plasma wave propagation along a magnetic field

The linearized response of a Vlasov plasma to the steady-state excitation of transverse plasma waves along an external magnetic field is examined. Assuming a delta-function excitation mechanism, and performing a detailed Vlasov-Maxwell equation analysis using Fourier-Laplace transforms, the plasma response is found to consist of three terms: a branch-cut term, a free-streaming term, and a dielectric-pole term. Also considered is the phenomenon of plasma wave echoes. The case of longitudinal electrostatic waves is extended to the case of transverse plasma waves that propagate along an external magnetic field. It is shown that a transverse echo results in lowest order only when one excitation is transverse and the other is longitudinal.

Olson, C. L.↗

Design for active flutter suppression and gust alleviation using state-space aeroelastic modeling

An analytical design technique for an active flutter-suppression and gust-alleviation control system is presented. It is based on a rational approximation of the unsteady aerodynamic loads in the entire Laplace domain, which yields matrix equations of motion with constant coefficients. Some existing rational approximation schemes are reviewed, and a new technique which yields a minimal number of augmented states for a desired accuracy is presented. The state-space aeroelastic model is used to design a constant gain, partial-feedback control system, which simultaneously assures stability and optimizes any desired combination of gust response parameters throughout the entire flight envelope.

Karpel, M.↗

Design for active and passive flutter suppression and gust alleviation

Analytical design techniques for active and passive control of aeroelastic systems are based on a rational approximation of the unsteady aerodynamic loads in the entire Laplace domain, which yields matrix equations of motion with constant coefficients. Some existing schemes are reviewed, the matrix Pade approximant is modified, and a technique which yields a minimal number of augmented states for a desired accuracy is presented. The state-space aeroelastic model is used to design an active control system for simultaneous flutter suppression and gust alleviation. The design target is for a continuous controller which transfers some measurements taken on the vehicle to a control command applied to a control surface. Structural modifications are formulated in a way which enables the treatment of passive flutter suppression system with the same procedures by which active control systems are designed.

Karpel, M.↗

Time-dependent shock acceleration of energetic electrons including synchrotron losses

The present investigation of the time-dependent particle acceleration problem in strong shocks, including synchrotron radiation losses, solves the transport equation analytically by means of Laplace transforms. The particle distribution thus obtained is then transformed numerically into real space for the cases of continuous and impulsive injections of particles at the shock. While in the continuous case the steady-state spectrum undergoes evolution, impulsive injection is noted to yield such unpredicted features as a pile-up of high-energy particles or a steep power-law with time-dependent spectral index. The time-dependent calculations reveal varying spectral shapes and more complex features for the higher energies which may be useful in the interpretation of outburst spectra.

Fritz, Klaus-Dieter↗

Drag Minimization for Wings and Bodies in Supersonic Flow

The minimization of inviscid fluid drag is studied for aerodynamic shapes satisfying the conditions of linearized theory, and subject to imposed constraints on lift, pitching moment, base area, or volume. The problem is transformed to one of determining two-dimensional potential flows satisfying either Laplace's or Poisson's equations with boundary values fixed by the imposed conditions. A general method for determining integral relations between perturbation velocity components is developed. This analysis is not restricted in application to optimum cases; it may be used for any supersonic wing problem.

Heaslet, Max A↗

Growth Angle: A Microscopic View

A microscopic continuum mechanical model of the growth angle is proposed. It is based on the van der Waals type framework that is used for surface force phenomena. The obtained augmented Laplace type integro-differential equations are, in general, difficult to analyze. Here we focused primarily on the particular case of equal melt and crystal surface energies. We derived an approximate equation for the meniscus shape, and obtained an analytical relationship between the contact and the growth angle. Interestingly, the same result can be obtained using the macroscopic model of Herring. The case of a macroscopically sharp corner is also considered. For this case, the macroscopic angle is not defined and it can be any angle between the contact angles of both flat surfaces. The microscopic model yields the smooth shape for the meniscus that also is not unique, but depends on the initial position of the meniscus.

Mazuruk, Konstantin↗

Thermodynamic analysis of journal bearings operating under steady state loading in laminar regime

Thermohydrodynamic effects in journal bearings operating under steady load in laminar regime are investigated. An analytical model for the finite and infinitely long journal bearings is formulated. The model includes correction factors for the cavitation effects in the unloaded region of the bearing and the mixing of the recirculating oil and supply oil at the oil inlet. A finite difference computer program is developed to numerically solve the governing equations of the continuity, Reynolds, energy, Laplace heat conduction, and a viscosity-temperature relation simultaneously. The program includes a numerical technique for obtaining an isothermal shaft temperature. The numerical results of temperature distribution and the heat effects on the bearing load carrying capacity agree closely with those of experimental findings. Several different sets of simpler boundary conditions for the energy equation are studied.

Khonsari, M. M.↗

Off-axis impact of unidirectional composites with cracks: Dynamic stress intensification

The dynamic response of unidirectional composites under off axis (angle loading) impact is analyzed by assuming that the composite contains an initial flaw in the matrix material. The analytical method utilizes Fourier transform for the space variable and Laplace transform for the time variable. The off axis impact is separated into two parts, one being symmetric and the other skew-symmetric with reference to the crack plane. Transient boundary conditions of normal and shear tractions are applied to a crack embedded in the matrix of the unidirectional composite. The two boundary conditions are solved independently and the results superimposed. Mathematically, these conditions reduce the problem to a system of dual integral equations which are solved in the Laplace transform plane for the transformation of the dynamic stress intensity factor. The time inversion is carried out numerically for various combinations of the material properties of the composite and the results are displayed graphically.

Sih, G. C.↗