OFFSET-AIM TARGET SEEKER TECHNIQUE FOR INTERPLANETARY BALLISTIC TRAJECTORIES
Iteration method using offset-aim technique to determine exact initial conditions for interplanetary ballistic trajectories
INTERPLANETARY TRAJECTORY↗
SEARCH · Engineering Papers
Results for “EXACT SOLUTIONS”
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.
The enigma of the Uranian satellites' orbital eccentricities
The eccentricity decay times for the Uranian satellites are calculated using recent observations (Brown et al., 1982) of the diameters and orbital elements of the satellites and assuming reasonable dissipation functions and rigidities for icy satellites. For the outer two satellites, Titania and Oberon, the decay times are found to be very long, whereas the inner three satellites, Miranda, Ariel, and Umbriel, have decay times on the order of 10 to the 7th to 10 to the 8th years and have a near-commensurability in their mean motions that cannot force their eccentricities. There are several possible solutions for the lack of resonant forcing: (1) the reported eccentricities are incorrect, and are very nearly zero, (2) the reported mean motions are incorrect, and an exact commensurability exists, (3) the physical properties assumed for the satellites are grossly in error, and (4) the system is evolving rapidly, perhaps from a previous state of higher eccentricity. A new lower bound of about 17,000 on the dissipation function of Uranus is calculated from the mass of Ariel and its proximity to Uranus.
Squyres, S. W.↗
Development of an integrated BEM approach for hot fluid structure interaction
The progress made toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-Orbit engine hot section components is reported. The convective viscous integral formulation was derived and implemented in the general purpose computer program GP-BEST. The new convective kernel functions, in turn, necessitated the development of refined integration techniques. As a result, however, since the physics of the problem is embedded in these kernels, boundary element solutions can now be obtained at very high Reynolds number. Flow around obstacles can be solved approximately with an efficient linearized boundary-only analysis or, more exactly, by including all of the nonlinearities present in the neighborhood of the obstacle. The other major accomplishment was the development of a comprehensive fluid-structure interaction capability within GP-BEST. This new facility is implemented in a completely general manner, so that quite arbitrary geometry, material properties and boundary conditions may be specified. Thus, a single analysis code (GP-BEST) can be used to run structures-only problems, fluids-only problems, or the combined fluid-structure problem. In all three cases, steady or transient conditions can be selected, with or without thermal effects. Nonlinear analyses can be solved via direct iteration or by employing a modified Newton-Raphson approach.
Dargush, G. F.↗
Magnetic Control of Solutal Buoyancy Driven Convection
Volumetric forces resulting from local density variations and gravitational acceleration cause buoyancy induced convective motion in melts and solutions. Solutal buoyancy is a result of concentration differences in an otherwise isothermal fluid. If the fluid also exhibits variations in magnetic susceptibility with concentration then convection control by external magnetic fields can be hypothesized. Magnetic control of thermal buoyancy induced convection in ferrofluids (dispersions of ferromagnetic particles in a carrier fluid) and paramagnetic fluids have been demonstrated. Here we show the nature of magnetic control of solutal buoyancy driven convection of a paramagnetic fluid, an aqueous solution of Manganese Chloride hydrate. We predict the critical magnetic field required for balancing gravitational solutal buoyancy driven convection and validate it through a simple experiment. We demonstrate that gravity driven flow can be completely reversed by a magnetic field but the exact cancellation of the flow is not possible. This is because the phenomenon is unstable. The technique can be applied to crystal growth processes in order to reduce convection and to heat exchanger devices for enhancing convection. The method can also be applied to impose a desired g-level in reduced gravity applications.
Ramachandran, N.↗
Error Analysis for Discontinuous Galerkin Method for Parabolic Problems
In the proposal, the following three objectives are stated: (1) A p-version of the discontinuous Galerkin method for a one dimensional parabolic problem will be established. It should be recalled that the h-version in space was used for the discontinuous Galerkin method. An a priori error estimate as well as a posteriori estimate of this p-finite element discontinuous Galerkin method will be given. (2) The parameter alpha that describes the behavior double vertical line u(sub t)(t) double vertical line 2 was computed exactly. This was made feasible because of the explicitly specified initial condition. For practical heat transfer problems, the initial condition may have to be approximated. Also, if the parabolic problem is proposed on a multi-dimensional region, the parameter alpha, for most cases, would be difficult to compute exactly even in the case that the initial condition is known exactly. The second objective of this proposed research is to establish a method to estimate this parameter. This will be done by computing two discontinuous Galerkin approximate solutions at two different time steps starting from the initial time and use them to derive alpha. (3) The third objective is to consider the heat transfer problem over a two dimensional thin plate. The technique developed by Vogelius and Babuska will be used to establish a discontinuous Galerkin method in which the p-element will be used for through thickness approximation. This h-p finite element approach, that results in a dimensional reduction method, was used for elliptic problems, but the application appears new for the parabolic problem. The dimension reduction method will be discussed together with the time discretization method.
Kaneko, Hideaki↗
Nuclear Cross Sections for Space Radiation Applications
The eikonal, partial wave (PW) Lippmann-Schwinger, and three-dimensional Lippmann-Schwinger (LS3D) methods are compared for nuclear reactions that are relevant for space radiation applications. Numerical convergence of the eikonal method is readily achieved when exact formulas of the optical potential are used for light nuclei (A = 16) and the momentum-space optical potential is used for heavier nuclei. The PW solution method is known to be numerically unstable for systems that require a large number of partial waves, and, as a result, the LS3D method is employed. The effect of relativistic kinematics is studied with the PW and LS3D methods and is compared to eikonal results. It is recommended that the LS3D method be used for high energy nucleon-nucleus reactions and nucleus-nucleus reactions at all energies because of its rapid numerical convergence and stability for both non-relativistic and relativistic kinematics.
Werneth, C. M.↗
A comparative study of an ABC and an artificial absorber for truncating finite element meshes
The type of mesh termination used in the context of finite element formulations plays a major role on the efficiency and accuracy of the field solution. The performance of an absorbing boundary condition (ABC) and an artificial absorber (a new concept) for terminating the finite element mesh was evaluated. This analysis is done in connection with the problem of scattering by a finite slot array in a thick ground plane. The two approximate mesh truncation schemes are compared with the exact finite element-boundary integral (FEM-BI) method in terms of accuracy and efficiency. It is demonstrated that both approximate truncation schemes yield reasonably accurate results even when the mesh is extended only 0.3 wavelengths away from the array aperture. However, the artificial absorber termination method leads to a substantially more efficient solution. Moreover, it is shown that the FEM-BI method remains quite competitive with the FEM-artificial absorber method when the FFT is used for computing the matrix-vector products in the iterative solution algorithm. These conclusions are indeed surprising and of major importance in electromagnetic simulations based on the finite element method.
Oezdemir, T.↗
Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation
Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.
ESG↗
Magnetic flux rope versus the spheromak as models for interplanetary magnetic clouds
Magnetic clouds form a subset of interplanetary ejecta with well-defined magnetic and thermodynamic properties. Observationally, it is well established that magnetic clouds expand as they propagate antisunward. The aim of this paper is to compare and contrast two models which have been proposed for the global magnetic field line topology of magnetic clouds: a magnetic flux tube geometry, on the one hand, and a spheromak geometry (including possible higher multiples), on the other. Traditionally, the magnetic structure of magnetic clouds has been modeled by force-free configurations. In a first step, we therefore analyze the ability of static force-free models to account for the asymmetries observed in the magnetic field profiles of magnetic clouds. For a cylindrical flux tube the magnetic field remains symmetric about closest approach to the magnetic axis on all spacecraft orbits intersecting it, whereas in a spheromak geometry one can have asymmetries in the magnetic field signatures along some spacecraft trajectories. The duration of typical magnetic cloud encounters at 1 AU (1 to 2 days) is comparable to their travel time from the Sun to 1 AU and thus magnetic clouds should be treated as strongly nonstationary objects. In a second step, therefore, we abandon the static approach and model magnetic clouds as self-similarly evolving MHD configurations. In our theory, the interaction of the expanding magnetic cloud with the ambient plasma is taken into account by a drag force proportional to the density and the velocity of expansion. Solving rigorously the full set of MHD equations, we demonstrate that the asymmetry in the magnetic signature may arise solely as a result of expansion. Using asymptotic solutions of the MHD equations, we least squares fit both theoretical models to interplanetary data. We find that while the central part of the magnetic cloud is adequately described by both models, the 'edges' of the cloud data are modeled better by the magnetic flux tube. Further comparisons of the two models necessarily involve thermodynamic properties, since real magnetic configurations are never exactly force-free and gas pressure plays an essential role. We consider a polytropic gas. Our theoretical analysis shows that the self-similar expansion of a magnetic flux tube requires the polytropic index gamma to be less than unity. For the spheromak, however, self-similar, radially expanding solutions are known only for gamma equal to 4/3. This difference, therefore, yields a good way of distinguishing between the two geometries. It has been shown recently that the polytropic relationship is applicable to magnetic clouds and that the corresponding polytropic index is approximately 0.5. This observational result is consistent with the self-similar model of the magnetic flux rope but is in conflict with the self-similar spheromak model.
Farrugia, C. J.↗
Numerical computation of three-dimensional blunt body flow fields with an impinging shock
A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.
Holst, T. L.↗
Response of linear dynamic systems with random coefficients
Numerous models of physical systems contain parameters whose values are not known exactly. The physical and mathematical complexities arising in the prediction of the statistical behavior of such systems are discussed. Although the discussions are far from providing a satisfactory solution to such problems, they perhaps, by utilization of simple examples, will create a greater awareness of the statistical effect of random parameters.
Dickerson, J.↗
An approximate analytic solution of a set of nonlinear model alpha-omega-dynamo equations for marginally unstable systems
An approximate analytic solution of a set of nonlinear model alpha-omega-dynamo equations is obtained. The reaction of the Lorentz force on the velocity shear which stretches and, hence, amplifies the magnetic field is incorporated into the model. To single out the effect of the Lorentz force on the omega-effect, the effect of the Lorentz force on the alpha-effect is neglected in this study. The solution represents a nonlinear oscillation with the amplitude and period determined by the dynamo number N. The amplitude is proportional to N - 1, while the period is almost exactly the same as the dissipation time of the unstable mode (proportional to N).
Hinata, S.↗
Electromagnetic scattering from two-dimensional thick material junctions
The problem of the plane wave diffraction is examined by an arbitrary symmetric two dimensional junction, where Generalized Impedance Boundary Conditions (GIBCs) and Generalized Sheet Transition Conditions (GSTCs) are employed to simulate the slabs. GIBCs and GSTCs are constructed for multilayer planar slabs of arbitrary thickness and the resulting GIBC/GSTC reflection coefficients are compared with exact counterparts to evaluate the GIBCs/GSTCs. The plane wave diffraction by a multilayer material slab recessed in a perfectly conducting ground plane is formulated and solved via the Generalized Scattering Matrix Formulation (GDMF) in conjunction with the dual integral equation approach. Various scattering patterns are computed and validated with exact results where possible. The diffraction by a material discontinuity in a thick dielectric/ferrite slab is considered by modelling the constituent slabs with GSTCs. A non-unique solution in terms of unknown constants is obtained, and these constants are evaluated for the recessed slab geometry by comparison with the solution obtained therein. Several other simplified cases are also presented and discussed. An eigenfunction expansion method is introduced to determine the unknown solution constants in the general case. This procedure is applied to the non-unique solution in terms of unknown constants; and scattering patterns are presented for various slab junctions and compared with alternative results where possible.
Ricoy, M. A.↗
Multilevel Error Estimation and Adaptive h-Refinement for Cartesian Meshes with Embedded Boundaries
This paper presents the development of a mesh adaptation module for a multilevel Cartesian solver. While the module allows mesh refinement to be driven by a variety of different refinement parameters, a central feature in its design is the incorporation of a multilevel error estimator based upon direct estimates of the local truncation error using tau-extrapolation. This error indicator exploits the fact that in regions of uniform Cartesian mesh, the spatial operator is exactly the same on the fine and coarse grids, and local truncation error estimates can be constructed by evaluating the residual on the coarse grid of the restricted solution from the fine grid. A new strategy for adaptive h-refinement is also developed to prevent errors in smooth regions of the flow from being masked by shocks and other discontinuous features. For certain classes of error histograms, this strategy is optimal for achieving equidistribution of the refinement parameters on hierarchical meshes, and therefore ensures grid converged solutions will be achieved for appropriately chosen refinement parameters. The robustness and accuracy of the adaptation module is demonstrated using both simple model problems and complex three dimensional examples using meshes with from 10(exp 6), to 10(exp 7) cells.
Aftosmis, M. J.↗
Modelling Sonic Boom Propagation through Planetary Boundary Layer Turbulence near the Lateral Extent of the Carpet
Nonlinear acoustic propagation of sonic booms in the atmospheric boundary layer is considered in the context of a one-way solution to a third order wave equation. A split-step approach is utilized to compute different physical effects efficiently. Unlike previous approaches, the diffraction effects are computed exactly in the forward direction with no restriction on the angle of propagation. This results in more accurate modelling of sonic boom near the carpet edge. Heterogenous flow effects are incorporated with a wide-angle parabolic approximation. Nonlinear propagation is computed with a Burgers-Hayes method. Turbulence in the medium is constructed by the method of random Fourier modes. The turbulence spectra are constructed using an altitude dependent 3-dimensional von Karman spectrum. The turbulence field is considered frozen, as the eddy turnover time in the atmospheric boundary layer is generally much larger than the propagation time of the acoustic wave. Comparisons with benchmark predictions from the PCBoom software are conducted. The benchmark cases provide insight into the accuracy of the current prediction code for non-turbulent atmospheres. Future work is discussed regarding the prediction of shaped booms, boundary conditions, and predictions near the lateral extent of the boom carpet where previous approaches fail to capture the diffraction effects.
sonic boom↗
Gravitational waves from rotating and precessing rigid bodies. 2: General solutions and computationally useful formulae
The classical mechanics results for free precession which are needed in order to calculate the weak field, slow-motion, quadrupole-moment gravitational waves are reviewed. Within that formalism, algorithms are given for computing the exact gravitational power radiated and waveforms produced by arbitrary rigid-body freely-precessing sources. The dominant terms are presented in series expansions of the waveforms for the case of an almost spherical object precessing with a small wobble angle. These series expansions, which retain the precise frequency dependence of the waves, may be useful for gravitational astronomers when freely-precessing sources begin to be observed.
Zimmerman, M.↗
Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction
High order accurate finite-volume schemes for solving the Euler equations of gasdynamics are developed. Central to the development of these methods are the construction of a k-exact reconstruction operator given cell-averaged quantities and the use of high order flux quadrature formulas. General polygonal control volumes (with curved boundary edges) are considered. The formulations presented make no explicit assumption as to complexity or convexity of control volumes. Numerical examples are presented for Ringleb flow to validate the methodology.
Barth, Timothy J.↗
Remarks on nonstandard least-squares problems
It is shown that after treating a least-squares problem in the standard fashion, the imposition of stochastic restrictions on certain parameters can be accomplished essentially by inverting a matrix of an order no larger than the total number of restricted parameters, and that generally there is an explicit term for correcting the unrestricted to the restricted solution. In addition, any known type of nonstandard least-squares adjustment problem can be reduced to an equivalent standard linear problem, that is, one in which there is exactly one observation per condition equation and in which the observations are not correlated.
Eichhorn, H.↗