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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗
NASA NTRS2010

The Atacama Cosmology Telescope: Cosmology from Galaxy Clusters Detected Via the Sunyaev-Zel'dovich Effect

We present constraints on cosmological parameters based on a sample of Sunyaev-Zel'dovich-selected galaxy clusters detected in a millimeter-wave survey by the Atacama Cosmology Telescope. The cluster sample used in this analysis consists of 9 optically-confirmed high-mass clusters comprising the high-significance end of the total cluster sample identified in 455 square degrees of sky surveyed during 2008 at 148 GHz. We focus on the most massive systems to reduce the degeneracy between unknown cluster astrophysics and cosmology derived from SZ surveys. We describe the scaling relation between cluster mass and SZ signal with a 4-parameter fit. Marginalizing over the values of the parameters in this fit with conservative priors gives (sigma)8 = 0.851 +/- 0.115 and w = -1.14 +/- 0.35 for a spatially-flat wCDM cosmological model with WMAP 7-year priors on cosmological parameters. This gives a modest improvement in statistical uncertainty over WMAP 7-year constraints alone. Fixing the scaling relation between cluster mass and SZ signal to a fiducial relation obtained from numerical simulations and calibrated by X-ray observations, we find (sigma)8 + 0.821 +/- 0.044 and w = -1.05 +/- 0.20. These results are consistent with constraints from WMAP 7 plus baryon acoustic oscillations plus type Ia supernova which give (sigma)8 = 0.802 +/- 0.038 and w = -0.98 +/- 0.053. A stacking analysis of the clusters in this sample compared to clusters simulated assuming the fiducial model also shows good agreement. These results suggest that, given the sample of clusters used here, both the astrophysics of massive clusters and the cosmological parameters derived from them are broadly consistent with current models.

Sehgal, Neelima↗
DOE OSTI2021

Search for charged-lepton-flavour violation in Z-boson decays with the ATLAS detector

Leptons with essentially the same properties apart from their mass are grouped into three families (or flavours). The number of leptons of each flavour is conserved in interactions, but this is not imposed by fundamental principles. Since the formulation of the standard model of particle physics, the observation of flavour oscillations among neutrinos has shown that lepton flavour is not conserved in neutrino weak interactions. So far, there has been no experimental evidence that this also occurs in interactions between charged leptons. Such an observation would be a sign of undiscovered particles or a yet unknown type of interaction. Here the ATLAS experiment at the Large Hadron Collider at CERN reports a constraint on lepton-flavour-violating effects in weak interactions, searching for Z -boson decays into a τ lepton and another lepton of different flavour with opposite electric charge. The branching fractions for these decays are measured to be less than 8.1 × 10 –6 ( eτ ) and 9.5 × 10 –6 ( μτ ) at the 95% confidence level using 139 fb –1 of proton–proton collision data at a centre-of-mass energy of $\sqrt{s}=13$ TeV and 20.3 fb –1 at $\sqrt{s}=8$ TeV. These results supersede the limits from the Large Electron–Positron Collider experiments conducted more than two decades ago.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
NASA NTRS1979

Numerical solution of transonic full stream function equations in conservation form

The stream function equation in conservation form is solved iteratively based on the artificial compressibility method. The density is not a unique function of the mass flux. In order to avoid the ambiguity near the sonic line, the density is updated in terms of the velocity, which is obtained through a simple integration of a first order equation step by step in the flow field. Iteration algorithms and finite difference approximations are discussed and numerical results of both conservative and nonconservative calculations are presented.

Hafez, M. M.↗
NASA NTRS1999

Monitoring Global Geophysical Fluids by Space Geodesy

Since its establishment on 1/1/1998 by the International Earth Rotation Service, the Coordinating Center for Monitoring Global Geophysical Fluids (MGGF) and its seven Special Bureaus have engaged in an effort to support and facilitate the understanding of the geophysical fluids in global geodynamics research. Mass transports in the atmosphere-hydrosphere-solid Earth-core system (the "global geophysical fluids") will cause the following geodynamic effects on a broad time scale: (1) variations in the solid Earth's rotation (in length-of-day and polar motion/nutation) via the conservation of angular momentum and effected by torques at the fluid-solid Earth interface; (2) changes in the global gravitational field according to Newton's gravitational law; and (3) motion in the center of mass of the solid Earth relative to that of the whole Earth ("geocenter") via the conservation of linear momentum. These minute signals have become observable by space geodetic techniques, primarily VLBI, SLR, GPS, and DORIS, with ever increasing precision/accuracy and temporal/spatial resolution. Each of the seven Special Bureaus within MGGF is responsible for calculations related to a specific Earth component or aspect -- Atmosphere, Ocean, Hydrology, Ocean Tides, Mantle, Core, and Gravity/Geocenter. Angular momenta and torques, gravitational coefficients, and geocenter shift will be computed for geophysical fluids based on global observational data, and from state-of-the-art models, some of which assimilate such data. The computed quantities, algorithm and data formats are standardized. The results are archived and made available to the scientific research community. This paper reports the status of the MGGF activities and current results.

Chao, Benjamin F.↗
NASA NTRS1997

Modelling the Centers of Galaxies

The key to studying central regions by means of nobody numerical experiments is to concentrate on the central few parsecs of a galaxy, replacing the remainder of the galaxy by a suitable boundary condition, rather after the manner in which stellar interiors can be studied without a detailed stellar atmosphere by replacing the atmosphere with a boundary condition. Replacements must be carefully designed because the long range gravitational force means that the core region is sensitive to mass outside that region and because particles can exchange between the outer galaxy and the core region. We use periodic boundary conditions, coupled with an iterative procedure to generate initial particle loads in isothermal equilibrium. Angular momentum conservation is ensured for problems including systematic rotation by a circular reflecting boundary and by integrating in a frame that rotates with the mean flow. Mass beyond the boundary contributes to the gravitational potential, but does not participate in the dynamics. A symplectic integration scheme has been developed for rotating coordinate systems. This combination works well, leading to robust configurations. Some preliminary results with this combination show that: (1) Rotating systems are extremely sensitive to non-axisymmetric external potentials, and (2) that a second core, orbiting near the main core (like the M31 second core system), shows extremely rapid orbital decay. The experimental setups will be discussed, along with preliminary results.

Smith, B. F.↗
DOE OSTI2023

Search for direct production of winos and higgsinos in events with two same-charge leptons or three leptons in $pp$ collision data at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for supersymmetry targeting the direct production of winos and higgsinos is conducted in final states with either two leptons (e or μ) with the same electric charge, or three leptons. The analysis uses 139 fb -1 of pp collision data at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector during Run 2 of the Large Hadron Collider. No significant excess over the Standard Model expectation is observed. Simplified and complete models with and without R-parity conservation are considered. In topologies with intermediate states including either Wh or WZ pairs, wino masses up to 525 GeV and 250 GeV are excluded, respectively, for a bino of vanishing mass. Higgsino masses smaller than 440 GeV are excluded in a natural R-parity-violating model with bilinear terms. Upper limits on the production cross section of generic events beyond the Standard Model as low as 40 ab are obtained in signal regions optimised for these models and also for an R-parity-violating scenario with baryon-number-violating higgsino decays into top quarks and jets. The analysis significantly improves sensitivity to supersymmetric models and other processes beyond the Standard Model that may contribute to the considered final states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Towards real-time forecasting of natural gas production by harnessing graph theory for stochastic discrete fracture networks

In this work, we compare hydrocarbon production curves obtained from a graph-based reduced-order model with the high-fidelity Discrete Fracture Network (DFN) predictions for a fracture network created using data from a real shale site. We observe that the bounds for the high fidelity DFN model lie within the bounds for the reduced order model, implying that the reduced-order model provides a conservative estimate. Moreover, we found that except for first-passage times and late arriving mass, the production curves from the reduced-order model predict transport accurately. However, it is to be noted that the results are inspite of trading a three-dimensional geometry for a reduced system in the form of a graph, one that is 500–1000 times faster in terms of computational efficiency (for this particular application). In addition, we also compare the production curves for large drawdown and small drawdown using our graph approach. The reduced-order model is successful in showing that the long term productivity is higher in case of small drawdown although the initial productivity is higher for large drawdown. Thus, this reduced-order model offers great potential in uncertainty quantification for production, as well as in providing operators with information to make real-time decisions for optimal production.

03 NATURAL GAS↗
DOE OSTI2022

Stability of exact solutions of the (2 + 1)-dimensional nonlinear Schrödinger equation with arbitrary nonlinearity parameter κ

In this work, we consider the nonlinear Schrödinger equation (NLSE) in 2+1 dimensions with arbitrary nonlinearity exponent κ in the presence of an external confining potential. Exact solutions to the system are constructed, and their stability as we increase the 'mass' (i.e., the L 2 norm) and the nonlinearity parameter κ is explored. Here we observe both theoretically and numerically that the presence of the confining potential leads to wider domains of stability over the parameter space compared to the unconfined case. Our analysis suggests the existence of a stable regime of solutions for all κ as long as their mass is less than a critical value M*(κ). Furthermore, we find that there are two different critical masses, one corresponding to width perturbations and the other one to translational perturbations. The results of Derrick's theorem are also obtained by studying the small amplitude regime of a four-parameter collective coordinate (4CC) approximation. A numerical stability analysis of the NLSE shows that the instability curve M*(κ) versus κ lies below the two curves found by Derrick's theorem and the 4CC approximation. In the absence of the external potential, κ = 1 demarcates the separation between the blowup regime and the stable regime. In this 4CC approximation, for κ < 1, when the mass is above the critical mass for the translational instability, quite complicated motions of the collective coordinates are possible. Energy conservation prevents the blowup of the solution as well as confines the center of the solution to a finite spatial domain. We call this regime the 'frustrated' blowup regime and give some illustrations. In an appendix, we show how to extend these results to arbitrary initial ground state solution data and arbitrary spatial dimension d.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Consistency of the pion form factor and unpolarized transverse momentum dependent parton distributions beyond leading twist in the light-front quark model

We investigate the interplay among the pion’s form factor, transverse momentum dependent distributions (TMDs), and parton distribution functions (PDFs) extending our light-front quark model (LFQM) computation based on the Bakamjian-Thomas construction for the two-point function [1,2] to the three-point and four-point functions. Ensuring the four-momentum conservation at the meson-quark vertex from the Bakamjian-Thomas construction, the meson mass is taken consistently as the corresponding invariant meson mass both in the matrix element and the Lorentz factor in our LFQM computation. We achieve the current-component independence in the physical observables such as the pion form factor and delve into the derivation of unpolarized TMDs and PDFs associated with the forward matrix element. We address the challenges posed by twist-4 TMDs and exhibit the fulfillment of the sum rule. Effectively, our LFQM successfully handles the light-front zero modes and offers insights for broader three-point and four-point functions and related observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

The modelling of latitudinal gradients in the solar wind in the outer solar system

A steady, axisymmetric, quasi-radial, global model is developed for thermally driven stellar winds with embedded magnetic fields. The asymptotic, linear results are presented for 0(1) latitudinal variations in the radial magnetic field, mass-loss rate, and radial velocity of the wind. The MHD equations are solved for the latitudinal dependence of the rotational velocity and magnetic field. They are driven by the meridional flows that develop naturally from internal magnetic stresses. Most flows open flux tubes in the stars equatorial plane, redistributing mass and magnetic flux as a function of stellar latitude. The plasma spins up to conserve angular momentum in fields and plasma.

Nerney, S.↗
NASA NTRS1985

Modeling the effects of latitudinal gradients in stellar winds, with application to the solar wind

A steady, axisymmetric, quasi-radial, global model previously developed for stellar winds with embedded magnetic fields has been extended to include latitudinal gradient effects on the azimuthal velocity and magnetic field. The linear results at large radii are presented for large-amplitude latitudinal variations in the radial magnetic field, mass loss rate, and radial velocity of the wind. The magnetohydrodynamic (MHD) equations predict meridional flows that develop naturally from internal magnetic stresses. The flows open flux tubes in the star's equatorial plane, redistributing mass and magnetic flux as a function of stellar latitude. The plasma spins up to conserve angular momentum in fields and plasma. The results are generally applicable to stellar winds (including radiatively driven winds), provided that the internal structure is not dominated by rotation. The asymptotic solutions do not explicitly depend on the form of the energy equation, although the assumed O(1) state which drives these solutions depends on the deposition of energy and momentum throughout the wind.

Nerney, S.↗
NASA NTRS1988

Molecular outflows in the Monoceros OB1 molecular cloud

Observations of J = 1-0 emission from CO in nine suspected molecular outflows in the Monoceros OB1 molecular cloud are presented. It is found that, if the five sources which are confirmed to be outflows conserve momentum as they evolve, they will sweep up at least 0.6 percent of the mass of the entire cloud before coming into pressure equilibrium with the ambient gas. This number indicates that it should take at most 160 episodes of similar outflow activity in order to sweep up the bulk of the Mon OB1 cloud to highly supersonic speeds.

Margulis, Michael↗
NASA NTRS1990

Implicit solution of three-dimensional internal turbulent flows

The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.

Michelassi, V.↗
NASA NTRS2007

Role of Coronal Mass Ejections in the Heliospheric Hale Cycle

The 11-year solar cycle variation in the heliospheric magnetic field strength can be explained by the temporary buildup of closed flux released by coronal mass ejections (CMEs). If this explanation is correct, and the total open magnetic flux is conserved, then the interplanetary-CME closed flux must eventually open via reconnection with open flux close to the Sun. In this case each CME will move the reconnected open flux by at least the CME footpoint separation distance. Since the polarity of CME footpoints tends to follow a pattern similar to the Hale cycle of sunspot polarity, repeated CME eruption and subsequent reconnection will naturally result in latitudinal transport of open solar flux. We demonstrate how this process can reverse the coronal and heliospheric fields, and we calculate that the amount of flux involved is sufficient to accomplish the reversal within the 11 years of the solar cycle.

Owens, M. J.↗
NASA NTRS2019

Using Power Diagrams to Build Optimal Unstructured Meshes for C-Grid Models

The Model for Prediction Across Scales (MPAS) for Ocean (-O), Sea-Ice (-SI) and Land-Ice (-LI), in addition to the Coastal Ocean Marine Prediction Across Scales (COMPAS) are two novel general circulation models designed to resolve coupled ocean-ice dynamics over variable spatial scales using non-uniform unstructured grids. Both models are based on a conservative mimetic finite-difference/volume formulation (TRiSK), in which staggered momentum, vorticity and mass-based degrees- of-freedom are distributed over an orthogonal 'primal-dual' mesh.

Engwirda, Darren↗