Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELLIPTIC INTEGRAL”

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 127 records · Page 7

Faster Algorithm For Computation Of Incompressible Flow

Improved algorithm yields faster numerical solutions of Navier-Stokes equations of steady or unsteady three-dimensional flow of incompressible fluid. In artificial-compressibility method, unsteady flow treated as incompressible in advancing from one time step to next, but at each time step (or in steady state), fluid treated as having variable compressibility enabling propagation of flow field, and subiterations performed in increments of pseudotime until effects of compressibility subside. Directly couples pressure and velocity fields at same time step and converts elliptic incompressible Navier-Stokes equations to hyperbolic form more amenable to numerical integration.

Rogers, S. E.↗

Homotopy Solutions of Kepler's Equations

Kepler's Equation is solved using an integrative algorithm developed using homotropy theory. The solution approach is applicable to both elliptic and hyperbolic forms of Kepler's Equation. The results from the proposed algorithm compare quite favorably with those from existing iterative schemes.

Fitz-Coy, Norman↗

A global multilevel atmospheric model using a vector semi-Lagrangian finite-difference scheme. I - Adiabatic formulation

An adiabatic global multilevel primitive equation model using a two time-level, semi-Lagrangian semi-implicit finite-difference integration scheme is presented. A Lorenz grid is used for vertical discretization and a C grid for the horizontal discretization. The momentum equation is discretized in vector form, thus avoiding problems near the poles. The 3D model equations are reduced by a linear transformation to a set of 2D elliptic equations, whose solution is found by means of an efficient direct solver. The model (with minimal physics) is integrated for 10 days starting from an initialized state derived from real data. A resolution of 16 levels in the vertical is used, with various horizontal resolutions. The model is found to be stable and efficient, and to give realistic output fields. Integrations with time steps of 10 min, 30 min, and 1 h are compared, and the differences are found to be acceptable.

Bates, J. R.↗

Comments on lunar origin.

Consideration of the tidal evolution of the earth-moon system in the light of recent published studies of lunar origin. The subjects include the factors of orbital evolution, the auxiliary Model Zero, orbital evolution along the outgoing branch, capture models, capture probability, capture with break-up, the tidal evolution of circular and elliptical fragmented rings, and precession and collisions in a ring of fragments. A numerical integration of precessions is performed for the actual nonzero mass of the moon and for a lunar distance of about 22 earth radii. A simplified Model Zero is found to compare well with Goldreich's solution (1968). It is concluded that the moon must have become freely yielding to tidal deformation at distances less than about 8 earth radii. Arguments are given in favor of the capture theory as a basis for the genesis of the moon.

Opik, E. J.↗

Meteoroid capture into earth orbit by atmospheric drag

A probabilistic analysis, based on orbital mechanics, is performed to evaluate the chances of meteoroids entering into earth orbit and the potential population of such objects. The problem is addressed in terms of meteoroids on earth collision courses, slowed by atmosphere entry/exit, entering elliptical orbits with apogees above the atmosphere. The overall capture probability is derived by integrating the capture fraction dependence on velocity and size over the probability frequency distribution of the collision courses. Account is taken of eventual orbit decay, lunar infall and ejection after encounters with the moon. The results indicate that the probability that a natural 10-100 m diam object has achieved earth orbit is negligibly small.

Friedlander, A. L.↗

Boundary integral techniques for multi-connected domains

Attention is given to iterative methods applicable to such changing domain cases of elliptic problems in multiconnected domains as those of incompressible inviscid flow with free surfaces. One such formulation is developed and tested on circular and elliptic annuli, using interpolated quadrature points to maintain accuracy when annuli regions are thin. The technique is noted to be applicable to the study of accelerating thin fluid shells.

Baker, G. R.↗

Experimental measurements for extracting nonlinear invariants

Nonlinear integrable optics are a promising alternative approach to lattice design. The integrable optics test accelerator (IOTA) at Fermilab has been constructed for dedicated studies of magnetostatic elliptical elements as described by Danilov and Nagaitsev. The most compelling verification of correct implementation of the NIO lattice is direct observation of the analytically expected invariants. This report outlines the experimental and analytical methods for extracting the nonlinear invariants of motion from data gathered in the last IOTA run.

43 PARTICLE ACCELERATORS↗

A numerical scheme to solve unstable boundary value problems

A new iterative scheme for solving boundary value problems is presented. It consists of the introduction of an artificial time dependence into a modified version of the system of equations. Then explicit forward integrations in time are followed by explicit integrations backwards in time. The method converges under much more general conditions than schemes based in forward time integrations (false transient schemes). In particular it can attain a steady state solution of an elliptical system of equations even if the solution is unstable, in which case other iterative schemes fail to converge. The simplicity of its use makes it attractive for solving large systems of nonlinear equations.

Kalnay Derivas, E.↗

Domain identification in impedance computed tomography by spline collocation method

A method for estimating an unknown domain in elliptic boundary value problems is considered. The problem is formulated as an inverse problem of integral equations of the second kind. A computational method is developed using a splice collocation scheme. The results can be applied to the inverse problem of impedance computed tomography (ICT) for image reconstruction.

Kojima, Fumio↗

The solar electric propulsion stage concept for high energy missions.

Definition of multimission and engine performance requirements for candidate solar electric propulsion stage configurations, considering launch vehicle compatibility, electric propulsion integration, payload requirements, and the effects of environmental extremes. Electric propulsion power options include two solar array power levels (15/22 kW), up to twelve electric thrustors of 30 cm diameter and 2.7 kW each, five to eight power conditioning units, and a maximum mercury propellant capacity of 1530 kg. In performance, the stage with a dry weight of 700 to 900 kg can deliver a net mass of 756 kg into Saturn orbit, 329 kg into a tight Mercury orbit, and 334 kg within 0.1 AU of the sun. The stage can also deliver a round trip payload of 3350 kg to geosynchronous orbit and return from an intermediate elliptical orbit using the Shuttle/Tug. Thus, a versatile stage is developed which competes effectively in performance with existing integrated spacecraft and promises considerable savings in total program costs.

Guttman, C. H.↗

Quantum properties of non-Dirichlet boundary conditions in gravity

The Euclidean path integral for gravity is enriched by the addition of boundaries, which provide useful probes of thermodynamic properties. Common boundary conditions include Dirichlet conditions on the boundary induced metric; microcanonical conditions, which refers to fixing some components of the Brown-York boundary stress tensor; and conformal conditions, in which the conformal structure of the induced metric and the trace of the extrinsic curvature are fixed. Boundaries also present interesting problems of consistency. The Dirichlet problem is known, under various (and generally different) conditions, to be inconsistent with perturbative quantization of graviton fluctuations, to exhibit thermodynamic instability, or to require infinite fine-tuning in the presence of matter fluctuations. We extend some of these results to other boundary conditions. We find that similarly to the Dirichlet problem, the graviton fluctuation operator is not elliptic with microcanonical boundaries, and the nonelliptic modes correspond to “boundary-moving diffeomorphisms.” However, we argue that microcanonical factorization of path integrals—essentially, the insertion of microcanonical constraints on two-sided surfaces in the bulk—is not affected by the same issues of ellipticity. We also show that for a variety of matter field boundary conditions, matter fluctuations renormalize the gravitational bulk and boundary terms differently, so that the classical microcanonical or conformal variational problems are not preserved unless an infinite fine-tuning is performed.

Draper, Patrick [Univ. of Illinois at Urbana-Champ↗

Prediction of drag at subsonic and transonic speeds using Euler methods

A technique for the evaluation of aerodynamic drag from flowfield solutions based on the Euler equations is discussed. The technique is limited to steady attached flows around three-dimensional configurations in the absence of active systems such as surface blowing/suction and propulsion. It allows the decomposition of the total drag into induced drag and wave drag and, consequently, it provides more information on the drag sources than the conventional surface-pressure integration technique. The induced drag is obtained from the integration of the kinetic energy (per unit distance) of the trailing vortex system on a wake plane and the wave drag is obtained from the integration of the entropy production on a plane just downstream of the shocks. The drag-evaluation technique is applied to three-dimensional flowfield solutions for the ONERA M6 wing as well as an aspect-ratio-7 wing with an elliptic spanwise chord distribution and an NACA-0012 section shape. Comparisons between the drag obtained with the present technique and the drag based on the integration of surface pressures are presented for two Euler codes.

Nikfetrat, K.↗

Elliptical cutouts in cylindrical shells

The stress concentrations due to an elliptical cutout in a shallow cylindrical shell are determined by means of a system of singular integral equations, solved numerically. Three loading conditions are considered: tension, and internal pressure. Results confirm those obtained by previous authors by other methods and extend the range of parameters over which results are available.

Tsai, C.-J.↗

The intermediate anomaly

Time transformations of the equation dt = cr to the n ds, where s is a variable called the intermediate anomaly, are known to reduce global error in the solution of gravitational systems obtained by numerical integration. Attention is given to the Sundman time transformation, and its relation to equations of Keplerian elliptical motion.

Nacozy, P.↗

Encke-Beta Predictor for Orion Burn Targeting and Guidance

The state vector prediction algorithm selected for Orion on-board targeting and guidance is known as the Encke-Beta method. Encke-Beta uses a universal anomaly (beta) as the independent variable, valid for circular, elliptical, parabolic, and hyperbolic orbits. The variable, related to the change in eccentric anomaly, results in integration steps that cover smaller arcs of the trajectory at or near perigee, when velocity is higher. Some burns in the EM-1 and EM-2 mission plans are much longer than burns executed with the Apollo and Space Shuttle vehicles. Burn length, as well as hyperbolic trajectories, has driven the use of the Encke-Beta numerical predictor by the predictor/corrector guidance algorithm in place of legacy analytic thrust and gravity integrals.

Robinson, Shane↗

Measurements at forward rapidity of elliptic flow of charged hadrons and open-heavy-flavor muons in Au + Au collisions at $\sqrt{𝑠_{𝑁⁢𝑁}}$ = 200 GeV

Here, we present the first forward-rapidity measurements of elliptic anisotropy of open-heavy-flavor muons at the Relativistic Heavy Ion Collider. The measurements are based on data samples of Au + Au collisions at $\sqrt{𝑠_{𝑁⁢𝑁}}$ = 200 GeV collected by the PHENIX experiment in 2014 and 2016 with integrated luminosity of 14.5 nb −1 . The measurements are performed in the pseudorapidity range 1.2 < |𝜂| < 2 and cover transverse momenta 1< 𝑝 𝑇 < 4 GeV/𝑐. The elliptic flow of charged hadrons as a function of transverse momentum is also measured in the same kinematic range. We observe significant elliptic flow for both charged hadrons and heavy-flavor muons. The results show clear mass ordering of elliptic flow of light- and heavy-flavor particles. The magnitude of the measured 𝑣 2 is comparable to that in the midrapidity region. This indicates that there is no strong longitudinal dependence in the quark-gluon-plasma evolution between midrapidity and the rapidity range of this measurement at $\sqrt{𝑠_{𝑁⁢𝑁}}$ = 200 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Downstream boundary conditions for viscous flow problems

The problem of the specification of artificial outflow conditions in flow problems is studied. It is shown that for transport type equations incorrect outflow conditions will adversely affect the solution only in a small region near the outflow boundary, while for elliptic equations, e.g. those governing the streamfunction or pressure, a correct boundary specification is essential. In addition, integral outflow boundary conditions for fluid dynamical problems are considered. It is shown that such conditions are well posed, and their effect on the solutions of the Navier-Stokes equations is also considered.

Fix, G.↗

The soaring kite: a tale of two punctured tori

We consider the 5-mass kite family of self-energy Feynman integrals and present a systematic approach for constructing an ε-form basis, along with its differential equation pulled back onto the moduli space of two tori. Each torus is associated with one of the two distinct elliptic curves this family depends on. We demonstrate how the locations of relevant punctures, which are required to parametrize the full image of the kinematic space onto this moduli space, can be extracted from integrals over maximal cuts. A boundary value is provided such that the differential equation is systematically solved in terms of iterated integrals over g-kernels and modular forms. Then, the numerical evaluation of the master integrals is discussed, and important challenges in that regard are emphasized. In an appendix, we introduce new relations between g-kernels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗