Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary-value problems”

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 55 records · Page 3

Numerical determination of the fundamental eigenvalue for the Laplace operator on a spherical domain

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (i.e., the membrane eignevalue problem for the vibration equation) on the unit spherical surface are developed. Two types of spherical surface domains are considered: the interior of a spherical triangle, and the exterior of a great circle arc extending for less than pi radians (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is employed. The fundamental eigenvalue is approximated by iteration utilizing the power method and point successive overrelaxation. Some numerical results are given and compared, in certain special cases, with analytical solutions to the eigenvalue problem. The significance of the numerical eigenvalue results is discussed in terms of the singularities in the solution of three-dimensional boundary-value problems near a polyhedral corner of the domain.

Walden, H.↗

Calculation of free-fall trajectories using numerical optimization methods.

An important problem in space flight is the calculation of trajectories for nonthrusting vehicles between fixed points in a given time. A new procedure based on Hamilton's principle for solving such two-point boundary-value problems is presented. It employs numerical optimization methods to perform the extremization required by Hamilton's principle. This procedure is applied to the calculation of an Earth-Moon trajectory. The results show that the initial guesses required to obtain an iteration procedure which converges are not critical and that convergence can be obtained to any predetermined degree of accuracy.

Hull, D. G.↗

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interplanetary Trajectory Optimization with Powerlimited Propulsion Systems

A trajectory-optimization process is described in which the optimum­ thrust equations are derived using the calculus of variations. The mag­nitude of the thrust is constrained within an upper and a lower bound, but the thrust direction is arbitrary. This formulation allows both the constant-thrust program and the variable-thrust program to be con­sidered. For the constant-thrust program, certain propulsion-system parameters are optimized for maximum final vehicle mass. This theory has been used to study interplanetary missions to Venus and Mars using a power-limited propulsion system. Both one-way and round­ trip rendezvous trajectories are considered. The analysis employs a two-body inverse-square force-field model of three dimensions. An iterative routine used to solve the two-point boundary-value problem is described in the Appendix.

TRAJECTORY↗

Approximate estimation for systems with quantized data.

Estimation of the state of a nonlinear discrete-time system using quantized data is considered. An exact solution for the maximum likelihood estimate is expressed as the solution of a nonlinear two-point boundary-value problem. Approximate recursive solutions for both the maximum likelihood and the conditional-mean estimates are obtained. The results of Monte-Carlo simulations are presented in which the performance of these two algorithms is compared with that of a Kalman filter in which the quantization error is approximated by white noise.-

Clements, K. A.↗

Unsteady, combined radiation and conduction in an absorbing, scattering, and emitting medium.

The transient cooldown of a gray, absorbing, isotropic scattering, emitting, and conducting medium bounded by gray, diffusely emitting and reflecting parallel plates is considered. Numerical solutions are obtained for the initial boundary-value problem with a discontinuous decrease in temperature at one boundary. The quasi-steady equation of radiative transfer is solved using Gaussian quadrature and a matrix eigenvector technique together with explicit numerical solution of the unsteady energy equation. Temperature and energy flux distributions are presented for variations of optical thickness, boundary emissivity, albedo, and conduction-radiation parameter.

Weston, K. C.↗

NBODY - a multipurpose trajectory optimization computer program

Documentation of the NBODY trajectory optimization program is presented in the form of a mathematical development plus a user's manual. Optimal multistage-launch ascent trajectories may be determined by variational thrust steering during the upper phase. Optimal low-thrust interplanetary spacecraft trajectories may also be calculated with solar power or constant power, all-propulsion or embedded coast arcs, fixed or optimal thrust angles, and a variety of terminal end conditions. A hybrid iteration scheme solves the boundary-value problem, while either transversality conditions or a univariate search scheme optimize vehicle or trajectory parameters.

Strack, W. C.↗

Lagrangian formulation of the one-dimensional Vlasov equation

A new formulation of the one-dimensional Vlasov equation is derived which is analogous to the Kalman-transformed cold-plasma equations. The equations are shown to yield nonsecular, nonlinear approximations to a source or boundary-value problem. It is suggested that the formulation may have other applications in nonlinear plasma theory.

Lewak, G. J.↗

Nonlinear theory of electron neutralization waves in ions beams with dissipation

An analytical theory of nonlinear neutralization waves generated by injection of electrons from a grid in the direction of a homogeneous ion beam of uniform velocity and infinite extension is presented. The electrons are assumed to interact with the ions through the self-consistent space charge field and by strong collective interactions, while diffusion in the pressure gradient is disregarded (zero-temperature approximation). The associated nonlinear boundary-value problem is solved in closed form by means of a von Mises transformation. It is shown that the electron gas moves into the ion space in the form of a discontinuous neutralization wave, which exhibits a periodic field structure (incomplete neutralization). This periodic wave structure is damped out by intercomponent momentum transfer - i.e., after a few relaxation lengths a quasi-neutral plasma results.

Wilhelm, H. E.↗

Nonlinear transient neutralization theory of ion beams with dissipation

An analytical theory of nonlinear neutralization waves generated by injection of electrons from a grid in the direction of a homogeneous ion beam of uniform velocity and infinite extension is presented. The electrons are assumed to interact with the ions through the self-consistent space charge field and by strong collective interactions. The associated nonlinear boundary-value problem is solved in closed form by means of a von Mises transformation. It is shown that the electron gas moves into the ion space in the form of a discontinuous neutralization wave. This periodic wave structure is damped out by intercomponent momentum transfer, i.e., after a few relaxation lengths a quasi-neutral beam results. The relaxation scale in space agrees with neutralization experiments of rarefied ion beams, if the collective momentum transfer between the electron and ion streams is assumed to be of the Buneman type.

Wilhelm, H. E.↗

Analysis of multicircuit shells of revolution by the field method

The reported investigation represents a continuation of studies conducted by Cohen (1974). The current analysis treats a much broader class of connected graphs which may contain multiple circuits. Attention is given to the definition of a boundary-value problem, the field relations, a plan of the field method, differential equations for field functions, initial values of field functions, and shells of revolution.

Cohen, G. A.↗

Plasma rotation by electric and magnetic fields in a discharge cylinder

A theoretical model for an electric discharge consisting of a spatially diverging plasma sustained electrically between a small ring cathode and a larger ring anode in a cylindrical chamber with an axial magnetic field is developed to study the rotation of the discharge plasma in the crossed electric and magnetic fields. The associated boundary-value problem for the coupled partial differential equations which describe the electric potential and the plasma velocity fields is solved in closed form. The electric field, current density, and velocity distributions are discussed in terms of the Hartmann number and the Hall coefficient. As a result of Lorentz forces, the plasma rotates with speeds as high as 1 million cm/sec around its axis of symmetry at typical conditions. As an application, it is noted that rotating discharges of this type could be used to develop a high-density plasma-ultracentrifuge driven by j x B forces, in which the lighter (heavier) ion and atom components would be enriched in (off) the center of the discharge cylinder.

Wilhelm, H. E.↗

A higher order panel method for general analysis and design applications in subsonic flow

A higher-order panel method is described for numerical solution of boundary-value problems relating to steady inviscid irrotational incompressible subsonic fluid flow in a domain. Both Neumann and Dirichlet boundary conditions are treated; two types of auxiliary conditions are used to remove the degrees of freedom that arise from specifying only the derivative of the perturbation velocity potential. Four general network types and two expansions of the induced potential kernel are employed in the numerical solution. Some results are presented which illustrate the modeling options and numerical characteristics of the method.

Johnson, F. T.↗

Nonlinear stress waves in a perfectly flexible string

This paper discusses nonlinear stress-wave propagation in a perfectly flexible string obeying a quasilinear (rate-dependent) constitutive equation. Wave speeds and compatibility relations valid along various families of characteristics were determined. It was shown that the compatibility relations associated with the transverse as well as the longitudinal waves readily yield a physical interpretation when they are expressed in suitable variables and in vector form. Coding based on the present information was completed for the machine solution of a class of mixed initial- and boundary-value problems of practical interest. Computer simulation of the stress-wave interaction in the 40-foot lanyard in the Arcas 'Rocoz' system during deployment was carried out using a stress-strain relation for nylon at the strain rate of 30/second. A method for estimating the maximum tension and strain in a string during the initial loading phase is proposed.

Fan, D.-N.↗

Generalized Theodorsen solution for singular integral equations of the airfoil class

A class of singular integral equations is considered which arise in various two-dimensional mixed boundary-value problems with simple harmonic time variation. A problem typical of this class is that of determining the lifting pressure distribution on an oscillating airfoil in an unbounded incompressible potential flow. It is shown that Theodorsen's (1935) solution to this problem, with some modification, is valid for a general class of unsteady kernel functions. The technique employed is to consider an equivalent steady problem and then show that the unsteady resolvent and unsteady homogeneous solution can be written directly in terms of the steady solutions and a single frequency-dependent function which reduces to the Theodorsen function for the steady kernel.

Williams, M. H.↗

Flutter of a uniform wing with an arbitrarily placed mass according to a differential-equation analysis and a comparison with experiment

A method is presented for the calculation of the flutter speed of a uniform wing carrying an arbitrarily placed concentrated mass. The method, an extension of recently published work by Goland and Luke, involves the solution of the differential equations of motion of the wing at flutter speed and therefore does not require the assumption of specific normal modes of vibration. The order of the flutter determinant to be solved by this method depends upon the order of the system of differential equations and not upon the number of modes of vibration involved. The differential equations are solved by operational methods, and a brief discussion of operational methods as applied to boundary-value problems is included in one of two appendixes. A comparison is made with experiment for a wing with a large eccentrically mounted weight and good agreement is obtained. Sample calculations are presented to illustrate the method; and curves of amplitudes of displacement, torque, and shear for a particular case are compared with corresponding curves computed from the first uncoupled normal modes.

Runyan, Harry L↗

Theoretical study of the transonic lift of a double-wedge profile with detached bow wave

A theoretical study is described of the aerodynamic characteristics at small angle of attack of a thin, double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis is carried out within the framework of the transonic (nonlinear) small-disturbance theory, and the effects of angle of attack are regarded as a small perturbation on the flow previously calculated at zero angle. The mixed flow about the front half of the profile is calculated by relaxation solution of a suitably defined boundary-value problem for transonic small-disturbance equation in the hodograph plane (i.e., the Tricomi equation). The purely supersonic flow about the rear half is found by an extension of the usual numerical method of characteristics. Analytical results are also obtained, within the framework of the same theory, for the range of speed in which the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗