Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass conservation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Advanced Small Perturbation Potential Flow Theory for Unsteady Aerodynamic and Aeroelastic Analyses

An advanced small perturbation (ASP) potential flow theory has been developed to improve upon the classical transonic small perturbation (TSP) theories that have been used in various computer codes. These computer codes are typically used for unsteady aerodynamic and aeroelastic analyses in the nonlinear transonic flight regime. The codes exploit the simplicity of stationary Cartesian meshes with the movement or deformation of the configuration under consideration incorporated into the solution algorithm through a planar surface boundary condition. The new ASP theory was developed methodically by first determining the essential elements required to produce full-potential-like solutions with a small perturbation approach on the requisite Cartesian grid. This level of accuracy required a higher-order streamwise mass flux and a mass conserving surface boundary condition. The ASP theory was further developed by determining the essential elements required to produce results that agreed well with Euler solutions. This level of accuracy required mass conserving entropy and vorticity effects, and second-order terms in the trailing wake boundary condition. Finally, an integral boundary layer procedure, applicable to both attached and shock-induced separated flows, was incorporated for viscous effects. The resulting ASP potential flow theory, including entropy, vorticity, and viscous effects, is shown to be mathematically more appropriate and computationally more accurate than the classical TSP theories. The formulaic details of the ASP theory are described fully and the improvements are demonstrated through careful comparisons with accepted alternative results and experimental data. The new theory has been used as the basis for a new computer code called ASP3D (Advanced Small Perturbation - 3D), which also is briefly described with representative results.

Batina, John T.↗
NASA NTRS1983

Flow in coronal loops with a mass source

This research studies the flow of plasma inside a coronal loop in which an injection of plasma through the lateral surface is permitted. The flow is assumed steady and polytropic. The problem covers two cases: (a) upflow at one footpoint, downflow at the other; (b) downflow at both footpoints. The first case can be shown to be quite similar to that of a mass-conserving flow with variable cross section; the second, instead, is characterized by solutions with a different type of topology; its main new feature is the obvious fact that all the solutions pass through a single point going from negative to positive velocities. In this second case the density ratio between footpoints and top can be much smaller than in a mass conserving flow. This can explain some properties of observed loops.

Noci, G.↗
DOE OSTI2024

Multifidelity methods for uncertainty quantification of a nonlocal model for phase changes in materials

This study is devoted to the construction of a multifidelity Monte Carlo (MFMC) method for the uncertainty quantification of a nonlocal, non-mass-conserving Cahn-Hilliard model for phase transitions with an obstacle potential. Here, we are interested in estimating the expected value of an output of interest (OoI) that depends on the solution of the nonlocal Cahn-Hilliard model. As opposed to its local counterpart, the nonlocal model captures sharp interfaces without the need for significant mesh refinement. However, the computational cost of the nonlocal Cahn-Hilliard model is higher than that of its local counterpart with similar mesh refinement, inhibiting its use for outer-loop applications such as uncertainty quantification. The MFMC method augments the desired high-fidelity, high-cost OoI with a set of lower-fidelity, lower-cost OoIs to alleviate the computational burden associated with nonlocality. Most of the computational budget is allocated to sampling the cheap surrogate models to achieve speedup, whereas the high-fidelity model is sparsely sampled to maintain accuracy. For the non-mass-conserving nonlocal Cahn-Hilliard model, the use of the MFMC method results in, for a given computational budget, about an order of magnitude reduction in the mean-squared error of the expected value of the OoI relative to that of the Monte Carlo method.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2005

Excitation of Earth Rotation Variations "Observed" by Time-Variable Gravity

Time variable gravity measurements have been made over the past two decades using the space geodetic technique of satellite laser ranging, and more recently by the GRACE satellite mission with improved spatial resolutions. The degree-2 harmonic components of the time-variable gravity contain important information about the Earth s length-of-day and polar motion excitation functions, in a way independent to the traditional "direct" Earth rotation measurements made by, for example, the very-long-baseline interferometry and GPS. In particular, the (degree=2, order= 1) components give the mass term of the polar motion excitation; the (2,O) component, under certain mass conservation conditions, gives the mass term of the length-of-day excitation. Combining these with yet another independent source of angular momentum estimation calculated from global geophysical fluid models (for example the atmospheric angular momentum, in both mass and motion terms), in principle can lead to new insights into the dynamics, particularly the role or the lack thereof of the cores, in the excitation processes of the Earth rotation variations.

Chao, Ben F.↗
DOE OSTI2022

A consistent and conservative Phase-Field method for multiphase incompressible flows

In the present study, a consistent and conservative Phase-Field method, including both the model and scheme, is developed for multiphase flows with an arbitrary number of immiscible and incompressible fluid phases. The consistency of mass conservation and the consistency of mass and momentum transport are implemented to address the issue of physically coupling the Phase-Field equation, which locates different phases, to the hydrodynamics. These two consistency conditions, as illustrated, provide the “optimal” coupling because (i) the new momentum equation resulting from them is Galilean invariant and implies the kinetic energy conservation, regardless of the details of the Phase-Field equation, and (ii) failures of satisfying the second law of thermodynamics or the consistency of reduction of the multiphase flow model only result from the same failures of the Phase-Field equation but are not due to the new momentum equation. Physical interpretation of the consistency conditions and their formulations are first provided, and general formulations that are obtained from the consistency conditions and independent of the interpretation of the velocity are summarized. Then, the present consistent and conservative multiphase flow model is completed by selecting a reduction consistent Phase-Field equation. Several novel techniques are developed to inherit the physical properties of the multiphase flows after discretization, including the gradient-based phase selection procedure, the momentum conservative method for the surface force, and the general theorems to preserve the consistency conditions on the discrete level. Equipped with those novel techniques, a consistent and conservative scheme for the present multiphase flow model is developed and analyzed. The scheme satisfies the consistency conditions, conserves the mass and momentum, and assures the summation of the volume fractions to be unity, on the fully discrete level and for an arbitrary number of phases. All those properties are numerically validated. Finally, numerical applications demonstrate that the present model and scheme are robust and effective in studying complicated multiphase dynamics, especially for those with large-density ratios.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2022

Ultraviolet Spectropolarimetry Conservative and Nonconservative Mass Transfer in OB Interacting Binaries

The current consensus is that at least half of the OB stars are formed in binary or multiple star systems. The evolution of OB stars is greatly influenced by whether the stars begin as close binaries, and the evolution of the binary systems depend on whether the mass transfer is conservative or nonconservative. FUV/NUV spectropolarimetry is poised to answer the latter question. This paper discusses how the Polstar spectropolarimetry mission can characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase, and quantify its amount. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 10 −3 precision or better. The spectroscopy will identify absorption by mass streams and other plasmas seen in projection against the stellar disk as a function of orbital phase, as well as scattering from extended splash structures, including jets. The polarimetry tracks the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Special attention will be given to the epochs of enhanced systemic mass loss that have been identified from IUE observations (pre-mass reversal and tangential gas stream impact). We show how the history of systemic mass and angular momentum loss/gain episodes can be inferred via ensemble evolution through the r–q diagram. Combining the above elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass and ultimate supernova of a large number of massive stars found in binaries at close enough separation to undergo interaction.

Geraldine J. Peters↗
NASA NTRS1997

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

In spite of the rapid advances in both scalar and parallel computational tools, the large number and breadth of variables involved in aerodynamic systems make the use of parabolized or even boundary layer fluid flow models impractical for both preliminary design and inverse design problems. Given this restriction, we have concluded that reduced or approximate models are an important family of tools for design purposes. This study of a combined perturbation/numerical modeling methodology with an application to ejector-mixer nozzles (shown schematically in the following figure) is nearing completion. The work is being funded by a grant from the NASA Lewis Research Center to Texas A&M University. These ejector-mixer nozzle models are designed to be of use to the High Speed Civil Transport Program and may be adopted by both NASA and industry. A computer code incorporating the ejector-mixer models is under development. This code, the Differential Reduced Ejector/Mixer Analysis (DREA), can be run fast enough to be used as a subroutine or to be called by a design optimization routine. Simplified conservation equations--x-momentum, energy, and mass conservation--are used to define the model. Unlike other preliminary design models, DREA requires minimal empirical input and includes vortical mixing and a fully compressible formulation among other features. DREA is being validated by comparing it with results obtained from open literature and proprietary industry data. Preliminary results for a subsonic ejector and a supersonic ejector are shown. In addition, dedicated experiments have been performed at Texas A&M. These experiments use a hydraulic/gas flow analog to provide information about the inviscid mixing interface structure. Final validation and documentation of this work is expected by May of 1997. However, preliminary versions of DREA can be expected in early 1997. In summary, DREA provides a sufficiently detailed and realistic ejector-mixer nozzle model at a computational cost compatible with preliminary design applications.

DeChant, Lawrence J.↗
DOE OSTI2026

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Transport of absolute angular momentum in quasi-axisymmetric equatorial jet streams

It is well known that prograde equatorial jet stresses cannot occur in an axisymmetric inviscid fluid, owing to the constraints of local angular momentum conservation. For a viscous fluid, the constraints of mass conservation prevent the formation of any local maximum of absolute angular momentum (m) without a means of transferring m against its gradient (delta m) in the meridional plane. The circumstances under which m can be diffused up-gradient by normal molecular viscosity are derived, and illustrated with reference to numerical simulations of axisymmetric flows in a cylindrical annulus. Viscosity is shown to act so as to tend to expel m from the interior outwards from the rotation axis. Such an effect can produce local super-rotation even in a mechanically isolated fluid. The tendency of viscosity to result in the expulsion of m is shown to be analogous in certain respects to a vorticity-mixing hypothesis for the effects of non-axisymmetric eddies of the zonally-averaged flow. It is shown how the advective and diffusive transport of m by non-axisymmetric eddies can be represented by the Transformed Eulerian Mean meridional circulation and the Eliassen-Palm (EP) flux of Andrews and McIntyre respectively, in the zonal mean. Constraints on the form and direction of the EP flux in an advective/diffusive flow for such eddies are derived, by analogy with similar constraints on the diffusive flux of m due to viscosity.

Read, P. L.↗
DOE OSTI2023

Operation Optimization Using Reinforcement Learning with Integrated Artificial Reasoning Framework

In large and complex systems, operational decision-making requires a systematic analysis with a vast amount of data from both process parameters and component status monitoring. In this paper, we present an integrated artificial reasoning approach for system state transition models that can help operational decision-making with explainable and traceable reasoning. The integrated artificial reasoning framework is a physics-based approach of defining the system structure in a Bayesian network, so we leveraged it in a Markov decision process (MDP) for finding optimal operational solutions. In our proposed framework, the MDP is implemented on a dynamic Bayesian network (DBN), which represents causalities in a system. The multilevel flow modeling was utilized in order to extract these causalities in a more efficient and objective manner. Since multilevel flow modeling is based on the fundamental energy and mass conservation laws, the target system is decomposed into several mass, energy, and information structures, which serve as the basis for a DBN. The MDP consists of the processes of finding a solution for the Bellman equation, which can be derived from the conditional probability equations of the constructed DBN. System operators can capture stochastic system dynamics as multiple subsystem state transitions based on their physical relations and uncertainties coming from the component degradation process or random failures. We analyzed a simplified example system to illustrate finding an optimal operational policy with this approach.

Kim, Junyung↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗
DOE OSTI2022

TYC 2990-127-1: An Algol-type SB2 binary system of subgiant and red giant with a probable ongoing mass-transfer

ABSTRACT We present a study of the spectroscopic binary TYC 2990-127-1 from the LAMOST survey. We use full-spectrum fitting to derive radial velocities and spectral parameters. The high mass ratio indicates that the system underwent mass transfer in the past. We compute the orbital solution and find that it is a very close sub-giant/red giant pair on circular orbit, slightly inclined to the sky-plane. Fitting of the TESS photometrical data confirms this and suggests an inclination of i ∼ 39.8°. The light curve and spectrum around Hα show signs of irregular variability, which supports ongoing mass transfer. The binary evolution simulations suggest that the binary may experience non-conservative mass transfer with accretion efficiency 0.3, and the binary will enter into common envelope (CE) phase in the subsequent evolution. The remnant product after the ejection of CE may be a detached double helium white dwarf (He WD) or a merger.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

Finite area method for nonlinear supersonic conical flows

A fully conservative numerical method for the computation of steady inviscid supersonic flow about general conical bodies at incidence is described. The procedure utilizes the potential approximation and implements a body conforming mesh generator. The conical potential is assumed to have its best linear variation inside each mesh cell; a secondary interlocking cell system is used to establish the flux balance required to conserve mass. In the supersonic regions the scheme is symmetrized by adding artificial viscosity in conservation form. The algorithm is nearly an order of a magnitude faster than present Euler methods and predicts known results accurately and qualitative features such as nodal point lift off correctly. Results are compared with those of other investigators.

Sritharan, S. S.↗
NASA NTRS1975

Numerical solution of the Navier-Stokes equations with topography

A finite difference scheme for solving the equations of fluid motion in a generalized coordinate system has been constructed. The scheme conserves mass and all the first integral moments of the motion. The scheme also advectively 'almost conserves' second moments, in that the magnitude of implicit numerical smoothing is typically about an order smaller than explicit viscosity and diffusion. Calculations with the model support the theoretical conjecture that the difference scheme is stable whenever the analogous Cartesian scheme is stable. The scheme has been used to calculate dry atmospheric convection due to differential heating between top and bottom of mountainous terrain. The general small-scale characteristics of mountain up-slope winds have been simulated. In addition, the results have demonstrated the crucial role played by the eddy diffusivities and the environmental stability, in determining both the quantitative and the qualitative features of the circulation.

Gal-Chen, T.↗
NASA NTRS1982

Finite area method for nonlinear conical flows

A fully conservative finite area method for the computation of steady inviscid flow about general conical bodies at incidence is described. The procedure utilizes the potential approximation and implements a body conforming mesh generator. The conical potential is assumed to have its best linear variation inside each mesh cell and a secondary interlocking cell system is used to establish the flux balance required to conserve mass. In the supersonic regions the scheme is desymmetrized by adding appropriate artificial viscosity in conservation form. The algorithm is nearly an order of a magnitude faster than present Euler methods and predicts known results accurately and qualitative features such as nodal point lift off correctly. Results are compared with those of other investigations.

Sritharan, S. S.↗
DOE OSTI2020

Flow Z-pinch plasma production on the FuZE experiment

The Fusion Z Pinch Experiment (FuZE) investigates sheared-flow stabilization of classic m = 0 and m = 1 instabilities in Z pinches with an embedded axial flow. FuZE consists of a 100 cm coaxial plasma accelerator, where neutral gas is ionized and accelerated in a pulsed electrical discharge, followed by a 50 cm assembly region, where pinches are formed. Maintaining the pinch requires continuous plasma injection provided by a deflagration mode in the coaxial accelerator. Two discharge modes, with and without deflagration, are investigated on FuZE. Pinch formation is observed with deflagration only. Plasma velocities in the assembly region are found to match the E × B velocity estimated in the accelerator based on a 1D circuit model, indicating that a 1D MHD approximation may offer a valid description of the plasma in the accelerator channel. The velocity of magnetic field propagation is found to agree with the snowplow model based on momentum conservation, and the lifetime of the pinch is shown to be in agreement with constraints imposed by mass conservation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗