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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2024

CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5

Abstract. This article presents the C-grid implementation of the CICE sea ice model, including the C-grid discretization of the momentum equation, the boundary conditions (BCs), and the modifications to the code required to use the incremental remapping transport scheme. To validate the new C-grid implementation, many numerical experiments were conducted and compared to the B-grid solutions. In idealized experiments, the standard advection method (incremental remapping with C-grid velocities interpolated to the cell corners) leads to a checkerboard pattern. A modal analysis demonstrates that this computational noise originates from the spatial averaging of C-grid velocities at corners. The checkerboard pattern can be eliminated by adjusting the departure regions to match the divergence obtained from the solution of the momentum equation. We refer to this novel approach as the edge flux adjustment (EFA) method. The C-grid discretization with edge flux adjustment allows for transport in channels that are one grid cell wide – a capability that is not possible with the B-grid discretization nor with the C-grid and standard remapping advection. Simulation results match the predicted values of a novel analytical solution for one-grid-cell-wide channels.

Lemieux, Jean-François (ORCID:0000000320845759)↗

NASA NTRS1976

A two-dimensional finite difference solution for the transient thermal behavior of tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique was analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

NASA NTRS2003

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

NASA NTRS1975

Improved solution for potential flow about arbitrary axisymmetric bodies by the use of a higher-order surface source method

An investigation is conducted of a case of axisymmetric bodies in which the application of main interest is an inlet, possibly with centerbody and ring vanes. The technique employed makes use of curved surface elements and a source density which varies over the element. Such an approach is designated a higher-order implementation. Questions of surface element geometry are discussed along with the computation of the induced velocity matrices and the organization of the calculation. The calculated results are compared with analytic solutions.

Hess, J. L.↗

NASA NTRS1980

A two-dimensional finite-difference solution for the transient thermal behavior of a tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique were analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗

NASA NTRS2013

Numerical Modeling of Ablation Heat Transfer

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Ewing, Mark E.↗

NASA NTRS1974

Diagnostics of seeded RF plasmas: An experimental study related to the gaseous core reactor

Measurements of the temperature profiles in an RF argon plasma were made over magnetic field intensities ranging from 20 amp turns/cm to 80 amp turns/cm. The results were compared with a one-dimensional numerical treatment of the governing equations and with an approximate closed form analytical solution that neglected radiation losses. The average measured temperatures in the plasma compared well with the numerical treatment, though the experimental profile showed less of an off center temperature peak than predicted by theory. This may be a result of the complex turbulent flow pattern present in the experimental torch and not modeled in the numerical treatment. The radiation term cannot be neglected for argon at the power levels investigated. The closed form analytical approximation that neglected radiation led to temperature predictions on the order of 1000 K to 2000 K higher than measured or predicted by the numerical treatment which considered radiation losses.

Thompson, S. D.↗

NASA NTRS1974

Solving multiple scattering problems in planetary atmospheres

Definitions are provided of the basic concepts occurring in the solution of multiple scattering problems involving planetary atmospheres and attention is given to aspects of problem characterization. Approaches are considered for finding the answer to a particular problem without the performance of detailed calculations. The characteristics of albedos are investigated, taking into account semiinfinite atmospheres and finite atmospheres. Questions of surface illumination are discussed along with aspects related to energy deposition in the atmosphere, intensity, and polarization. Precise numerical methods are examined and analytical solutions are presented.

Irvine, W. M.↗

NASA NTRS1982

Rime ice accretion and its effect on airfoil performance

A methodology was developed to predict the growth of rime ice, and the resulting aerodynamic penalty on unprotected, subcritical, airfoil surfaces. The system of equations governing the trajectory of a water droplet in the airfoil flowfield is developed and a numerical solution is obtained to predict the mass flux of super cooled water droplets freezing on impact. A rime ice shape is predicted. The effect of time on the ice growth is modeled by a time-stepping procedure where the flowfield and droplet mass flux are updated periodically through the ice accretion process. Two similarity parameters, the trajectory similarity parameter and accumulation parameter, are found to govern the accretion of rime ice. In addition, an analytical solution is presented for Langmuir's classical modified inertia parameter. The aerodynamic evaluation of the effect of the ice accretion on airfoil performance is determined using an existing airfoil analysis code with empirical corrections. The change in maximum lift coefficient is found from an analysis of the new iced airfoil shape. The drag correction needed due to the severe surface roughness is formulated from existing iced airfoil and rough airfoil data. A small scale wind tunnel test was conducted to determine the change in airfoil performance due to a simulated rime ice shape.

Bragg, M. B.↗

NASA NTRS2005

Application of an Extended Parabolic Equation to the Calculation of the Mean Field and the Transverse and Longitudinal Mutual Coherence Functions Within Atmospheric Turbulence

Solutions are derived for the generalized mutual coherence function (MCF), i.e., the second order moment, of a random wave field propagating through a random medium within the context of the extended parabolic equation. Here, "generalized" connotes the consideration of both the transverse as well as the longitudinal second order moments (with respect to the direction of propagation). Such solutions will afford a comparison between the results of the parabolic equation within the pararaxial approximation and those of the wide-angle extended theory. To this end, a statistical operator method is developed which gives a general equation for an arbitrary spatial statistical moment of the wave field. The generality of the operator method allows one to obtain an expression for the second order field moment in the direction longitudinal to the direction of propagation. Analytical solutions to these equations are derived for the Kolmogorov and Tatarskii spectra of atmospheric permittivity fluctuations within the Markov approximation.

Manning, Robert M.↗

NASA NTRS1983

A new explicit method for the numerical solution of parabolic differential equations

A new method is derived for solving parabolic partial differential equations arising in transient heat conduction or in boundary-layer flows. The method is based on a combination of the modified differential quadrature (MDQ) method with the rational Runge-Kutta time-integration scheme. It is fully explicit, requires no matrix inversion, and is stable for any time-step for the heat equations. Burgers equation and the one- and two-dimensional heat equations are solved to demonstrate the accuracy and efficiency of the proposed algorithm. The present method is found to be very accurate and efficient when results are compared with analytic solutions.

Satofuka, N.↗

DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗

DOE OSTI2024

Kubo-Anderson theory of polariton line shape

We apply the Kubo-Anderson stochastic theory of molecular spectral line shapes to the case of polaritons formed in the collective strong-coupling regime. We investigate both the fast and slow limits of the random frequency modulation of the emitter as well as the intermediate regime and show how the interplay between the characteristic timescales of the cavity and the molecular disorder is expressed in the observed polaritons line shapes. The analytical solution obtained for the slow limit is valid for any ratio between the inhomogeneous broadening of the molecules and the Rabi splitting, which is especially relevant for molecular polaritons where these two quantities can be of the same order of magnitude.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1971

Comparison of a discrete steepest ascent method with the continuous steepest ascent method for optimal programing

A discrete steepest ascent method which allows controls which are not piecewise constant (for example, it allows all continuous piecewise linear controls) was derived for the solution of optimal programming problems. This method is based on the continuous steepest ascent method of Bryson and Denham and new concepts introduced by Kelley and Denham in their development of compatible adjoints for taking into account the effects of numerical integration. The method is a generalization of the algorithm suggested by Canon, Cullum, and Polak with the details of the gradient computation given. The discrete method was compared with the continuous method for an aerodynamics problem for which an analytic solution is given by Pontryagin's maximum principle, and numerical results are presented. The discrete method converges more rapidly than the continuous method at first, but then for some undetermined reason, loses its exponential convergence rate. A comparsion was also made for the algorithm of Canon, Cullum, and Polak using piecewise constant controls. This algorithm is very competitive with the continuous algorithm.

Childs, A. G.↗

NASA NTRS1990

Comparison of solution approaches for minimum-fuel, low-thrust, power-limited orbital transfers

An initial assessment of the feasibility of a function space gradient method for computing solutions to minimum-fuel power-limited transfers encompassing a wide range of thrust to weight ratios is conducted. Three transfers between coplanar ellipses are used as test cases. The gradient method performs best at the high end of the thrust to weight ratio range. At the lower end, there is reduced sensitivity of the fuel consumption to the control profiles. The minimum fuel consumption and the trajectory are computed quite accurately but the control profiles are in error. An approximate analytical solution, obtained by Edelbaum using the method of averaging, is discussed.

Mease, Kenneth D.↗

NASA NTRS1993

Analytical stress intensity solution for the Stable Poisson Loaded specimen

An analytical calibration of the Stable Poisson Loaded (SPL) specimen is presented. The specimen configuration is similar to the ASTM E-561 compact-tension specimen with displacement controlled wedge loading used for R-curve determination. The crack mouth opening displacements (CMODs) are produced by the diametral expansion of an axially compressed cylindrical pin located in the wake of a machined notch. Due to the unusual loading configuration, a three-dimensional finite element analysis was performed with gap elements simulating the contact between the pin and specimen. In this report, stress intensity factors, CMODs, and crack displacement profiles, are reported for different crack lengths and different contacting conditions. It was concluded that the computed stress intensity factor decreases sharply with increasing crack length thus making the SPL specimen configuration attractive for fracture testing of brittle, high modulus materials.

Ghosn, Louis J.↗

NASA NTRS1993

Thermocapillary convection in two immiscible liquid layers with free surface

Thermocapillary convection is studied in two immiscible liquid layers with one free surface, one liquid/liquid interface, and differential heating applied parallel to the interfaces. An analytical solution is introduced for infinite horizontal layers. The defining parameter for the flow pattern is lambda, the ratio of the temperature coefficient of the interfacial tension to that of the surface tension. Four different flow patterns exist under zero gravity conditions. 'Halt' conditions which halt the fluid motion in the lower encapsulated liquid layer have been found. A numerical experiment is carried out to study effects of vertical end walls on the double layer convection in a 2D cavity. The halt condition obtained from the analytical study is found to be valid in the limit of small Reynolds numbers. The flow in the encapsulated liquid layer can be suppressed substantially.

Doi, Takao↗