SEARCH · Engineering Papers

Results for “Differentiable predictive control”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Towards an Enhanced Droplet Activation Scheme for Multi-Moment Bulk Microphysics Schemes

Initial droplet spectra produced upon activation impact the ensuing chain of microphysical processes andtherefore play a crucial role in cloud evolution. This work re-examines dependencies of newly formed clouddroplet size distribution (CDSD) characteristics on environmental and aerosol properties via parcel model simulationsthat serve as the basis for a multi-moment bulk microphysics droplet activation scheme suitable for acloud-resolving model (CRM). It is found that applying a fixed size threshold to define activated droplets versusemploying physical considerations can lead to erroneous activation and overly broad CDSDs for high aerosolconcentration and weak updraft conditions. Aerosol distributions characterized by larger median sizes and/orincreased solubility can result in greater activated droplet numbers, whereas impacts of these parameters onCDSD spectral width depend on both aerosol number concentration and updraft velocity. An expansion of theactivation scheme to include CDSD spectral width is proposed to aid efforts to extend high-order momentprediction to cloud droplet categories in CRMs as well as better represent variability in the activation process onthe cloud scale.simulations to investigate the regime dependence of the relative dispersion(d)1 of newly activated CDSDs, where d is the ratio of dropletradius standard deviation (σ) to the mean radius (r ). C16 demonstratedthat increasing Na resulted in increasing (decreasing) d values via reducedcondensational narrowing (spectral broadening) rates within theAL (UL) regime, with d values peaking in the TR regime. Their findingssuggest a similar regime dependence for d as R09 noted for Nc and helpexplain reportedly conflicting relationships between Na and CDSDspectral characteristics (cf. Hudson and Noble, 2014; Liu et al., 2014),although the applicability of these results within bulk microphysicalschemes was not addressed.Simulating aerosol-cloud interactions with CRMs employing bulkmicrophysics requires that the scheme minimally predict two CDSDparameters, namely mass and number concentrations, and represent thedroplet activation process. Various activation schemes aim to determineNc from aerosol and environmental properties and include analyticalexpressions (e.g., Abdul-Razzak et al., 1998; Morrison et al., 2005) aswell as lookup tables (LUTs) based on detailed parcel model calculations(e.g., Saleeby and Cotton, 2004, hereafter SC04; Segal and Khain,2006; Thompson and Eidhammer, 2014). Expressions to diagnose CDSDspectral width from Nc (Grabowski, 1998; Liu et al., 2006; Morrison andGrabowski, 2007) or cloud water content (Geoffroy et al., 2010) havealso been developed, although more robust methods to obtain CDSDspectral width upon activation are presently lacking. This latter point isrelevant for triple-moment (3 M) bulk microphysics that aim to predictdistribution spectral width alongside number and mass concentrations(e.g., Loftus et al., 2014; Milbrandt and Yau, 2005).The current work extends the findings of C16 to the current LUTbasedaerosol activation scheme used in the Regional AtmosphericModeling System (RAMS) (Cotton et al., 2003; SC04; Saleeby and vanden Heever, 2013, hereafter SvdH13) and additionally examinesaerosol size and solubility impacts on newly activated CDSD properties.Because early cloud development processes such as condensationalgrowth, evaporation, and droplet self-collection depend on and impactCDSD spectral width (Hudson and Yum, 1997; Seifert and Beheng 2001;Lu and Seinfeld, 2006; Igel and van den Heever, 2017), an expansion ofthe activation LUTs to include CDSD spectral width is proposed as apreliminary step for extending 3M prediction to CDSDs in CRMs forimproved simulations of aerosol-cloud interactions.2. MethodologyThe current RAMS two-moment microphysics module determinesthe fractional number of aerosol particles that activate to cloud dropletsfrom five-dimensional LUTs based on model predicted air temperature(T), w, Na, and the geometric median radius (rg) and soluble fraction (ε)of the aerosol size distribution (SvdH13). These LUTs are created offlineusing a one-dimensional Lagrangian adiabatic parcel model (Feingoldand Heymsfield, 1992; Heymsfield and Sabin, 1989; SC04) to simulateexplicit droplet activation and initial CDSD growth for a range of ambientatmospheric conditions [T, w] and binned lognormal aerosol sizedistributions given by= ⎡⎣ ⎢− ⎤⎦ ⎥N r Nr π σr rσ( )2 lnexp[ln( / )]2(ln )aggg22 (1)where r is the dry aerosol particle bin radius and σg is the geometricstandard deviation of the distribution. As the parcel model simulationsfocus on the activation process, other processes such as coalescence,sedimentation, and mixing are not considered. Details of the parcelmodel can be found in SC04 and SvdH13, and only a brief description isprovided here. At the onset of parcel model calculations, the initiallydry aerosol particles in all bins first deliquesce and reach theirequilibrium diameters in a sub-saturated environment based on theKöhler equation for solution droplets. The parcel is then lifted at a fixedupward velocity w, and particle growth by vapor diffusion, along withconcurrent changes in the ambient environment, are iteratively computedusing the Variable-coefficient Ordinary Differential Equation(VODE) solver (Brown et al., 1989). The time resolution of these calculationsis determined within the VODE solver, and the frequency atwhich the solver is called is controlled by a longer model time stepbased on fixed upward parcel displacement increments (Δz) at thespecified w (Δt=Δz/w). Model calculations proceed until the parcelreaches a height 50m beyond the level of maximum saturation ratio(Smax) or total parcel displacement exceeds 2 km. Upon model termination,Smax and the fractional number of aerosols (factv) resulting innewly formed cloud droplets, defined as particles having diameters of atleast 2 μm, are cataloged in the LUTs according to the specified T, w, Na,rg, and ε parameter values.A critical point regarding the creation of these LUTs is the use of afixed minimum diameter (Dmin) to define cloud droplets in the parcelmodel, which can produce erroneous CDSD characteristics, particularlywithin the UL regime. For aerosol distributions with large rg valuesunder low SS conditions, for example, deliquesced aerosols within thelarge tail of the distribution can exceed 2 μm in diameter yet remainunactivated as ‘haze’ particles (Levin and Cotton, 2009; McFigganset al., 2006). For this study, aerosol particles activate to cloud dropletsbased on the critical diameter Dcrit as a function of parcel supersaturationratio (Sr) as in R09:D = σ MS RTρ83 ln( ) critsol wr w (2)where σsol is the surface tension of a solution droplet, Mw and ρw are themolar mass and density of liquid water, respectively, and R is theuniversal gas constant. Additionally, at relatively large w values withinthe AL regime, Nc stabilizes shortly after reaching supersaturation.However, parcel ascent and condensational growth continue beyondthe level of Smax, potentially causing additional narrowing of the CDSD.In the current work, model calculations terminate upon reaching Smaxas changes in Nc are negligible with continued ascent (Peng et al., 2007;R09).Parcel model simulations are performed to examine the sensitivitiesof CDSD characteristics to w, Na, rg, and ε, with the ranges for theseparameters listed in Table 1. Aerosols are assumed to be a mix of solubleand insoluble material of equal density, specified by ε, where fullysoluble aerosols correspond to ammonium sulfate with hygroscopicityparameter κ=0.61 (Petters and Kreidenweis, 2007). FollowingSvdH13, aerosol geometric standard deviation is fixed at σg=1.8, andaerosol distributions (Eq. 1) are partitioned into 100 logarithmicallyspacedbins spanning a size range specific to each rg value. For all simulations,Δz=1 m, and initial values of relative humidity, air temperatureand pressure are set to RH=0.99, T=10 °C and p=900 hPa,respectively.

Loftus, Adrian M.↗

NASA NTRS1994

Differential GPS and system integration of the Low Visibility Landing and Surface Operations (LVLASO) demonstration

The LVLASO Flight Demonstration of ASTA concepts (FDAC) integrates NASA-Langley's electronic moving map display and Transport Systems Research Vehicle (TSRV) (a modified Boeing 737 aircraft); ARINC's VHF data link, GPS ground station, and automated controller workstation; and Norden's surface radar/airport movement safety system. Aircraft location is shown on the electronic map display in the cockpit. An approved taxi route as well as other aircraft and surface traffic are also displayed. An Ashtech Z12 Global Positioning System (GPS) receiver on the TSRV estimates the aircraft's position. In Differential mode (DSPS), the Ashtech receiver accepts differential C/A code pseudorange corrections from a GPS ground station. The GPS ground station provides corrections up to ten satellites. The corrections are transmitted on a VHF data link at a 1 Hz. rate using the RTCM-104 format. DGPS position estimates will be within 5 meters of actual aircraft position. DGPS position estimates are blended with position, velocity, acceleration, and heading data from the TSRV Air Data/Inertial Reference System (ADIRS). The ADIRS data is accurate in the short-term, but drifts over time. The DGPS data is used to keep the ADIRS position accurate. Ownship position, velocity, heading, and turn rate are sent at a 20 Hz. rate to the electronic map display. Airport traffic is detected by the airport surface radar system. Aircraft and vehicles such as fuel trucks and baggage carts are detected. The traffic's location, velocity, and heading are sent to the TSRV. To prevent traffic symbology from jumping each second when a location update arrives, velocity and heading are used to predict a new traffic location for each display update. Possible runway incursions and collisions can be shown on the electronic map. Integrating the different systems used in the FDAC requires attention to the underlying coordinate systems. The airport diagram displayed on the electronic map is obtained from published navigational charts. The charts reference the North American Datum of 1927 (NAD27) or a local state-plane coordinate system. GPS uses the World Geodetic Standard of 1984 (WGS84). Both NAD27 and WGS84 model the Earth as an ellipsoid, however, they use a different origin and different size ellipsoids. Latitudes and longitudes given in these systems can be converted to a Cartesian system with the origin at the Earth's center. The surface radar detects traffic in a locally-level, rho-theta coordiante system. The electronic airport diagram is stored using a flat XY coordinate system. The map origin is at the tower and is referenced as True North up. All ownship and other traffic positions must be converted to the electronic map's frame of reference for display.

Rankin, James M.↗

DOE OSTI2023

Depletion-driven thermochemistry of molten salt reactors: review, method, and analysis

Molten salt reactors (MSRs) are innovative advanced nuclear reactors that utilize nuclear fuel by dissolving it in a high-temperature liquid salt. This unique feature differentiates MSRs from other types of reactors and allows for enhanced safety and economic performance. The liquid fuel also entails several multiphysics effects that can complicate reactor design and operation. One primary effect termed here as depletion-driven thermochemistry is a driving force in altering the multiphysics behavior of the reactor. Essentially, depletion-driven thermochemistry is the effect that fuel depletion has on changing the chemical redox potential of the fuel salt over time. As the fuel is consumed, the redox potential shifts toward a more oxidizing state. Without active control, the changing chemistry due to depletion increases corrosion thereby limiting reactor component lifetimes. Additionally, the changing redox potential of the fuel salt alters the vapor pressures of chemical species dissolved in the fuel salt. Changing vapor pressures of species in the fuel salt is an important parameter to understand when off-gassing volatile species during normal reactor operation, and for source term characterization during accident scenario transients. The present work represents a fundamental step toward modeling and coupling the driving physics (i.e., neutronics and chemistry) involved in altering the redox potential in an MSR. Here, the neutronic code Griffin models the depletion of the fuel-salt system, while the chemical equilibrium code Thermochimica calculates the thermochemical state of the isotopic inventory, using the Molten Salt Thermodynamic Database - Thermochemical (MSTDB-TC). These two codes are tightly coupled to predict the impact of fuel depletion in altering the chemistry in MSR systems. Redox potential control methods are discussed and can be modeled using this multiphysics approach. The vapor pressures of chemical species that could be extracted to an off-gas system, as determined by the reactor’s thermochemical state, are examined. The neutronics-chemistry coupling developed in this work is expected to have potential application for analyzing corrosion, source term evolution, and material safeguards in MSR systems. Lastly, suggestions for areas of further improvements of the models to expand these capabilities by incorporating other coupled physics effects is provided.

Walker, Samuel A.↗

DOE OSTI2023

Physics-informed machine-learning model of temperature evolution under solid phase processes

We model temperature dynamics during Shear Assisted Proccess Extrusion (ShAPE), a solid phase process that plasticizes feedstock with a rotating tool and subsequently extrudes it into a consolidated tube, rod, or wire. Control of temperature is critical during ShAPE processing to avoid liquefaction, ensure smooth extrusion, and develop desired material properties in the extruded products. Accurate modeling of the complicated thermo-mechanical feedbacks between process inputs, material temperature, and heat generation presents a significant barrier to predictive modeling and process design. In particular, connecting micro-structural scale mechanisms of heat generation to macro-scale predictions of temperature can become computationally intractable. In this work we use a neural network (NN) model of heat generation to bridge this gap, by combining it with a simplified model of the temperature dynamics due to conduction and convection to capture the macro scale evolution of temperature. We inform the construction of the NN heat generation model using crystal plasticity simulations at the micro-structural scale to model the effects of process inputs on generation of heat. We achieved close fits of the temperature dynamics model to a diverse experimental data-set. Further, the relationships learned by the NN model between process inputs and heat generation showed qualitative agreement with those predicted by crystal plasticity simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Review of Wave Energy Converter Power Take-Off Systems, Testing Practices, and Evaluation Metrics: Preprint

While the field of wave energy has been the subject of numerical simulation, scale model testing, and precommercial project testing for decades, wave energy technologies remain in the early stages of development and must continuing proving themselves as a promising modern renewable energy field. One of the difficulties that wave energy systems have been struggling to overcome is the design of highly efficient energy conversion systems that can convert the mechanical power, derived from the oscillation of wave activated bodies, into another useful product. Often the power take-off (PTO) is defined as the single unit responsible for converting mechanical power into another usable form such as electricity, pressurized fluid, compressed air, and others. The PTO, and the entire power conversion chain (PCC), is of great importance as it affects not only how efficient wave power is converted into electricity, but also contributes to the mass, size, structural dynamics, and levelized cost of energy (LCOE) of the wave energy converter (WEC). Unlike wind and solar, there is no industrial standard device, or devices, for wave energy conversion and this diversity is transferred to the PTO system. The majority of current WEC PTO systems incorporate a mechanical or hydraulic drive train, power generator, and an electrical control system. The challenge of WEC PTO designs is designing a mechanical-to-electrical component that can efficiently convert irregular, bi-directional, low frequency and low alternating velocity wave motions. While gross average power levels can be predicted in advance, the variable wave elevation input has to be converted into smooth electrical output and hence usually necessitates some type of energy storage system, such as battery storage, accumulator super capacitors, etc., or other means of compensation such as an array of devices. One of the primary challenges for wave energy converter systems is the fluctuating nature of wave resources, which require WEC components to be designed to handle loads (i.e. torques, forces, and powers) that are many times greater than the average load. This approach requires a much greater PTO capacity than the average power output and indicates a higher cost. In addition, supporting mechanical coupling and or gearing can be added to the PCC to help alleviate the difficulties with transmission and control of fluctuating large loads with low frequencies (indicative of wave forcing) into smaller loads at higher frequencies (optimum for conventional electrical machine design) can quickly increase the complexity of the PCC which could result in a greater number of failure modes and increased maintenance costs. All of the previous points demonstrate how the PTO influences WEC dynamics, reliability, performance and cost which are critical design factors. This paper further explores these topics by providing a review of the state-of-the-art PTO systems currently under development, how these novel PTO systems are tested and derisked prior to precommercial deployment, and the evaluation metrics historically used to differentiate between PTO designs and how they can be improved to support control co-design focused development of wave energy systems.

laboratory testing↗

DOE OSTI2024

PPI DataHub Project Data Package: S. elongatus PCC 7942 Circadian Control Bioproduction Transcriptomics (PB-DP3)

The purpose of this experiment was to evaluate how circadian clock regulation impacts carbon partitioning between storage, growth, and product synthesis in Synechococcus elongatus PCC 7942 in providing insights to strategies for enhanced bioproduction. Sample data was acquired using a Illumina HiSeq sequencer system and processed for RNA sequencing (RNA-Seq) expression analysis. Transcriptomic differential expression analysis revealed coordinated circadian clock-driven adjustment of the cell cycle and rewiring of energy and carbon metabolism. Processed RNA-Seq datasets are openly accessible from the PNNL DataHub project dataset download page and contain secondary processed RNA-seq results files and supporting metadata materials linked to relevant source code information supporting data transparency and reuse.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Defining Lipidomic Responses to Coronavirus Infection

Highly pathogenic human coronavirus infection can cause a severe atypical, rapid onset pneumonia with a mortality rate of up to 10% for severe acute respiratory syndrome coronavirus 1 (SARS-CoV 1), 35% for Middle East respiratory syndrome coronavirus (MERS-CoV), or 1% for severe acute respiratory syndrome coronavirus 2 (SARS-CoV 2 causative agent of COVID 19). While medical countermeasures successfully controlled the worldwide SARS-CoV epidemic, the MERS-CoV epidemic is still ongoing and continues to be a concern for travelers in the Middle East and the multi-year SARS-CoV 2 pandemic underscore the importance of defining biomarkers that are diagnostic and/or predictive of severe disease outcomes for respiratory viruses. Systems biology approaches provide global snapshots of infection induced changes in host cells/tissues and provide extremely rich datasets for understanding host pathogen interactions. Metabolites, especially lipids, are critical for viral replication but less is understood about infection induced changes to lipids due to limits in lipid species detection and identification. To characterize how individual lipid species and proteins with lipid associated functions contribute to highly pathogenic human coronavirus replication and disease severity, existing datasets were probed to determine cell type specific lipid responses to MERS-CoV infection and verification studies were performed to determine if modification of the host lipid signature (by inhibiting enzymatic functions that produce specific lipid species) can perturb CoV replication in human lung cells. MERS-CoV infects human lung epithelial, endothelial and fibroblast cells. All three cell types were infected with MERS-CoV and samples collected to analyze lipids, proteins, metabolites, and transcripts from 12 to 48 hours post infection. Time matched mock-infected cells were collected in parallel for each cell type. Following sample and statistical analysis, functional enrichment and bioinformatic analysis was performed to identify differentially expressed lipids and proteins. Two lipid species were found to be significantly upregulated following MERS-CoV infection, ceramides, and triacylglycerol both of which are indicative of cells undergoing apoptotic cell death. In contrast, sphingomyelins (lipid molecules that can serve as a precursor for one pathway for ceramide synthesis) had significantly decreased expression in MERS-CoV infected cells. An inhibitor of acid sphingomyelinase (that converts sphingomyelin to ceramides and phosphorylcholine) reduced MERS-CoV replication suggesting that production of ceramides is key for successful viral replication and transmission. Acyl-CoA-synthetase 3 (ACSL3), the only protein whose function is lipid associated and had increased differential expression in our dataset, regulates the synthesis of triacylglycerol (increased expression). Inhibitors that directly block ACSL3 (Triacsin C) but not steps later in the triacylglycerol synthesis pathway (Etomoxir) inhibit MERS-CoV replication suggesting that triacylglycerol production and/or ACSL3 activity is also key for viral replication. ACSL3 expression is also upregulated in lung cancer cells and is predicted to promote continued cell viability which would also enhance viral replication. The differentially expressed lipid and lipid-associated protein species identified in our studies suggest that MERS-CoV infection is activating cellular death pathways to limit the number of cells that become infected but also stimulating the production of lipid-associated enzymes that can prolong host cell viability and the amount of time progeny virions can be produced and released. As the inhibitors that worked against MERS-CoV infection were also efficacious against SARS-CoV 2 infection, countermeasures that target these host pathways may provide novel ways to block highly pathogenic human coronavirus infection and/or prevent severe disease outcomes.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2014

Analysis and Prediction of Weather Impacted Ground Stop Operations

When the air traffic demand is expected to exceed the available airport's capacity for a short period of time, Ground Stop (GS) operations are implemented by Federal Aviation Administration (FAA) Traffic Flow Management (TFM). The GS requires departing aircraft meeting specific criteria to remain on the ground to achieve reduced demands at the constrained destination airport until the end of the GS. This paper provides a high-level overview of the statistical distributions as well as causal factors for the GSs at the major airports in the United States. The GS's character, the weather impact on GSs, GS variations with delays, and the interaction between GSs and Ground Delay Programs (GDPs) at Newark Liberty International Airport (EWR) are investigated. The machine learning methods are used to generate classification models that map the historical airport weather forecast, schedule traffic, and other airport conditions to implemented GS/GDP operations and the models are evaluated using the cross-validations. This modeling approach produced promising results as it yielded an 85% overall classification accuracy to distinguish the implemented GS days from the normal days without GS and GDP operations and a 71% accuracy to differentiate the GS and GDP implemented days from the GDP only days.

Analysis↗

NASA NTRS2016

Real-time Identification and Control of Satellite Signal Impairments Solution and Application of the Stratonovich Equation Part 1. Theoretical Development

As satellite communications systems become both more complex and reliant with respect to their operating environment, it has become imperative to be able to identify, during real-time operation, the onset of one or more impairments to the quality of overall communications system integrity. One of the most important aspects to monitor of a satellite link operating within the Earth's atmosphere is the signal fading due to the occurrence of rain and/or phase scintillations. This, of course, must be done in the presence of the associated measurement uncertainty or potentially faulty measurement equipment such as in the Advanced Communication Technology Satellite (ACTS) experiment. In the present work, an approach originally suggested in 1991, and apparently still considered iconoclastic, will be significantly developed and applied to the satellite communications link on which the deleterious composite signal fade is the result of one or many component fade mechanisms. Through the measurement (with the attendant uncertainty or 'error' in the measurement) of such a composite fading satellite signal, it is desired to extract the level of each of the individual fading mechanisms so they can be appropriately mitigated before they impact the overall performance of the communications network. Rather than employing simple-minded deterministic filtering to the real-time fading, the present approach is built around all the models and/or descriptions used to describe the individual fade components, including their dynamic evolution. The latter is usually given by a first-order Langevin equation. This circumstance allows the description of the associated temporal transition probability densities of each of the component processes. By using this description, along with the real-time measurements of the composite fade (along with the measurement errors), one can obtain statistical estimates of the levels of each of the component fading mechanisms as well as their predicted values into the future. This is all accomplished by the use of the well-known Stratonovich integro-differential equation that results from the model of the measured signal fade that is also tailored to adaptively adjust the values of the parameters used in the statistical models of the individual fade mechanisms. Three examples of increasing complexity are addressed and solved for the iterative determination of fade component levels from the measured composite signal fade in the presence of measurement error and, in the last case, with uncertainty in the model parameters.

Countermeasure↗

NASA NTRS2005

New Directions in Space Operations Services in Support of Interplanetary Exploration

To gain access to the necessary operational processes and data in support of NASA's Lunar/Mars Exploration Initiative, new services, adequate levels of computing cycles and access to myriad forms of data must be provided to onboard spacecraft and ground based personnel/systems (earth, lunar and Martian) to enable interplanetary exploration by humans. These systems, cycles and access to vast amounts of development, test and operational data will be required to provide a new level of services not currently available to existing spacecraft, on board crews and other operational personnel. Although current voice, video and data systems in support of current space based operations has been adequate, new highly reliable and autonomous processes and services will be necessary for future space exploration activities. These services will range from the more mundane voice in LEO to voice in interplanetary travel which because of the high latencies will require new voice processes and standards. New services, like component failure predictions based on data mining of significant quantities of data, located at disparate locations, will be required. 3D or holographic representation of onboard components, systems or family members will greatly improve maintenance, operations and service restoration not to mention crew morale. Current operational systems and standards, like the Internet Protocol, will not able to provide the level of service required end to end from an end point on the Martian surface like a scientific instrument to a researcher at a university. Ground operations whether earth, lunar or Martian and in flight operations to the moon and especially to Mars will require significant autonomy that will require access to highly reliable processing capabilities, data storage based on network storage technologies. Significant processing cycles will be needed onboard but could be borrowed from other locations either ground based or onboard other spacecraft. Reliability will be a key factor with onboard and distributed backup processing an absolutely necessary requirement. Current cluster processing/Grid technologies may provide the basis for providing these services. An overview of existing services, future services that will be required and the technologies and standards required to be developed will be presented. The purpose of this paper will be to initiate a technological roadmap, albeit at a high level, of current voice, video, data and network technologies and standards (which show promise for adaptation or evolution) to what technologies and standards need to be redefined, adjusted or areas where new ones require development. The roadmap should begin the differentiation between non manned and manned processes/services where applicable. The paper will be based in part on the activities of the CCSDS Monitor and Control working group which is beginning the process of standardization of the these processes. Another element of the paper will be based on an analysis of current technologies supporting space flight processes and services at JSC, MSFC, GSFC and to a lesser extent at KSC. Work being accomplished in areas such as Grid computing, data mining and network storage at ARC, IBM and the University of Alabama at Huntsville will be researched and analyzed.

Bradford, Robert N.↗

DOE OSTI2021

Above-ground net primary productivity in regenerating seasonally dry tropical forest: Contributions of rainfall, forest age and soil

Identifying factors controlling forest productivity is critical to understanding forest-climate change feedbacks, modelling vegetation dynamics and carbon finance schemes. However, little research has focused on productivity in regenerating tropical forests which are expanding in their fraction of global area have an order of magnitude larger carbon uptake rates relative to older forest. Here We examined above-ground net primary productivity (ANPP) and its components (wood production and litterfall) over 10 years in forest plots that vary in successional age, soil characteristics and species composition using band dendrometers and litterfall traps in regenerating seasonally dry tropical forests in northwestern Costa Rica. We show that the components of ANPP are differentially driven by age and annual rainfall and that local soil variation is important. Total ANPP was explained by a combination of age, annual rainfall and soil variation. Wood production comprised 35% of ANPP on average across sites and years, and was explained by annual rainfall but not forest age. Conversely, litterfall increased with forest age and soil fertility yet was not affected by annual rainfall. In this region, edaphic variability is highly correlated with plant community composition. Thus, variation in ecosystem processes explained by soil may also be partially explained by species composition. These results suggest that future changes in annual rainfall can alter the secondary forest carbon sink, but this effect will be buffered by the litterfall flux which varies little among years. In determining the long-term strength of the secondary forest carbon sink, both rainfall and forest age will be critical variables to track. We also conclude that detailed understanding of local site variation in soils and plant community may be required to accurately predict the impact of changing rainfall on forest carbon uptake. Synthesis. We show that in seasonally dry tropical forest, annual rainfall has a positive relationship with the growth of above-ground woody tissues of trees and that droughts lead to significant reductions in above-ground productivity. These results provide evidence for climate change—carbon cycle feedbacks in the seasonal tropics and highlight the value of longitudinal data on forest regeneration.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2011

Assessment of an Automated Touchdown Detection Algorithm for the Orion Crew Module

Orion Crew Module (CM) touchdown detection is critical to activating the post-landing sequence that safe?s the Reaction Control Jets (RCS), ensures that the vehicle remains upright, and establishes communication with recovery forces. In order to accommodate safe landing of an unmanned vehicle or incapacitated crew, an onboard automated detection system is required. An Orion-specific touchdown detection algorithm was developed and evaluated to differentiate landing events from in-flight events. The proposed method will be used to initiate post-landing cutting of the parachute riser lines, to prevent CM rollover, and to terminate RCS jet firing prior to submersion. The RCS jets continue to fire until touchdown to maintain proper CM orientation with respect to the flight path and to limit impact loads, but have potentially hazardous consequences if submerged while firing. The time available after impact to cut risers and initiate the CM Up-righting System (CMUS) is measured in minutes, whereas the time from touchdown to RCS jet submersion is a function of descent velocity, sea state conditions, and is often less than one second. Evaluation of the detection algorithms was performed for in-flight events (e.g. descent under chutes) using hi-fidelity rigid body analyses in the Decelerator Systems Simulation (DSS), whereas water impacts were simulated using a rigid finite element model of the Orion CM in LS-DYNA. Two touchdown detection algorithms were evaluated with various thresholds: Acceleration magnitude spike detection, and Accumulated velocity changed (over a given time window) spike detection. Data for both detection methods is acquired from an onboard Inertial Measurement Unit (IMU) sensor. The detection algorithms were tested with analytically generated in-flight and landing IMU data simulations. The acceleration spike detection proved to be faster while maintaining desired safety margin. Time to RCS jet submersion was predicted analytically across a series of simulated Orion landing conditions. This paper details the touchdown detection method chosen and the analysis used to support the decision.

Gay, Robert S.↗

NASA NTRS2012

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

DOE OSTI2023

Carbon Organisms Rhizosphere and Protection in Soil Environment model script and input data for soil moisture-respiration responses in tropical forests

Objectives: Climatic drying is predicted for many tropical forests, yet models remain poorly parameterized for tropical forests, hampering predictions of forest-climate feedbacks. We applied an integrated model–experiment approach, parameterizing an ecosystem model Carbon Organisms Rhizosphere and Protection in the Soil Environment (CORPSE) with tropical forest observational data, and comparing model predictions with a field drying manipulation. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We used the field data to parameterize and run tests in the model.Results: Measured CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. We used this data to parameterize the model, which then predicted increased soil CO2 fluxes in wetter and fertile forests with drying, and decreased fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment in the forests initially suppressed soil CO2 fluxes across forests, with sustained suppression after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season), as predicted by the model. The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Code files:CORPSE_array.py: Defines the equations of the CORPSE modelCORPSE_solvers: Functions for running the CORPSE model using either iterative or ordinary differential equation (ODE) solversrun_Panama_sims.py: Read in datasets and run the model simulations for this studyInput data:PanamaGradientEcosystemChem_BT_CPools_20152016CO2_DC_20190615.xlsx: Plot characteristics used in running model simulationsLiCor compiled surface flux only to 2020_03 DC_20200825.xlsx: Surface gas exchange fluxes used in model-data comparisonsPARCHED litterfall data for Ben Sulman LD 20200902.xlsx: Litterfall data used to drive model simulationsInitialization data:state_500y_20190823.csv: Initial state of model pools based on previous spinup runsOutput data:Outputs/prev_moisture_response.csv: Simulations of multiple sites using original model moisture response function.Outputs/updated_moisture_response.csv: Simulations of multiple sites using updated model moisture response function.Outputs/dry15_prev_moisture_response.csv: Simulations with soil moisture reduced by 15%, using original moisture response function.Outputs/dry15_updated_moisture_response.csv: Simulations with soil moisture reduced by 15%, using updated moisture response function.Outputs/dry30_prev_moisture_response.csv: Simulations with soil moisture reduced by 30%, using original moisture response function.Outputs/dry30_updated_moisture_response.csv: Simulations with soil moisture reduced by 30%, using updated moisture response function.Outputs/latestart_prev_moisture_response.csv: Simulations with extended dry season, using original moisture response function.Outputs/latestart_updated_moisture_response.csv: Simulations with extended dry season, using updated moisture response function.Outputs/[site name]_oneyear.csv: One-year simulation for each site in expanded site list using original moisture response function.Outputs/[site name]_oneyear_dried.csv: One-year simulation for each site in expanded site list using original moisture response function, with soil moisture reduced by 25%.Outputs/[site name]_oneyear_updated_moisture_response.csv: One-year simulation for each site in expanded site list using updated moisture response function.Outputs/[site name]_oneyear_updated_moisture_response_dried.csv: One-year simulation for each site in expanded site list using updated moisture response function, with soil moisture reduced by 25%.Field plot location data:There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site).

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2002

Trends in lignin modification: a comprehensive analysis of the effects of genetic manipulations/mutations on lignification and vascular integrity

A comprehensive assessment of lignin configuration in transgenic and mutant plants is long overdue. This review thus undertook the systematic analysis of trends manifested through genetic and mutational manipulations of the various steps associated with monolignol biosynthesis; this included consideration of the downstream effects on organized lignin assembly in the various cell types, on vascular function/integrity, and on plant growth and development. As previously noted for dirigent protein (homologs), distinct and sophisticated monolignol forming metabolic networks were operative in various cell types, tissues and organs, and form the cell-specific guaiacyl (G) and guaiacyl-syringyl (G-S) enriched lignin biopolymers, respectively. Regardless of cell type undergoing lignification, carbon allocation to the different monolignol pools is apparently determined by a combination of phenylalanine availability and cinnamate-4-hydroxylase/"p-coumarate-3-hydroxylase" (C4H/C3H) activities, as revealed by transcriptional and metabolic profiling. Downregulation of either phenylalanine ammonia lyase or cinnamate-4-hydroxylase thus predictably results in reduced lignin levels and impaired vascular integrity, as well as affecting related (phenylpropanoid-dependent) metabolism. Depletion of C3H activity also results in reduced lignin deposition, albeit with the latter being derived only from hydroxyphenyl (H) units, due to both the guaiacyl (G) and syringyl (S) pathways being blocked. Apparently the cells affected are unable to compensate for reduced G/S levels by increasing the amounts of H-components. The downstream metabolic networks for G-lignin enriched formation in both angiosperms and gymnosperms utilize specific cinnamoyl CoA O-methyltransferase (CCOMT), 4-coumarate:CoA ligase (4CL), cinnamoyl CoA reductase (CCR) and cinnamyl alcohol dehydrogenase (CAD) isoforms: however, these steps neither affect carbon allocation nor H/G designations, this being determined by C4H/C3H activities. Such enzymes thus fulfill subsidiary processing roles, with all (except CCOMT) apparently being bifunctional for both H and G substrates. Their severe downregulation does, however, predictably result in impaired monolignol biosynthesis, reduced lignin deposition/vascular integrity, (upstream) metabolite build-up and/or shunt pathway metabolism. There was no evidence for an alternative acid/ester O-methyltransferase (AEOMT) being involved in lignin biosynthesis.The G/S lignin pathway networks are operative in specific cell types in angiosperms and employ two additional biosynthetic steps to afford the corresponding S components, i.e. through introduction of an hydroxyl group at C-5 and its subsequent O-methylation. [These enzymes were originally classified as ferulate-5-hydroxylase (F5H) and caffeate O-methyltransferase (COMT), respectively.] As before, neither step has apparently any role in carbon allocation to the pathway; hence their individual downregulation/manipulation, respectively, gives either a G enriched lignin or formation of the well-known S-deficient bm3 "lignin" mutant, with cell walls of impaired vascular integrity. In the latter case, COMT downregulation/mutation apparently results in utilization of the isoelectronic 5-hydroxyconiferyl alcohol species albeit in an unsuccessful attempt to form G-S lignin proper. However, there is apparently no effect on overall G content, thereby indicating that deposition of both G and S moieties in the G/S lignin forming cells are kept spatially, and presumably temporally, fully separate. Downregulation/mutation of further downstream steps in the G/S network [i.e. utilizing 4CL, CCR and CAD isoforms] gives predictable effects in terms of their subsidiary processing roles: while severe downregulation of 4CL gave phenotypes with impaired vascular integrity due to reduced monolignol supply, there was no evidence in support of increased growth and/or enhanced cellulose biosynthesis. CCR and CAD downregulation/mutations also established that a depletion in monolignol supply reduced both lignin contents supply reduced both lignin contents and vascular integrity, with a concomitant shift towards (upstream) metabolite build-up and/or shunting.The extraordinary claims of involvement of surrogate monomers (2-methoxybenzaldehyde, feruloyl tyramine, vanillic acid, etc.) in lignification were fully disproven and put to rest, with the investigators themselves having largely retracted former claims. Furthermore analysis of the well-known bm1 mutation, a presumed CAD disrupted system, apparently revealed that both G and S lignin components were reduced. This seems to imply that there is no monolignol specific dehydrogenase, such as the recently described sinapyl alcohol dehydrogenase (SAD) for sinapyl alcohol formation. Nevertheless, different CAD isoforms of differing homology seem to be operative in different lignifying cell types, thereby giving the G-enriched and G/S-enriched lignin biopolymers, respectively. For the G-lignin forming network, however, the CAD isoform is apparently catalytically less efficient with all three monolignols than that additionally associated with the corresponding G/S lignin forming network(s), which can more efficiently use all three monolignols. However, since CAD does not determine either H, G, or S designation, it again serves in a subsidiary role-albeit using different isoforms for different cell wall developmental and cell wall type responses.The results from this analysis contrasts further with speculations of some early investigators, who had viewed lignin assembly as resulting from non-specific oxidative coupling of monolignols and subsequent random polymerization. At that time, though, the study of the complex biological (biochemical) process of lignin assembly had begun without any of the (bio)chemical tools to either address or answer the questions posed as to how its formation might actually occur. Today, by contrast, there is growing recognition of both sophisticated and differential control of monolignol biosynthetic networks in different cell types, which serve to underscore the fact that complexity of assembly need not be confused any further with random formation. Moreover, this analysis revealed another factor which continues to cloud interpretations of lignin downregulation/mutational analyses, namely the serious technical problems associated with all aspects of lignin characterization, whether for lignin quantification, isolation of lignin-enriched preparations and/or in determining monomeric compositions. For example, in the latter analyses, some 50-90% of the lignin components still cannot be detected using current methodologies, e.g. by thioacidolysis cleavage and nitrobenzene oxidative cleavage. This deficiency in lignin characterization thus represents one of the major hurdles remaining in delineating how lignin assembly (in distinct cell types) and their configuration actually occurs.

Review, Academic↗

NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗

NASA NTRSDate not supplied

Thermo-Poro-Mechanical Modeling of RTV Intumescence

silicone intumescence↗

DOE OSTI2020

Prediction for Pressure Differential Across HEPA Filter Media Based on Media Characteristics and Particle Size Distribution

A new method for predicting the pressure drop across High Efficiency Particulate Air (HEPA) filter media is proposed based upon mass deposited onto the filter and known physical characteristics of the filter media. Detailed are the methods used in conjunction with current filter loading models to predict the pressure drop, as well as tests conducted to validate the methods. The benefit of a prediction model for practical use lies in the manufacturing and service life of nuclear grade HEPA filters. HEPA filters for use in nuclear facilities have a prescribed expiration date and a maximum allowable operating pressure drop. Therefore, the ability to predict the pressure drop across a filter and relate it to an expected length of service time can enable a reduction in wasted filters. This will allow for more informed decisions to be made based upon the dictated life cycle of the filters. Also, understanding and predicting how the pressure drop of HEPA filter media behaves as a function of physical characteristics and loaded mass can assist in future design and manufacturing of filter media. Currently, most existing pressure drop models are either computationally based or analytical methods relying on data gathered during tests. Neither are practical for prediction; the computational methods are difficult to implement, and the current analytical methods are more useful as tools for analysis. The current analytical model developed by Bergman et al. is based upon the pressure drop from each media fiber and modeling the deposited particles as newly formed fibers. Using Bergman's model and the media properties, air properties, and known aerosol particle size distribution, a pressure drop curve as a function of loaded mass can be created. This curve implies initial loading in the depth of the filter media with a transition to pure surface loading on the filter face, and a sensitivity to an evolving particle size distribution as mass continues to load onto the filter. To implement the predictive model, knowledge of the mean filter fiber diameter and porosity of the filter media is required. Traditionally, the mean fiber diameter is calculated as an effective diameter from prior media testing, however, in this study a Scanning Electron Microscope (SEM) technique is used to acquire this variable. Validation of the predictive model is provided by flat sheet media tests under a controlled environment with a measured particle size distribution of the challenge aerosol and shows reasonable preliminary agreement.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗