Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytical and numerical 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.

At least 19 records

Finite-analytic numerical solution of heat transfer in two-dimensional cavity flow

Heat transfer in cavity flow is numerically analyzed by a new numerical method called the finite-analytic method. The basic idea of the finite-analytic method is the incorporation of local analytic solutions in the numerical solutions of linear or nonlinear partial differential equations. In the present investigation, the local analytic solutions for temperature, stream function, and vorticity distributions are derived. When the local analytic solution is evaluated at a given nodal point, it gives an algebraic relationship between a nodal value in a subregion and its neighboring nodal points. A system of algebraic equations is solved to provide the numerical solution of the problem. The finite-analytic method is used to solve heat transfer in the cavity flow at high Reynolds number (1000) for Prandtl numbers of 0.1, 1, and 10.

Chen, C.-J.↗

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗

Finite analytic numerical solution of two-dimensional channel flow over a backward-facing step

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.

Ho, K.-S.↗

Finite analytic numerical solution of heat transfer and flow past a square channel cavity

A numerical solution of flow and heat transfer characteristics is obtained by the finite analytic method for a two dimensional laminar channel flow over a two-dimensional square cavity. The finite analytic method utilizes the local analytic solution in a small element of the problem region to form the algebraic equation relating an interior nodal value with its surrounding nodal values. Stable and rapidly converged solutions were obtained for Reynolds numbers ranging to 1000 and Prandtl number to 10. Streamfunction, vorticity and temperature profiles are solved. Local and mean Nusselt number are given. It is found that the separation streamlines between the cavity and channel flow are concave into the cavity at low Reynolds number and convex at high Reynolds number (Re greater than 100) and for square cavity the mean Nusselt number may be approximately correlated with Peclet number as Nu(m) = 0.365 Pe exp 0.2.

Chen, C.-J.↗

Analytic numerical solutions for shock waves

Study of weak solutions of simple wave equation, inviscid Burgers equation, and Euler equations has resulted in technique for accurate prediction of shock waves occurring in inviscid supersonic flows.

Maccormack, R. W.↗

Radiative transfer in highly scattering materials - numerical solution and evaluation of approximate analytic solutions

Numerical solutions for radiative transport in a class of anisotropically scattering materials are presented. Conditions for convergence and divergence of the iterative method are given and supported by computed results. The relation of two flux theories to the equation of radiative transfer for isotropic scattering is discussed. The adequacy of the two flux approach for the reflectance, radiative flux and radiative flux divergence of highly scattering media is evaluated with respect to solutions of the radiative transfer equation.

Weston, K. C.↗

Analytical and Numerical Solution for a Solidifying Liquid Alloy Slab

Numerical and analytical solutions are presented for the temperature and concentration distributions during the solidification of a binary liquid alloy slab. The slab is taken to be of a finite depth but infinite in the horizontal direction. The solidification process is started by withdrawing a fixed amount of heat from the lower surface of the slab. The upper surface of the slab is subjected to both radiation and convective conditions. The solution gives the concentration and temperature profiles and the interface position as a function of time. Due to the smallness of the mass diffusion coefficient in the solid, the numerical solution method breaks down whenever the ratio of the diffusivities in the solid and the liquid falls below a certain value. An analytical method is developed which gives accurate solution for any value of the diffusivity ratio.

Antar, B. N.↗

The Finite Analytic Method for steady and unsteady heat transfer problems

A new numerical method called the Finite Analytical Method for solving partial differential equations is introduced. The basic idea of the finite analytic method is the incorporation of the local analytic solution in obtaining the numerical solution of the problem. The finite analytical method first divides the total region of the problem into small subregions in which local analytic solutions are obtained. Then an algebraic equation is derived from the local analytic solution for each subregion relating an interior nodal value at a point P in the subregion to its neighboring nodal values. The assembly of all the local analytic solutions thus provides the finite-analytic numerical solution of the problem. In this paper the finite analytic method is illustrated in solving steady and unsteady heat transfer problems.

Chen, C.-J.↗

Albedos of homogeneous semi-infinite canopies - Comparison of two-stream analytic and numerical solutions

The albedo of plant canopies is treated as a problem in radiative transfer. Albedos calcualted from an iterative multistream numerical model are compared with those calculated with an analytic two-stream solution. With the assumption of a randomly homogeneous distribution of leaf positions and orientations and isotropic scattering by individual leaves, the single-scattering albedo of the canopy can be found analytically. This single-scattering solution is incorporated into the two-stream solution and used to benchmark the multistream numerical model in the single-scattering limit. Relative errors so established in the multistream model are O(0.3 percent) or less. The two-stream model is also found to be remarkably accurate, with the error in multiply scattered radiation O(5 percent) or less, corresponding to absolute errors in visible albedo of less than 0.001 and near-infrared albedo of less than or equal to 0.01. Thus the two-stream model should be adequate for many purposes, such as climate modeling, provided the assumptions of homogeneous canopy and isotropic scattering are not too unrealistic.

Dickinson, Robert E.↗

A comparison of numerical and analytical solution of the creeping flame spread over thermally thin material

The present numerical solution for the de Ris (1969) problem of flame-spread over thin condensed fuel in an opposed-flow environment is obtained by reformulating the problem in terms of four nondimensional parameters, while retaining all assumptions of the original theory. While the de Ris theory sees the location of the leading edge and the eigenlocation of the onset of evaporation as identical, this analysis treats the leading edge as part of the solution; the location of the flame leading edge is in this way established to be upstream of the eigenlocation, with significant consequences for the spread rate formula.

Bhattacharjee, Subrata↗

Numerical and Analytical Solutions of Hypersonic Interactions Involving Surface Property Discontinuities

The local viscous-inviscid interaction field generated by a wall temperature jump on a flat plate in supersonic flow and on the windside of a Reusable Launch Vehicle in hypersonic flow is studied in detail by both a Navier-Stokes numerical code and an analytical triple-deck model. Treatment of the rapid heat transfer changes both upstream and downstream of the jump is included. Closed form relationships derived from the triple-deck theory are presented. The analytically predicted pressure and heating variations including upstream influence are found to be in generally good agreement with the Computational Fluid Dynamic (CFD) predictions. These analyses not only clarify the interactive physics involved but also are useful in preliminary design of thermal protection systems and as an insertable module to improve CFD code efficiency when applied to such small-scale interaction problems. The analyses only require conditions at the wall and boundary-layer edge which are easily extracted from a baseline, constant wall temperature, CFD solution.

Gnoffo, Peter A.↗

Satellite attitude motion models for capture and retrieval investigations

The primary purpose of this research is to provide mathematical models which may be used in the investigation of various aspects of the remote capture and retrieval of uncontrolled satellites. Emphasis has been placed on analytical models; however, to verify analytical solutions, numerical integration must be used. Also, for satellites of certain types, numerical integration may be the only practical or perhaps the only possible method of solution. First, to provide a basis for analytical and numerical work, uncontrolled satellites were categorized using criteria based on: (1) orbital motions, (2) external angular momenta, (3) internal angular momenta, (4) physical characteristics, and (5) the stability of their equilibrium states. Several analytical solutions for the attitude motions of satellite models were compiled, checked, corrected in some minor respects and their short-term prediction capabilities were investigated. Single-rigid-body, dual-spin and multi-rotor configurations are treated. To verify the analytical models and to see how the true motion of a satellite which is acted upon by environmental torques differs from its corresponding torque-free motion, a numerical simulation code was developed. This code contains a relatively general satellite model and models for gravity-gradient and aerodynamic torques. The spacecraft physical model for the code and the equations of motion are given. The two environmental torque models are described.

Cochran, John E., Jr.↗

The stability of motion of satellites with cavities partially filled with liquid

The stability and time dependent motion of a spinning satellite, simulated by a rigid body with a cavity partially filled with liquid is examined. The problem formulation, consisting of the boundary-value problem for the liquid and moment equations for the entire system is presented. Because of large Reynold's numbers involved, viscosity effects are negligible everywhere except for a thin boundary layer near the wetted surface. Using a boundary-layer analysis, the effect of the boundary layer is replaced by modified boundary conditions for the liquid. The solution of the differential equations for the inviscid problem is solved in closed form. A semi-analytical numerical solution of the inviscid equations subject to the viscous boundary condition has proved unsucessful.

Nayfeh, A. H.↗

Gasdynamics of explosions today.

A brief review is given of blast and detonation wave phenomena and some of their uses in war and peace. It is concluded that great strides have been made over the last three decades toward the physical understanding, the analytical-numerical solution, and the measurement of dynamic and thermodynamic quantities, also taking into consideration severe environments and extremely short durations. Questions of internal ballistics are discussed together with hypervelocity launchers and shock tubes, collapsing cylindrical drivers, spherical implosions, explosive weapons, dynamic response, and equation of state data.

Brode, H. L.↗

Conflict Probability Estimation for Free Flight

The safety and efficiency of free flight will benefit from automated conflict prediction and resolution advisories. Conflict prediction is based on trajectory prediction and is less certain the farther in advance the prediction, however. An estimate is therefore needed of the probability that a conflict will occur, given a pair of predicted trajectories and their levels of uncertainty. A method is developed in this paper to estimate that conflict probability. The trajectory prediction errors are modeled as normally distributed, and the two error covariances for an aircraft pair are combined into a single equivalent covariance of the relative position. A coordinate transformation is then used to derive an analytical solution. Numerical examples and Monte Carlo validation are presented.

Paielli, Russell A.↗

Three-dimensional convection in horizontal cylinders - Numerical solutions and comparison with experimental and analytical results

The results of a three-dimensional numerical simulation of Boussinesq free convection in a horizontal differentially heated cylinder are presented. The computation was based on a Samarskii-Andreyev scheme (described by Leong, 1981) and a false-transient advancement in time, with vorticity, velocity, and temperature as dependent variables. Solutions for velocity and temperature distributions were obtained for Rayleigh numbers (based on the radius) Ra = 74-18,700, thus covering the core- and boundary-layer-driven regimes. Numerical solutions are compared with asymptotic analytical solutions and experimental data. The numerical results well represent the complex three-dimensional flows found experimentally.

Smutek, C.↗

Dynamic Beam Solutions for Real-Time Simulation and Control Development of Flexible Rockets

In this study, flexible rockets are structurally represented by linear beams. Both direct and indirect solutions of beam dynamic equations are sought to facilitate real-time simulation and control development for flexible rockets. The direct solution is completed by numerically integrate the beam structural dynamic equation using an explicit Newmark-based scheme, which allows for stable and fast transient solutions to the dynamics of flexile rockets. Furthermore, in the real-time operation, the bending strain of the beam is measured by fiber optical sensors (FOS) at intermittent locations along the span, while both angular velocity and translational acceleration are measured at a single point by the inertial measurement unit (IMU). Another study in this paper is to find the analytical and numerical solutions of the beam dynamics based on the limited measurement data to facilitate the real-time control development. Numerical studies demonstrate the accuracy of these real-time solutions to the beam dynamics. Such analytical and numerical solutions, when integrated with data processing and control algorithms and mechanisms, have the potential to increase launch availability by processing flight data into the flexible launch vehicle's control system.

flexible body dyanmics fiber optic strain sensor p↗