Global solutions by analytic parameter continuation for a class of non-linear singular integral equations
Solution of non-linear singular integral equations from hydroelasticity and electroelasticity by analytic parameter continuation
SEARCH · Engineering Papers
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 non-linear singular integral equations from hydroelasticity and electroelasticity by analytic parameter continuation
Direct analytical solutions can be useful in validating the Core formulation of numerical systems. In this document an exact analytical solution to the nonlinear Navier-Stokes equation is compared to the numerical results from the three-dimensional Terminal Area Simulation System (TASS). This exact solution, of which the derivation is included, is for Beltrami type flow. Direct comparison of TASS to the analytical Beltrami solution is then used in evaluating the accuracy of TASS.
Semi-analytic reconstruction uses the analytic solution to a second-order, steady, ordinary differential equation (ODE) to simultaneously evaluate the convective and diffusive flux at all interfaces of a finite volume formulation. The second-order ODE is itself a linearized approximation to the governing first- and second- order partial differential equation conservation laws. Thus, semi-analytic reconstruction defines a family of formulations for finite volume interface fluxes using analytic solutions to approximating equations. Limiters are not applied in a conventional sense; rather, diffusivity is adjusted in the vicinity of changes in sign of eigenvalues in order to achieve a sufficiently small cell Reynolds number in the analytic formulation across critical points. Several approaches for application of semi-analytic reconstruction for the solution of one-dimensional scalar equations are introduced. Results are compared with exact analytic solutions to Burger s Equation as well as a conventional, upwind discretization using Roe s method. One approach, the end-point wave speed (EPWS) approximation, is further developed for more complex applications. One-dimensional vector equations are tested on a quasi one-dimensional nozzle application. The EPWS algorithm has a more compact difference stencil than Roe s algorithm but reconstruction time is approximately a factor of four larger than for Roe. Though both are second-order accurate schemes, Roe s method approaches a grid converged solution with fewer grid points. Reconstruction of flux in the context of multi-dimensional, vector conservation laws including effects of thermochemical nonequilibrium in the Navier-Stokes equations is developed.
A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.
Knowledge of thermal diffusivity data for high-temperature materials (solids and liquids) is very important in analyzing a number of processes, among them solidification, crystal growth, and welding. However, reliable thermal diffusivity versus temperature data, particularly those for high-temperature liquids, are still far from complete. The main measurement difficulties are due to the presence of convection and the requirement for a container. Fortunately, the availability of levitation techniques has made it possible to solve the containment problem. Based on the feasibility of the levitation technology, a new laser flash technique which is applicable to both levitated liquid and solid samples is being developed. At this point, the analysis for solid samples is near completion and highlights of the technique are presented here. The levitated solid sample which is assumed to be a sphere is subjected to a very short burst of high power radiant energy. The temperature of the irradiated surface area is elevated and a transient heat transfer process takes place within the sample. This containerless process is a two-dimensional unsteady heat conduction problem. Due to the nonlinearity of the radiative plus convective boundary condition, an analytic solution cannot be obtained. Two options are available at this point. Firstly, the radiation boundary condition can be linearized, which then accommodates a closed-form analytic solution. Comparison of the analytic curves for the temperature rise at different points to the experimentally-measured values will then provide the thermal diffusivity values. Secondly, one may set up an inverse conduction problem whereby experimentally obtained surface temperature history is used as the boundary conditions. The thermal diffusivity can then be elevated by minimizing the difference between the real heat flux boundary condition (radiation plus convection) and the measurements. Status of an experimental study directed at measuring the thermal diffusivity of high-temperature solid samples of pure Nickel and Inconel 718 superalloys are presented. Preliminary measurements showing surface temperature histories are discussed.
A technique for solving n-dimensional games is developed and applied to two pursuit-evasion games. The first is a two-dimensional game similar to the homicidal chauffeur but modified to resemble an airplane-helicopter engagement. The second is a five-dimensional game of two airplanes at constant altitude and with thrust and turning controls. The performance function to be optimized by the pursuer and evader was the distance between the evader and a given target point in front of the pursuer. The analytic solution to the first game reveals that both unique and nonunique solutions exist. A comparison between the gradient results and the analytic solution shows a dependence on the nominal controls in regions where nonunique solutions exist. In the unique solution region, the results from the two methods agree closely. The results for the five-dimensional two-airplane game are also shown to be dependent on the nominal controls selected and indicate that initial conditions are in a region of nonunique solutions.
A defect stream function formulation for nonadiabatic flow with small crossflow is developed. The first-integral property of this formulation provides for two removal of the streamline curvature term in the governing equation so that the form of the reduced equation for small crossflow is the same as that for two dimensional flow. The combined law of the wall and wake is used in place of the no-slip boundary condition. The tangential velocity equation for law-of-the-wall flow is shown to be the same for three-dimensions as for two when the Boussinesq approximation applies, and a closed form solution for the crossflow angle in the inner region is obtained. Analytic solutions for nonadiabatic, compressible, equilibrium flow with a Clauser outer-region eddy-viscosity model are obtained, and excellent agreement with experimental skin friction and velocity profile data for nonadiabatic, compressible flat-plate flow is achieved. An analytic solution for a linear inner-region eddy-viscosity model is also obtained; the wake function part of this solution is found to be inconsistent with the empirically established law of the wake.
Approximate analytic solutions for transient and steady-state 180 deg domain wall motion in bulk magnetic material are obtained from the dynamic torque equations with a Gilbert damping term. The results for the Walker region in which the transient solution approaches the familiar Walker steady-state solution are presented in a slightly new form for completeness. An analytic solution corresponding to larger drive fields predicts an oscillatory motion with an average value which decreases with drive field for reasonable values of the damping parameter. These results agree with those obtained by a computer solution of the torque equation and those obtained with the assumption of a very large anisotropy field.
Approximate analytic solutions for transient and steady-state 180 deg domain-wall motion in bulk magnetic material are obtained from the dynamic torque equations with a Gilbert damping term. The results for the Walker region in which the transient solution approaches the familiar Walker steady-state solution are presented in a slightly new form for completeness. An analytic solution corresponding to larger drive fields predicts an oscillatory motion with an average value of the velocity which decreases with drive field for reasonable values of the damping parameter. These results agree with those obtained by others from a computer solution of the torque equation and those obtained by others with the assumption of a very large anisotropy field.
The intensity of sunlight multiply scattered in model atmospheres is derived from the equation of radiative transfer by an analytical small-angle approximation. The approximate analytical solutions are compared to rigorous numerical solutions of the same problem. Results obtained from an aerosol-laden model atmosphere are presented. Agreement between the rigorous and the approximate solutions is found to be within a few per cent. The analytical solution to the problem which considers an aerosol-laden atmosphere is then inverted to yield a phase function which describes a single scattering event at small angles. The effect of noisy data on the derived phase function is discussed.
This work summarizes theoretical and experimental concepts used to design the flight experiment mission for SHIVA - Spaceflight Holography Investigation in a Virtual Apparatus. SHIVA is a NASA project that exploits a unique, holography-based, diagnostics tool to understand the behavior of small particles subjected to transient accelerations. The flight experiments are designed for testing model equations, measuring g, g-jitter, and other microgravity phenomena. Data collection will also include experiments lying outside of the realm of existing theory. The regime under scrutiny is the low Reynolds number, Stokes regime or creeping flow, which covers particles and bubbles moving at very low velocity. The equations describing this important regime have been under development and investigation for over 100 years and yet a complete analytical solution of the general equation had remained elusive yielding only approximations and numerical solutions. In the course of the ongoing NASA NRA, the first analytical solution of the general equation was produced by members of the investigator team using the mathematics of fractional derivatives. This opened the way to an even more insightful and important investigation of the phenomena in microgravity. Recent results include interacting particles, particle-wall interactions, bubbles, and Reynolds numbers larger than unity. The Space Station provides an ideal environment for SHIVA. Limited ground experiments have already confirmed some aspects of the theory. In general the space environment is required for the overall experiment, especially for cases containing very heavy particles, very light particles, bubbles, collections of particles and for characterization of the space environment and its effect on particle experiments. Lightweight particles and bubbles typically rise too fast in a gravitational field and heavy particles sink too fast. In a microgravity environment, heavy and light particles can be studied side-by-side for long periods of time.
SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus) will expand our understanding of the fundamental physics of particle movement in fluids by exploiting the power of holography in a spaceflight experiment'. In addition, the study will exploit the movement of particles in fluids to observe and quantify microgravity phenomena that are extremely important in materials sciences with applications both in space and on earth. The regime under scrutiny is the low Reynolds number, Stokes regime or creeping flow, which covers particles and bubbles moving at very low velocity. The equations describing this important regime have been under development and investigation for over 100 years and yet a complete analytical solution of the general equation had remained elusive yielding only approximations and numerical solutions. In the course of the ongoing NASA NRA, the first analytical solution of the general equation was produced by members of the investigator team using the mathematics of fractional derivatives. This opened the way to an even more insightful and important investigation of the phenomena in microgravity.
The present paper describes the verification and validation of a quasi one-dimensional pressure based finite volume algorithm, implemented in Generalized Fluid System Simulation Program (GFSSP), for predicting compressible flow with friction, heat transfer and area change. The numerical predictions were compared with two classical solutions of compressible flow, i.e. Fanno and Rayleigh flow. Fanno flow provides an analytical solution of compressible flow in a long slender pipe where incoming subsonic flow can be choked due to friction. On the other hand, Raleigh flow provides analytical solution of frictionless compressible flow with heat transfer where incoming subsonic flow can be choked at the outlet boundary with heat addition to the control volume. Nonuniform grid distribution improves the accuracy of numerical prediction. A benchmark numerical solution of compressible flow in a converging-diverging nozzle with friction and heat transfer has been developed to verify GFSSP's numerical predictions. The numerical predictions compare favorably in all cases.
The role that distributed parameter model can play in CSI is demonstrated, in particular in combined structures controls optimization problems of importance in preliminary design. Closed form solutions can be obtained for performance criteria such as rms attitude error, making possible analytical solutions of the optimization problem. This is in contrast to the need for numerical computer solution involving the inversion of large matrices in traditional finite element model (FEM) use. Another advantage of the analytic solution is that it can provide much needed insight into phenomena that can otherwise be obscured or difficult to discern from numerical computer results. As a compromise in level of complexity between a toy lab model and a real space structure, the lattice truss used in the EPS (Earth Pointing Satellite) was chosen. The optimization problem chosen is a generic one: of minimizing the structure mass subject to a specified stability margin and to a specified upper bond on the rms attitude error, using a co-located controller and sensors. Standard FEM treating each bar as a truss element is used, while the continuum model is anisotropic Timoshenko beam model. Performance criteria are derived for each model, except that for the distributed parameter model, explicit closed form solutions was obtained. Numerical results obtained by the two model show complete agreement.
This report discusses the progress accomplished during the first year of a two-year contract to study Propulsive Re-Entry Aerodynamics. Analytical methods were developed during the first year to predict the flow field in the immediate vicinity of a planetary entry system composed of an aeroshell with a single retrorocket which exhausts into a subsonic or supersonic counterflowing planetary atmosphere. In the second year a digital computer program will be developed for the analytical solutions, and the analytical model will be tested by performing an experimental test program and comparing the experimental results with the analytical predictions. The analytical investigation considers the flow field due to a single supersonic retrorocket exhausting from an aeroshell into oncoming subsonic, transonic, or supersonic streams in the large thrusting coefficient regime, 0.5 < C(sub T) = T/q(sub ∞) A(sub m) < 15, which is the range of thrusting coefficients encountered during the terminal landing phase for planetary landing on Mars or Venus. Analyses are provided for all the geometrical features of the flow field such as the location of the terminal shock, jet boundary, and the profile for the interface between the jet and atmospheric gases. In addition, analyses are developed which describe the growth of the mixing layers along the interface and at the edge of the jet boundary and dead-air region and the reattachment of the shear layer to the aeroshell or the recompression of the shear layer in wake type flow. A comprehensive discussion of the experimental wind tunnel test program to be accomplished in Phase II is presented.
Laminar channel flow over a backward-facing step is investigated. The finite analytic (FA) method is used to obtain the numerical solution. The FA solutions predict the recirculation zone lengths and the recirculated mass flow rates for Reynolds numbers, Re, of 25, 50, 73, 125, 191 and 229 which correlate well with experimental measurements. The general flow patterns of the recirculation region flows for the Reynolds numbers considered in this study are similar to each other.
An approach using a two-dimensional thin aerofoil, allied with the theory of viscous bluff-body separation, is used to study the initial cross-over from massive separation to an attached flow in a single-row unstaggered cascade. Analytic solutions are developed for the limit of small cascade-spacing. From the analytic solutions several interesting features of the cascade are examined, including multiple-solution branches and multiple regions of hysteresis. In addition, numerical results are presented for several selected aerofoils. Some of the aerofoils are found to contain markedly enlarged regions of hysteresis for certain critical cascade spacings.
The problem of retrieving ligntning, ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered, It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (e.g., minimization). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A summary of a quasi-analytic extension of the spherical Earth solution to an oblate spheroid Earth geometry is also provided.