SEARCH · Engineering Papers
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.
LANSCE is primarily focused on three science questions for NNSA and contributes to several other programs. We will primarily consider the NNSA questions: (1) How can we advance our understanding of dynamic material behavior using focused experiments? (2) What tools are needed to address Advanced Manufacturing and Aging? (3) How do we constrain the nuclear reaction networks involved in weapons?
Soils store more carbon than the atmosphere and vegetation combined, yet Earth system models still struggle to predict how this vast reservoir will respond to environmental change. A central limitation is that most soil biogeochemical models represent organic matter using bulk conceptual pools or chemically homogeneous fractions, preventing direct use of rapidly expanding molecular-scale datasets. Here we develop and test a new soil decomposition framework that explicitly integrates high-resolution information on organic matter composition. First, we construct a molecularly informed litter decomposition module in which plant inputs are partitioned into five functional compound classes—carbohydrates, proteins, lignin-like aromatics, lipids, and carbonyls—using a molecular mixing model calibrated to solid-state 13 C Nuclear Magnetic Resonance (NMR) spectra. Class-specific kinetics, lignin-dependent physical protection, and substrate-driven microbial carbon use efficiency allow the module to capture metabolic tradeoffs associated with enzyme production and nutrient limitation. We then embed this litter module within a microbially explicit whole-soil model that tracks the transformation of these compound classes through particulate organic matter, dissolved organic matter, mineral-associated organic matter, and microbial biomass. High-resolution Fourier Transform Ion Cyclotron Resonance mass spectrometry (FTICR-MS) data are used to link internal pools to measurable soil organic matter fractions and to constrain key process parameters. Applications at soil-core and ecosystem scales demonstrate that the new model reproduces observed soil respiration dynamics while providing mechanistic attribution of CO 2 fluxes to specific chemical classes and pools. Compared to existing frameworks such as the Community Land Model soil biogeochemistry module and the Millennial model, our approach maintains competitive predictive skill while substantially improving interpretability and opportunities for data–model integration. This work illustrates a viable pathway for leveraging molecular-scale observations to reduce structural uncertainty in soil carbon–climate feedback projections.
The Heliostat Consortium (HelioCon) Wind Load Subtask was initiated with the aim of bringing research work pertaining to wind load measurement, characterization, and prediction taking place across several tasks, including Advanced Manufacturing, Components and Controls, and Field Deployment. The cross-cutting wind load subtopic in the HelioCon roadmap report [1] highlighted standardized methods and tools that are needed for a more detailed understanding of the static and dynamic loads on a heliostat. This will enable cost reduction of wind-dependent heliostat components to avoid unnecessarily conservative, overly constrained designs and increase field efficiency/reliability, to reduce the risk of component failures due to high-wind events (>15 m/s). Gaps related to heliostat wind load include site characterization for wind measurements, critical load cases for heliostat design, turbulence impacts on heliostat tracking error, testing of heliostat array configurations, understanding spatial variation of maximum loads across the solar field, and heliostat field layout and operating strategies. The recommended highest-priority pathway as first steps taken by HelioCon to address these gaps are to develop site characterization guidelines for heliostat design and field load measurements.
Mathematical model, designed for zero gravity conditions, analyzes dynamic effects of large amplitude fluid motion interior to a rigid body. It has two advantages over other mathematical models: (1) constrains slosh motion to given region in natural manner, and (2) allows equilibrium position of slosh mass to be anywhere on slosh surface.
This observational program (in conjunction with R. B. Tully (IfA, Honolulu), and J. Bland (Rice U., Houston)) aims to constrain the kinematic organization and dominant excitation mechanisms of ionized gas in active galaxies. More generally, researchers are interested in the dynamics of radiative, supersonic flows in the Interstellar Medium (ISM). Imaging Fabry-Perot interferometers and low-noise Charge Coupled Devices (CCDs) are used for complete spatial coverage of the complex gas distribution in circumnuclear narrow-line regions (NLRs). Extranuclear emission line widths in NLRs can exceed 3000 km s(-1), so to avoid inter-order confusion researchers use an etalon of 4000 km s(-1) free spectral range to map (N II) lambda lambda 6548, 6583 and H alpha. To maximize spatial resolution, researchers select nearby active systems independent of luminosity but known to possess interesting morphologies and/or high-velocity extranuclear ionized gas. Monochromatic images Full Width Half Maximum (FWHM) approx. 65 km s(-1) have thus far been obtained in 1 second or better seeing at the U. Hawaii 2.2m, CFH 3.6m, and CTIO 4.0m telescopes. These are stacked into grids of line profiles, of spectrophotometric quality, at sub-arcsecond increments across a 3 second field. To handle the approx. 20,000 to 300,000 useful spectra that arise from a typical night's work, researchers have developed a complete analysis and reduction package for VAX and Sun image workstations. Reduction involves parametrization of approx. 10 to the 8th raw data points to a few maps (e.g., velocities of each kinematic subsystem, continuum-free line fluxes) containing approx. 10 to the 5th pixels. Researchers identify kinematic and structural symmetries by examining these maps and the point to point variations of the synthesized line profiles. The combination of monochromatic images and full spatial sampling of line profiles has allowed them to isolate such symmetries and has led to reliable kinematic deprojections.
The terrestrial biosphere absorbs about 20% of fossil-fuel CO2 emissions. The overall magnitude of this sink is constrained by the difference between emissions, the rate of increase in atmospheric CO2 concentrations, and the ocean sink. However, the land sink is actually composed of two largely counteracting fluxes that are poorly quantified: fluxes from land-use change andCO2 uptake by terrestrial ecosystems. Dynamic global vegetation model simulations suggest that CO2 emissions from land-use change have been substantially underestimated because processes such as tree harvesting and land clearing from shifting cultivation have not been considered. As the overall terrestrial sink is constrained, a larger net flux as a result of land-use change implies that terrestrial uptake of CO2 is also larger, and that terrestrial ecosystems might have greater potential to sequester carbon in the future. Consequently, reforestation projects and efforts to avoid further deforestation could represent important mitigation pathways, with co-benefits for biodiversity. It is unclear whether a larger land carbon sink can be reconciled with our current understanding of terrestrial carbon cycling. Our possible underestimation of the historical residual terrestrial carbon sink adds further uncertainty to our capacity to predict the future of terrestrial carbon uptake and losses.
The Aeronautics Research Mission Directorate (ARMD) sent out an NASA Research Announcement (NRA) for proposals soliciting research and technical development. The proposed research program was aimed at addressing the desired milestones and outcomes of ROA (ROA-2006) Subtopic A.4.1.1 Advanced Computational Methods. The second milestone, SUP.1.06.02 Robust, validated mesh adaptation and error quantification for near field Computational Fluid Dynamics (CFD), was addressed by the proposed research. Additional research utilizing the direct links to geometry through a CAD interface enabled by this work will allow for geometric constraints to be applied and address the final milestone, SUP2.07.06 Constrained low-drag supersonic aerodynamic design capability. The original product of the proposed research program was an integrated system of tools that can be used for the mesh mechanics required for rapid high fidelity analysis and for design of supersonic cruise vehicles. These Euler and Navier-Stokes volume grid manipulation tools were proposed to efficiently use parallel processing. The mesh adaptation provides a systematic approach for achieving demonstrated levels of accuracy in the solutions. NASA chose to fund only the mesh generation/adaptation portion of the proposal. So this report describes the completion of the proposed tasks for mesh creation, manipulation and adaptation as it pertains to sonic boom prediction of supersonic configurations.
Clay minerals and salts in the Gale crater stratigraphy investigated by the Mars Science Laboratory (MSL) Curiosity rover preserve information about sediment sources, depositional environment, and post-depositional alteration history. The compositions of amorphous components in samples have bearing on the redox potential, pH, and chemistry of alteration fluids, and on the heterogeneity of these conditions. Interlayer clay hydration in present-day Gale crater could be an indication of high water/rock ratios of late-stage fluid alteration events, which increases the potential habitability of these post-depositional environments. The MSL Dynamic Albedo of Neutrons instrument (DAN) measures bulk hydration and here we attempt to determine the hydration states of both clay minerals and amorphous components throughout the Curiosity traverse to constrain the abundance of individual amorphous phases and determine if clay interlayer water is present today. These implications are also more broadly applicable to similar sedimentary successions outside Gale crater. In this study, we compare hydration results from DAN to mineralogical results from the MSL Chemistry and Mineralogy instrument (CheMin) for 16 different fluvio-lacustrine sample locations. These samples typically contain one or more hydrated crystalline phases, a hydrated X-ray amorphous fraction composed of multiple amorphous phases, and one or more clay mineral phases which CheMin has consistently measured as collapsed (i.e., no interlayer water). However, Gale samples are known to dehydrate during sample acquisition and handling, and because DAN measures hydration to a depth of tens of cm, we can test for the presence of clay interlayer water in situ by comparing DAN bulk hydration to CheMin mineralogy. We can also constrain the hydration of the samples’ amorphous fractions, and with geochemical data from the ChemCam and APXS instruments, this will allow us to better define the phases which compose this fraction.
ABSTRACT Algal‐fungal symbionts share water, nutrients, and gases via an architecture unique to lichens. Because lichen activity is controlled by moisture dynamics, understanding water transport is prerequisite to understand their fundamental biology. We propose a model of water distributions within foliose lichens governed by laws of fluid motion. Our model differentiates between water stored in symbionts, on extracellular surfaces, and in distinct morphological layers. We parameterize our model with hydraulic properties inverted from laboratory measurements of Flavoparmelia caperata and validate for wetting and drying. We ask: (1) Where is the bottleneck to water transport? (2) How do hydration and dehydration dynamics differ? and (3) What causes these differences? Resistance to vapor flow is concentrated at thallus surfaces and acts as the bottleneck for equilibrium, while internal resistances are small. The model captures hysteresis in hydration and desiccation, which are shown to be controlled by nonlinearities in hydraulic capacitance. Muting existing nonlinearities slowed drying and accelerated wetting, while exaggerating nonlinearities accelerated drying and slowed wetting. The hydraulic nonlinearity of F. caperata is considerable, which may reflect its preference for humid and stable environments. The model establishes the physical foundation for future investigations of transport of water, gas, and sugar between symbionts.
An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.
The paper presents a modified model reference adaptive control (M-MRAC) for multi-input multi-output nonlinear dynamical systems with time varying parametric uncertainties and bounded external disturbances. It uses a prediction model to rapidly generate adaptive estimates of the system's uncertainties with adjustable errors that converge to a small neighborhood of the origin. A sufficient condition is derived to specify the region of attraction in the space of initialization errors, design parameters and external commends. The approach is applied to thrust controlled multi-rotor air vehicles operating in an urban environment. It is shown that the designed controller can provide a good tracking of a given trajectory in the unknown urban wind field, assuming that the maximum thrust generated by the rotors is known. The performance of the algorithms are demonstrated in simulations.
Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.
Due to the high computational cost associated with unsteady aeroelastic analysis, state-of-the-art aeroelastic optimizations based on computational fluid dynamics typically ignore critical constraints like flutter and aeroelastic gust response. The linearized frequency-domain method offers an approach for adding high-fidelity flutter constraints to multidisciplinary optimizations at relatively low cost compared to other unsteady computational fluid dynamics methods. In recent work, sensitivities have been implemented for the linearized frequency-domain method in FUN3D. In this work, the linearized frequency-domain method and associated sensitivities are applied to gradient-based aeroelastic optimization with flutter constraints based on computational fluid dynamics. An overview of the flutter constraint formulation and implementation is provided, and then two optimization problems are presented. The first optimization increases the flutter speed of a pitch and plunge airfoil at transonic conditions using the minimal amount of geometric changes. The second optimization minimizes the mass of the AGARD 445.6 wing subject to a flutter constraint.
Dynamic optimization problems arise in many applications including flow control, full waveform inversion, and medical imaging. These problems are plagued by significant computational challenges. One such challenge — and the focus of this work — is the memory limitation induced by the size of the underlying dynamical system. In particular, the entire dynamic trajectory is required for derivative computation and therefore must be stored or recomputed using, e.g., checkpointing. Although recent work demonstrated the use of adaptive randomized sketching to overcome the memory challenge, that work only applies to smooth unconstrained problems, prohibiting its use for nonsmooth regularized and constrained problems. The inclusion of nonsmooth regularizers and constraints is critical as they often arise in an attempt to preserve certain physical properties or to promote sparsity. To solve these problems, we introduce a trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function that leverages randomized sketching to compress the dynamical system trajectories and adaptively adjust the sketch rank to satisfy a gradient inexactness condition. We prove convergence of this algorithm and demonstrate that it achieves substantial memory reduction on three discretized PDE-constrained optimization applications.
Abstract. Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.
The primary purpose of the first mission utilizing a Mariner Mark II spacecraft will be to rendezvous with the short-period comet Wild 2 and conduct a proximal examination during the comet's quiescent and active phases. The purpose of the present study is to elucidate, using an accurate dynamic model, some of the characteristics of the interplanetary trajectory that are relevant to the mission design. Optimal transfer trajectories are generated, constrained by the Shuttle/Centaur launch vehicle capability and the total spacecraft mass and performance. Asteroid flyby opportunities are examined, and the effects of selected asteroid encounters on mission performance are discussed.
This paper presents an extension of a FEM formulation based on variational principles. A general formulation for handling internal boundary conditions and discontinuities in the state equations is presented, and the general formulation is modified for optimal control problems subject to state-variable inequality constraints. Solutions which only touch the state constraint and solutions which have a boundary arc of finite length are considered. Suitable shape and test functions are chosen for a FEM discretization. All element quadrature (equivalent to one-point Gaussian quadrature over each element) may be done in closed form. The final form of the algebraic equations is then derived. A simple state-constrained problem is solved. Then, for a practical application of the use of the FEM formulation, a launch vehicle subject to a dynamic pressure constraint (a first-order state inequality constraint) is solved. The results presented for the launch-vehicle trajectory have some interesting features, including a touch-point solution.
Advantages of flight at high angles of attack include increased maneuverability and lift capabilities. These are beneficial not only for fighter aircraft, but also for future supersonic and hypersonic transport aircraft during take-off and landing. At high angles of attack the aerodynamics of the vehicle are dominated by separation, vortex shedding and possibly vortex breakdown. These phenomena severely compromise the effectiveness of conventional control surfaces. As a result, controlled flight at high angles of attack is not feasible for current aircraft configurations. Alternate means to augment the control of the vehicle at these flight regimes are therefore necessary. The present work investigates the augmentation of an aircraft flight control system by the injection of a thin sheet of air tangentially to the forebody of the vehicle. This method, known as Forebody Tangential Blowing (FTB), has been proposed as an effective means of increasing the controllability of aircraft at high angles of attack. The idea is based on the fact that a small amount of air is sufficient to change the separation lines on the forebody. As a consequence, the strength and position of the vortices are altered causing a change on the aerodynamic loads. Although a very effective actuator, forebody tangential blowing is also highly non-linear which makes its use for aircraft control very difficult. In this work, the feasibility of using FTB to control the roll-yaw motion of a wind tunnel model was demonstrated both through simulations and experimentally. The wind tunnel model used in the experiments consists of a wing-body configuration incorporating a delta wing with 70-degree sweep angle and a cone-cylinder fuselage. The model is equipped with forebody slots through which blowing is applied. There are no movable control surfaces, therefore blowing is the only form of actuation. Experiments were conducted at a nominal angle of attack of 45 degrees. A unique apparatus that constrains the model to two degrees-of-freedom, roll and yaw, was designed and built. The apparatus was used to conduct dynamic experiments which showed that the system was unstable, its natural motion divergent. A model for the unsteady aerodynamic loads was developed based on the basic physics of the flow and results from flow visualization experiments. Parameters of the aerodynamic model were identified from experimental data. The model was validated using data from dynamic experiments. The aerodynamic model completes the equations of motion of the system which were used in the design of control laws using blowing as the only actuator. The unsteady aerodynamic model was implemented as part of the real-time vehicle control system. A control strategy using asymmetric blowing was demonstrated experimentally. A discrete vortex method was developed to help understand the main physics of the flow. The method correctly captures the interactions between forebody and wing vortices. Moreover, the trends in static loads and flow structure are correctly represented. Flow visualization results revealed the vortical structure of the flow to be asymmetric even for symmetric flight conditions. The effects of blowing, and roll and yaw angles on the flow structure were determined. It is shown that superimposing symmetric and asymmetric blowing has a linearizing effect on the actuator characteristics. Transient responses of roll and yaw moments to step input blowing were characterized, and their differences were explained based on the physical mechanisms through which these loads are generated.