On integration of Hamilton-Jacobi partial differential equation.
Hamilton-Jacobi partial differential equation integration method, using variational principle and transversality condition
SEARCH · Engineering Papers
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.
Hamilton-Jacobi partial differential equation integration method, using variational principle and transversality condition
An analytical model is discussed that can be used to examine the effects of wind tunnel walls on helicopter rotational noise. A complete physical model of an acoustic source in a wind tunnel is described and a simplified version is then developed. This simplified model retains the important physical processes involved, yet it is more amenable to analysis. The simplified physical model is then modeled as a mathematical problem. An inhomogeneous partial differential equation with mixed boundary conditions is set up and then transformed into an integral equation. Details of generating a suitable Green's function and integral equation are included and the equation is discussed and also given for a two-dimensional case.
A previously developed local-global stiffness matrix methodology for the response of a composite half plane, arbitrarily layered with isotropic, orthotropic or monoclinic plies, to indentation by a rigid parabolic punch is further extended to accommodate the presence of layers with complex eigenvalues (e.g., honeycomb or piezoelectric layers). First, a generalized plane deformation solution for the displacement field in an orthotropic layer or half plane characterized by complex eigenvalues is obtained using Fourier transforms. A local stiffness matrix in the transform domain is subsequently constructed for this class of layers and half planes, which is then assembled into a global stiffness matrix for the entire multilayered half plane by enforcing continuity conditions along the interfaces. Application of the mixed boundary condition on the top surface of the half plane indented by a rigid punch results in an integral equation for the unknown pressure in the contact region. The integral possesses a divergent kernel which is decomposed into Cauchy-type and regular parts using the asymptotic properties of the local stiffness matrix and a relationship between Fourier and finite Hilbert transform of the contact pressure. The solution of the resulting singular integral equation is obtained using a collocation technique based on the properties of orthogonal polynomials developed by Erdogan and Gupta. Examples are presented that illustrate the important influence of low transverse properties of layers with complex eigenvalues, such as those exhibited by honeycomb, on the load versus contact length response and contact pressure distributions for half planes containing typical composite materials.
Class of functional integral equations in space of continuous functions
Relaxation schemes based on approximate and incomplete factorization technique (AF) are described. The AF schemes allow construction of a fast multigrid method for solving integral equations of the second and first kind. The smoothing factors for integral equations of the first kind, and comparison with similar results from the second kind of equations are a novel item. Application of the MD algorithm shows convergence to the level of truncation error of a second order accurate panel method.
A solution to the problem of plane wave scattering by a rectangular microstrip patch on a grounded dielectric substrate is presented. The model does not include the microstrip feed, and thus doe not include the so-called 'antenna mode' component of the scattering. The solution begins by formulating an electric field integral equation for the surface current density on the microstrip patch. The integral equation is solved using the method of moments. Computed data for the patch radar cross section (RCS) is found to be in close agreement with measurements over a broad frequency range. The microstrip RCS versus frequency consists of a number of large peaks which are identified as impedance or pattern factor resonance peaks.
The third phase of the development of the computer codes for scattering by coated bodies that has been part of an ongoing effort in the Electromagnetics Laboratory of the Electrical Engineering and Computer Science Department at the University of Illinois at Chicago is described. The work reported discusses the analytical and numerical results for the scattering of an obliquely incident plane wave by impedance bodies of revolution with phi variation of the surface impedance. Integral equation formulation of the problem is considered. All three types of integral equations, electric field, magnetic field, and combined field, are considered. These equations are solved numerically via the method of moments with parametric elements. Both TE and TM polarization of the incident plane wave are considered. The surface impedance is allowed to vary along both the profile of the scatterer and in the phi direction. Computer code developed for this purpose determines the electric surface current as well as the bistatic radar cross section. The results obtained with this code were validated by comparing the results with available results for specific scatterers such as the perfectly conducting sphere. Results for the cone-sphere and cone-cylinder-sphere for the case of an axially incident plane were validated by comparing the results with the results with those obtained in the first phase of this project. Results for body of revolution scatterers with an abrupt change in the surface impedance along both the profile of the scatterer and the phi direction are presented.
The scattering effect of harmonic flexural waves at a through crack in an elastic plate carrying electrical current is investigated. In this context, the Kirchhoffean bending plate theory is extended as to include magnetoelastic interactions. An incident wave giving rise to bending moments symmetric about the longitudinal z-axis of the crack is applied. Fourier transform technique reduces the problem to dual integral equations, which are then cast to a system of two singular integral equations. Efficient numerical computation is implemented to get the bending moment intensity factor for arbitrary frequency of the incident wave and of arbitrary electrical current intensity. The asymptotic behaviour of the bending moment intensity factor is analysed and parametric studies are conducted.
The results are reported of two unrelated studies. The first was an investigation of the formulation of the equations for non-uniform unsteady flows, by perturbation of an irrotational flow to obtain the linear Green's equation. The resulting integral equation was found to contain a kernel which could be expressed as the solution of the adjoint flow equation, a linear equation for small perturbations, but with non-constant coefficients determined by the steady flow conditions. It is believed that the non-uniform flow effects may prove important in transonic flutter, and that in such cases, the use of doublet type solutions of the wave equation would then prove to be erroneous. The second task covered an initial investigation into the use of the Monte Carlo method for solution of acoustical field problems. Computed results are given for a rectangular room problem, and for a problem involving a circular duct with a source located at the closed end.
The general relativistic correction to the position of a satellite is found by retaining Newtonian physics for an observer on the satellite and introducing a potential. The potential is expanded in terms of the Keplerian elements of the orbit and substituted in Lagrange's equations. Integration of the equations shows that a typical earth satellite with small orbital eccentricity is displaced by about 17 cm. from its unperturbed position after a single orbit, while the periodic displacement over the orbit reaches a maximum of about 3 cm. The moon is displaced by about the same amounts. Application of the equations to Mercury gives a total displacement of about 58 km. after one orbit and a maximum periodic displacement of about 12 km.
An orthotropic laminate composite containing a completely broken layer is considered. The problem is formulated in terms of integral transforms and then reduced to a singular integral equation which is solved numerically. The strength of stress singularity at the crack tip is determined from a characteristic equation which is obtained by studying the dominant part of the singular integral equation near the end points. The stress intensity factors are given for various material properties.
An orthotropic laminate composite containing a completely broken layer is considered. The problem is formulated in terms of integral transforms and then reduced to a singular integral equation which is solved numerically. The strength of stress singularity at the crack tip is determined from a characteristic equation which is obtained by studying the dominant part of the singular integral equation near the end points. The stress intensity factors are given for various material properties.
A new approach to the solution of matrix equations resulting from integral equations is presented and applied to the solution of two-dimensional Neumann problems describing the inviscid, incompressible flow past an airfoil. The problem is reformulated in terms of a preselected set of mode functions giving an equivalent matrix equation to be solved for the mode-function expansion coefficients. Because of the inherent smoothness of the original problem, the coefficient problem can be solved approximately without significantly affecting the accuracy of the final solution. Very promising two-dimensional results are obtained and the extension of the method to three-dimensional problems is investigated. On the basis of these results it is shown that the computing time for the matrix solution for a large three-dimensional panel method calculation could be reduced by an order of magnitude compared with that required for a direct solution.
The plane elasticity problem for a nonhomogeneous medium containing a crack is considered. It is assumed that the Poisson's ratio of the medium is constant and the Young's modulus E varies exponentially with the coordinate parallel to the crack. First the half plane problem is formulated and the solution is given for arbitrary tractions along the boundary. Then the integral equation for the crack problem is derived. It is shown that the integral equation having the derivative of the crack surface displacement as the density function has a simple Cauchy type kernel. Hence, its solution and the stresses around the crack tips have the conventional square root singularity. The solution is given for various loading conditions. The results show that the effect of the Poisson's ratio and consequently that of the thickness constraint on the stress intensity factors are rather negligible.
A half space of combustible material is subjected to an arbitrary energy flux at the boundary where convection heat loss is also allowed. An asymptotic analysis of the temperature growth reveals two conditions necessary for ignition to occur. Cases of both large and order unity Lewis number are shown to lead to a nonlinear integral equation governing the thermal runaway. Some global and asymptotic properties of the integral equation are obtained.
Indirect methods are employed to estimate the wintertime and summertime mean vertical velocity fields of the extratropical Northern Hemisphere and intercomparisons are made, together with comparisons with mean seasonal patterns of cloudiness and precipitation. Twice-daily NMC operational analyses produced general circulation statistics for 11 winters and 12 summers, permitting calculation of the seasonal NMC averages for 6 hr forecasts, solution of the omega equation, integration of continuity equation downward from 100 mb, and solution of the thermodynamic energy equation in the absence of diabatic heating. The methods all yielded similar vertical velocity patterns; however, the magnitude of the vertical velocities could not be calculated with great accuracy. Orography was concluded to have less of an effect in summer than in winter, when winds are stronger.
The plane elasticity problem for a nonhomogeneous medium containing a crack is considered. It is assumed that the Poisson's ratio of the medium is constant and the Young's modulus E varies exponentially with the coordinate parallel to the crack. First the half plane problem is formulated and the solution is given for arbitrary tractions along the boundary. Then the integral equation for the crack problem is derived. It is shown that the integral equation having the derivative of the crack surface displacement as the density function has a simple Cauchy type kernel. Hence, its solution and the stresses around the crack tips have the conventional square root singularity. The solution is given for various loading conditions. The results show that the effect of the Poisson's ratio and consequently that of the thickness constraint on the stress intensity factors are rather negligible.
A numerical procedure is presented for computing the unsteady transonic flow field about three dimensional swept wings undergoing general time dependent motion. The outer inviscid portion of the flow is assumed to be governed by the modified unsteady transonic small disturbance potential equation which is integrated in the time domain by means of an efficient alternating direction implicit approximate factorization algorithm. Gross dominant effects of the shock boundary layer interaction are accounted for by a simple empirically defined model. Viscous flow regions adjacent to the wing surface and in the trailing wake are described by a set of integral equations appropriate for compressible turbulent shear layers. The two dimensional boundary layer equations are applied quasi-statically stripwise across the span. Coupling with the outer inviscid flow is implemented through use of the displacement thickness concept within the limitations of small disturbance theory. Validity of the assumptions underlying the method is established by comparison with experimental data for the flow about a high aspect ratio transport wing having an advanced airfoil section.