Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical calculations: variational”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1923

General biplane theory

This report deals with the air forces on biplane cellule. The first part of the report deals with the two-dimensional problem neglecting viscosity. The variation of the section, chord, gap, stagger, and decalage are investigated, a great number of examples are calculated, and all numerical results are given in tables. For the biplane without stagger it is found that the loss of lift in consequence of the mutual influence of the two wing sections is only half as much if the lift is produced by the curvature of the section as it is when the lift is produced by the inclination of the chord to the direction of motion. The second part deals with the influence of the lateral dimensions. It is found that the loss of lift due to induction is almost unchanged, whether the biplane is staggered or not. In the third part conclusions from previous investigations are drawn, viscosity and experimental experience are brought in, and the method is simplified for practical application. Simple formulas give the drag, lift, and moment. In order to make use of the simple formulas more convenient, tables for the dynamical pressure, induced drag, and angle of attack are added so that practically no computation is needed for the application of the results.

Munk, Max M↗
NASA NTRS2003

What is the Time Scale for Orbital Forcing of the Martian Water Cycle?

Calculation of the periodic variations in the martian orbital parameters by Ward and subsequent refinements to the theory have inspired numerous models of variation of the martian water cycle. Most of these models have focused on variations in planetary obliquity on a both a short-term (110 kyr) time scale as well as larger oscillations occuring over millions of years. To a lesser extent, variations in planetary eccentricity have also been considered. The third and fastest mode of variation, the precession of the longitude of perihelion, has generally been deemphasized because, among the three parameters, it is the only one that does not change the integrated annual insolation. But as a result of this precession, the asymmetry in peak summer insolation between the poles exceeds 50%, with the maximum cycling between poles every 25.5 kyrs. The relative contribution of these different elements to orbital forcing of climate takes on particular importance in the context of apparently recent waterrelated features such as gullies or polar layered deposits (PLD). Christensen, for example, recently indentified mantling of heavily gullied crater walls as residual dust-covered snow deposits that were responsible for the formation of the gullies in a previous epoch. Christensen assumed that the snow was originally deposited at a period of high obliquity which was stabilized against sublimation by a lag deposit of dust. It is suggested here that not obliquity, but the shortterm oscillations associated with precession of the perihelion may play the dominant role in the formation of gullies, major strata in the polar layered deposits (PLD), and other water-related features.

M H Hecht↗
NASA NTRS2003

Speculations on Orbital Forcing of Sublimation from the Polar Caps

Calculation of the periodic variations in the martian orbital parameters by Ward and subsequent refinements to the theory have inspired numerous models of variation of the martian water cycle. The limitations of models such as this can not be overstated. Albedo, for example, has a much greater influence on temperature than orbital forcing, and seems to be linked to insolation, possibly via its influence on dust transport. To make further progress, either new constraints must be identified, or additional clues must be obtained from remote sensing and in situ exploration.

Hecht, M. H.↗
NASA NTRS2003

Speculations on orbital forcing of sublimination from the polar caps

Calculation of the periodic variations in the martian orbital parameters by Ward and subsequent refinements to the theory, have inspired numerous models of variation of the martian water cycle. The limitations of models such as this can not be overstated. Albedo, for example, has a much greater influence on temperature than orbital forcing, and seems to be linked to insolation, possible via its influence on dust transport. To make further progress, either new constraints must be identified, or additional clues must be obtained from remote sensing and in situ exploration

Hecht, Michael H.↗
NASA NTRS2001

Free Vibration of Uncertain Unsymmetrically Laminated Beams

Monte Carlo Simulation and Stochastic FEA are used to predict randomness in the free vibration response of thin unsymmetrically laminated beams. For the present study, it is assumed that randomness in the response is only caused by uncertainties in the ply orientations. The ply orientations may become random or uncertain during the manufacturing process. A new 16-dof beam element, based on the first-order shear deformation beam theory, is used to study the stochastic nature of the natural frequencies. Using variational principles, the element stiffness matrix and mass matrix are obtained through analytical integration. Using a random sequence a large data set is generated, containing possible random ply-orientations. This data is assumed to be symmetric. The stochastic-based finite element model for free vibrations predicts the relation between the randomness in fundamental natural frequencies and the randomness in ply-orientation. The sensitivity derivatives are calculated numerically through an exact formulation. The squared fundamental natural frequencies are expressed in terms of deterministic and probabilistic quantities, allowing to determine how sensitive they are to variations in ply angles. The predicted mean-valued fundamental natural frequency squared and the variance of the present model are in good agreement with Monte Carlo Simulation. Results, also, show that variations between plus or minus 5 degrees in ply-angles can affect free vibration response of unsymmetrically and symmetrically laminated beams.

Kapania, Rakesh K.↗
NASA NTRS1992

Multi-dimensional high order essentially non-oscillatory finite difference methods in generalized coordinates

The nonlinear stability of compact schemes for shock calculations is investigated. In recent years compact schemes were used in various numerical simulations including direct numerical simulation of turbulence. However to apply them to problems containing shocks, one has to resolve the problem of spurious numerical oscillation and nonlinear instability. A framework to apply nonlinear limiting to a local mean is introduced. The resulting scheme can be proven total variation (1D) or maximum norm (multi D) stable and produces nice numerical results in the test cases. The result is summarized in the preprint entitled 'Nonlinearly Stable Compact Schemes for Shock Calculations', which was submitted to SIAM Journal on Numerical Analysis. Research was continued on issues related to two and three dimensional essentially non-oscillatory (ENO) schemes. The main research topics include: parallel implementation of ENO schemes on Connection Machines; boundary conditions; shock interaction with hydrogen bubbles, a preparation for the full combustion simulation; and direct numerical simulation of compressible sheared turbulence.

Shu, Chi-Wang↗
DOE OSTI2023

Variational Monte Carlo calculations of n + H 3 scattering

A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1955

Analysis of behavior of simply supported flat plates compressed beyond the buckling load into the plastic range

An analysis is presented of the postbuckling behavior of a simply supported square flat plate with straight edges compressed beyond the buckling load into the plastic range. The method of analysis involves the application of a variational principle of the deformation theory of plasticity in conjunction with computations carried out on a high-speed calculating machine. Numerical results are obtained for several plate proportions and for one material. The results indicate plate strengths greater than those that have been found experimentally on plates that do not satisfy straight-edge conditions. (author)

Mayers, J↗
NASA NTRS1993

Plasma-screening effects on the electron-impact excitation of hydrogenic ions in dense plasmas

Plasma-screening effects are investigated on electron-impact excitation of hydrogenic ions in dense plasmas. Scaled cross sections Z(exp 4) sigma for 1s yields 2s and 1s yields 2p are obtained for a Debye-Hueckel model of the screened Coulomb interaction. Ground and excited bound wave functions are modified in the screened Coulomb potential (Debye-Hueckel model) using the Ritz variation method. The resulting atomic wave functions and their eigenenergies agree well with the numerical and high-order perturbation theory calculations for the interesting domain of the Debye length not less than 10. The Born approximation is used to describe the continuum states of the projectile electron. Plasma screening effects on the atomic electrons cannot be neglected in the high-density cases. Including these effects, the cross sections are appreciably increased for 1s yields 2s transitions and decreased for 1s yields 2p transitions.

Jung, Young-Dae↗
DOE OSTI2022

Quantum algorithm for electronic band structures with local tight-binding orbitals

While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum “toolbox”. With this in mind, we consider the quantum calculation of a periodic system’s single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

The Passive Microwave Remote Sensing of Soil Moisture: the Effect of Tilled Row Structure

The tilled rowstructure is known to be one of the important factors affecting the observations of the microwave emission from a natural surface. Measurements of this effect were carried out with both I and X band radiometers mounted on a mobile truck on a bare 40 m x 45 m row tilled field. The soil moisture content during the measurements ranged from approximately 10 percent to approximately 30 percent by dry weight. The results of these measurements showed that the variations of the antenna temperatures with incident angle theta changed with the azimuthal angle a measured from the row direction. A numerical calculation based on a composite surface roughness was made and found to predict the observed features within the model's limit of accuracy. It was concluded that the difference between the horizontally and vertically polarized temperatures was due to the change in the local angle of field emission within the antenna field of view caused by the large scale row structure.

Wang, J. R.↗
NASA NTRS1989

The Hard X-rays and Gamma-rays from Solar Flares

Radiation of energies from 10 KeV to greater than 10 MeV has been observed during solar flares, and is interpreted to be due to bremsstrahlung by relativistic electrons. A complete treatment of this problem requires solution of the kinetic equation for relativistic electrons and inclusion of synchrotron energy losses. Using the electron distributions obtained from numerical solutions of this equation the bremsstrahlung spectra in the impulsive x ray and gamma-ray regimes are calculated, and the variation of these spectral indices and directivities with energy and observation angle are described. The dependences of these characteristics of the radiation of changes in the solar atmospheric model, including the convergence of the magnetic field, the injected electron spectral index, and most importantly, in the anisotropy of the injected electrons and of the convergence of the magnetic field are also described. The model results are compared with stereoscopic observations of individual flares and the constraints that this data sets on the models are discussed.

Mctiernan, James M.↗
NASA NTRS1990

The hard X-rays and gamma rays from solar flares

Radiation of energies from 10 KeV to greater than 10 MeV has been observed during solar flares, and is interpreted to be due to bremsstrahlung by relativistic electrons. A complete treatment of this problem requires solution of the kinetic equation for relativistic electrons and inclusion of synchrotron energy losses. Using the electron distributions obtained from numerical solutions of this equation the bremsstrahlung spectra in the impulsive X-ray and gamma ray regimes are calculated, and the variation of these spectral indices and directivities with energy and observation angle are described. The dependences of these characteristics of the radiation of changes in the solar atmospheric model, including the convergence of the magnetic field, the injected electron spectral index, and most importantly, in the anisotropy of the injected electrons and of the convergence of the magnetic field are also described. The model results are compared with stereoscopic observations of individual flares and the constraints that this data sets on the models are discussed.

Mctiernan, James M.↗
NASA NTRS2003

High Temperature Permeability of Carbon Cloth Phenolic Composite

The carbon fiber phenolic resin composite material used for the RSRM nozzle insulator occasionally experiences problems during operation from pocketing or spalling-like erosion and lifting of plies into the char layer. This phenomenon can be better understood if the permeability of the material at elevated temperatures is well defined. This paper describes an experimental approach to determining high temperature permeability of the carbon phenolic material used as the RSRM nozzle liner material. Two different approaches were conducted independently using disk and bar type specimens with the designed permeability apparatus. The principle of the apparatus was to subject a test specimen to a high pressure differential and a heat supply and to monitor both the pressure and temperature variations resulting from gas penetration through the permeable wall between the two chambers. The bar types, especially designed to eliminate sealing difficulties at a high temperature environment, were directly exposed to real time temperature elevation from 22 C to 260 C during the test period. The disk types were pre-heat treated up to 300 C for 8 hours and cooled to room temperature before testing. Nonlinear variation of downstream pressure at a certain temperature range implied moisture release and matrix pyrolysis. Permeability was calculated using a semi-numerical model of quasi-steady state. The test results and the numerical model are discussed in the paper.

Park, O. Y.↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
NASA NTRS1991

Biological production models as elements of coupled, atmosphere-ocean models for climate research

Process models of phytoplankton production are discussed with respect to their suitability for incorporation into global-scale numerical ocean circulation models. Exact solutions are given for integrals over the mixed layer and the day of analytic, wavelength-independent models of primary production. Within this class of model, the bias incurred by using a triangular approximation (rather than a sinusoidal one) to the variation of surface irradiance through the day is computed. Efficient computation algorithms are given for the nonspectral models. More exact calculations require a spectrally sensitive treatment. Such models exist but must be integrated numerically over depth and time. For these integrations, resolution in wavelength, depth, and time are considered and recommendations made for efficient computation. The extrapolation of the one-(spatial)-dimension treatment to large horizontal scale is discussed.

Platt, Trevor↗
NASA NTRS1985

On a class of TVD schemes for gas dynamic calculations

The purpose of this paper is to review a class of explicit and implicit second-order accurate Total Variation Diminishing (TVD) schemes and to show by numerical experiments, the performance of these schemes to the Euler equations of gas dynamics. The method of constructing these second-order accurate TVD schemes is sometimes known as the modified flux approach.

Yee, H. C.↗