A computational method for two-impulse orbital rendezvous and transfer problems
Computational method for two impulse orbital transfer and rendezvous between two elliptical orbits
RENDEZVOUS↗
SEARCH · Engineering Papers
Results for “Elliptic 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.
A semi-direct procedure using a local relaxation factor and its application to an internal flow problem
Generally, fast direct solvers are not directly applicable to a nonseparable elliptic partial differential equation. This limitation, however, is circumvented by a semi-direct procedure, i.e., an iterative procedure using fast direct solvers. An efficient semi-direct procedure which is easy to implement and applicable to a variety of boundary conditions is presented. The current procedure also possesses other highly desirable properties, i.e.: (1) the convergence rate does not decrease with an increase of grid cell aspect ratio, and (2) the convergence rate is estimated using the coefficients of the partial differential equation being solved.
Chang, S. C.↗
Key design issues for user terminals operating with HEO satellites
Coverage of high-latitude areas in the northern hemisphere using geostationary satellites can not be achieved with good elevation angles, and highly inclined elliptical orbits (HEO) present one possible way to overcome this problem. The high elevation angles they achieve and the relative motion of the satellites across the sky, including satellite hand-over, present a specific set of requirements for the design of user terminals, be they for communication or broadcast-reception purposed. This paper reviews key design issues for user terminals operating with satellites in highly inclined elliptical orbits. While all aspects of user terminal design are broadly covered, special emphasis is given to the antenna and RF front-end design. The types of terminals discussed include communications terminals: the hand-held, the notebook, and the vehicle-mounted types, as well as the vehicle-mounted digital audio broadcasting (DAB) receiver terminal. A number of solutions are outlined, based on tradeoffs with respect to the most significant system interface parameters such as G/T and equivalent isotropically radiated power (EIRP). Candidate antenna and front-end configurations were studied, in particular patch and helical antennas with the principal aim of deriving low cost, yet effective solutions. Issues related to achieving best EIRP for the communications terminals are presented, taking into account radiation safety limits. The differences in the design of receive only versus receive/transmit RF front-ends for HEO system terminals are highlighted and solutions for both proposed. This paper summarizes the results of internal ESTEC studies and external ESA-sponsored work concerning major user terminal parameters and design aspects performed for Archimedes -- a HEO satellite constellation for multimedia broadcasting and communications.
Stojkovic, I.↗
A fast semi-implicit algorithm for problems of mixed type
Certain physical processes are modeled by partial differential equations which are parabolic over part of the domain and elliptic over the remainder. A family of semi-implicit algorithms which are well suited to initial-boundary value problems of this mixed type is discussed. One important feature of these algorithms is the use of an approximate inverse for the solution of the implicit linear system. A strong error analysis results in an estimate of the total error as a function of approximate inverse error e and time step h.
Frederickson, P. O.↗
Sum-Frequency Generation from a Thin Cylindrical Layer
In the Rayleigh–Gans–Debye approximation, we have solved the problem of the sum-frequency generation by two plane elliptically polarized electromagnetic waves from the surface of a dielectric particle of a cylindrical shape that is coated by a thin layer possessing nonlinear optical properties. The formulas that describe the sum-frequency field have been presented in the tensor and vector forms for the second-order nonlinear dielectric susceptibility tensor, which was chosen in the general form, containing chiral components. Expressions describing the sum-frequency field from the cylindrical particle ends have been obtained for the case of a nonlinear layer possessing chiral properties. Three-dimensional directivity patterns of the sum-frequency radiation have been analyzed for different combinations of parameters (angles of incidence, degrees of ellipticity, orientations of polarization ellipses, cylindrical particle dimensions). The mathematical properties of the spatial distribution functions of the sum-frequency field, which characterize the symmetry of directivity patterns, have been revealed.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Benchmark Problems for Space Mission Formation Flying
To provide a high-level focus to distributed space system flight dynamics and control research, several benchmark problems are suggested for space mission formation flying. The problems cover formation flying in low altitude, near-circular Earth orbit, high altitude, highly elliptical Earth orbits, and large amplitude lissajous trajectories about co-linear libration points of the Sun-Earth/Moon system. These problems are not specific to any current or proposed mission, but instead are intended to capture high-level features that would be generic to many similar missions that are of interest to various agencies.
Carpenter, J. Russell↗
On the Orbital Elements of the Two-body Problem with Slowly Decreasing Mass: The Gyldén–Mestchersky Cases
We study the dynamics of two-body problems with gravitational parameters μ = μ(t) decreasing with the time. In particular we focus our attention on the behavior of the orbital elements in the first two Gyldén–Mestchersky cases: in the first case (GM-I) μ(t)=μ{sub 0}(1+αt){sup −1}, where μ{sub 0} is the constant initial value of the parameter and α a small positive constant. In the second case (GM-II) μ(t)=μ{sub 0}(1+2αt){sup −1/2}. In the GM-I problem, using a suitable transformation of the radius vector and the physical time, which reduces the original problem to a pure Kepler problem, we are able to obtain explicit analytical expressions of the orbital elements that hold true for all t⩾0 and therefore derive their asymptotic behavior. Moreover using the alternative scaled time these expressions show clearly their periodic and secular terms. For the GM-II problem a similar transformation to the GM-I problem leads to a perturbed Kepler problem with a radial direction perturbation. Although exact solutions in terms of elliptic functions may be obtained, the use of some integrals leads to analytical expressions of orbital elements that allow us to study their behavior. Finally, some comparisons between the time evolution of elements in both problems are also included.
79 ASTRONOMY AND ASTROPHYSICS↗
On Integral of a Semi-Markov Diffusion Process
Let (X(t)) (t ≥ 0) be a semi-Markov diffusion process. The process (J(T )) (T ≥ 0) equal to the integral of (X(t)) on interval [0, T ) is studied. The relation between one-dimensional differential equation of the second order of elliptical type and asymptotics of a solution to Dirichlet problem on an interval with length tending to zero is established. This relation is used to derive a differential equation for the Laplace transform of the semi-Markov generating function of the process (J(t)).
97 MATHEMATICS AND COMPUTING↗
Three-dimensional viscous flows with large secondary velocity
A new system of approximation equations is derived for three-dimensional steady viscous compressible flows in which a primary-flow direction is present, but in which both transverse velocity components can be large. Previous approaches which address simplification of the steady Navier-Stokes equations are discussed, and a new approach is proposed. The transverse velocity vector which corrects a given potential flow has been decomposed into potential and rotational components. It is found that the potential-velocity vector may be assumed small, whereas the rotational-velocity vector may be assumed small, whereas the rotational velocity vector and hence the composite secondary flow can be of order unity. This assumption leads to a system of governing equations whose characteristic polynomial has a non-elliptic form for arbitrary Mach numbers. The resulting non-elliptic approximation equations can be solved as an initial/boundary-value problem. Computed results confirm the small scalar-potential approximation.
Briley, W. R.↗
Best estimate of luminal cross-sectional area of coronary arteries from angiograms
We have reexamined the problem of estimating the luminal area of an elliptically-shaped coronary artery cross section from two or more radiographic diameter measurements. The expected error is found to be much smaller than the maximum potential error. In the cae of two orthogonal views, closed form expressions have been derived for calculating the area and the uncertainty. Assuming that the underlying ellipse has limited ellipticity (major/minor axis ratio less than five), it is shown that the average uncertainty in the area is less than 14 percent. When more than two views are available, we suggest using a least-squares fit method to extract all available information from the data.
Lee, P. L.↗
Algorithms for the automatic generation of 2-D structured multi-block grids
Two different approaches to the fully automatic generation of structured multi-block grids in two dimensions are presented. The work aims to simplify the user interactivity necessary for the definition of a multiple block grid topology. The first approach is based on an advancing front method commonly used for the generation of unstructured grids. The original algorithm has been modified toward the generation of large quadrilateral elements. The second method is based on the divide-and-conquer paradigm with the global domain recursively partitioned into sub-domains. For either method each of the resulting blocks is then meshed using transfinite interpolation and elliptic smoothing. The applicability of these methods to practical problems is demonstrated for typical geometries of fluid dynamics.
Schoenfeld, Thilo↗
Building a better leapfrog
In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.
Hut, Piet↗
Optimizing Power Generation of a Bottom-Raised Oscillating Surge Wave Energy Converter Using a Theoretical Model: Preprint
Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.
ACE↗
Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model
Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.
ACE↗
The Restricted Three-Body-Problem as a Perturbation of Euler's Problem of Two Fixed Centers and Its Application to Lunar Trajectories
The restricted Three-Body-Problem considers the motion of an infinitesimal mass under the gravitational attraction of two finite masses, which revolve about their common center of gravity in coplanar circles. It is well known that Euler's problem of two fixed centers, consisting of the motion of an infinitesimal mass under the gravitational attraction of two finite masses fixed in space, can be solved by elliptic functions. The idea presented here is to take the solution of Euler's problem as the solution of the restricted Three-Body-Problem by allowing the initial values to be functions of time now. Differential equations for the perturbed initial values are established. These equations can be given in closed form by using the fact that the transformation to the perturbed initial values of Euler's problem is canonical. Thus, an approximation can be obtained for the solution of the restricted Three-Body-Problem. The method can also be used to represent classes of neighboring trajectories for guidance purposes.
Euler equation↗
Systems of Functions Orthogonal Over the Domain and Their Application in Boundary-Value Problems of Mathematical Physics
We formulate a boundary-value problem for eigenvalues and eigenfunctions of the Helmholtz equation in a complex domain with the use of mutually conjugated complex variables. The obtained systems of functions are orthogonal in this domain and constructed by using the Bessel functions and powers of the conformal mappings of the analyzed domains onto a circle. The solutions of boundary-value problems for the principal equations of mathematical physics (hyperbolic, parabolic and elliptic types) are obtained in the form of the sums of series in the systems of functions orthogonal over the domain.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A regularization method for extrapolation of solar potential magnetic fields
The mathematical basis of a Tikhonov regularization method for extrapolating the chromospheric-coronal magnetic field using photospheric vector magnetograms is discussed. The basic techniques show that the Cauchy initial value problem can be formulated for potential magnetic fields. The potential field analysis considers a set of linear, elliptic partial differential equations. It is found that, by introducing an appropriate smoothing of the initial data of the Cauchy potential problem, an approximate Fourier integral solution is found, and an upper bound to the error in the solution is derived. This specific regularization technique, which is a function of magnetograph measurement sensitivities, provides a method to extrapolate the potential magnetic field above an active region into the chromosphere and low corona.
Gary, G. A.↗
Tearing instability in a tokamak with a noncircular cross section
Using a straight-column model to describe tokamak plasma with a noncircular cross section, it is shown how to (i) find the boundary of tearing instability from the condition of existence of a magnetohydrodynamic plasma equilibrium different from that of a straight cylinder by solving a two-dimensional linear boundary-value problem with a second-order equation with respect to the flux coordinate and (ii) find the spatial structure of the tearing mode and the corresponding effective Δ' when there is only one resonance magnetic surface in the plasma for a given axial wavenumber by solving some kind of a boundary-value problem for the perturbation. The proposed approach is illustrated by numerical calculations for the case of an elliptical cross section as an example.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗