Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗
NASA NTRS2003

Flame Synthesis Of Single-Walled Carbon Nanotubes And Nanofibers

Carbon nanotubes are widely sought for a variety of applications including gas storage, intercalation media, catalyst support and composite reinforcing material [1]. Each of these applications will require large scale quantities of CNTs. A second consideration is that some of these applications may require redispersal of the collected CNTs and attachment to a support structure. If the CNTs could be synthesized directly upon the support to be used in the end application, a tremendous savings in post-synthesis processing could be realized. Therein we have pursued both aerosol and supported catalyst synthesis of CNTs. Given space limitations, only the aerosol portion of the work is outlined here though results from both thrusts will be presented during the talk. Aerosol methods of SWNT, MWNT or nanofiber synthesis hold promise of large-scale production to supply the tonnage quantities these applications will require. Aerosol methods may potentially permit control of the catalyst particle size, offer continuous processing, provide highest product purity and most importantly, are scaleable. Only via economy of scale will the cost of CNTs be sufficient to realize the large-scale structural and power applications on both earth and in space. Present aerosol methods for SWNT synthesis include laser ablation of composite metalgraphite targets or thermal decomposition/pyrolysis of a sublimed or vaporized organometallic [2]. Both approaches, conducted within a high temperature furnace, have produced single-walled nanotubes (SWNTs). The former method requires sophisticated hardware and is inherently limited by the energy deposition that can be realized using pulsed laser light. The latter method, using expensive organometallics is difficult to control for SWNT synthesis given a range of gasparticle mixing conditions along variable temperature gradients; multi-walled nanotubes (MWNTs) are a far more likely end products. Both approaches require large energy expenditures and produce CNTs at prohibitive costs, around $500 per gram. Moreover these approaches do not possess demonstrated scalability. In contrast to these approaches, flame synthesis can be a very energy efficient, low-cost process [3]; a portion of the fuel serves as the heating source while the remainder serves as reactant. Moreover, flame systems are geometrically versatile as illustrated by innumerable boiler and furnace designs. Addressing scalability, flame systems are commercially used for producing megatonnage quantities of carbon black [4]. Although it presents a complex chemically reacting flow, a flame also offers many variables for control, e.g. temperature, chemical environment and residence times [5]. Despite these advantages, there are challenges to scaling flame synthesis as well.

Wal, Randy L. Vander↗
NASA NTRS2012

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.↗
NASA NTRS2011

A Starshade Petal Error Budget for Exo-Earth Detection and Characterization

We present a starshade error budget with engineering requirements that are well within the current manufacturing and metrology capabilities. The error budget is based on an observational scenario in which the starshade spins about its axis on timescales short relative to the zodi-limited integration time, typically several hours. The scatter from localized petal errors is smoothed into annuli around the center of the image plane, resulting in a large reduction in the background flux variation while reducing thermal gradients caused by structural shadowing. Having identified the performance sensitivity to petal shape errors with spatial periods of 3-4 cycles/petal as the most challenging aspect of the design, we have adopted and modeled a manufacturing approach that mitigates these perturbations with 1-meter-long precision edge segments positioned using commercial metrology that readily meets assembly requirements. We have performed detailed thermal modeling and show that the expected thermal deformations are well within the requirements as well. We compare the requirements for four cases: a 32 meter diameter starshade with a 1.5 meter telescope, analyzed at 75 and 90 milliarcseconds, and a 40 meter diameter starshade with a 4 meter telescope, analyzed at 60 and 75 milliarcseconds.

error budget↗
NASA NTRS2019

A Mixed Integer Efficient Global Optimization Algorithm with Multiple Infill Strategy - Applied to a Wing Topology Optimization Problem

With the advancement in high performance computing and numerical optimization techniques,engineering design optimization problems are becoming more complex, larger scale,higher fidelity, and computationally more demanding, requiring longer run times than ever before. There exists methodologies and techniques that can address some of these challenges but very few can address all, and most are limited in the extent that these concerns can be addressed. With the goal of addressing such challenging engineering problems, we developed anew optimization framework, named AMIEGO, that combines concepts from surrogate-based optimization approaches, gradient-based numerical methods, Partial Least Squares, evolutionary algorithms, and Branch-and-Bound, providing newer capabilities that were not previouslyperceived. However, the original version of this framework, in the process of adaptive samplingto explore and exploit the design space, finds only a single sample point per iteration. The efforthere builds upon this previously developed optimization framework to include multiple infillsampling capability that combines the concept of generalized expected improvement function,unsupervised learning, and multi-objective evolutionary technique. To demonstrate, AMIEGOwith the multiple infill capability (called AMIEGO-MIMOS) solves a series of increasingly difficultengineering design optimization problems. The results reveal the performance of the newapproach is problem dependent. When applied to a ten-bar truss problem, the newly proposedmultiple infill strategy consistently leads to a better design solutions when compared to theexisting CPTV method (implemented with the context of the AMIEGO framework). On theother hand, when applied to a mixed-integer high fidelity wing topology optimization problem- MIMOS, despite showing a steeper convergence at the start, eventually leads to an inferiorsolution as compared to CPTV approach. These results also reveal that a small number ofstarting points, in general, are sufficient to lead to a good overall solution.

Mixed-integer optimization↗
DOE OSTI2024

A Straightforward Model for Quantifying Local pH Gradients Governing the Oxygen Evolution Reaction

The production and consumption of protons by an electrocatalyst will, under certain conditions, generate localized microenvironments with properties distinct from those of the bulk solution. These local properties are particularly impactful for reactions involving proton-coupled electron transfer, where the generation of locally basic or acidic environments may significantly influence the energy efficiency and reaction selectivity of the electrocatalyst. Whereas local pH environments have been observed and characterized in reductive half-reactions, including the CO 2 reduction and hydrogen evolution reactions, the incompatibility of conventional techniques and materials has limited studies in oxidative half-reactions, including the oxygen evolution reaction (OER), which provides the reducing equivalents for solar-to-fuels electrolysis. With the straightforward parameters bulk pH, buffer composition and pK a , and mass transport, we develop a model for describing local pH as a function of current density regardless of the microscopic details of the mechanism. Using an acid-stable PbO x OER catalyst, we observe the formation and dissipation of pH gradients during the OER and validate the model with voltammetric and potentiometric studies. Here, the model predicts how local acidic environments can develop over a narrow OER current density window, thus providing further motivation for the development of OER catalysts that are stable to acid, even when operating in basic aqueous conditions. More generally, the model is not restricted to the OER and is useful for determining the onset of local pH gradients for other electrocatalytic reactions that involve the consumption or generation of protons in energy conversion reactions.

Anions↗
DOE OSTI2019

Optimization of a Method To Quantify Soil Bacterial Abundance by Flow Cytometry

ABSTRACT Bacterial abundance is a fundamental metric for understanding the population dynamics of soil bacteria and their role in biogeochemical cycles. Despite its importance, methodological constraints hamper our ability to assess bacterial abundance in terrestrial environments. Here, we aimed to optimize the use of flow cytometry (FCM) to assay bacterial abundances in soil while providing a rigorous quantification of its limitations. Soil samples were spiked with Escherichia coli to evaluate the levels of recovery efficiency among three extraction approaches. The optimized method added a surfactant (a tetrasodium pyrophosphate [TSP] buffer) to 0.1 g of soil, applied an intermediate degree of agitation through shaking, and used a Nycodenz density gradient to separate the cells from background debris. This procedure resulted in a high (average, 89%) level of cell recovery. Recovery efficiencies did not differ significantly among sites across an elevation gradient but were positively correlated with percent carbon in the soil samples. Estimated abundances were also highly repeatable between technical replicates. The method was applied to samples from two field studies and, in both cases, was sensitive enough to detect treatment and site differences in bacterial abundances. We conclude that FCM offers a fast and sensitive method to assay soil bacterial abundance from relatively small amounts of soil. Further work is needed to assay differential biases of the method across a wider range of soil types. IMPORTANCE The ability to quantify bacterial abundance is important for understanding the contributions of microbial communities in soils, but such assays remain difficult and time-consuming. Flow cytometry offers a fast and direct way to count bacterial cells, but several concerns remain in applying the technique to soils. This study aimed to improve the efficiency of the method for soil while quantifying its limitations. We demonstrated that an optimized procedure was sensitive enough to capture differences in bacterial abundances among treatments and ecosystems in two field studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Results for aliovalent doping of CeBr{sub 3} with Ca{sup 2+}

Despite the outstanding scintillation performance characteristics of cerium tribromide (CeBr{sub 3}) and cerium-activated lanthanum tribromide, their commercial availability and application are limited due to the difficulties of growing large, crack-free single crystals from these fragile materials. This investigation employed aliovalent doping to increase crystal strength while maintaining the optical properties of the crystal. One divalent dopant (Ca{sup 2+}) was used as a dopant to strengthen CeBr{sub 3} without negatively impacting scintillation performance. Ingots containing nominal concentrations of 1.9% of the Ca{sup 2+} dopant were grown, i.e., 1.9% of the CeBr{sub 3} molecules were replaced by CaBr{sub 2} molecules, to match our target replacement of 1 out of 54 cerium atoms be replaced by a calcium atom. Precisely the mixture was composed of 2.26 g of CaBr{sub 2} added to 222.14 g of CeBr{sub 3}. Preliminary scintillation measurements are presented for this aliovalently doped scintillator. Ca{sup 2+}-doped CeBr{sub 3} exhibited little or no change in the peak fluorescence emission for 371 nm optical excitation for CeBr{sub 3}. The structural, electronic, and optical properties of CeBr{sub 3} crystals were studied using the density functional theory within the generalized gradient approximation. Calculated lattice parameters are in agreement with the experimental data. The energy band structures and density of states were obtained. The optical properties of CeBr{sub 3}, including the dielectric function, were calculated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Origin of Unusual Acidity and Li+ Diffusivity in a Series of Water-in-Salt Electrolytes.

Superconcentrated aqueous electrolytes ("water-in-salt" electrolytes, or WiSEs) enable various aqueous battery chemistries beyond the voltage limits imposed by the Pourbaix diagram of water. However, their detailed structural and transport properties remain unexplored and could be better understood through added studies. Here, we report on our observations of strong acidity (pH 2.4) induced by lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) at superconcentration (at 20 mol/kg). Multiple nuclear magnetic resonance (NMR) and pulsed-field gradient (PFG) diffusion NMR experiments, density functional theory (DFT) calculations, and molecular dynamics (MD) simulations reveal that such acidity originates from the formation of nanometric ion-rich structures. The experimental and simulation results indicate the separation of water-rich and ion-rich domains at salt concentrations >= 5 m and the acidity arising therefrom is due to deprotonation of water molecules in the ion-rich domains. As such, the ion-rich domain is composed of hydrophobic -CF3 (of TFSI-) and hydrophilic hydroxyl (OH-) groups. At 20 m concentration, the tortuosity and radius of water diffusion channels are estimated to be similar to 10 and similar to 1 nm, respectively, which are close to values obtained from hydrated Nafion membranes that also have hydrophobic polytetrafluoroethylene (PTFE) backbones and hydrophilic channels consisting of SO3- ion cluster networks providing for the transport of ions and water. Thus, we have discovered the structural similarity between WiSE and hydrated Nafion membranes on the nanometer scale.

Han, Kee Sung↗
DOE OSTI2023

Dynamic subcanopy leaf traits drive resistance of net primary production across a disturbance severity gradient

Across the globe, the forest carbon sink is increasingly vulnerable to an expanding array of low- to moderate-severity disturbances. However, some forest ecosystems exhibit functional resistance (i.e., the capacity of ecosystems to continue functioning as usual) following disturbances such as extreme weather events and insect or fungal pathogen outbreaks. Unlike severe disturbances (e.g., stand-replacing wildfires), moderate severity disturbances do not always result in near-term declines in forest production because of the potential for compensatory growth, including enhanced subcanopy production. Community-wide shifts in subcanopy plant functional traits, prompted by disturbance-driven environmental change, may play a key mechanistic role in resisting declines in net primary production (NPP) up to thresholds of canopy loss. However, the temporal dynamics of these shifts, as well as the upper limits of disturbance for which subcanopy production can compensate, remain poorly characterized. In this study, we leverage a 4-year dataset from an experimental forest disturbance in northern Michigan to assess subcanopy community trait shifts as well as their utility in predicting ecosystem NPP resistance across a wide range of implemented disturbance severities. Through mechanical girdling of stems, we achieved a gradient of severity from 0% (i.e., control) to 45, 65, and 85% targeted gross canopy defoliation, replicated across four landscape ecosystems broadly representative of the Upper Great Lakes ecoregion. We found that three of four examined subcanopy community weighted mean (CWM) traits including leaf photosynthetic rate (p = 0.04), stomatal conductance (p = 0.07), and the red edge normalized difference vegetation index (p < 0.0001) shifted rapidly following disturbance but before widespread changes in subcanopy light environment triggered by canopy tree mortality. Surprisingly, stimulated subcanopy production fully compensated for upper canopy losses across our gradient of experimental severities, achieving complete resistance (i.e., no significant interannual differences from control) of whole ecosystem NPP even in the 85% disturbance treatment. Additionally, we identified a probable mechanistic switch from nutrient-driven to light-driven trait shifts as disturbance progressed. Our findings suggest that remotely sensed traits such as the red edge normalized difference vegetation index (reNDVI) could be particularly sensitive and robust predictors of production response to disturbance, even across compositionally diverse forests. The potential of leaf spectral indices to predict post-disturbance functional resistance is promising given the capabilities of airborne to satellite remote sensing. We conclude that dynamic functional trait shifts following disturbance can be used to predict production response across a wide range of disturbance severities.

resistance↗
DOE OSTI2025

Double-Crucible Vertical Bridgman Technique for Stoichiometry-Controlled Chalcogenide Crystal Growth

Precise stoichiometry control in single-crystal growth is essential for both technological applications and fundamental research. However, conventional growth methods often face challenges such as non-stoichiometry, compositional gradients, and phase impurities, particularly in non-congruent melting systems. Even in congruent melting systems like Bi₂Se₃, deviations from the ideal stoichiometric composition can lead to significant property degradation, such as excessive bulk conductivity, which limits its topological applications. In this study, we introduce the double-crucible vertical Bridgman (DCVB) method, a novel approach that enhances stoichiometry control through the combined use of continuous source material feeding, traveling-solvent growth, and liquid encapsulation, which suppresses volatile element loss under high pressure. Using Bi₂Se₃ as a model system, we demonstrate that crystals grown via DCVB exhibit enhanced stoichiometric control, significantly reducing defect density and achieving much lower carrier concentrations compared to those produced by conventional Bridgman techniques. Moreover, the continuous feeding of source material enables the growth of large crystals. As a result, this approach presents a promising strategy for synthesizing high-quality, large-scale crystals, particularly for metal chalcogenides and pnictides that exhibit challenging non-congruent melting behaviors.

Crystal structure↗
DOE OSTI2012

Imaging x-ray Thomson scattering spectrometer design and demonstration (invited)

In many laboratory astrophysics experiments, intense laser irradiation creates novel material conditions with large, one-dimensional gradients in the temperature, density, and ionization state. X-ray Thomson scattering is a powerful technique for measuring these plasma parameters. However, the scattered signal has previously been measured with little or no spatial resolution, which limits the ability to diagnose inhomogeneous plasmas. We report on the development of a new imaging x-ray Thomson spectrometer (IXTS) for the Omega laser facility. The diffraction of x-rays from a toroidally curved crystal creates high-resolution images that are spatially resolved along a one-dimensional profile while spectrally dispersing the radiation. This focusing geometry allows for high brightness while localizing noise sources and improving the linearity of the dispersion. Preliminary results are presented from a scattering experiment that used the IXTS to measure the temperature profile of a shocked carbon foam.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A Trajectory-Driven Algorithm for Differentiating SRB Measures on Unstable Manifolds

Sinai-Ruelle-Bowen (SRB) measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory- driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. Here, we numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Thermoelastic analysis of solar window of a vacuum chamber

It has been noted during space simulation test programs that the chamber penetration window becomes the limiting factor as the solar intensity increases. The temperature at the center of the window normally attains a level of 550 C, while the edge temperature remains at about 35 C. The stresses produced by this thermal gradient combined with the atmospheric pressure are studied within the framework of classical theory of elasticity. The solution for the problem is obtained by solving the static uncoupled thermoelastic equations for a circular blank, subjected to nonuniform heat on one face, by introduction of two auxiliary variables. The displacements are expressed in terms of Bessel functions and parameters of interest are computed at various mesh points of the window. The maximum principal stresses thus obtained show that even with steep temperature gradient, the safety factor is acceptable.

Amba-Rao, C. L.↗
NASA NTRS1990

Experimental and numerical study of the British Experimental Rotor Programme blade

Wind-tunnel tests on the British Experimental Rotor Programme (BERP) tip are described, and the results are compared with computational fluid dynamics (CFD) results. The test model was molded using the Lynx-BERP blade tooling to provide a semispan, cantilever wing comprising the outboard 30 percent of the rotor blade. The tests included both surface-pressure measurements and flow visualization to obtain detailed information of the flow over the BERP tip for a range of angles of attack. It was observed that, outboard of the notch, favorable pressure gradients exist which ensure attached flow, and that the tip vortex also remains stable to large angles of attack. On the rotor, these features yield a very gradual break in control loads when the retreating-blade limit is eventually reached. Computational and experimental results were generally found to be in good agreement.

Brocklehurst, Alan↗
NASA NTRS1991

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗
NASA NTRS1994

Ecosystem and physiological controls over methane production in northern wetlands

Peat chemistry appears to exert primary control over methane production rates in the Canadian Northern Wetlands Study (NOWES) area. We determined laboratory methane production rate potentials in anaerobic slurries of samples collected from a transect of sites through the NOWES study area. We related methane production rates to indicators of resistance to microbial decay (peat C: N and lignin: N ratios) and experimentally manipulated substrate availability for methanogenesis using ethanol (EtOH) and plant litter. We also determined responses of methane production to pH and temperature. Methane production potentials declined along the gradient of sites from high rates in the coastal fens to low rates in the interior bogs and were generally highest in surface layers. Strong relationships between CH4 production potentials and peat chemistry suggested that methanogenesis was limited by fermentation rates. Methane production at ambient pH responded strongly to substrate additions in the circumneutral fens with narrow lignin: N and C: N ratios (delta CH4/delta EtOH = 0.9-2.3 mg/g) and weakly in the acidic bogs with wide C: N and lignin: N ratios (delta CH4/delta EtOH = -0.04-0.02 mg/g). Observed Q(sub 10) values ranged from 1.7 to 4.7 and generally increased with increasing substrate availability, suggesting that fermentation rates were limiting. Titration experiments generally demonstrated inhibition of methanogenesis by low pH. Our results suggest that the low rates of methane emission observed in interior bogs during NOWES likely resulted from pH and substrate quality limitation of the fermentation step in methane production and thus reflect intrinsically low methane production potentials. Low methane emission rates observed during NOWES will likely be observed in other northern wetland regions with similar vegetation chemistry.

Valentine, David W.↗
DOE OSTI2016

Revisiting Fenton Hill Phase I reservoir creation and stimulation mechanisms through the GTO code comparison effort

A series of experiments performed at the Fenton Hill hot dry rock site after stage 2 drilling of Phase I reservoir provided intriguing field observations on the reservoir’s responses to injection and venting under various conditions. Two teams participating in the US DOE Geothermal Technologies Office (GTO)’s Code Comparison Study (CCS) used different numerical codes to model these five experiments with the objective of inferring the hydraulic stimulation mechanism involved. The codes used by the two teams are based on different numerical principles, and the assumptions made were also different, due to intrinsic limitations in the codes and the modelers’ personal interpretations of the field observations. Both sets of models were able to produce the most important field observations and both found that it was the combination of the vertical gradient of the fracture opening pressure, injection volume, and the use/absence of proppant that yielded the different outcomes of the five experiments.

Fu, Pengcheng↗