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 1,243 records · Page 69

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

Pythagorean coupling: Complete population transfer in a four-state system

Complete population transfer in a four-coupled-modes system is analyzed from a geometrical point of view. An analytical solution of the dynamics is written by the use of two distinct frequencies, the generalization of the single Rabi frequency of the two-state dynamics. We also present its visualization on two separate Bloch spheres with two independent torque equations. With this scheme we analytically derive the requirements for complete population transfer in a four-state quantum system. Interestingly, the solutions are found to be linked to fundamental number theory, whereas complete population transfer occurs only if the ratios between coupling coefficients exactly match a set of Pythagorean triples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Asymptotic solution for solar electric low thrust orbit raising with eclipse penalty

This paper describes an analytic solution for the solar electric low thrust orbit raising problem derived by application of the two-variable asymptotic expansion method. Included in the formulation are the effects of (1) high initial eccentricity, (2) large change of inclination, (3) variable mass, and (4) earth shadow penalty including the rotation of the shadow axis relative to the osculating orbit. An advantage of this mathematical strategy is that the effects of various steering policies can be evaluated readily from a physical point of view.

Flandro, G. A.↗

Simplified solution of diffraction from a Lyot system

This paper presents a derivation of a simplified analytical solution of diffraction from a Lyot (1939) system designed for observation of the solar corona outside of eclipses. Applying the theorem of Papoulis (1986) to simplify the calculations, a simplified solution is derived which is found to be in reasonable agreement with the exact solution. The simplified solution suffices for the preliminary evaluation of the amount of diffraction reduction needed to meet certain system requirements and also serves as a guideline for further apodization.

Wang, Yaujen↗

The MCNP6 Analytic Criticality Benchmark Suite

Analytical benchmarks provide an invaluable tool for verifying computer codes used to simulate neutron transport. Several collections of analytical benchmark problems [1-4] are used routinely in the verification of production Monte Carlo codes such as MCNP® [5,6]. Verification of a computer code is a necessary prerequisite to the more complex validation process. The verification process confirms that a code performs its intended functions correctly. The validation process involves determining the absolute accuracy of code results vs. nature. In typical validations, results are computed for a set of benchmark experiments using a particular methodology (code, cross-section data with uncertainties, and modeling) and compared to the measured results from the set of benchmark experiments. The validation process determines bias, bias uncertainty, and possibly additional margins. Verification is generally performed by the code developers, while validation is generally performed by code users for a particular application space. The VERIFICATION_KEFF suite of criticality problems [1,2] was originally a set of 75 criticality problems found in the literature for which exact analytical solutions are available. Even though the spatial and energy detail is necessarily limited in analytical benchmarks, typically to a few regions or energy groups, the exact solutions obtained can be used to verify that the basic algorithms, mathematics, and methods used in complex production codes perform correctly. The present work has focused on revisiting this benchmark suite. A thorough review of the problems resulted in discarding some of them as not suitable for MCNP benchmarking. For the remaining problems, many of them were reformulated to permit execution in either multigroup mode or in the normal continuous-energy mode for MCNP. Execution of the benchmarks in continuous-energy mode provides a significant advance to MCNP verification methods.

97 MATHEMATICS AND COMPUTING↗

Soliton solutions in two-dimensional Lorentz-violating higher derivative scalar theory

This paper shows a new approach to obtain analytical topological defects of a 2D Myers–Pospelov Lagrangian for two scalar fields. Such a Lagrangian presents higher-order kinetic terms, which lead us to equations of motion which are non-trivial to be integrated. Here we describe three possible scenarios for the equations of motion, named by timelike, spacelike and lightlike respectively. We started our investigation with a kink-like traveling wave Ansatz for the free theory, which led us to constraints for the dispersion relations of each scenario. We also introduced a procedure to obtain analytical solutions for the general theory in the three mentioned scenarios. We exemplified the procedure and discussed the behavior of the defect solutions carefully. It is remarkable that the methodology presented in this study led to analytical models, despite the complexity of the equations of motion derived from the 2D Myers–Pospelov Lagrangian. The methodology here tailored can be applied to several Lagrangians with higher-order derivative terms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Efficient co-solution of time step size and independent state in simulations of fluid-driven fracture propagation with embedded meshes

Here we present an efficient time-continuation scheme for fluid-driven fracture propagation problems in the extended finite element method framework. The approach applies a monolithic solution strategy to a fully coupled and implicit approximation of hydro-mechanical systems in conjunction with simultaneous linear elastic propagation of multiple fractures. At the end of each time step, the process ensures that the weakest fracture tip is in an equilibrium propagation regime. Furthermore, the solution process provides an initialization procedure for the newly created fracture spaces and an a priori estimate of the stress intensity factor growth rate, improving simulation robustness, and efficiency. The solution process is validated using the Kristianovich-Geertsma-de Klerk analytical solution under the toughness- and viscosity-dominated regimes. It is also extended to and demonstrated on problems with multiple fractures undergoing simultaneous propagation with stress shadow interactions. Numerical examples demonstrate that the solution process can reduce the required computational cost by one order of magnitude compared to other existing methods.

42 ENGINEERING↗

Electric fields near undulating dielectric membranes

Dielectric interfaces are crucial to the behavior of charged membranes, from graphene to synthetic and biological lipid bilayers. Understanding electrolyte behavior near these interfaces remains a challenge, especially in the case of rough dielectric surfaces. A lack of analytical solutions consigns this problem to numerical treatments. Here, we report an analytic method for determining electrostatic potentials near curved dielectric membranes in a two-dimensional periodic “slab” geometry using a periodic summation of Green’s functions. This method is amenable to simulating arbitrary groups of charges near surfaces with two-dimensional deformations. We concentrate on one-dimensional undulations. We show that increasing membrane undulation increases the asymmetry of interfacial charge distributions due to preferential ionic repulsion from troughs. In the limit of thick membranes, we recover results mimicking those for electrolytes near a single interface. Our work demonstrates that rough surfaces generate charge patterns in electrolytes of charged molecules or mixed-valence ions.

Chemistry↗

Verification Test of the SURF and SURFplus Models in xRage: Part III Affect of Mesh Alignment

The previous studies used an underdriven detonation wave in 1-dimension (steady ZND reaction zone profile followed by a scale-invariant rarefaction wave) for PBX 9502 as a verification test of the implementation of the SURF and SURFplus models in the xRage code. Since the SURF rate is a function of the lead shock pressure, the question arises as to the effect on accuracy of variations in the detected shock pressure due to the alignment of the shock front with the mesh. To study the effect of mesh alignment we simulate a cylindrically diverging detonation wave using a planar 2-D mesh. The leading issue is the magnitude of azimuthal asymmetries in the numerical solution. The 2-D test case does not have an exact analytic solution. To quantify the accuracy, the 2-D solution along rays through the origin are compared to a highly resolved 1-D simulation in cylindrical geometry.

97 MATHEMATICS AND COMPUTING↗

On the statistical uncertainties associated with line profile fitting

The paper considers the statistical uncertainties in the best-fit parameters associated with the fitting of observational data by model line profiles of specific form. The parameters involved are amplitude, center, and width. A general procedure is described for any profile type of parametric form. The case of Gaussian profiles is treated in detail, and numerical solutions are presented for equally spaced data of equal weight. Approximate analytical solutions are derived that give good general agreement with the numerical solutions over the whole range of interest. Approximate solutions are also obtained for the case of Lorentzian profiles. It is shown that the Lorentzian results can be used in conjunction with the Gaussian results to estimate the corresponding parameter uncertainties in the case of the general Voigt profile.

Landman, D. A.↗

Numerical solutions for heat flow in adhesive lap joints

The present formulation for the modeling of heat transfer in thin, adhesively bonded lap joints precludes difficulties associated with large aspect ratio grids required by standard FEM formulations. This quasi-static formulation also reduces the problem dimensionality (by one), thereby minimizing computational requirements. The solutions obtained are found to be in good agreement with both analytical solutions and solutions from standard FEM programs. The approach is noted to yield a more accurate representation of heat-flux changes between layers due to a disbond.

Howell, P. A.↗

Analysis of gas absorption to a thin liquid film in the presence of a zero-order chemical reaction

The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Rajagopalan, S.↗

COMPARISON OF THE ACCURACY OF VARIOUS SPATIAL DISCRETIZATION SCHEMES OF THE DISCRETE ORDINATES EQUATIONS IN 2D CARTESIAN GEOMETRY

We present a comprehensive error estimation of four spatial discretization schemes of the two-dimensional Discrete Ordinates (SN) equations on Cartesian grids utilizing a Method of Manufactured Solution (MMS) benchmark suite based on variants of Larsen’s benchmark featuring different orders of smoothness of the underlying exact solution. The considered spatial discretization schemes include the arbitrarily high order transport methods of the nodal (AHOTN) and characteristic (AHOTC) types, the discontinuous Galerkin Finite Element method (DGFEM) and the recently proposed higher order diamond difference method (HODD) of spatial expansion orders 0 through 3. While AHOTN and AHOTC rely on approximate analytical solutions of the transport equation within a mesh cell, DGFEM and HODD utilize a polynomial expansion to mimick the angular flux profile across each mesh cell. Intuitively, due to the higher degree of analyticity, we expect AHOTN and AHOTC to feature superior accuracy compared with DGFEM and HODD, but at the price of potentially longer grind times and numerical instabilities. The latter disadvantages can result from the presence of exponential terms evaluated at the cell optical thickness that arise from the semianalytical solution process. This work quantifies the order of accuracy and the magnitude of the error of all four discretization methods for different optical thicknesses, scattering ratios and degrees of smoothness of the underlying exact solutions in order to verify or contradict the aforementioned intuitive expectation.

99 GENERAL AND MISCELLANEOUS↗

Comparison of analytic models of instability of rarefied gas flow in a channel

Numerical and analytical results are compared concerning the limit properties of the trajectories, attractors and bifurcations of rarefied gas flows in channels. The cascade of bifurcations obtained in our previous analytical and numerical investigations is simulated numerically for different scattering functions V generalizing the ray-diffuse reflection of gas particles from the surface. The main purpose of numerical simulation by Monte Carlo method is the investigation of the properties of different analytic nonlinear dynamic systems corresponding to rarefied gas flow in a channel. The results are compared as well for the models suggested originally by R. N. Miroshin, as for the approximations considered for the first time or for studied in our subsequent papers. Analytical solutions we obtained earlier for the ray reflection which means only one determined velocity of scattered from the walls gas atoms, generally different from the specular reflection. The nonlinear iterative equation describing a rarefied gas flow in a long channel becomes unstable in some regions of corresponding parameters of V (it means the sensitivity to boundary conditions). The values of the parameters are found from analytical approximations in these regions. Numerical results show that the chaotic behavior of the nonlinear dynamic system corresponds to strange attractors and distinguishes clearly from Maxwellian distribution and from the equilibrium on the whole. In the regions of instability (as the dimension of the attractor increases) the search for a corresponding state requires a lot more computation time and a lot of data (the amount of data required increases exponentially with embedding dimension). Therefore the main complication in the computation is reducing as well the computing time as the amount of data to find a suitably close solution. To reduce the computing time our analytical results are applied. Flow conditions satisfying the requirements to the experiment are indicated where the instability of considered type can be detected. The advantages and the drawbacks of considered approximations are detected and the features of obtained solutions are indicated. The recommendations are given for applying the results in practical applications and in numerical calculations of rarefied gas flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Continuum calculations of continental deformation in transcurrent environments

A thin viscous sheet approximation is used to investigate continental deformation near a strike-slip boundary. The vertically averaged velocity field is calculated for a medium characterized by a power law rheology with stress exponent n. Driving stresses include those applied along boundaries of the sheet and those arising from buoyancy forces related to lateral differences in crustal thickness. Exact and approximate analytic solutions for a region with a sinusoidal strike-slip boundary condition are compared with solutions for more geologically relevant boundary conditions obtained using a finite element technique. The across-strike length scale of the deformation is approximately 1/4pi x sq rt n times the dominant wavelength of the imposed strike-slip boundary condition for both the analytic and the numerical solutions; this result is consistent with length scales observed in continental regions of large-scale transcurrent faulting. An approximate, linear relationship between displacement and rotation is found that depends only on the deformation length scale and the rheology. Calculated displacements, finite rotations, and distribution of crustal thicknesses are consistent with those observed in the region of the Pacific-North America plate boundary in California.

Sonder, L. J.↗

Performance and Feature Improvements in Parareal-based Power System Dynamic Simulation

In recent years, a novel Parareal-based approach has been developed for fast transient simulations of large power system interconnections. Parareal belongs to the class of Parallel-in-time algorithms for solution of systems of differential-algebraic equations in parallel over an interval of time. The selection of a reasonably fast and accurate coarse solution is crucial to improve the performance of Parareal algorithm. Semi-analytical solution methods are one promising approach to achieve this goal. They have been investigated, and some preliminary results are presented here. In addition, Parareal-based simulator has been expanded to enable co-simulation with OpenDSS, a widely used open-source distribution system simulator. Preserving the parallel nature of the Parareal approach and taking advantage of the parallel capabilities of the latest versions of OpenDSS, each distribution system can be solved in their entirety on different processors in parallel within the main Parareal simulator. This paper also presents the structure of the transmission and distribution co-simulation and some results with different dynamic models of inverter-based resources in the distribution systems.

Park, Byungkwon↗