Virtual mass technique for computing N-body solutions
N body numerical integrator using virtual mass technique for computing trajectories with speed and accuracy suitable for spacecraft onboard computer applications
Novak, D. H.↗
SEARCH · Search NASA
Results for “Integral solution”
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.
Solution of elastoplastic torsion problem by boundary integral method
The boundary integral method was applied to the elastoplastic analysis of the torsion of prismatic bars, and the results are compared with those obtained by the finite difference method. Although fewer unknowns were used, very good accuracy was obtained with the boundary integral method. Both simply and multiply connected bodies can be handled with equal ease.
Mendelson, A.↗
Significance of Strain in Formulation in Theory of Solid Mechanics
The basic theory of solid mechanics was deemed complete circa 1860 when St. Venant provided the strain formulation or the field compatibility condition. The strain formulation was incomplete. The missing portion has been formulated and identified as the boundary compatibility condition (BCC). The BCC, derived through a variational formulation, has been verified through integral theorem and solution of problems. The BCC, unlike the field counterpart, do not trivialize when expressed in displacements. Navier s method and the stiffness formulation have to account for the extra conditions especially at the inter-element boundaries in a finite element model. Completion of the strain formulation has led to the revival of the direct force calculation methods: the Integrated Force Method (IFM) and its dual (IFMD) for finite element analysis, and the completed Beltrami-Michell formulation (CBMF) in elasticity. The benefits from the new methods in elasticity, in finite element analysis, and in design optimization are discussed. Existing solutions and computer codes may have to be adjusted for the compliance of the new conditions. Complacency because the discipline is over a century old and computer codes have been developed for half a century can lead to stagnation of the discipline.
Patnaik, Surya N.↗
The Emergence of Agent-Based Technology as an Architectural Component of Serious Games
The evolution of games as an alternative to traditional simulations in the military context has been gathering momentum over the past five years, even though the exploration of their use in the serious sense has been ongoing since the mid-nineties. Much of the focus has been on the aesthetics of the visuals provided by the core game engine as well as the artistry provided by talented development teams to produce not only breathtaking artwork, but highly immersive game play. Consideration of game technology is now so much a part of the modeling and simulation landscape that it is becoming difficult to distinguish traditional simulation solutions from game-based approaches. But games have yet to provide the much needed interactive free play that has been the domain of semi-autonomous forces (SAF). The component-based middleware architecture that game engines provide promises a great deal in terms of options for the integration of agent solutions to support the development of non-player characters that engage the human player without the deterministic nature of scripted behaviors. However, there are a number of hard-learned lessons on the modeling and simulation side of the equation that game developers have yet to learn, such as: correlation of heterogeneous systems, scalability of both terrain and numbers of non-player entities, and the bi-directional nature of simulation to game interaction provided by Distributed Interactive Simulation (DIS) and High Level Architecture (HLA).
Phillips, Mark↗
The Importance of Data Visualization: Incorporating Storytelling into the Scientific Presentation
From its inception in 2000, one of the primary tasks of the Biomedical Data Reduction Analysis (BDRA) group has been translation of large amounts of data into information that is relevant to the audience receiving it. BDRA helps translate data into an integrated model that supports both operational and research activities. This data integrated model and subsequent visual data presentations have contributed to BDRA's success in delivering the message (i.e., the story) that its customers have needed to communicate. This success has led to additional collaborations among groups that had previously not felt they had much in common until they worked together to develop solutions in an integrated fashion. As more emphasis is placed on working with "big data" and on showing how NASA's efforts contribute to the greater good of the American people and of the world, it becomes imperative to visualize the story of our data to communicate the greater message we need to share. METHODS To create and expand its data integrated model, BDRA has incorporated data from many different collaborating partner labs and other sources. Data are compiled from the repositories of the Lifetime Surveillance of Astronaut Health and the Life Sciences Data Archive, and from the individual laboratories at Johnson Space Center that support collection of data from medical testing, environmental monitoring, and countermeasures, as designated in the Medical Requirements Integration Documents. Ongoing communication with the participating collaborators is maintained to ensure that the message and story of the data are retained as data are translated into information and visual data presentations are delivered in different venues and to different audiences. RESULTS We will describe the importance of storytelling through an integrated model and of subsequent data visualizations in today's scientific presentations and discuss the collaborative methods used. We will illustrate the discussion with examples of graphs from BDRA's past work supporting operations and/or research efforts.
Babiak-Vazquez, A.↗
Driving Force Analysis in an Infinite Anisotropic Plate with Multiple Crack Interactions
The methodology and a rigorous solution formulation are presented for stress intensity factors (SIF's, k) and total strain energy release rates (SERR, G(sub T)) of a multicracked plate, that has fully interacting cracks and is subjected to a far-field arbitrary stress state. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. Parametric studies are conducted for two, three, and four crack problems. The sensitivity and characteristics of the model is demonstrated.
Binienda, Wieslaw K.↗
Viscous wing theory development. Volume 1: Analysis, method and results
Viscous transonic flows at large Reynolds numbers over 3-D wings were analyzed using a zonal viscid-inviscid interaction approach. A new numerical AFZ scheme was developed in conjunction with the finite volume formulation for the solution of the inviscid full-potential equation. A special far-field asymptotic boundary condition was developed and a second-order artificial viscosity included for an improved inviscid solution methodology. The integral method was used for the laminar/turbulent boundary layer and 3-D viscous wake calculation. The interaction calculation included the coupling conditions of the source flux due to the wing surface boundary layer, the flux jump due to the viscous wake, and the wake curvature effect. A method was also devised incorporating the 2-D trailing edge strong interaction solution for the normal pressure correction near the trailing edge region. A fully automated computer program was developed to perform the proposed method with one scalar version to be used on an IBM-3081 and two vectorized versions on Cray-1 and Cyber-205 computers.
Chow, R. R.↗
An exterior Poisson solver using fast direct methods and boundary integral equations with applications to nonlinear potential flow
A general method is developed combining fast direct methods and boundary integral equation methods to solve Poisson's equation on irregular exterior regions. The method requires O(N log N) operations where N is the number of grid points. Error estimates are given that hold for regions with corners and other boundary irregularities. Computational results are given in the context of computational aerodynamics for a two-dimensional lifting airfoil. Solutions of boundary integral equations for lifting and nonlifting aerodynamic configurations using preconditioned conjugate gradient are examined for varying degrees of thinness.
Young, D. P.↗
On the theory of the synthesis of single and dual offset shaped reflector antennas
Two reflectors shaped so as to transforn a geometrical optics (GO) feed power pattern into a uniform amplitude and phase aperture distribution are indicated schematically. Single- and dual-shaped geometrical optics problems are solved via exact solutions obtained by progressive integration. It is concluded that an infinite set of such solutions exists.
Galindo-Israel, Victor↗
Large-scale structure of the solar corona and inner heliosphere
The large-scale structure of the solar corona influences solar activity particularly coronal mass ejections (CMEs). The helmet streamers that are observed to dominate the structure of the inner corona are formed by the interaction of the solar wind with coronal magnetic fields. We have simulated this interaction in three dimensions using the magnetohydrodynamic (MHD) equations. In order to create a realistic model, we use the magnetic field that is observed at the Sun's surface (deduced from daily Wilcox Solar Observatory magnetograms) as input, in combination with specified density and temperature profiles at the surface. A self-consistent 3D solar-wind solution is developed by integrating the MHD equations in time to steady state. Such solutions can reproduce the observed structures that are seen in coronagraph images and eclipse photographs of the corona. This model allows us to accurately determine the position of the heliospheric current sheet. We will compare the results obtained from our model with Ulysses observations during the period May-June 1993, and with an eclipse photograph of the corona on November 3, 1994.
Mikic, Z.↗
Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements
This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.
Gould, Dana C.↗
Quasi 1D Modeling of Mixed Compression Supersonic Inlets
The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.
Kopasakis, George↗
Second-order small-disturbance solutions for hypersonic flow over power-law bodies
Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.
Townsend, J. C.↗
Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations
Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.
CFD↗
Photovoltaic array reliability optimization
An overview of the photovoltaic array reliability problem is presented, and a high reliability/minimum cost approach to this problem is presented. Design areas covered are cell failure, interconnect fatigue, and electrical insulation breakdown, and three solution strategies are discussed. The first involves controlling component failures in the solar cell (cell cracking, cell interconnects) and at the module level (must be statistically treated). Second, a fault tolerant circuit is designed which reduces array degradation, improves module yield losses, and controls hot-spot heating. Third, cost optimum module replacement strategies are also effective in reducing array degradation. This can be achieved by minimizing the life-cycle energy cost of the photovoltaic system. The integration of these solutions is aimed at reducing the 0.01% failure rate.
Ross, R. G., Jr.↗
Control of Combustion-Instabilities Through Various Passive Devices
It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.
Frendi, Kader↗
Fast two-stream method for computing diurnal-mean actinic flux in vertically inhomogeneous atmospheres
Based on a derivation of the two-stream daytime-mean equations of radiative flux transfer, a method for computing the daytime-mean actinic fluxes in the absorbing and scattering vertically inhomogeneous atmosphere is suggested. The method applies direct daytime integration of the particular solutions of the two-stream approximations or the source functions. It is valid for any duration of period of averaging. The merit of the method is that the multiple scattering computation is carried out only once for the whole averaging period. It can be implemented with a number of widely used two-stream approximations. The method agrees with the results obtained with 200-point multiple scattering calculations. The method was also tested in runs with a 1-km cloud layer with optical depth of 10, as well as with aerosol background. Comparison of the results obtained for a cloud subdivided into 20 layers with those obtained for a one-layer cloud with the same optical parameters showed that direct integration of particular solutions possesses an 'analytical' accuracy. In the case of the source function interpolation, the actinic fluxes calculated above the one-layer and 20-layer clouds agreed within 1%-1.5%, while below the cloud they may differ up to 5% (in the worst case). The ways of enhancing the accuracy (in a 'two-stream sense') and computational efficiency of the method are discussed.
Filyushkin, V. V.↗
Vibration suppression using a constrained rate-feedback Threshold control strategy
Quasi-closed form solutions are derived for the finite time, minimum force rate-feedback threshold controller to bring a system with or without known external disturbances back into an 'allowable' state manifold in finite time. The disturbances are assumed to be expandable in terms of Fourier series. The quasi-closed form solutions replace the solution of the two-point boundary value problem and definite integral constraints with the solution of algebraic equations and the calculation of matrix exponentials. Examples demonstrate the threshold control technique and compare the quasi-closed form solutions with MACSYMA generated exact solutions (for small system order) and with the numerical solution of the two-point boundary value problem.
Zimmerman, D. C.↗