Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Number Theory”

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 145 records · Page 8

Continuous-variable quantum computation of the O(3) model in 1+1 dimensions

We formulate the $O(3)$ non-linear sigma model in $1+1$ dimensions as a limit of a three-component scalar field theory restricted to the unit sphere in the large squeezing limit. This allows us to describe the model in terms of the continuous variable (CV) approach to quantum computing. Here we construct the ground state and excited states using the coupled cluster ansatz and find excellent agreement with the exact diagonalization results for a small number of lattice sites. We then present the simulation protocol for the time evolution of the model using CV gates, estimate the discretization error, and present numerical results obtained from a photonic quantum simulator. We expect that the methods developed in this work will be useful for exploring interesting dynamics for a wide class of sigma models and gauge theories, as well as for simulating scattering events on quantum hardware in the coming decade.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using a finite element structural model and an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗

Good Practices in Free-energy Calculations

As access to computational resources continues to increase, free-energy calculations have emerged as a powerful tool that can play a predictive role in drug design. Yet, in a number of instances, the reliability of these calculations can be improved significantly if a number of precepts, or good practices are followed. For the most part, the theory upon which these good practices rely has been known for many years, but often overlooked, or simply ignored. In other cases, the theoretical developments are too recent for their potential to be fully grasped and merged into popular platforms for the computation of free-energy differences. The current best practices for carrying out free-energy calculations will be reviewed demonstrating that, at little to no additional cost, free-energy estimates could be markedly improved and bounded by meaningful error estimates. In energy perturbation and nonequilibrium work methods, monitoring the probability distributions that underlie the transformation between the states of interest, performing the calculation bidirectionally, stratifying the reaction pathway and choosing the most appropriate paradigms and algorithms for transforming between states offer significant gains in both accuracy and precision. In thermodynamic integration and probability distribution (histogramming) methods, properly designed adaptive techniques yield nearly uniform sampling of the relevant degrees of freedom and, by doing so, could markedly improve efficiency and accuracy of free energy calculations without incurring any additional computational expense.

Pohorille, Andrew↗

Three-nucleon lepton-number-violating potentials in chiral effective field theory and their matrix elements in light nuclei

Here, we derive the three-nucleon neutrinoless double-𝛽 decay potential in a Δ-full chiral effective field theory through next-to-next-to-next-to leading order in Weinberg's power counting. The matrix elements of the resulting operators are computed in light nuclei using variational Monte Carlo with wave functions constructed from the Norfolk family of nuclear interactions. We find that three-nucleon corrections induce a modest quenching of the total nuclear matrix elements. We discuss model dependencies and the potential impact of these corrections on the sensitivity of experimental programs to probe lepton number violating parameters. These results provide a benchmark for many-body methods capable of reaching heavier nuclei of experimental interest.

Chambers-Wall, Graham [Washington University, St. ↗

Experimental Investigation of Composite Pressure Vessel Performance and Joint Stiffness for Pyramid and Inverted Pyramid Joints

The focus of this study is on the suitability in the application of classical laminate theory analysis tools for filament wound pressure vessels with adhesive laminated joints in particular: pressure vessel wall performance, joint stiffness and failure prediction. Two 18-inch diameter 12-ply filament wound pressure vessels were fabricated. One vessel was fabricated with a 24-ply pyramid laminated adhesive double strap butt joint. The second vessel was fabricated with the same number of plies in an inverted pyramid joint. Results from hydrostatic tests are presented. Experimental results were used as input to the computer programs GENLAM and Laminate, and the output compared to test. By using the axial stress resultant, the classical laminate theory results show a correlation within 1% to the experimental results in predicting the pressure vessel wall pressure performance. The prediction of joint stiffness for the two adhesive joints in the axial direction is within 1% of the experimental results. The calculated hoop direction joint stress resultant is 25% less than the measured resultant for both joint configurations. A correction factor is derived and used in the joint analysis. The correction factor is derived from the hoop stress resultant from the tank wall performance investigation. The vessel with the pyramid joint is determined to have failed in the joint area at a hydrostatic pressure 33% value below predicted failure. The vessel with the inverted pyramid joint failed in the wall acreage at a hydrostatic pressure within 10% of the actual failure pressure.

Verhage, Joseph M.↗

The prediction of three-dimensional liquid-propellant rocket nozzle admittances

Crocco's three-dimensional nozzle admittance theory is extended to be applicable when the amplitudes of the combustor and nozzle oscillations increase or decrease with time. An analytical procedure and a computer program for determining nozzle admittance values from the extended theory are presented and used to compute the admittances of a family of liquid-propellant rocket nozzles. The calculated results indicate that the nozzle geometry entrance Mach number and temporal decay coefficient significantly affect the nozzle admittance values. The theoretical predictions are shown to be in good agreement with available experimental data.

Bell, W. A.↗

Calculation of unsteady transonic aerodynamics for oscillating wings with thickness.

An analytical approach is presented to account for some of the nonlinear characteristics of the transonic flow equation for finite thickness wings undergoing harmonic oscillation at sonic flight speed in an inviscid, shock-free fluid. The thickness effect is accounted for in the analysis through use of the steady local Mach number distribution over the wing at its mean position by employing the local linearization concept and a coordinate transformation. Computed results are compared with that of the linearized theory and experiments. Application to a flutter problem is shown.

Ruo, S. Y.↗

Calculation of unsteady transonic aerodynamics for oscillating wings with thickness

An analytical approach is presented to account for some of the nonlinear characteristics of the transonic flow equation for finite thickness wings undergoing harmonic oscillation at sonic flight speed in an inviscid, shock-free fluid. The thickness effect is accounted for in the analysis through use of the steady local Mach number distribution over the wing at its mean position by employing the local linearization concept and a coordinate transformation. Computed results are compared with that of the linearized theory and experiments. Based on the local linearization concept, an alternate formulation avoiding the limitations of the coordinate transformation method is presented.

Ruo, S. Y.↗

Multifrequency and multipolarization radar scatterometry of sand dunes and comparison with spaceborne and airborne radar images

Airborne radar scatterometer data on sand dunes, acquired at multiple frequencies and polarizations, are reported. Radar backscatter from sand dunes is very sensitive to the imaging geometry. At small incidence angles the radar return is mainly due to quasi-specular reflection from dune slopes favorably oriented toward the radar. A peak return usually occurs at the incidence angle equal to the angle of repose for the dunes. The peak angle is the same at all frequencies as computed from specular reflection theory. At larger angles the return is significantly weaker. The scatterometer measurements verified observations made with airborne and spaceborne radar images acquired over a number of dune fields in the U.S., central Africa, and the Arabian peninsula. The imaging geometry constraints indicate that possible dunes on other planets, such as Venus, will probably not be detected in radar images unless the incidence angle is less than the angles of repose of such dunes and the radar look direction is approximately orthogonal to the dune trends.

Blom, Ronald↗

Navier-Stokes simulations of axisymmetric vortex rings

The solution of the incompressible Navier-Stokes equations by means of a spectral method is used in the present numerical investigation of viscous axisymmetric vortex rings, in the cases of a single ring over the 0.001-1000 Re number range, and a pair of interacting rings. The propagation speeds of vortex rings of varying Re are computed and compared with the results of Saffman's (1970) theory; the present error estimate is found to be smaller than Saffman's, and to decrease with increasing Re at fixed core/ring radius ratio. The 'leapfrogging' of two vortex rings of equal sign and Re=1000 each is observed, indicating severe strain of the vortex core through the first pass and merging during the second pass.

Stanaway, S. K.↗

Inclined Bodies of Various Cross Sections at Supersonic Speeds

To aid in assessing effects of cross-sectional shape on body aerodynamics, the forces and moments have been measured for bodies with circular, elliptic, square, and triangular cross sections at Mach numbers 1.98 and 3.88. Results for bodies with noncircular cross sections have been compared with results for bodies of revolution having the same axial distribution of cross-sectional area (and, thus, the same equivalent fineness ratio). Comparisons have been made for bodies of fineness ratios 6 and 10 at angles of attack from 0 deg to about 20 deg and for Reynolds numbers, based on body length, of 4.0 x 10(exp 6) and 6.7 x 10(exp 6). The results of this investigation show that distinct aerodynamic advantages can be obtained by using bodies with noncircular cross sections. At certain angles of bank, bodies with elliptic, square, and triangular cross sections develop considerably greater lift and lift-drag ratios than equivalent bodies of revolution. For bodies with elliptic cross sections, lift and pitching-moment coefficients can be correlated with corresponding coefficients for equivalent circular bodies. It has been found that the ratios of lift and pitching-moment coefficients for an elliptic body to those for an equivalent circular body are practically constant with change in both angle of attack and Mach number. These lift and moment ratios are given very accurately by slender-body theory. As a result of this agreement, the method of NACA Rep. 1048 for computing forces and moments for bodies of revolution has been simply extended to bodies with elliptic cross sections. For the cases considered (elliptic bodies of fineness ratios 6 and 10 having cross-sectional axis ratios of 1.5 and 2), agreement of theory with experiment is very good. As a supplement to the force and moment results, visual studies of the flow over bodies have been made by use of the vapor-screen, sublimation, and white-lead techniques. Photographs from these studies are included in the report.

Jorgensen, Leland H.↗

A computational procedure for automated flutter analysis

A direct solution procedure for computing the flutter Mach number and the flutter frequency is applied to the aeroelastic analysis of propfans using a finite element structural model and an unsteady aerodynamic model based on a three-dimensional subsonic compressible lifting surface theory. An approximation to the Jacobian matrix that improves the efficiency of the iterative process is presented. The Jacobian matrix is indirectly approximated from approximate derivatives of the flutter matrix, which are updated only in the direction of the last move. Examples are used to illustrate the convergence properties. The direct solution procedure facilitates the automated flutter analysis in addition to contributing to the efficient use of computer time as well as the analyst's time.

Murthy, Durbha V.↗

A multimedia adult literacy program: Combining NASA technology, instructional design theory, and authentic literacy concepts

For a number of years, the Software Technology Branch of the Information Systems Directorate has been involved in the application of cutting edge hardware and software technologies to instructional tasks related to NASA projects. The branch has developed intelligent computer aided training shells, instructional applications of virtual reality and multimedia, and computer-based instructional packages that use fuzzy logic for both instructional and diagnostic decision making. One outcome of the work on space-related technology-supported instruction has been the creation of a significant pool of human talent in the branch with current expertise on the cutting edges of instructional technologies. When the human talent is combined with advanced technologies for graphics, sound, video, CD-ROM, and high speed computing, the result is a powerful research and development group that both contributes to the applied foundations of instructional technology and creates effective instructional packages that take advantage of a range of advanced technologies. Several branch projects are currently underway that combine NASA-developed expertise to significant instructional problems in public education. The branch, for example, has developed intelligent computer aided software to help high school students learn physics and staff are currently working on a project to produce educational software for young children with language deficits. This report deals with another project, the adult literacy tutor. Unfortunately, while there are a number of computer-based instructional packages available for adult literacy instruction, most of them are based on the same instructional models that failed these students when they were in school. The teacher-centered, discrete skill and drill-oriented, instructional strategies, even when they are supported by color computer graphics and animation, that form the foundation for most of the computer-based literacy packages currently on the market may not be the most effective or most desirable way to use computer technology in literacy programs. This project is developing a series of instructional packages that are based on a different instructional model - authentic instruction. The instructional development model used to create these packages is also different. Instead of using the traditional five stage linear, sequential model based on behavioral learning theory, the project uses the recursive, reflective design and development model (R2D2) that is based on cognitive learning theory, particularly the social constructivism of Vygotsky, and an epistemology based on critical theory. Using alternative instructional and instructional development theories, the result of the summer faculty fellowship is LiteraCity, a multimedia adult literacy instructional package that is a simulation of finding and applying for a job. The program, which is about 120 megabytes, is distributed on CD-ROM.

Willis, Jerry W.↗

Theory and Computation of Optimal Low- and Medium- Thrust Orbit Transfers

This report presents new theoretical results which lead to new algorithms for the computation of fuel-optimal multiple-burn orbit transfers of low and medium thrust. Theoretical results introduced herein show how to add burns to an optimal trajectory and show that the traditional set of necessary conditions may be replaced with a much simpler set of equations. Numerical results are presented to demonstrate the utility of the theoretical results and the new algorithms. Two indirect methods from the literature are shown to be effective for the optimal orbit transfer problem with relatively small numbers of burns. These methods are the Minimizing Boundary Condition Method (MBCM) and BOUNDSCO. Both of these methods make use of the first-order necessary conditions exactly as derived by optimal control theory. Perturbations due to Earth's oblateness and atmospheric drag are considered. These perturbations are of greatest interest for transfers that take place between low Earth orbit altitudes and geosynchronous orbit altitudes. Example extremal solutions including these effects and computed by the aforementioned methods are presented. An investigation is also made into a suboptimal multiple-burn guidance scheme. The FORTRAN code developed for this study has been collected together in a package named ORBPACK. ORBPACK's user manual is provided as an appendix to this report.

Goodson, Troy D.↗

Approximate and exact numerical computation of supersonic flow over an airfoil

The computation of inviscid supersonic flow over a two-dimensional airfoil is considered. There are two main nonlinear approaches which lead to approximate solutions. Small-amplitude theory gives valid solutions provided the airfoil thickness is not too great and the Mach number is not too high. The second type of approximation, shock expansion theory, employs the fact that even for flows with strong shocks the effect of the down-running characteristics remains small. This leads to an analytic solution at the airfoil. It is pointed out that the approximate solution is accurate enough for many cases of interest. The numerical method furnishes a rapid correction to the solution in those cases where it is not. The characteristic-streamline coordinate system is useful both for the computation of the approximate solution and the corrections, and is also convenient for displaying and interpreting the results.

Lewis, T. S.↗

Analytic and Computational Perspectives of Multi-Scale Theory for Homogeneous, Laminated Composite, and Sandwich Beams and Plates

This paper reviews the theoretical foundation and computational mechanics aspects of the recently developed shear-deformation theory, called the Refined Zigzag Theory (RZT). The theory is based on a multi-scale formalism in which an equivalent single-layer plate theory is refined with a robust set of zigzag local layer displacements that are free of the usual deficiencies found in common plate theories with zigzag kinematics. In the RZT, first-order shear-deformation plate theory is used as the equivalent single-layer plate theory, which represents the overall response characteristics. Local piecewise-linear zigzag displacements are used to provide corrections to these overall response characteristics that are associated with the plate heterogeneity and the relative stiffnesses of the layers. The theory does not rely on shear correction factors and is equally accurate for homogeneous, laminated composite, and sandwich beams and plates. Regardless of the number of material layers, the theory maintains only seven kinematic unknowns that describe the membrane, bending, and transverse shear plate-deformation modes. Derived from the virtual work principle, RZT is well-suited for developing computationally efficient, C0-continuous finite elements; formulations of several RZT-based elements are highlighted. The theory and its finite elements provide a unified and reliable computational platform for the analysis and design of high-performance load-bearing aerospace structures.

Tessler, Alexander↗

Analytic and Computational Perspectives of Multi-Scale Theory for Homogeneous, Laminated Composite, and Sandwich Beams and Plates

This paper reviews the theoretical foundation and computational mechanics aspects of the recently developed shear-deformation theory, called the Refined Zigzag Theory (RZT). The theory is based on a multi-scale formalism in which an equivalent single-layer plate theory is refined with a robust set of zigzag local layer displacements that are free of the usual deficiencies found in common plate theories with zigzag kinematics. In the RZT, first-order shear-deformation plate theory is used as the equivalent single-layer plate theory, which represents the overall response characteristics. Local piecewise-linear zigzag displacements are used to provide corrections to these overall response characteristics that are associated with the plate heterogeneity and the relative stiffnesses of the layers. The theory does not rely on shear correction factors and is equally accurate for homogeneous, laminated composite, and sandwich beams and plates. Regardless of the number of material layers, the theory maintains only seven kinematic unknowns that describe the membrane, bending, and transverse shear plate-deformation modes. Derived from the virtual work principle, RZT is well-suited for developing computationally efficient, C(sup 0)-continuous finite elements; formulations of several RZT-based elements are highlighted. The theory and its finite element approximations thus provide a unified and reliable computational platform for the analysis and design of high-performance load-bearing aerospace structures.

Tessler, Alexander↗