Search NASA⌕ Search

SEARCH · Search NASA

Results for “Closed form 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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

An analytic guidance technique for planetary and lunar approach trajectories.

An explicit, analytic guidance technique is developed for the hyperbolic approach phases of interplanetary and lunar spacecraft trajectories. The guidance technique is based upon a first-order analytic solution for the perturbed planet-centered (or moon-centered) trajectory. This trajectory is represented as the sum of two components: (1) the unperturbed osculating hyperbola at pericenter, and (2) first-order position and velocity perturbations due to gravitational effects of the sun and other planets. A closed-form analytic approximation for these perturbations valid for the entire approach trajectory is derived, thereby eliminating the need for numerical integration of the equations of motion. By means of this analytic trajectory model, the approach guidance problem is reduced to an equivalent two-body problem. The guidance objectives are specified in terms of actual, attainable conditions at pericenter, and the required corrective velocity is determined explicitly for both fixed and variable times of arrival.

Carlson, N. A.↗
NASA NTRS1974

Lectures on the scattering of light

The exact (Mie) theory for the scattering of a plane wave by a dielectric sphere is presented. Since this infinite series solution is computationally impractical for large spheres, another formulation is given in terms of an integral equation valid for a bounded, but otherwise general array of scatterers. This equation is applied to the scattering by a single sphere, and several methods are suggested for approximating the scattering cross section in closed form. A tensor scattering matrix is introduced, in terms of which some general scattering theorems are derived. The application of the formalism to multiple scattering is briefly considered.

Saxon, D. S.↗
NASA NTRS1989

Real-time approximate optimal guidance laws for the advanced launch system

An approach to optimal ascent guidance for a launch vehicle is developed using an expansion technique. The problem is to maximize the payload put into orbit subject to the equations of motion of a rocket over a rotating spherical earth. It is assumed that the thrust and gravitational forces dominate over the aerodynamic forces. It is shown that these forces can be separated by a small parameter epsilon, where epsilon is the ratio of the atmospheric scale height to the radius of the earth. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in a series where the zeroth-order term (epsilon = 0) can be obtained in closed form. The zeroth-order problem is that of putting maximum payload into orbit subject to the equations of motion of a rocket in a vacuum over a flat earth. The neglected inertial and aerodynamic terms are included in higher order terms of the expansion, which are determined from the solution of first-order linear partial differential equations requiring only quadrature integrations. These quadrature integrations can be performed rapidly, so that real-time approximate optimization can be used to construct the launch guidance law.

Speyer, Jason L.↗
NASA NTRS1951

Flow through cascades in tandem

An exact treatment of the problem of finding the incompressible, inviscid two-dimensional flow around two cascades in tandem is presented. The analysis includes solutions of both the direct and the inverse problems. These problems are solved by conformally mapping the tandem cascade onto the region between two concentric circles in which region there are suitably placed flow singularities. Formulas for the velocity and the potential in the annular region are presented in a closed form by means of elliptic functions. The equations are presented in a form suitable for computation.

Spraglin, William E↗
NASA NTRS1984

Effect of sidewall suction on flow in two-dimensional wind tunnels

A closed-form analysis of flow in a two-dimensional subsonic wind tunnel which uses sidewall suction around the model to reduce sidewall boundary-layer effects is presented. The model problem which is treated involves a flat plate airfoil in a tunnel with a suction window shaped to permit an analytic solution. This solution shows that the lift coefficient depends explicitly on the porosity parameter of the suction window and implicitly on the suction pressure differential. For a given sidewall displacement thickness, the lift coefficient increases as the suction-window porosity decreases.

Barnwell, R. W.↗
NASA NTRS1993

Effect of sidewall suction on flow in two-dimensional wind tunnels

A closed-form analysis of flow in a two-dimensional subsonic wind tunnel that uses sidewall suction around the model to reduce sidewall boundary-layer effects is presented. The model problem that is treated involves a flat plate airfoil in a tunnel with a suction window shaped to permit an analytic solution. This solution shows that the lift coefficient depends explicitly on the porosity parameter of the suction window and implicitly on the suction pressure differential. For a given sidewall displacement thickness, the lift coefficient increases as the suction-window porosity decreases.

Barnwell, Richard W.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
NASA NTRS1994

Unsteady panel method for flows with multiple bodies moving along various paths

A potential flow based three-dimensional panel method was modified to treat time-dependent conditions in which several submerged bodies can move within the fluid along different trajectories. This modification was accomplished by formulating the momentary solution in an inertial frame of reference, attached to the undisturbed stationary fluid. Consequently, the numerical interpretation of the multiple-body, solid-surface boundary condition and the viscous wake rollup was considerably simplified. The usteady capability of this code was calibrated and validated by comparing computed results with closed-form analytical results available for an airfoil, which was impulsively set into a constant speed forward motion. To demonstrate the multicomponent capability, computations were made for two wings following closely intersecting paths (i.e., simulations aimed at avoiding mid-air collisions) and for a flowfield with relative rotation (i.e., the case of a helicopter rotor rotating relative to the fuselage). Computed results for the cases were compared to experimental data, when such data was available.

Richardson, Thomas F.↗
NASA NTRSDate not supplied

Thermal Modelling of the in-Situ Consolidation of Automated Fiber Placement of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Thermoplastics↗
NASA NTRSDate not supplied

Thermal Modelling of the In-Situ Consolidation of Automated Fiber Placement Of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Christopher J. Stelter↗
DOE OSTI2023

Boson star normal modes

Boson stars are gravitationally bound objects that arise in ultralight dark matter models and form in the centers of galactic halos or axion miniclusters. We systematically study the excitations of a boson star, taking into account the mixing between positive and negative frequencies introduced by gravity. We show that the spectrum contains zero-energy modes in the monopole and dipole sectors resulting from spontaneous symmetry breaking by the boson star background. We analyze the general properties of the eigenmodes and derive their orthogonality and completeness conditions which have non-standard form due to the positive-negative frequency mixing. The eigenvalue problem is solved numerically for the first few energy levels in different multipole sectors and the results are compared to the solutions of the Schrödinger equation in fixed boson star gravitational potential. The two solutions differ significantly for the lowest modes, but get close for higher levels. We further confirm the normal mode spectrum in 3D wave simulations where we inject perturbations with different multipoles. As an application of the normal mode solutions, we compute the matrix element entering the evaporation rate of a boson star immersed in a hot axion gas. The computation combines the use of exact wavefunctions for the low-lying bound states and of the Schrödinger approximation for the high-energy excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2004

Containerless Measurement of Thermophysical Properties of Ti-Zr-Ni Alloys

The surface tension, viscosity, density, and thermal expansion of Ti-Zr-Ni alloys were measured for a number of compositions by electrostatic levitation methods. Containerless methods greatly reduce heterogeneous nucleation, increasing access to the undercooled liquid regime at finite cooling rates. The density and thermal expansion are measured optically, while the surface tension and viscosity are measured by the oscillating drop method. The measured alloys include compositions which form a metastable quasicrystal phase from the undercooled liquid, and alloys close to the composition of several multi-component bulk metallic glass-forming alloys. Measurements of surface tension show behavior typical of transition metals at high temperature, but a sudden decrease in the deeply undercooled liquid for alloys near the quasicrystal-forming composition range, but not for compositions which form the solid-solution phase first.

Hyers, Robert↗
DOE OSTI2022

Large-Area Periodic Arrays of Atomically Flat Single-Crystal Gold Nanotriangles Formed Directly on Substrate Surfaces

The advancement of nanoenabled wafer-based devices requires the establishment of core competencies related to the deterministic positioning of nanometric building blocks over large areas. Within this realm, plasmonic single-crystal gold nanotriangles represent one of the most attractive nanoscale components but where the formation of addressable arrays at scale has heretofore proven impracticable. In this work, a benchtop process is presented for the formation of large-area periodic arrays of gold nanotriangles. The devised growth pathway sees the formation of an array of defect-laden seeds using lithographic and vapor-phase assembly processes followed by their placement in a growth solution promoting planar growth and threefold symmetric side-faceting. The nanotriangles formed in this high-yield synthesis distinguish themselves in that they are epitaxially aligned with the underlying substrate, grown to thicknesses that are not readily obtainable in colloidal syntheses, and present atomically flat pristine surfaces exhibiting gold atoms with a close-packed structure. As such, they express crisp and unambiguous plasmonic modes and form photoactive surfaces with highly tunable and readily modeled plasmon resonances. The devised methods, hence, advance the integration of single-crystal gold nanotriangles into device platforms and provide an overall fabrication strategy that is adaptable to other nanomaterials.

42 ENGINEERING↗
DOE OSTI2025

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

Use of Polarization Lidar for Investigation of Meteorological Formations

This paper presents the results of theoretical and experimental, investigations of depolarization characteristics of different meteorological formations. Experimental investigations are carried out with a monostatic lidar. The ruby laser radiation is polarized in a vertical plane. The radiation reflected is accepted by a lens system of 150 rom in diameter and a viewing angle of 12' and further it is divided by Wollaston prism into the components polarized orthogonally. In this case the principal plane of the prism is exposed parallel with the laser polarization plane. Investigations show the degree of radiation polarization, reflected from water clouds, to be changed within 1/0.7 (seldom up to 0.6) depending on their density. In most cases a signal reflected from the cloud leading edge is polarized completely. The time shift is observed between polarized and crosspolarized components of a Fignal, reflected from a cloud, depending on the density of a meteorological object. While penetrating into the cloud depth a degree of polarization decreases up to 0.8-0.7, and the character of this decrease is different for various types of clouds. For crystal clouds the shift between the components of the reflected signal is not observed and the magnitude of polarization degree amounts to 0.1/0.3 in comparison with water clouds. The polarization degree of radiation reflected by fog is not less than 0.6, and that in the rains of average intensity (about 5 mm/h) is always about 1. The authors have suggested an algorithm of numerical solution of nonstationary transfer equation in the vector form to forecast the influence of multiple scattering effects on polarization characteristics of the lidar light signal. The method of statistical simulation (Monte-Carlo technique) forms the basis of the algorithm. Numerical estimates obtained for a model of stratocumulus at lambda = 0.6943 microns under boundary conditions close to the conditions of natural experiment being discussed proved to be in a good agreement with the results of observation. Specifically, Fig. 1 shows the profiles of polarization (p) versus depth (L) of the following drop formations: fog (curve 1) with horizontal meteorological visibility of 4 km two stratocumulus at a height of 1100 m with the attenuation factors delta = 0.01m(exp -1) (curve 2) and delta = 0.05m(exp -1) (curve 3). Curve 3 shows the results of numerical estimates and the value of their statistical error.

Balin, Yu. S.↗
DOE OSTI2020

Iodate interactions with calcite: implications for natural attenuation

Solid-phase interactions and speciation are important to radioiodine transport in groundwater. At the Hanford Site in Southeastern Washington State, iodate (IO 3 - ) is the main aqueous species in dilute radioiodine groundwater plumes. Like other oxyanions, IO 3 - may be incorporated into and/or adsorbed onto calcite, a common mineral at Hanford, decreasing its mobility in the environment. A series of macroscale batch experiments combined with solid phase characterization were conducted to identify variables impacting time-dependent aqueous IO 3 - removal via calcite precipitation and determine the location of IO 3 - within the calcite crystal structure. Results demonstrated 11.5-97% aqueous IO 3 - removal during initial rapid calcite precipitation. Incorporation was apparently the main removal mechanism, although later slower precipitation and/or adsorption may have also contributed to IO 3 - removal. Using a higher concentration of the calcite forming solutions (i.e., using 1M vs. 0.1M concentrations) resulted in an increase in the amount of precipitated calcite and a greater percentage of IO 3 - removed; however, calcite formed with lower molarity solutions resulted in higher IO 3 - mass (µg/g) removal. Solubility testing of laboratory produced calcites showed only small differences in solubility for calcite with and without IO 3 - incorporated into the structure. Evidence collected from SEM/FIB and TEM/SAED suggested that the IO 3 - incorporated into calcite was present in regions close to surface (implying potential easy release upon calcite dissolution).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Incremental Interval Assignment by Integer Linear Algebra with Improvements

Interval Assignment (IA) is the problem of selecting the number of mesh edges (intervals) for each curve for conforming quad and hex meshing. The intervals x is fundamentally integer-valued. Many other approaches perform numerical optimization then convert a floating-point solution into an integer solution, which is slow and error prone. We avoid such steps: we start integer, and stay integer. Incremental Interval Assignment (IIA) uses integer linear algebra (Hermite normal form) to find an initial solution to the meshing constraints, satisfying the integer matrix equation Solving for reduced row echelon form provides integer vectors spanning the nullspace of A. Here we add vectors from the nullspace to improve the initial solution, maintaining Ax = b Heuristics find good integer linear combinations of nullspace vectors that provide strict improvement towards variable bounds or goals. IIA always produces an integer solution if one exists. In practice we usually achieve solutions close to the user goals, but there is no guarantee that the solution is optimal, nor even satisfies variable bounds, e.g. has positive intervals. We describe several algorithmic changes since first publication that tend to improve the final solution. The software is freely available.

97 MATHEMATICS AND COMPUTING↗