Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical balances”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

autoGEMM: Pushing the Limits of Irregular Matrix Multiplication on Arm Architectures

This paper presents an open-source library that pushes the limits of performance portability for irregular General Matrix Multiplication (GEMM) on the widely-used Arm architectures. Our library, autoGEMM, is designed to support a wide range of Arm processors: from edge devices to HPC-grade CPUs. autoGEMM generates optimized kernels for various hardware configurations by auto-combining fragments of autogenerated micro-kernels that employ hand-written optimizations to maximize computational efficiency. We optimize the kernel pipeline by tuning the register reuse and the data load/store overlapping. In addition, we use a dynamic tiling scheme to generate balanced tile shapes. Finally, we position autoGEMM on top of the TVM framework where our dynamic tiling scheme prunes the search space for TVM to identify the optimal combination of parameters for code optimization. Evaluations on five different classes of Arm chips demonstrate the advantages of autoGEMM. For small matrices, autoGEMM achieves 98% of peak and up to 2.0x speedup over state-of-the-art libraries such as LIBXSMM and LibShalom. For irregular matrices (i.e. tall skinny and long rectangles), autoGEMM is 1.3-2.0x faster than widely-used libraries such as OpenBLAS and Eigen. autoGEMM is publicly available at: https://github.com/wudu98/autoGEMM.

Wu, Du↗
NASA NTRS1993

Simulations of the general circulation of the Martian atmosphere. II - Seasonal pressure variations

The CO2 seasonal cycle of the Martian atmosphere and surface is simulated with a hybrid energy balance model that incorporates dynamical and radiation information from a large number of general circulation model runs. This information includes: heating due to atmospheric heat advection, the seasonally varying ratio of the surface pressure at the two Viking landing sites to the globally averaged pressure, the rate of CO2 condensation in the atmosphere, and solar heating of the atmosphere and surface. The predictions of the energy balance model are compared with the seasonal pressure variations measured at the two Viking landing sites and the springtime retreat of the seasonal polar cap boundaries. The following quantities are found to have a strong influence on the seasonal pressures at the Viking landing sites: albedo of the seasonal CO2 ice deposits, emissivity of this deposit, atmospheric heat advection, and the pressure ratio.

Pollack, James B.↗
DOE OSTI2021

Arctic tundra shrubification: a review of mechanisms and impacts on ecosystem carbon balance

Abstract Vegetation composition shifts, and in particular, shrub expansion across the Arctic tundra are some of the most important and widely observed responses of high-latitude ecosystems to rapid climate warming. These changes in vegetation potentially alter ecosystem carbon balances by affecting a complex set of soil–plant–atmosphere interactions. In this review, we synthesize the literature on (a) observed shrub expansion, (b) key climatic and environmental controls and mechanisms that affect shrub expansion, (c) impacts of shrub expansion on ecosystem carbon balance, and (d) research gaps and future directions to improve process representations in land models. A broad range of evidence, including in-situ observations, warming experiments, and remotely sensed vegetation indices have shown increases in growth and abundance of woody plants, particularly tall deciduous shrubs, and advancing shrublines across the circumpolar Arctic. This recent shrub expansion is affected by several interacting factors including climate warming, accelerated nutrient cycling, changing disturbance regimes, and local variation in topography and hydrology. Under warmer conditions, tall deciduous shrubs can be more competitive than other plant functional types in tundra ecosystems because of their taller maximum canopy heights and often dense canopy structure. Competitive abilities of tall deciduous shrubs vs herbaceous plants are also controlled by variation in traits that affect carbon and nutrient investments and retention strategies in leaves, stems, and roots. Overall, shrub expansion may affect tundra carbon balances by enhancing ecosystem carbon uptake and altering ecosystem respiration, and through complex feedback mechanisms that affect snowpack dynamics, permafrost degradation, surface energy balance, and litter inputs. Observed and projected tall deciduous shrub expansion and the subsequent effects on surface energy and carbon balances may alter feedbacks to the climate system. Land models, including those integrated in Earth System Models, need to account for differences in plant traits that control competitive interactions to accurately predict decadal- to centennial-scale tundra vegetation and carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Allometric relationships and trade‐offs in 11 common M editerranean‐climate grasses

Abstract Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant–environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co‐occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process‐based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above‐ and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade‐offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast‐growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf‐to‐stem ratio. However, contrary to our hypothesis, there were no differences in root‐to‐leaf ratio between perennial and annual or C 4 and C 3 plants. Our study provides species‐specific and functional‐type‐specific allometry equations for both above‐ and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade‐offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Wealth inequality in the prehispanic northern US Southwest: from Malthus to Tyche

Persistent differences in wealth and power among prehispanic Pueblo societies are visible from the late AD 800s through the late 1200s, after which large portions of the northern US Southwest were depopulated. In this paper we measure these differences in wealth using Gini coefficients based on house size, and show that high Ginis (large wealth differences) are positively related to persistence in settlements and inversely related to an annual measure of the size of the unoccupied dry-farming niche. We argue that wealth inequality in this record is due first to processes inherent in village life which have internally different distributions of the most productive maize fields, exacerbated by the dynamics of systems of balanced reciprocity; and second to decreasing ability to escape village life owing to shrinking availability of unoccupied places within the maize dry-farming niche as villages get enmeshed in regional systems of tribute or taxation. We embed this analytical reconstruction in the model of an ‘Abrupt imposition of Malthusian equilibrium in a natural-fertility, agrarian society’ proposed by Puleston et al. (Puleston C, Tuljapurkar S, Winterhalder B. 2014 PLoS ONE 9, e87541 (doi:10.1371/journal.pone.0087541)), but show that the transition to Malthusian dynamics in this area is not abrupt but extends over centuries. This article is part of the theme issue ‘Evolutionary ecology of inequality’.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1979

Measurement of dynamic direct and cross-coupling derivatives due to oscillatory roll

An oscillatory roll apparatus is described with which the complete set of moment and force derivatives due to an oscillation in roll can be obtained in a wind tunnel. A newly developed technique to reduce the dynamic data originating from the four-component balance used with the above apparatus is briefly presented. The technique is intended to account for the dynamic behavior of the model-balance subsystem in the presence of motions in the various degrees of freedom allowed by the balance webs. Finally, results obtained from both dynamic calibration as well as wind-tunnel tests are shown.

Hanff, E. S.↗
NASA NTRS2018

Shake Test Results and Dynamic Characterization of the Large Rotor Test Apparatus in the NFAC 40- by 80-Foot Wind Tunnel

Prior to the full-scale wind tunnel test of the UH-60A Airloads rotor, a shake test was completed on the Large Rotor Test Apparatus. The goal of the shake test was to characterize the oscillatory response of the test rig and provide a dynamic calibration of the balance to permit accurate measurement of vibratory hub loads. This paper provides a summary of the shake test results, including balance, shaft bending gauge, and accelerometer measurements. Sensitivity to hub mass and angle of attack were investigated during the shake test. Hub mass was found to have an important impact on the vibratory forces and moments measured at the balance, especially near the UH-60A 4/rev frequency. Comparisons were made between the accelerometer data and an existing finite-element model, showing agreement on mode shapes, but not on natural frequencies. The results suggest that the shake test data can be used to correct in-plane loads measurements up to 10 Hz and normal loads up to 30 Hz.

Shake Test↗
DOE OSTI2025

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Facile Measurement of the Rotation of a Single Optically Trapped Nanoparticle Using the Diagonal Ratio of a Quadrant Photodiode

Optical tweezers are a powerful tool for exploring physical properties of particles in various environments through analysis of their dynamics in a trapping potential. Analysis of the trapped particle's position is often done using interferometric back-focal-plane detection microscopy in which interference between the trapping laser and the light interacting with the trapped particle is projected on a quadrant photo-diode (QPD). The QPD measurements are almost universally the voltage ratios between the transverse components of the detector (left/right; up/down) that are calibrated to ascertain particle positions and orientations (i.e., rotation). In this paper, we demonstrate how measuring the diagonal ratio of the QPD allows tracking the spinning and orbital motion of optically trapped objects with enhanced sensitivity compared to the conventional transverse ratio and reveals additional information about particle dynamics. The diagonal ratio "balances" the signal from opposite sides of the QPD, allowing common mode noise reduction and increased measurement sensitivity of the rotating nanospheres at frequencies from below 200 Hz to above 15 kHz. We show that the variability in particle rotation rates is due to slight differences in the aspect ratio of the nanoparticles. We also measure in situ a gradual acceleration in the spinning frequency of trapped nanoparticles due to photothermal shedding of surface ligands and the resulting decrease of rotational friction in solution. The accelerating particles spin at rates exceeding 30 kHz and are eventually "optically printed" onto the glass coverslip on top of the liquid cell enclosure. The diagonal ratio is a simple, readily available, and powerful enabler for measuring rotational dynamics with a 5-fold sensitivity increase over conventional QPD measurements. This affords opportunities for controlling the orientational dynamics of trapped nanoparticles, which is crucial for various applications, including nanoscale mechanical motors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗
NASA NTRS2001

Parallel Processing of Adaptive Meshes with Load Balancing

Many scientific applications involve grids that lack a uniform underlying structure. These applications are often also dynamic in nature in that the grid structure significantly changes between successive phases of execution. In parallel computing environments, mesh adaptation of unstructured grids through selective refinement/coarsening has proven to be an effective approach. However, achieving load balance while minimizing interprocessor communication and redistribution costs is a difficult problem. Traditional dynamic load balancers are mostly inadequate because they lack a global view of system loads across processors. In this paper, we propose a novel and general-purpose load balancer that utilizes symmetric broadcast networks (SBN) as the underlying communication topology, and compare its performance with a successful global load balancing environment, called PLUM, specifically created to handle adaptive unstructured applications. Our experimental results on an IBM SP2 demonstrate that the SBN-based load balancer achieves lower redistribution costs than that under PLUM by overlapping processing and data migration.

Das, Sajal K.↗
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
NASA NTRS1983

Arctic Sea ice by passive microwave observations from the Nimbus-5 Satellite

The results of a dynamic/thermodynamic numerical model of Arctic sea ice are compared with satellite images from the Nimbus 5 electrically scanning microwave radiometer. The model combines aspects of two previous sea ice models those of Parkinson and Washington and Ling, Rasmussen, and Campbell. A solid/fluid model basically follows the formulation of the Parkinson and Washington model with the addition of the constitutive equation and equation of state from the Ling model. The Parkinson and Washington model simulates the seasonal cycle of sea ice thicknesses and concentrations with a horizontal resolution of roughly 200 km and a timestep of 8 hours. The thermodynamics are calculated through energy balances at the interfaces between ice and air, water and ice, and water and air. The ice dynamics are calculated through a momentum equation balancing air stress, water stress, dynamic topography, and Coriolis force, with an adjustment for internal ice resistance.

Campbell, W. J.↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
DOE OSTI2025

Distinctive features of structural evolution and thermodynamic response in wide-bandgap semiconductors driven by intense electronic excitation

Radiation-tolerant material selection requires balancing lattice rigidity, defect dynamics, and electronic stability, as shown by covalent SiC outperforming ionic Ga 2 O 3 and GaN under extreme environments. Responding to intense electronic excitation, irradiation-driven phase segregation (β → δ/κ in Ga 2 O 3 ) and core–shell track (disordered structure in GaN), accompanied by elemental redistribution, contrastingly, exceptional radiation tolerance manifested by comparatively minimal lattice distortion (0.17 % strain variation) was demonstrated in SiC. These differential responses are primarily attributed to two fundamental mechanisms: (i) thermodynamic driving forces governing defect migration and phase separation, and (ii) the synergistic effects of robust covalent bonding composition coupled with efficient defect recombination processes. Here, the stronger electron–phonon (e-ph) coupling in Ga 2 O 3 (4.34 × 1018 W m −3 K −1 ) and GaN (3.55 × 10 18 W m −3 K −1 ) enhances lattice energy deposition, triggering thermal spikes (ΔT ≫ T m ) and structural transition behaviors, whereas weaker e-ph coupling in SiC (3.69 × 10 18 W m −3 K −1 ), relatively high thermodynamic parameters and efficient energy dissipation suppress thermal spikes to maintaining lattice integrity. The photoresponse degradation driven by enhanced radiative recombination is dominant in N-doped SiC, while V-doped systems achieve defect-mediated photoconduction optimization characterized by abrupt current transitions, matching fluorescence yield evolutions, and directly connecting defect engineering to optoelectronic performance.

Intense electronic excitation↗
DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluating Methods for Measuring Grid Frequency in Low-Inertia Power Systems: Preprint

Accurate measurement of grid frequency is a critical component of reliable grid control. Traditionally, methods such as phase locked loops (PLLs) and discrete Fourier transforms (DFTs) have been used in inverters and phasor measurement units (PMUs) to measure frequency. However, as the percentage of inverter-based resources (IBRs) such as solar and wind has increased, these conventional frequency measurement methods are proving unable to guarantee reliable control in some cases. One challenge is measuring frequency during transient events, where there is a disruption in the steady state sinusoidal voltage. During these events, the underlying frequency of the grid may barely change, but measurement methods report a large spike in frequency due to the disrupted waveform. New methods must balance between suppressing spikes in frequency during faults, and providing fast, accurate, measurements in all other grid operation conditions, especially during events with high rate-of-change-of frequency (ROCOF), which are more prevalent in high-IBR power systems. This paper first surveys frequency measurement methods that have been proposed to reduce measurement errors during transient events. Then, both conventional and more novel frequency measurement methods are tested against an IEEE standard and industry recommendations, and their performance is evaluated for events simulated in PSCAD. Results quantify the trade-offs in performance during different grid conditions and lead to suggestions for the most appropriate frequency and ROCOF measurement methods for low inertia grids.

frequency↗
DOE OSTI2024

Tropical root responses to global changes: A synthesis

Tropical ecosystems face escalating global change. These shifts can disrupt tropical forests' carbon (C) balance and impact root dynamics. Since roots perform essential functions such as resource acquisition and tissue protection, root responses can inform about the strategies and vulnerabilities of ecosystems facing present and future global changes. However, root trait dynamics are poorly understood, especially in tropical ecosystems. We analyzed existing research on tropical root responses to key global change drivers: warming, drought, flooding, cyclones, nitrogen (N) deposition, elevated (e) CO 2 , and fires. Based on tree species- and community-level literature, we obtained 266 root trait observations from 93 studies across 24 tropical countries. We found differences in the proportion of root responsiveness to global change among different global change drivers but not among root categories. In particular, we observed that tropical root systems responded to warming and eCO 2 by increasing root biomass in species-scale studies. Drought increased the root: shoot ratio with no change in root biomass, indicating a decline in aboveground biomass. Despite N deposition being the most studied global change driver, it had some of the most variable effects on root characteristics, with few predictable responses. Episodic disturbances such as cyclones, fires, and flooding consistently resulted in a change in root trait expressions, with cyclones and fires increasing root production, potentially due to shifts in plant community and nutrient inputs, while flooding changed plant regulatory metabolisms due to low oxygen conditions. The data available to date clearly show that tropical forest root characteristics and dynamics are responding to global change, although in ways that are not always predictable. This synthesis indicates the need for replicated studies across root characteristics at species and community scales under different global change factors.

54 ENVIRONMENTAL SCIENCES↗