Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytic 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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Small body mass determination from a flyby

This paper investigates the accuracy with which the mass of a small body can be determined from a flyby. Due to the nearly linear nature of hyperbolic flybys of small bodies, analytic solutions can be derived. Four analytic solutions are developed, each of increasing complexity with different applications. These solutions are compared to a complex computer simulation for accuracy. When these methods are applied to nearly linear trajectories, the analytic solutions agree with the computer simulation results to within 4.0 percent except in one boundary case which is discussed in the text. Also described is the linear range for which these equations are applicable. These solutions are then applied to the proposed asteroid flyby of the CRAF mission to predict the mass estimate accuracy.

Sybert, Catherine B.↗
DOE OSTI2013

Leptogenesis from first principles in the resonant regime

The lepton asymmetry generated by the out-of-equilibrium decays of heavy Majorana neutrinos with a quasi-degenerate mass spectrum is resonantly enhanced. In this work, we study this scenario within a first-principle approach. The quantum field theoretical treatment is applicable for mass splittings of the order of the width of the Majorana neutrinos, for which the enhancement is maximally large. The non-equilibrium evolution of the mixing Majorana neutrino fields is described by a formal analytical solution of the Kadanoff-Baym equations, that is obtained by neglecting the back-reaction. Based on this solution, we derive approximate analytical expressions for the generated asymmetry and compare them to the Boltzmann result. We find that the resonant enhancement obtained from the Kadanoff-Baym approach is smaller compared to the Boltzmann approach, due to additional contributions that describe coherent transitions between the Majorana neutrino species. We also discuss corrections to the masses and widths of the degenerate pair of Majorana neutrinos that are relevant for very small mass splitting, and compare the approximate analytical result for the lepton asymmetry with numerical results. - Highlights: Black-Right-Pointing-Pointer Treatment of leptogenesis in the resonant regime by a first-principle approach. Black-Right-Pointing-Pointer Applicable for mass differences of the order of the width. Black-Right-Pointing-Pointer Non-equilibrium evolution described by Kadanoff-Baym equations. Black-Right-Pointing-Pointer Analytical solution can be obtained neglecting backreaction. Black-Right-Pointing-Pointer Resonant enhancement is smaller compared to the classical Boltzmann approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Study on 3D Internal Magnetic Field Distribution and Dynamic Mechanics of a Giant Magnetostrictive Actuator

The three-dimensional (3D) coupled magneto-elastic dynamic model of a giant magnetostrictive actuator (GMA) was built based on Maxwell’s equations, elasticity, and vibration theories. The general analytical expressions for the 3D internal magnetic field distribution of GMA were derived, and the corresponding analytical solutions were obtained. This analytical model is able to predict the internal magnetic field distribution and the output displacement of GMA system with different physical parameters under external excitations. The dynamic behaviors, such as the amplitude vs. frequency responses and the time domain waveforms of the GMA, were also studied. Analytical solutions showed that the internal magnetic field has a prominent end effect along the axial direction of a giant magnetostrictive material (GMM) rod, whereas a minimal skin effect along its radial direction, which are matched very well with the FEM results. The amplitude-frequency relationship exhibited prominent resonance and anti-resonance characteristic, which is a typical hysteresis feature of the GMA system. In order to testify the validity of the model prediction, a GMA testing system was also established. Excellent agreement was found between the experiment measurements and model predictions in the range of our research (frequency is lower than 300 Hz). The analytical results of the present study have potential applications in the fields of active vibration control and other precision machining of the GMA system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Kinematic Design of Serial Link Manipulators From Task Specifications

The Reconfigurable Modular Manipulator System (RMMS) consists of modular links and joints that can be assembled into many manipulator configurations. This capability allows the RMMS to be rapidly reconfigured to custom tailor it to specific tasks. An important issue related to the RMMS is the determination of the optimal manipulator configuration for a specific task. This article addresses the problem of mapping kinematic task specifications into a kinematic manipulator configuration. For the design of two-degrees-of-freedom (2- DOF) planar manipulators, an analytical solution is derived. Because analytical solutions become impractical for problems with more than two design parameters, we have also developed a numerical approach for the design of 6-DOF manipulators. The numerical procedure determines the Denavit-Hartenberg (D-H) parameters of a nonredundant manipulator with joint limits that can reach a set of specified positions/orientations in an environment that may include parallelepiped-shaped obstacles. Finally, this approach is demonstrated with a three- dimensional example for a 6-DOF manipulator.

42 ENGINEERING↗
DOE OSTI2016

An added-mass partition algorithm for fluid–structure interactions of compressible fluids and nonlinear solids

We describe an added-mass partitioned (AMP) algorithm for solving fluid–structure interaction (FSI) problems involving inviscid compressible fluids interacting with nonlinear solids that undergo large rotations and displacements. The computational approach is a mixed Eulerian–Lagrangian scheme that makes use of deforming composite grids (DCG) to treat large changes in the geometry in an accurate, flexible, and robust manner. The current work extends the AMP algorithm developed in Banks et al. [1] for linearly elasticity to the case of nonlinear solids. To ensure stability for the case of light solids, the new AMP algorithm embeds an approximate solution of a nonlinear fluid–solid Riemann (FSR) problem into the interface treatment. The solution to the FSR problem is derived and shown to be of a similar form to that derived for linear solids: the state on the interface being fundamentally an impedance-weighted average of the fluid and solid states. Numerical simulations demonstrate that the AMP algorithm is stable even for light solids when added-mass effects are large. The accuracy and stability of the AMP scheme is verified by comparison to an exact solution using the method of analytical solutions and to a semi-analytical solution that is obtained for a rotating solid disk immersed in a fluid. The scheme is applied to the simulation of a planar shock impacting a light elliptical-shaped solid, and comparisons are made between solutions of the FSI problem for a neo-Hookean solid, a linearly elastic solid, and a rigid solid. The ability of the approach to handle large deformations is demonstrated for a problem of a high-speed flow past a light, thin, and flexible solid beam.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Smoothed particle hydrodynamics model for Landau-Lifshitz-Navier-Stokes and advection-diffusion equations

We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1981

Comments on the Leimanis solution of self-excited rigid body

An attempt is made to apply the solutions contained in a book by Leimanis (1965), which includes the Bodewadt solution, to the performance assessment of the Galileo spacecraft during spin up and spin down maneuvers. Because the Galileo is not precisely symmetric, however, the Bodewadt solution fails to achieve the desired accuracy for useful analysis. It is noted that in the case of the Leimanis solution, the analytic results are incorrect for theoretical reasons. In order to remedy this situation, analytic solutions are developed which satisfy the criterion of high accuracy and provide a useful analytic tool for the performance assessment and maneuver analysis of the Galileo spacecraft and many similar near symmetric spinning spacecraft. Simulation results suggest the relative accuracies of the Bodewadt solution and the author's solution, while the restricted regions of validity of the Leimanis solution are clearly indicated.

Longuski, J. M.↗
DOE OSTI2023

Semi-Analytical Leakage Solutions for Aquifers (SALSA) v1

SALSA is a semi-analytical modeling tool that can provide assessments of pressure perturbations and brine leakage because of fluid injection and extraction activities in sedimentary basins. Sedimentary basins typically contain vertical sequences of near-horizontal aquifers separated by less-permeable aquitards. In the context of geologic carbon sequestration (GCS), while aquitards are relied upon to contain injected buoyant free-phase CO2 and limit inter-aquifer brine flow, the large number of leaky wells (e.g., unused/unsuccessful exploration wells, or improperly plugged water or oil/gas wells) that are present in sedimentary basins creates concerns for potential leakage and groundwater contamination. The mathematical theory and solution method included in SALSA is presented in our recent manuscript under review. SALSA is the first to account for the coupled leakage through aquitards and leaky wells, with geologic pressurization sources. SALSA based on semi-analytical solutions can be computationally very efficient for problems containing many injection and leaky wells in multilayered sedimentary systems, compared to any numerical simulation software requiring local mesh refinement around each well to obtain accurate results. Therefore, stakeholders (e.g., regulators, operators) for site screening, injection and post-injection pressure behavior and assessing leakage risks can use SALSA as a fast-predictive tool in GCS applications in multilayered aquifer systems.

Cihan, Abdullah↗
DOE OSTI2010

Analyte detection using an active assay

Analytes using an active assay may be detected by introducing an analyte solution containing a plurality of analytes to a lacquered membrane. The lacquered membrane may be a membrane having at least one surface treated with a layer of polymers. The lacquered membrane may be semi-permeable to nonanalytes. The layer of polymers may include cross-linked polymers. A plurality of probe molecules may be arrayed and immobilized on the lacquered membrane. An external force may be applied to the analyte solution to move the analytes towards the lacquered membrane. Movement may cause some or all of the analytes to bind to the lacquered membrane. In cases where probe molecules are presented, some or all of the analytes may bind to probe molecules. The direction of the external force may be reversed to remove unbound or weakly bound analytes. Bound analytes may be detected using known detection types.

Morozov, Victor↗
NASA NTRS1984

Comparison of an analytical approximation for particle motion in a current sheet with precise numerical calculations

Approximate analytic solutions exist for particle motion in a one-dimensional current sheet with a constant normal magnetic field component. These solutions are tested against precise numerical calculations, and a range of validity of the analytic solutions is inferred. For example, in the geomagnetic tail neutral sheet, for a dawn-dusk electric field of 0.1-1 mV/m, lobe field of 10-40 nT, and sheet thickness of 1000 km, the analytic solutions serve as a good predictor of particle motion when the normal magnetic field component is less than 3 or 4 nT. By using the analytic solutions, initial distribution functions are mapped into final (accelerated) distributions, and the analytic mappings are compared with numerical mappings.

Speiser, T. W.↗
DOE OSTI2012

NUMERICAL SIMULATIONS OF OPTICALLY THICK ACCRETION ONTO A BLACK HOLE. I. SPHERICAL CASE

Modeling the radiation generated by accreting matter is an important step toward realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al., have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but it is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

Rotating Rake Mode Measurements Over Passive Treatment in a Ducted Fan

The NASA Glenn Research Center s Rotating Rake mode measurement system has been successful in measuring the modal content propagating in hardwall ducts. This paper proposes an extension of the Rotating Rake measurement and analysis technique to treated sections by developing basis functions based on wall impedance boundary conditions for flow conditions (i.e., constant duct area and Mach number) where the closed form analytical solution exists. Analytical equations developed to estimate mode power are incorporated. This method is verified by decomposing and analyzing radial pressure profiles generated numerically by the Eversman propagation code. Several modes, frequencies and impedances are evaluated. Data from a low-speed ducted fan with several different impedance conditions was acquired and reduced to determine the best fit to the data. Using the impedance boundary conditions result in better mode measurement solutions.

Sutliff, Daniel L.↗
DOE OSTI2011

A TWO-DIMENSIONAL METHOD OF MANUFACTURED SOLUTIONS BENCHMARK SUITE BASED ON VARIATIONS OF LARSEN'S BENCHMARK WITH ESCALATING ORDER OF SMOOTHNESS OF THE EXACT SOLUTION

The quantification of the discretization error associated with the spatial discretization of the Discrete Ordinate(DO) equations in multidimensional Cartesian geometries is the central problem in error estimation of spatial discretization schemes for transport theory as well as computer code verification. Traditionally fine mesh solutions are employed as reference, because analytical solutions only exist in the absence of scattering. This approach, however, is inadequate when the discretization error associated with the reference solution is not small compared to the discretization error associated with the mesh under scrutiny. Typically this situation occurs if the mesh of interest is only a couple of refinement levels away from the reference solution or if the order of accuracy of the numerical method (and hence the reference as well) is lower than expected. In this work we present a Method of Manufactured Solutions (MMS) benchmark suite with variable order of smoothness of the underlying exact solution for two-dimensional Cartesian geometries which provides analytical solutions aver- aged over arbitrary orthogonal meshes for scattering and non-scattering media. It should be emphasized that the developed MMS benchmark suite first eliminates the aforementioned limitation of fine mesh reference solutions since it secures knowledge of the underlying true solution and second that it allows for an arbitrary order of smoothness of the underlying ex- act solution. The latter is of importance because even for smooth parameters and boundary conditions the DO equations can feature exact solution with limited smoothness. Moreover, the degree of smoothness is crucial for both the order of accuracy and the magnitude of the discretization error for any spatial discretization scheme.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1980

Radiatively driven winds for different power law spectra

The analytic solutions for radiatively driven winds are given for the case in which the winds are driven by absorption of line and continuum radiation. The wind solutions are analytically estimated for different parameters of the central source and for different power law spectra. For flat spectra, three sonic points can exist; it is shown, however, that only one of these sonic points is physically realistic. Parameters of the central source are given which generate winds of further interest for explaining the narrow and broad absorption lines in quasars. For the quasar model presented here, winds which could give rise to the narrow absorption lines are generated by central sources with parameters which are not realistic for quasars.

Beltrametti, M.↗
NASA NTRS1992

A boundary element alternating method for two-dimensional mixed-mode fracture problems

A boundary element alternating method, denoted herein as BEAM, is presented for two dimensional fracture problems. This is an iterative method which alternates between two solutions. An analytical solution for arbitrary polynomial normal and tangential pressure distributions applied to the crack faces of an embedded crack in an infinite plate is used as the fundamental solution in the alternating method. A boundary element method for an uncracked finite plate is the second solution. For problems of edge cracks a technique of utilizing finite elements with BEAM is presented to overcome the inherent singularity in boundary element stress calculation near the boundaries. Several computational aspects that make the algorithm efficient are presented. Finally, the BEAM is applied to a variety of two dimensional crack problems with different configurations and loadings to assess the validity of the method. The method gives accurate stress intensity factors with minimal computing effort.

Raju, I. S.↗
NASA NTRS1953

Impingement of water droplets on wedges and diamond airfoils at supersonic speeds

An analytical solution has been obtained for the equations of motion of water droplets impinging on a wedge in a two-dimensional supersonic flow field with a shock wave attached to the wedge. The closed-form solution yields analytical expressions for the equation of the droplet trajectory, the local rate of impingement and the impingement velocity at any point on the wedge surface, and the total rate of impingement. The analytical expressions are utilized to determine the impingement on the forward surfaces of diamond airfoils in supersonic flow fields with attached shock waves. The results presented include the following conditions: droplet diameters from 2 to 100 microns, pressure altitudes from sea level to 30,000 feet, free-stream static temperatures from 420 degrees to 460 degrees R. Also, free-stream Mach numbers from 1.1 to 2.0, semi-apex angles for the wedge from 1.14 degrees to 7.97 degrees, thickness-to-chord ratios for the diamond airfoil from 0.02 to 0.14, chord lengths from 1 to 20 feet, and angles of attack from zero to the inverse tangent of the airfoil thickness-to-chord ratio.

SURFACE CONDITIONS - WING SECTIONS↗
NASA NTRS1954

Impingement of water droplets on wedges and double-wedge airfoils at supersonic speeds

An analytical solution has been obtained for the equations of motion of water droplets impinging on a wedge in a two-dimensional supersonic flow field with a shock wave attached to the wedge. The closed-form solution yields analytical expressions for the equation of the droplet trajectory, the local rate of impingement and the impingement velocity at any point on the wedge surface, and the total rate of impingement. The analytical expressions are utilized to determine the impingement on the forward surfaces of diamond airfoils in supersonic flow fields with attached shock waves. The results presented include the following conditions: droplet diameters from 2 to 100 microns, pressure altitudes from sea level to 30,000 feet, free-stream static temperatures from 420 degrees r, free stream Mach numbers from 1.1 to 2.0, semiapex angles for the wedge from 1.14 degrees to 7.97 degrees, thickness-to-chord ratios for the diamond airfoil from 0.02 to 0.14, chord lengths from 1 to 20 feet, and angles of attack from zero to the inverse tangent of the airfoil thickness-to-chord ratio.

Serafini, John S↗
DOE OSTI2021

Effective Permeability of a Nuclear Fuel Assembly

This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗