Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conservative variables”

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 127 records · Page 7

Application of a finite element algorithm for the solution of steady transonic Euler equations

A finite element algorithm for the solution of two-dimensional, steady Euler equations is presented which, through a Clebsch-type transformation for the velocity vector, solves the conservation of mass equation with one primary variable and two additional equations for the convection of two new variables. The accuracy and efficiency of this scheme is discussed, and the second-order accuracy attained in the analysis of the convection of the vorticity is demonstrated together with the efficient treatment of rotational and irrotational flow subregions. A sample problem is used to show the accuracy of the numerical scheme and its convergence characteristics.

Akay, H. U.↗

Modelling and measurement of crack closure and crack growth following overloads and underloads

Ignoring crack growth retardation following overloads can result in overly conservative life predictions in structures subjected to variable amplitude fatigue loading. Crack closure is believed to contribute to the crack growth retardation, although the specific closure mechanism is dabatable. The delay period and corresponding crack growth rate transients following overload and overload/underload cycles were systematically measured as a function of load ratio and overload magnitude. These responses are correlated in terms of the local 'driving force' for crack growth, i.e. the effective stress intensity factor range. Experimental results are compared with the predictions of a Dugdale-type (1960) crack closure model, and improvements in the model are suggested.

Dexter, R. J.↗

A petabyte size electronic library using the N-Gram memory engine

A model library containing petabytes of data is proposed by Triada, Ltd., Ann Arbor, Michigan. The library uses the newly patented N-Gram Memory Engine (Neurex), for storage, compression, and retrieval. Neurex splits data into two parts: a hierarchical network of associative memories that store 'information' from data and a permutation operator that preserves sequence. Neurex is expected to offer four advantages in mass storage systems. Neurex representations are dense, fully reversible, hence less expensive to store. Neurex becomes exponentially more stable with increasing data flow; thus its contents and the inverting algorithm may be mass produced for low cost distribution. Only a small permutation operator would be recalled from the library to recover data. Neurex may be enhanced to recall patterns using a partial pattern. Neurex nodes are measures of their pattern. Researchers might use nodes in statistical models to avoid costly sorting and counting procedures. Neurex subsumes a theory of learning and memory that the author believes extends information theory. Its first axiom is a symmetry principle: learning creates memory and memory evidences learning. The theory treats an information store that evolves from a null state to stationarity. A Neurex extracts information data without a priori knowledge; i.e., unlike neural networks, neither feedback nor training is required. The model consists of an energetically conservative field of uniformly distributed events with variable spatial and temporal scale, and an observer walking randomly through this field. A bank of band limited transducers (an 'eye'), each transducer in a bank being tuned to a sub-band, outputs signals upon registering events. Output signals are 'observed' by another transducer bank (a mid-brain), except the band limit of the second bank is narrower than the band limit of the first bank. The banks are arrayed as n 'levels' or 'time domains, td.' The banks are the hierarchical network (a cortex) and transducers are (associative) memories. A model Neurex was built and studied. Data were 50 MB to 10 GB samples of text, data base, and images: black/white, grey scale, and high resolution in several spectral bands. Memories at td, S(m(sub td)), were plotted against outputs of memories at td-1. S(m(sub td)) was Boltzman distributed, and memory frequencies exhibited self-organized criticality (SOC); i.e., 'l/f(sup beta)' after long exposures to data. Whereas output signals from level n may be encoded with B(sub output) = O(-log(2)f(sup beta)) bits, and input data encoded with B(sub input) = O((S(td)/S(td-1))(sup n)), B(sup output)/B(sub input) is much less than 1 always, the Neurex determines a canonical code for data and it is a lossless data compressor. Further tests are underway to confirm these results with more data types and larger samples.

Bugajski, Joseph M.↗

Blue Carbon Stocks Along the Pacific Coast of North America Are Mainly Driven by Local Rather Than Regional Factors

Coastal wetlands, including seagrass meadows, emergent marshes, mangroves, and temperate tidal swamps, can efficiently sequester and store large quantities of sediment organic carbon (SOC). However, SOC stocks may vary by ecosystem type and along environmental or climate gradients at different scales. Quantifying such variability is needed to improve blue carbon accounting, conservation effectiveness, and restoration planning. We analyzed SOC stocks in 1,284 sediment cores along >6,500 km of the Pacific coast of North America that included large environmental gradients and multiple ecosystem types. Tidal wetlands with woody vegetation (mangroves and swamps) had the highest mean stocks to 1 m depth (357 and 355 Mg ha −1 , respectively), 45% higher than marshes (245 Mg ha −1 ), and more than 500% higher than seagrass (68 Mg ha −1 ). Unvegetated tideflats, though not often considered a blue carbon ecosystem, had noteworthy stocks (148 Mg ha −1 ). Stocks increased with tidal elevation and with fine (<63 μm) sediment content in several ecosystems. Stocks also varied by dominant plant species within individual ecosystem types. At larger scales, marsh stocks were lowest in the Sonoran Desert region of Mexico, and swamp stocks differed among climate zones; otherwise stocks showed little correlation with ecoregion or latitude. More variability in SOC occurred among ecosystem types, and at smaller spatial scales (such as individual estuaries), than across regional climate gradients. These patterns can inform coastal conservation and restoration priorities across scales where preserving stored carbon and enhancing sequestration helps avert greenhouse gas emissions and maintains other vital ecosystem services.

emergent marsh↗

The MOOSE fluid properties module

The fluid properties module in MOOSE~\cite{lindsay2022moose} serves a variety of fluid simulation applications based on MOOSE, including the MOOSE Navier Stokes module~\cite{moose_ns}, Pronghorn~\cite{pgh}, the MOOSE thermal hydraulics module, SAM~\cite{sam}, RELAP-7~\cite{relap7}, Sockeye~\cite{sockeye}, Pronghorn-subchannel~\cite{subchannel} and the MOOSE porous flow module~\cite{porous}. These applications are used to solve coarse mesh multi-dimensional thermal-hydraulics~\cite{pgh}, 1D systems analysis~\cite{sam,relap7} in nuclear reactor analysis, heat pipe modeling~\cite{sockeye} and porous flow simulations~\cite{porous} for underground gas storage and water seepage. The use of consistent fluid properties across fluid flow applications facilitates coupled flow simulations~\cite{anl_sam_pgh,osti_1889653}. Each application has historically driven the implementation of several fluid properties, which were later extended to be compatible with other applications. The unique diversity of applications of the module, due to its presence in MOOSE, has driven its expansion to new fluids, such as advanced nuclear reactor coolants and, more recently, arbitrary functions or tables-based property definitions, as detailed in section~\ref{content}, as well as numerous thermophysical properties and variable sets, as detailed in subsection~\ref{sec:prop}. The need for different discretizations of flow equations based on the compressibility of the fluid has motivated support for both primitive (pressure- and temperature-based) and conservative (internal energy- and specific volume-based) flow variables in the module. Thermodynamic relations are used to convert between these two formulations, as needed. The module serves a dual purpose of providing fluid properties for direct use in flow simulations and facilitating the implementation of user-specific fluid properties. Contributions of new properties for existing fluids or new fluids are strongly encouraged.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Adaptive control with variable dead-zone nonlinearities

It has been found that fixed error dead-zones as defined in the existing literature result in serious degradation of performance, due to the conservativeness which characterizes the determination of their width. In the present paper, variable width dead-zones are derived for the adaptive control of plants with unmodeled dynamics. The derivation makes use of information available about the unmodeled dynamics both a priori as well as during the adaptation process, so as to stabilize the adaptive loop and at the same time overcome the conservativeness and performance limitations of fixed-dead zone adaptive or fixed gain controllers.

Orlicki, D.↗

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗

Critical habitat identification of peripheral Sage Thrashers under climate change

Abstract The Sage Thrasher ( Oreoscoptes montanus ) has been assessed as “Endangered” in Canada since 1992. Like other species with a geographic range that barely extends into Canada, Sage Thrashers are rare. Thirty‐one percent of Canada's bird species listed for recovery under the Canadian Species at Risk Act (SARA) are at the periphery of their range. A listing of “endangered” under SARA requires identification of critical habitat for the species. With anticipated climate change, recovery of species requires a more proactive intervention than relying on historical occurrence to locate suitable habitat. We synthesized 19 years of Sage Thrasher occurrence and related habitat data across the species' northern range in British Columbia (BC) and Washington (WA) to define critical habitat characteristics. We found Sage Thrashers selected less leaf litter and less grass cover in flat or low‐slope regions farther from anthropogenic or natural habitat breaks; habitat sensitive to the expected climate change impacts of fire, changes in precipitation, and invasive species establishment. By augmenting the BC data collected in the species' peripheral range with data from their core distribution in WA, we identified key habitat elements of an otherwise data‐poor species that do not breed in sufficient numbers in Canada to reliably characterize their habitat. These methods improve the identification of “critical habitat” for peripheral species like Sage Thrashers in preparation for climate‐induced range expansion northward. The framework developed demonstrates a useful template for conservation strategies for data‐limited peripheral populations in other regions. Focusing on the landscape‐level variables that indicate good habitat, and not the locations of habitat, can identify suitable future areas for conservation.

Millikin, Rhonda L.↗

Vibration Reduction of Helicopter Blade Using Variable Dampers: A Feasibility Study

In the report, the investigation of controlling helicopter-blade lead-lag vibration is described. Current practice of adding passive damping may be improved to handle large dynamic range of the blade with several peaks of vibration resonance. To minimize extra-large damping forces that may damage the control system of blade, passive dampers should have relatively small damping coefficients, which in turn limit the effectiveness. By providing variable damping, a much larger damping coefficient to suppress the vibration can be realized. If the damping force reaches the maximum allowed threshold, the damper will be automatically switched into the mode with smaller damping coefficient to maintain near-constant damping force. Furthermore, the proposed control system will also have a fail-safe feature to guarantee the basic performation of a typical passive damper. The proposed control strategy to avoid resonant regions in the frequency domain is to generate variable damping force in combination with the supporting stiffness to manipulate the restoring force and conservative energy of the controlled blade system. Two control algorithms are developed and verified by a prototype variable damper, a digital controller and corresponding algorithms. Primary experiments show good potentials for the proposed variable damper: about 66% and 82% reductions in displacement at 1/3 length and the root of the blade respectively.

Lee, George C.↗

The impact of wind, solar, and other factors on the decline in wholesale power prices in the United States

Across multiple organized wholesale power markets in the United States, annual average prices declined by 19–64 dollars per MWh between 2008 and 2017 while retirements of thermal power plants accelerated. Several prominent changes over the last decade are often discussed as contributors to this decline in prices. These include growth in wind and solar, a reduction in the price of natural gas, and weakened load growth. Here we construct a fundamental supply curve model for each of seven organized wholesale market regions and use counterfactual simulations to assess the degree to which wind and solar—among other factors—have influenced wholesale electricity prices. We find that growth in wind and solar since 2008 reduced average annual wholesale electricity prices by less than 3 dollars per MWh. In contrast the decline in natural gas prices reduced wholesale prices by 7–53 dollars per MWh, depending on the region. This suggests that recent thermal-plant retirements in the U.S. are primarily due to low natural gas prices, not growth in wind and solar. Fully isolating the impact of individual factors, however, is limited by non-linear interactions between factors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Structural and kinetic characterization of an acetoacetyl-Coenzyme A: acetate Coenzyme A transferase from the extreme thermophile Thermosipho melanesiensis

Family 1 Coenzyme A transferases (CtfAB) from the extremely thermophilic bacterium, Thermosipho melanesiensis, has been used for in vivo acetone production up to 70°C. This enzyme has tentatively been identified as the rate-limiting step, due to its relatively low-binding affinity for acetate. However, existing kinetic and mechanistic studies on this enzyme are insufficient to evaluate this hypothesis. Here, kinetic analysis of purified recombinant T. melanesiensis CtfAB showed that it has a ping-pong bi-bi mechanism typical of Coenzyme A (CoA) transferases with Km values for acetate and acetoacetyl-CoA of 85 mM and 135 μM, respectively. Product inhibition by acetyl-CoA was competitive with respect to acetoacetyl-CoA and non-competitive with respect to acetate. Crystal structures of wild-type and mutant T. melanesiensis CtfAB were solved in the presence of acetate and in the presence or absence of acetyl-CoA. These structures led to a proposed structural basis for the competitive and non-competitive inhibition of acetyl-CoA: acetate binds independently of acetyl-CoA in an apparent low-affinity binding pocket in CtfA that is directly adjacent to a catalytic glutamate in CtfB. Similar to other CoA transferases, acetyl-CoA is bound in an apparent high-affinity binding site in CtfB with most interactions occurring between the phospho-ADP of CoA and CtfB residues far from the acetate binding pocket. This structural-based mechanism also explains the organic acid promiscuity of CtfAB. High-affinity interactions are predominantly between the conserved phospho-ADP of CoA, and the variable organic acid binding site is a low-affinity binding site with few specific interactions.

Biochemistry & Molecular Biology↗

Multiphysics Modeling Framework for Soft PVC Gel Sensors with Experimental Comparisons

Polyvinyl chloride (PVC) gels have recently been found to exhibit mechanoelectrical transduction or sensing capabilities under compressive loading applications. This phenomenon is not wholly understood but has been characterized as an adsorption-like phenomena under varying amounts and types of plasticizers. A different polymer lattice structure has also been tested, thermoplastic polyurethane, which showed similar sensing characteristics. This study examines mechanical and electrical properties of these gel sensors and proposes a mathematical framework of the underlying mechanisms of mechanoelectrical transduction. COMSOL Multiphysics is used to show solid mechanics characteristics, electrostatic properties, and transport of interstitial plasticizer under compressive loading applications. The solid mechanics takes a continuum mechanics approach and includes a highly compressive Storakers material model for compressive loading applications. The electrostatics and transport properties include charge conservation and a Langmuir adsorption migration model with variable diffusion properties based on plasticizer properties. Results show both plasticizer concentration gradient as well as expected voltage response under varying amounts and types of plasticizers. Experimental work is also completed to show agreeance with the modeling results.

42 ENGINEERING↗

A rapid implicit-explicit solution to the two-dimensional time dependent incompressible Navier-Stokes equations

A second-order time-accurate and spatially factored algorithm was used in a finite difference scheme for the numerical solution of the time-dependent, incompressible, two dimensional Navier-Stokes equations in conservation-law form using vorticity and stream function variables. The systems of equations are solved at each time step by an iterative technique. Numerical results were obtained for a circular cylinder at a Reynolds number of 15, and an NACA 0012 airfoil at zero angle of attack at Reynolds numbers of 10 to the third and 10 to the fourth powers. The results are in agreement with another numerical technique, and the computing time required to obtain the steady state solution at the Reynolds number of 10 to the 4th power was 49.7 sec on CDC 7600 computer using a 65 x 84 computational grind.

Davis, J. E.↗

The evolution of AAOE observed constituents with the polar vortex

One of the difficulties in determining constituent trends from the ER-2 flight data is the large amount of day to day variability generated by the motion of the polar vortex. To reduce this variability, the observations have been transformed into the conservative (Lagrangian) reference frames consisting of the coordinate pairs, potential temperature (PT) and potential vorticity (PV), or PT and N2O. The requirement of only two independent coordinates rests on the assumption that constituent distributions and their chemical processes are nearly zonal in that coordinate system. Flight data is used everywhere for these transformation except for potential vorticity. Potential vorticity is determined from level flight segments, and NMC PV values during flight dives and takeoffs are combined with flight data in a smooth fashion.

Schoeberl, Mark R.↗

A guide to the use of theoretical models of the solar nebula for the interpretation of the meteoritic record

Attempts to derive a theoretical framework for the interpretation of the meteoritic record have been frustrated by our incomplete understanding of the fundamental processes that controlled the evolution of the primitive solar nebula. Nevertheless, it is possible to develop qualitative models of the nebula that illuminate its dynamic character, as well as the roles of some key parameters. These models draw on the growing body of observational data on the properties of disks around young, solar-type stars, and are constructed by applying the results of known solutions of protostellar collapse problems; making simple assumptions about the radial variations of nebular variables; and imposing the integral constraints demanded by conservation of mass, angular momentum, and energy. The models so constructed are heuristic, rather than predictive; they are intended to help us think about the nebula in realistic ways, but they cannot provide a definitive description of conditions in the nebula.

Cassen, Pat↗

Volume 2: Explicit, multistage upwind schemes for Euler and Navier-Stokes equations

The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using monotonic upstream schemes for conservation laws (MUSCL)-type differencing to obtain state variables at cell interface. Variable interpolations were written in the k-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator. Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-state predictor-corrector schemes, and multistage time-stepping schemes. The coefficients of the multistage time-stepping schemes have been modified successfully to achieve better performance with upwind differencing. A technique was developed to optimize the coefficients for good high-frequency damping at relatively high CFL numbers. Local time-stepping, implicit residual smoothing, and multigrid procedure were added to the explicit time stepping scheme to accelerate convergence to steady-state. The developed algorithm was implemented successfully in a multi-block code, which provides complete topological and geometric flexibility. The only requirement is C degree continuity of the grid across the block interface. The algorithm has been validated on a diverse set of three-dimensional test cases of increasing complexity. The cases studied were: (1) supersonic corner flow; (2) supersonic plume flow; (3) laminar and turbulent flow over a flat plate; (4) transonic flow over an ONERA M6 wing; and (5) unsteady flow of a compressible jet impinging on a ground plane (with and without cross flow). The emphasis of the test cases was validation of code, and assessment of performance, as well as demonstration of flexibility.

Elmiligui, Alaa↗

Heat Transfer from Radiatively Heated Material in a Low Reynolds Number Microgravity Environment

A mathematical model of the transient three-dimensional heat transfer between a slowly moving ambient gas stream and a thermally thick or thin flat surface heated by external radiation in a microgravity environment is presented. The problem is motivated in part by fire safety issues in spacecraft. The gas phase is represented by variable property convection-diffusion energy and mass conservation equations valid at low Reynolds numbers. The absence of gravity and low Reynolds number together permit the flow to be represented by a self-consistent velocity potential determined by the ambient velocity and the thermal expansion in the gas. The solid exchanges energy with the gas by conduction/convection and with the surroundings by surface absorption and re-emission of radiation. Heat conduction in the solid is assumed to be one dimensional at each point on the surface as a consequence of the limited times (of order of 10 seconds) of interest in these simulations. Despite the apparent simplicity of the model, the results show a complex thermally induced flow near the heated surface. The thermal exchange between the gas and solid produces an outward sourcelike flow upstream of the center of the irradiated area and a sinklike flow downstream. The responses of the temperature fields and the associated flows to changes in the intensity of the external radiation and the ambient velocity are discussed.

Yamashita, H.↗

Three Dimensional Flow and Pressure Patterns in a Single Pocket of a Hydrostatic Journal Bearing

The flow in a hydrostatic pocket is described by a mathematical model that uses the three dimensional Navier-Stokes equations written in terms of the primary variables, u, v, w, and p. Using a conservative formulation, a finite volume multi-block method is applied through a collocated, body fitted grid. The flow is simulated in a shallow pocket with a depth/length ratio of 0.02. The flow structures obtained and described by the authors in their previous two dimensional models are made visible in their three dimensional aspect for the Couette flow. It has been found that the flow regimes formed central and secondary vortical cells with three dimensional corkscrew-like structures that lead the fluid on an outward bound path in the axial direction of the pocket. The position of the central vortical cell center is at the exit region of the capillary restrictor feedline. It has also been determined that a fluid turn around zone occupies all the upstream space between the floor of the pocket and the runner, thus preventing any flow exit through the upstream port. The corresponding pressure distribution under the shaft presented as well. It was clearly established that for the Couette dominated case the pressure varies significantly in the pocket in the circumferential direction, while its variation is less pronounced axially.

Braun, M. Jack↗