Engineering Papers⌕ Search

SEARCH · Engineering Papers

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

At least 541 records · Page 30

The heating of a thermally conducting stratified medium. II - A simple plane model of an atmosphere

Exact solutions of the following theoretical problem are presented: A plane atmosphere is in hydrostatic equilibrium with a uniform gravity. The ideal gas law is assumed. Heat is generated everywhere at a rate proportional to the local density. The atmosphere is maintained in a steady state through cooling by thermal conduction and radiation. This problem is reducible to quadratures for a thermal conductivity which is an arbitrary, but prescribed, function of the temperature, and for a radiative loss which is expressible as the product of the density and an arbitrary, but prescribed, function of the pressure. The analysis is carried out for the case of power law thermal conductivity, and a radiative loss proportional to the square of the density and to the first power of the temperature. The radiative cooling function adopted here has the basic mathematical form for an optically thin medium. The solutions reproduce the macroscopic ordering of a hot 'corona' separated from a 'photosphere' by a layer of temperature minimum. The analytic solutions allow direct illustration of the interplay between steady energy transport and the requirements of hydrostatic equilibrium.

Lerche, I.↗

Vacillations induced by interference of stationary and traveling planetary waves

The interference pattern produced when a traveling planetary wave propagates over a stationary forced wave is explored, examining the interference signature in a variety of diagnostics. The wave field is first restricted to a diatomic spectrum consisting of two components: a single stationary wave and a single monochromatic traveling wave. A simple barotropic normal mode propagating over a simple stationary plane wave is considered, and closed form solutions are obtained. The wave fields are then restricted spatially, providing more realistic structures without sacrificing the advantages of an analytical solution. Both stationary and traveling wave fields are calculated numerically with the linearized Primitive Equations in a realistic basic state. The mean flow reaction to the fluctuating eddy forcing which results from interference is derived. Synoptic geopotential behavior corresponding to the combined wave and mean flow fields is presented, and the synoptic signature in potential vorticity on isentropic surfaces is examined.

Salby, Murry L.↗

Adaptive Shape Functions and Internal Mesh Adaptation for Modelling Progressive Failure in Adhesively Bonded Joints

Enhanced finite elements are elements with an embedded analytical solution that can capture detailed local fields, enabling more efficient, mesh independent finite element analysis. The shape functions are determined based on the analytical model rather than prescribed. This method was applied to adhesively bonded joints to model joint behavior with one element through the thickness. This study demonstrates two methods of maintaining the fidelity of such elements during adhesive non-linearity and cracking without increasing the mesh needed for an accurate solution. The first method uses adaptive shape functions, where the shape functions are recalculated at each load step based on the softening of the adhesive. The second method is internal mesh adaption, where cracking of the adhesive within an element is captured by further discretizing the element internally to represent the partially cracked geometry. By keeping mesh adaptations within an element, a finer mesh can be used during the analysis without affecting the global finite element model mesh. Examples are shown which highlight when each method is most effective in reducing the number of elements needed to capture adhesive nonlinearity and cracking. These methods are validated against analogous finite element models utilizing cohesive zone elements.

computational mechanics↗

Study of recirculating liquid fuel in a 1D critical stationary system

Several studies have been conducted to investigate the physics of liquid fuel reactors, showing also applications with molten salts. A liquid nuclear fuel implies changes in the neutron balance equation to take into account the precursors' displacement and the emission of delayed neutrons in a different position than at the original fission. This requires to upgrade the computer codes normally used to calculate nuclear reactors using only solid fuel. In this work, we revisit a simple problem with liquid fuel, which is proposed for the verification of the numerical solutions obtained by advanced computer codes. This problem studies criticality with constant coefficients, thus neglecting thermal feedback. We elaborate on the analytical solution of the problem, deriving also a generalized eigenvalue problem by finite-volume integration over the cells of a discretized mesh to study the evolution of the dominance ratio with fuel velocity. Finally, we investigate the influence of the fuel velocity on the reactivity of the system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Derivation of a three dimensional numerical water quality model for estuary and continental shelf application

A derivation is given for a three dimensional mass transport equation which is appropriate for numerical modeling of estuary and continental shelf water quality variations for both the time dependent and steady state cases. A finite difference approximation to the derived equation is presented and a solution scheme for the resulting equations outlined. Preliminary results are obtained using the model for the extremely simple problems which have analytical solutions. The numerical model, as presented, will provide a scheme to study water quality problems in coastal waters for both steady state and time dependent cases.

Spaulding, M.↗

Bidirectional reflectance spectroscopy. I - Theory

An approximate analytic solution is derived for the radiative transfer equation describing particulate surface light scattering, taking into account multiple scattering and mutual shadowing. Analytical expressions for the following quantities are found: bidirectional reflectance, radiance coefficient and factor, the normal, Bond, hemispherical, and physical albedos, integral phase function and phase integral, and limb-darkening profile. Scattering functions for mixtures can be calculated, as well as corrections for comparisons of experimental transmission or reflection spectra with observational planetary spectra. The theory should be useful for the interpretation of reflectance spectroscopy of laboratory surfaces and the photometry of solar system objects.

Hapke, B.↗

Optimal two-impulse transfer between specified terminal states of Keplerian orbits

The integrals of motion for an inverse-square force field and a necessary condition for optimality are used in a transfer orbit co-ordinate system to formulate the solution of the optimal two-impulse transfer between fixed position and velocity vectors on Keplerian orbits. In this co-ordinate system, the equations reveal two asymptotes that are useful in analysing the solution. It is shown that there are only two real extremals (both minimums) which are separated by an asymptote. An approximate analytic solution is also obtained by ignoring a quadratic term whose coefficient is approximately zero for a large class of orbit transfer problems.

Jezewski, D. J.↗

Confinement of thermocapillary floating zone flow by uniform rotation

Smith and Greenspan (1984) examined theoretically the idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow in crystal growth experiments to the melt sidewall, leaving the interior of the melt passive. Here, that model is extended to a full zone with a more realistic temperature distribution imposed on the sidewall, and both linear and nonlinear thermocapillary flows are theoretically studied. Linearized, analytical solutions are found using singular perturbation theory and the various sidewall boundary layers described by Greenspan (1969) for rotating fluids. The analytical and linearized numerical results are compared, and the linear and nonlinear flows are discussed. The results demonstrate that the thermocapillary flow is strong and that rotation cannot confine the flow. Temperature advection by strongly nonlinear flow is significant even for the small Prandtl number of silicon.

Fowlis, W. W.↗

Application of finite-element-based solution technologies for viscoplastic structural analyses

Finite-element solution technology developed for use in conjunction with advanced viscoplastic models is described. The development of such solution technology is necessary for performing stress/life analyses of engineering structural problems where the complex geometries and loadings make the conventional analytical solutions difficult. The versatility of the solution technology is demonstrated by applying it to viscoplastic models possessing different mathematical structures and encompassing isotropic and anisotropic material. The computational results qualitatively replicate deformation behavior observed in experiments on prototypical structural components.

Arya, V. K.↗

Application of finite-element-based solution technologies for viscoplastic structural analyses

Finite-element solution technology developed for use in conjunction with advanced viscoplastic models is described. The development of such solution technology is necessary for performing stress/life analyses of engineering structural problems where the complex geometries and loadings make the conventional analytical solutions difficult. The versatility of the solution technology is demonstrated by applying it to viscoplastic models possessing different mathematical structures and encompassing isotropic and anisotropic materials. The computational results qualitatively replicate deformation behavior observed in experiments on prototypical structural components.

Arya, Vinod K.↗

Torque Transient of Magnetically Drive Flow for Viscosity Measurement

Viscosity is a good indicator of structural changes for complex liquids, such as semiconductor melts with chain or ring structures. This paper discusses the theoretical and experimental results of the transient torque technique for non-intrusive viscosity measurement. Such a technique is essential for the high temperature viscosity measurement of high pressure and toxic semiconductor melts. In this paper, our previous work on oscillating cup technique was expanded to the transient process of a magnetically driven melt flow in a damped oscillation system. Based on the analytical solution for the fluid flow and cup oscillation, a semi-empirical model was established to extract the fluid viscosity. The analytical and experimental results indicated that such a technique has the advantage of short measurement time and straight forward data analysis procedures

Ban, Heng↗

Solving time-dependent two-dimensional eddy current problems

Results of transient eddy current calculations are reported. For simplicity, a two-dimensional transverse magnetic field which is incident on an infinitely long conductor is considered. The conductor is assumed to be a good but not perfect conductor. The resulting problem is an interface initial boundary value problem with the boundary of the conductor being the interface. A finite difference method is used to march the solution explicitly in time. The method is shown. Treatment of appropriate radiation conditions is given special consideration. Results are validated with approximate analytic solutions. Two stringent test cases of high and low frequency incident waves are considered to validate the results.

Lee, Min Eig↗

Preconditioned upwind methods to solve 3-D incompressible Navier-Stokes equations for viscous flows

A computational method for calculating low-speed viscous flowfields is developed. The method uses the implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. An extensive study of optimizing the preconditioned system is carried out for two viscous flow problems. Computed results are compared with analytical solutions and experimental data.

Hsu, C.-H.↗

A Markov model reduction technique for fault tolerant processor reliability analysis

A fault tolerant processor (FTP) plays a key role in many high performance, safety-critical control system applications. Realistic modeling of an FTP is crucial to gaining a high degree of confidence in the reliability and safety analysis of such a system. While fidelity is clearly a major consideration, a practical model must also be kept to a moderate size to allow its incorporation into the overall system model. This paper presents a systematic reduction technique that starts from a complex, detailed model of a triple, redundant FTP and produces a low order approximation of very high accuracy. The existence of two distinct time scales represents the key to the success of the technique. No eigenvalue solution or coordinate transformation are needed. The reduced model captures all the important features of the detailed model, is amenable to an analytical solution and provides insight into the reconfiguration behavior of an FTP.

Schor, Andrei L.↗

The double universal joint wrist on a manipulator: Solution of inverse position kinematics and singularity analysis

This paper presents three methods to solve the inverse position kinematics position problem of the double universal joint attached to a manipulator: (1) an analytical solution for two specific cases; (2) an approximate closed form solution based on ignoring the wrist offset; and (3) an iterative method which repeats closed form position and orientation calculations until the solution is achieved. Several manipulators are used to demonstrate the solution methods: cartesian, cylindrical, spherical, and an anthropomorphic articulated arm, based on the Flight Telerobotic Servicer (FTS) arm. A singularity analysis is presented for the double universal joint wrist attached to the above manipulator arms. While the double universal joint wrist standing alone is singularity-free in orientation, the singularity analysis indicates the presence of coupled position/orientation singularities of the spherical and articulated manipulators with the wrist. The cartesian and cylindrical manipulators with the double universal joint wrist were found to be singularity-free. The methods of this paper can be implemented in a real-time controller for manipulators with the double universal joint wrist. Such mechanically dextrous systems could be used in telerobotic and industrial applications, but further work is required to avoid the singularities.

Williams, Robert L., III↗

A uniformly valid solution for motion about the interior libration point of the perturbed elliptic-restricted problem

Bounded motion about Lagrange collinear libration points is considered for a perturbed elliptic-restricted problem. A practical application is the motion of a satellite near a libration point collinear with the sun and the earth-moon barycenter. Such a system is treated here as an earth-sun-satellite elliptic restricted problem with lunar perturbations. The method of dual time scales is utilized to develop a uniformly valid three-dimensional analytical solution to the satellite's equations of motion. The analytical development applies somewhat generally to that class of four-body problems where the second primary mass is much greater than the first, and the third primary mass is much greater than the second.

Richardson, D. L.↗

Verification and Validation Activities for the Multi-Fidelity Toolkit

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

42 ENGINEERING↗

Diffusion-convection function of cosmic rays

The fundamental properties and some numerical results of the solution of the diffusion equation of an impulsive cosmic-ray point source in an uniform, unbounded and spherically symmetrical moving medium is presented. The diffusion-convection(D-C) function is an elementary composite function of the solution of the D-C equation for the particles injected impulsively from a diffusive point source into the medium. It is the analytic solution derived by the dimensional method for the propagation equation of solar cosmic rays in the heliosphere, i.e. the interplanetary space. Because of the introduction of convection effect of solar wind, a nonhomogeneous term appears in the propagation equation, it is difficult to express its solution in terms of the ordinary special functions. The research made so far has led to a solution containing only the first order approximation of the convection effect.

Zhang, G.↗