Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Primitive equations”

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.

387 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Circular-Orbit Maintenance Strategies for Primitive Body Orbiters

For missions to smaller primitive bodies, solar radiation pressure (SRP) is a significant perturbation to Keplerian dynamics. For most orbits, SRP drives large oscillations in orbit eccentricity, which leads to large perturbations from the irregular gravity field at periapsis. Ultimately, chaotic motion results that often escapes or impacts that body. This paper presents an orbit maintenance strategy to keep the orbit eccentricity small, thus avoiding the destabilizing secondary interaction with the gravity field. An estimate of the frequency and magnitude of the required maneuvers as a function of the orbit and body parameters is derived from the analytic perturbation equations.

comet↗
NASA NTRS1988

Developing fluid flow in a curved duct of square cross-section and its fully developed dual solutions

Developing fluid flow in a curved duct of square cross-section is studied numerically by a factored ADI finite-difference method on a staggered grid. A central-difference scheme with primitive variables is used inside the computational domain to reduce numerical diffusion. Two Reynolds numbers, 574 and 790, based upon a bulk velocity and hydraulic diameter are chosen for curvature ratios of 1/6.45 and 1/2.3, respectively. It is found that the secondary flow is far more complicated than expected, with the appearance of at least two pairs of vortices. Main-flow separation is also observed for the higher curvature ratio. Furthermore, it is observed that the flow develops into two quite different states downstream, depending upon the inlet conditions. Solutions of the fully developed Navier-Stokes equations is shown to be not unique beyond a certain critical Reynolds number. Developing flow seems to evolve into the fully developed state along a particular branch into which the fully developed solution bifurcates.

Soh, W. Y.↗
DOE OSTI2023

symPACK: A GPU-Capable Fan-Out Sparse Cholesky Solver

Sparse symmetric positive definite systems of equations are ubiquitous in scientific workloads and applications. Parallel sparse Cholesky factorization is the method of choice for solving such linear systems. Therefore, the development of parallel sparse Cholesky codes that can efficiently run on today’s large-scale heterogeneous distributed-memory platforms is of vital importance. Modern supercomputers offer nodes that contain a mix of CPUs and GPUs. To fully utilize the computing power of these nodes, scientific codes must be adapted to offload expensive computations to GPUs. We present symPACK, a GPU-capable parallel sparse Cholesky solver that uses one-sided communication primitives and remote procedure calls provided by the UPC++ library. We also utilize the UPC++ "memory kinds" feature to enable efficient communication of GPU-resident data. We show that on a number of large problems, symPACK outperforms comparable state-of-the-art GPU-capable Cholesky factorization codes by up to 14x on the NERSC Perlmutter supercomputer.

Bellavita, Julian↗
DOE OSTI2025

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neuromorphic Computing is Turing-Complete

Neuromorphic computing is a non-von Neumann computing paradigm that performs computation by emulating the human brain. Neuromorphic systems are extremely energy-efficient and known to consume thousands of times less power than CPUs and GPUs. They have the potential to drive critical use cases such as autonomous vehicles, edge computing and internet of things in the future. For this reason, they are sought to be an indispensable part of the future computing landscape. Neuromorphic systems are mainly used for spike-based machine learning applications, although there are some non-machine learning applications in graph theory, differential equations, and spike-based simulations. These applications suggest that neuromorphic computing might be capable of general-purpose computing. However, general-purpose computability of neuromorphic computing has not been established yet. In this work, we prove that neuromorphic computing is Turing-complete and therefore capable of general-purpose computing. Specifically, we present a model of neuromorphic computing, with just two neuron parameters (threshold and leak), and two synaptic parameters (weight and delay). We devise neuromorphic circuits for computing all the μ-recursive functions (i.e., constant, successor and projection functions) and all the μ-recursive operators (i.e., composition, primitive recursion and minimization operators). Given that the μ-recursive functions and operators are precisely the ones that can be computed using a Turing machine, this work establishes the Turing-completeness of neuromorphic computing.

Date, Prasanna↗
NASA NTRS1989

Numerical solutions of forward-flight rotor flow using an upwind method

A finite-volume upwind algorithm for solving the 3-D Euler equations with a moving grid has been developed for computing helicopter forward-flight rotor flows. The computed pressure distributions and shock positions of high-speed rotor flow are compared with various experimental data as well as with other numerical results, and the agreement is encouraging. A comparison of quasi-steady solutions with unsteady solutions reveals that when a shock occurs in the flowfield, the assumption of quasi-steady flow may fail due to the time-lag of the shock motion. Similarly, three-dimensional effects cannot be neglected. Sufficient subiterations for each time step are required to avoid numerical lag effects in using the present method. The redistribution of the residual due to the coordinate transformation is discussed. For high-order MUSCL-type schemes, a coordinate-independent solution can be obtained by interpolating primitive variables.

Chen, C. L.↗
NASA NTRS1984

Exponential-fitted methods for integrating stiff systems of ordinary differential equations: Applications to homogeneous gas-phase chemical kinetics

Conventional algorithms for the numerical integration of ordinary differential equations (ODEs) are based on the use of polynomial functions as interpolants. However, the exact solutions of stiff ODEs behave like decaying exponential functions, which are poorly approximated by polynomials. An obvious choice of interpolant are the exponential functions themselves, or their low-order diagonal Pade (rational function) approximants. A number of explicit, A-stable, integration algorithms were derived from the use of a three-parameter exponential function as interpolant, and their relationship to low-order, polynomial-based and rational-function-based implicit and explicit methods were shown by examining their low-order diagonal Pade approximants. A robust implicit formula was derived by exponential fitting the trapezoidal rule. Application of these algorithms to integration of the ODEs governing homogenous, gas-phase chemical kinetics was demonstrated in a developmental code CREK1D, which compares favorably with the Gear-Hindmarsh code LSODE in spite of the use of a primitive stepsize control strategy.

Pratt, D. T.↗
NASA NTRS1993

Particle-Gas Dynamics in the Midplane of a Protoplanetary Nebula

In this paper, we study the stage of planetary formation during which the particulate material is still dispersed as centimeter-to-meter sized primordial aggregates. During this stage, particles are able to settle toward the midplane into a layer of mass density comparable to or much greater than that of the gas. Nonlinear, coupled interactions between the particles and the nebula gas become significant and ultimately determine the vertical profiles of the particle density and the mean velocities of the particles and the nearby gas. This is the environment in which the earliest planetesimals probably form. Our numerical models rely on the Reynolds averaged Navier-Stokes equations for the gas and particles and are fully viscous and turbulent (and, for the particles. compressible). Our turbulence modeling uses a Prandtl local shear parametrization, validated by laboratory experiments. We have developed a new model for particle diffusion involving the profile of the nebula gas turbulence and the particle Schmidt number, which is a function of particle size and density. We model a cool, quiescent nebula at 1 AU (280 K) and 10 AU (90 K) and a higher temperature stage (1000 K) at 1 AU. Our main results include: (a) rapid accretion of planetesimals by gravitationally unstable fragmentation on an orbital time scale (the 'Goldreich-Ward instability') is unlikely to occur until objects have already accreted by some other process to the mass of the largest known meteorite samples, if at all; (b) from 'seeds' as small as 10 m, growth of 10- to 100-km planetesimals can proceed rapidly by drift-augmented accretion in the particle-rich midplane with orbital decay of about 1% for the growing planetesimals; (c) outward transport of vapor and small entrained chips can account for significant radial compositional and mineralogical mixing in primitive meteorite parent bodies. Other implications may be drawn concerning particle random collisional velocities and the time scales of midplane turbulence.

Cuzzi, Jeffrey N.↗
NASA NTRS1966

Fully Deuterated Aliphatic Hydrocarbons Obtained From Iron Carbide Treated with DCl and D2O

According to Oparin, Mendeleev thought that the origin of petroleum was the result of the hydrolysis of iron carbides by superheated steam under pressure from the deep interior of the Earth through geological formations where the metal carbides exist. As early as 1877, Mendeleev described the reaction leading to the synthesis of hydrocarbons according to the general equation 3Fe(sub m)C(sub n) + mH2O yields mFe3O4 + C(sub 3n)H(sub 8m). Other experimental studies on the production of hydrocarbons from cast iron have been reported. Because of the possibility that hydrocarbons may have been trapped within the carbon matrix of the cast iron, which usually has a high content of carbon, we have studied the reaction of pure iron carbide with deuterium chloride and deuterated water. This was done in order to distinguish any newly formed deuterated hydrocarbons from any possible impurities of trapped hydrocarbons. The experiments were carried out by simply allowing iron carbide to react with concentrated deuterium chloride in D2O. The volatile hydrocarbon fraction examined by gas chromatography-mass spectrometry (GC/MS), using a Finnigan 1020/OWA instrument. contained low molecular weight hydrocarbons in a range C3 to C7. Lower molecular weight hydrocarbons were not detected by GC/MS because the MS scanning mode was preset above mass 40 to exclude components of air. The identified hydrocarbons are similar to those obtained under prebiotic conditions using high frequency discharge. The hydrocarbons found in common were propane, butane, pentane, 3-methylpentane, hexane, and heptane. The percent yields decline with increasing carbon number (propane 11%, n-heptane 1%). Similar results have been obtained by the direct treatment of metal carbides by pulse laser vaporization mass spectrometry. These results show that the hydrolysis of iron carbides may have been a significant source of hydrocarbons on the primitive Earth. There appears to be a predominance of straight chain isomers, which may have been important in the prebiotic synthesis of membrane-forming compounds.

Marquez, C.↗