Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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 847 records · Page 47

Experimental Verification of a Progressive Damage Model for IM7/5260 Laminates Subjected to Tension-Tension Fatigue

The durability and damage tolerance of laminated composites are critical design considerations for airframe composite structures. Therefore, the ability to model damage initiation and growth and predict the life of laminated composites is necessary to achieve structurally efficient and economical designs. The purpose of this research is to experimentally verify the application of a continuum damage model to predict progressive damage development in a toughened material system. Damage due to monotonic and tension-tension fatigue was documented for IM7/5260 graphite/bismaleimide laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables to predict stiffness loss in unnotched laminates. A damage dependent finite element code predicted the stiffness loss for notched laminates with good agreement to experimental data. It was concluded that the continuum damage model can adequately predict matrix damage progression in notched and unnotched laminates as a function of loading history and laminate stacking sequence.

Coats, Timothy W.↗

Parallelization and Visual Analysis of Multidimensional Fields: Application to Ozone Production, Destruction, and Transport in Three Dimensions

The three-dimensional, spectral transport model used in the current project was first successfully integrated over climatological time scales by Dr. Guang Ping Lou for the simulation of atmospheric N2O using the United Kingdom Meteorological Office (UKMO) 4-dimensional, assimilated wind and temperature data set. A non-parallel, FORTRAN version of this integration using a fairly simple N2O chemistry package containing only photo-chemical reactions was used to verify our initial parallel model results. The integrations reproduced the gross features of the observed stratospheric climatological N2O distributions but also simulated the structure of the stratospheric Antarctic vortex and its evolution. Subsequently, Dr. Thomas Kindler, who produced much of the parallel version of our model, enlarged the N2O model chemistry package to include N2O reactions involving O(D-1) and also introduced assimilated wind data from NASA as well as UKMO. Initially, transport calculations without chemistry were run using Carbon-14 as a non-reactive tracer gas with the result that large differences in the transport properties of the two assimilated wind data sets were apparent from the resultant Carbon-14 distributions. Subsequent calculations for N2O, including its chemistry, with the two input winds data sets with verification from UARS satellite observations have refined the transport differences between the two such that the model's steering capabilities could be used to infer the correct climatological vertical velocity fields required to support the N2O observations. During this process, it was also discovered that both the NASA and the UKMO data contained spurious values in some of the higher frequency wave components, leading to incorrect local transport calculations and ultimately affecting the large scale properties of the model's N2O distributions, particularly at tropical latitudes. Subsequent model runs with wind data that had been filtered to remove some of the high frequency components produced much more realistic N2O distributions. During the past few months, the UKMO wind data base for a complete two-year period was processed into spectral form for model use. This new version of the input transport data base now includes complete temperature fields as well as the necessary wind data. This was done to facilitate advanced chemical calculations in the parallel model which often depend upon temperature. Additional UKMO data is being added as it becomes available.

Schwan, Karsten↗

Use of an adjoint model for finding triggers for Alpine lee cyclogenesis

The authors propose a new procedure, designated the adjoint-based genesis diagnostic (AGD) procedure, for studying triggering mechanisms and the subsequent genesis of the synoptic phenomena of interest. This procedure makes use of a numerical model sensitivity to initial conditions and the nonlinear evolution of the initial perturbations that are designed using this sensitivity. The model sensitivity is evaluated using the associated adjoint model. This study uses the dry version of the National Center for Atmospheric Research Mesoscale Adjoint Modeling System (MAMS) for the numerical experiments. The authors apply the AGD procedure to two cases of Alpine lee cyclogenesis that were observed during the Alpine Experiment special observations period. The results show that the sensitivity fields that are produced by the adjoint model and the associated initial perturbations are readily related to the probable triggering mechanisms for these cyclones. Additionally, the nonlinear evolution of these initial perturbations points toward the physical processes involved in the lee cyclone formation. The AGD experiments for a weak cyclone case indicate that the MAMS forecast model has an underrepresented topographic forcing due to the sigma vertical coordinate and that this model error can be compensated by adjustments in the initial conditions that are related to the triggering mechanisms, which is not associated with the topographic blocking mechanism.

Vukicevic, Tomislava↗

Intermountain West Energy Sustainability & Transitions Initiative: CO2 Transport and Geologic Storage Modeling Results

These resources provide the full set of cost modeling results and methods as part of the I-WEST Roadmap Initiative that were used to compile figures as part of the "Pathways to CO2 Utilization and Storage for the Intermountain West Region" chapter. The data were generated from a series of National Energy Technology Laboratory (NETL) cost models and relate to carbon dioxide (CO2) transport costs, CO2 enhanced oil recovery (CO2-EOR) economics, and saline storage economics. These models were used to analyze various business cases given changes in technical and financial assumptions for the I-WEST region as a means to explore how these assumptions influence CO2 transport and storage costs, as well as to evaluate the effect of changing oil prices on the viability of CO2-EOR and the mass of CO2 stored via CO2-EOR. The accompanying report titled "Intermountain West Energy Sustainability & Transitions Initiative: CO2 Transport and Geologic Storage Modeling Results" provides a detailed overview on the models, assumptions, and parameters used in the modeling, as well as example results..

CO2 EOR costs,CO2 Storage Costs,CO2 transportation↗

Representation and evolution of urban weather boundary conditions in downtown Chicago

Our study presents a novel computing technique for data exchange and coupling between a high-resolution weather simulation model and a building energy model, with a goal of evaluating the impact of urban weather boundary conditions on energy performance of urban buildings. The Weather Research and Forecasting (WRF) model is initialized with the operational High-Resolution Rapid Refresh (HRRR) dataset to provide hourly weather conditions over the Chicago region. We utilize the building footprint, land use, and building stock datasets to generate building energy models using EnergyPlus. We mapped the building exterior surfaces to local air nodes to import simulated microclimate data and to export buildings' heat emissions to their local environment. Preliminary experiments for a test area in Chicago show that predicted building cooling energy use differs by about 4.7% for the selected date when compared with simulations using TMY weather data and without considering the urban microclimate boundary conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DCMIP2016: the splitting supercell test case

Abstract. This paper describes the splitting supercell idealized test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). These storms are useful test beds for global atmospheric models because the horizontal scale of convective plumes is O(1 km), emphasizing non-hydrostatic dynamics. The test case simulates a supercell on a reduced-radius sphere with nominal resolutions ranging from 4 to 0.5 km and is based on the work of Klemp et al. (2015). Models are initialized with an atmospheric environment conducive to supercell formation and forced with a small thermal perturbation. A simplified Kessler microphysics scheme is coupled to the dynamical core to represent moist processes. Reference solutions for DCMIP2016 models are presented. Storm evolution is broadly similar between models, although differences in the final solution exist. These differences are hypothesized to result from different numerical discretizations, physics–dynamics coupling, and numerical diffusion. Intramodel solutions generally converge as models approach 0.5 km resolution, although exploratory simulations at 0.25 km imply some dynamical cores require more refinement to fully converge. These results can be used as a reference for future dynamical core evaluation, particularly with the development of non-hydrostatic global models intended to be used in convective-permitting regimes.

Zarzycki, Colin M.↗

COUPLED GRIFFIN AND PRONGHORN SIMULATION OF THE MOLTEN SALT FAST REACTOR (MSFR) FOR THE VIRTUAL TEST BED

The Virtual Test Bed (VTB) repository hosts a wide range of challenge problems to showcase modeling and simulation capabilities to support advanced reactor demonstration. An overview of a coupled multiphysics model for Molten Salt Reactors (MSR) is presented here. The analysis leverages MOOSE-based tools (notably Griffin and Pronghorn) for neutronic and thermal hydraulic simulation. Neutron precursor drifting capability is showcased, along with some initial turbulence models. Both steady-state and transient coupled multiphysics results for an MSR concept are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Observation Denial and Performance of a Local Mesoscale Model

Forecasters at the 45th Weather Squadron (45 WS) use observations from the Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) wind tower network and the CCAFS (XMR) daily rawinsonde observations (RAOB) to issue and verify wind advisories and warnings for operations. These observations are also used by the Spaceflight Meteorology Group (SMG) in Houston, Texas and the Melbourne, Florida National Weather Service office to initialize their locally run mesoscale models. SMG also uses the observations to support shuttle landings at the KSC Shuttie Landing Facility. Due to impending budget cuts, some or all of the KSC/CCAFS wind towers on the east-central Florida mainland and the XMR RAOBs may be eliminated. The loss of these data may impact the forecast capability of the 45 WS and SMG. The Applied Meteorology Unit (AMU) was tasked to conduct a modeling study to determine how important these observations are to the accuracy of the model output used by the forecasters as input to their forecasts. To accomplish this, the AMU performed a sensitivity study using the Weather Research and Forecasting (NRF) model initialized with and without KSC/CCAFS wind tower and XMR RAOB observations. The AMU assessed the accuracy of model output by comparing peak wind forecasts with operationally significant wind advisory and warning criteria forecast by the 45 WS. To assess model performance when initialized with and without some of the wind tower and XMR RAOB observations, the AMU conducted a subjective analysis by displaying model wind forecasts graphically with the observations overlaid for comparison and they conducted an objective analysis by comparing the maximum peak wind forecast to the maximum peak wind observed within the KSC/CCAFS wind tower network. Data were collected for twelve warm season cases and eight cool season cases from June - September 2007 and November - January 2008, respectively. For each case chosen, the 45 WS must have issued a wind advisory or warning for the KSC/CCAFS area and the KSC/CCAFS wind towers must have recorded significant wind events, or winds greater than 18 kt.

Watson, Leela R.↗

Application of the metabolic modeling pipeline in KBase to categorize reactions, predict essential genes, and predict pathways in an isolate genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

DOE knowledgebase↗

Application of the Metabolic Modeling Pipeline in KBase to Categorize Reactions, Predict Essential Genes, and Predict Pathways in an Isolate Genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

Allen, Benjamin↗

Atomic protein structure refinement using all-atom graph representations and SE(3)-equivariant graph transformer

Abstract Motivation The state-of-art protein structure prediction methods such as AlphaFold are being widely used to predict structures of uncharacterized proteins in biomedical research. There is a significant need to further improve the quality and nativeness of the predicted structures to enhance their usability. In this work, we develop ATOMRefine, a deep learning-based, end-to-end, all-atom protein structural model refinement method. It uses a SE(3)-equivariant graph transformer network to directly refine protein atomic coordinates in a predicted tertiary structure represented as a molecular graph. Results The method is first trained and tested on the structural models in AlphaFoldDB whose experimental structures are known, and then blindly tested on 69 CASP14 regular targets and 7 CASP14 refinement targets. ATOMRefine improves the quality of both backbone atoms and all-atom conformation of the initial structural models generated by AlphaFold. It also performs better than two state-of-the-art refinement methods in multiple evaluation metrics including an all-atom model quality score—the MolProbity score based on the analysis of all-atom contacts, bond length, atom clashes, torsion angles, and side-chain rotamers. As ATOMRefine can refine a protein structure quickly, it provides a viable, fast solution for improving protein geometry and fixing structural errors of predicted structures through direct coordinate refinement. Availability and implementation The source code of ATOMRefine is available in the GitHub repository (https://github.com/BioinfoMachineLearning/ATOMRefine). All the required data for training and testing are available at https://doi.org/10.5281/zenodo.6944368.

59 BASIC BIOLOGICAL SCIENCES↗

The future of the El Niño–Southern Oscillation: using large ensembles to illuminate time-varying responses and inter-model differences

Future changes in the El Niño–Southern Oscillation (ENSO) are uncertain, both because future projections differ between climate models and because the large internal variability of ENSO clouds the diagnosis of forced changes in observations and individual climate model simulations. By leveraging 14 single model initial-condition large ensembles (SMILEs), we robustly isolate the time-evolving response of ENSO sea surface temperature (SST) variability to anthropogenic forcing from internal variability in each SMILE. We find nonlinear changes in time in many models and considerable inter-model differences in projected changes in ENSO and the mean-state tropical Pacific zonal SST gradient. We demonstrate a linear relationship between the change in ENSO SST variability and the tropical Pacific zonal SST gradient, although forced changes in the tropical Pacific SST gradient often occur later in the 21st century than changes in ENSO SST variability, which can lead to departures from the linear relationship. Single-forcing SMILEs show a potential contribution of anthropogenic forcing (aerosols and greenhouse gases) to historical changes in ENSO SST variability, while the observed historical strengthening of the tropical Pacific SST gradient sits on the edge of the model spread for those models for which single-forcing SMILEs are available. Our results highlight the value of SMILEs for investigating time-dependent forced responses and inter-model differences in ENSO projections. The nonlinear changes in ENSO SST variability found in many models demonstrate the importance of characterizing this time-dependent behavior, as it implies that ENSO impacts may vary dramatically throughout the 21st century.

54 ENVIRONMENTAL SCIENCES↗

The dynamics of the Venus ionosphere. I - A simulation of the solar wind compression of the upper dayside ionosphere

One of the most exciting discoveries by the Pioneer Venus mission is the extreme variability in the structure of the Venus atmosphere. The solar wind plays a major, although as yet not well-defined, role in the dynamics of the ionosphere. An investigation is being conducted regarding the response of the dayside Venus ionosphere to changing solar wind conditions, in particular to the varying solar wind dynamic pressure. In the present study the dynamics of the upper (h equal to or greater than 200 km) ionosphere are simulated numerically using a one-dimensional, spherically symmetric, Lagrangian hydrodynamic code developed by Stein and Schwartz (1972). The ionosphere is assumed to be unmagnetized and is represented with a two-fluid model. The initial ionosphere is chosen to be in pressure equilibrium with the solar wind at the ionopause. It is shown how some of the time-dependent features of the Venus ionosphere may be simulated with the considered model.

Wolff, R. S.↗

Dynamics of the transition zone in coastal zone color scanner-sensed ocean color in the North Pacific during oceanographic spring

A transition zone in phytoplankton concentration running across the North Pacific basin at 30 deg to 40 deg north latitude corresponds to a basin-wide front in surface chlorophyll observed in a composite of coastal zone color scanner (CZCS) images for May, June, and July 1979-1986. This transition zone with low chlorophyll to the south and higher chlorophyll to the north can be simulated by a simple model of the concentration of phytoplankton, zooplankton, and dissolved nutrient (nitrate) in the surface mixed layer of the ocean applied to the North Pacific basin for the climatological conditions during oceanographic springtime (May, June, and July). The model is initialized with a 1 deg x 1 deg gridded estimate of wintertime (February, March, and April) mixed layer nitrate concentrations calculated from an extensive nutrient database and a similarly gridded mixed layer depth data set. Comparison of model predictions with CZCS data provides a means of evaluating the dynamics of the transition zone. We conclude that in the North Pacific, away from major boundary currents and coastal upwelling zones, wintertime vertical mixing determines the total nutrient available to the plankton ecosystem in the spring. The transition zone seen in basin-scale CZCS images is a reflection of the geographic variation in the wintertime mixed layer depth and the nitracline, leading to a latitudinal gradient in phytoplankton chlorophyll.

Glover, David M.↗

Answering the Call for Model-Relevant Observations of Aerosols and Clouds

We describe a technique for combining multiple A-Train aerosol data sets, namely MODIS spectral AOD (aerosol optical depth), OMI AAOD (absorption aerosol optical depth) and CALIOP aerosol backscatter retrievals (hereafter referred to as MOC retrievals) to estimate full spectral sets of aerosol radiative properties, and ultimately to calculate the 3-D distribution of direct aerosol radiative effects (DARE). We present MOC results using almost two years of data collected in 2007 and 2008, and show comparisons of the aerosol radiative property estimates to collocated AERONET retrievals. We compare the spatio-temporal distribution of the MOC retrievals and MOC-based calculations of seasonal clear-sky DARE to values derived from four models that participated in the Phase II AeroCom model intercomparison initiative. Comparisons of seasonal aerosol property to AeroCom Phase II results show generally good agreement best agreement with forcing results at TOA is found with GMI-MerraV3.We discuss the challenges in making observations that really address deficiencies in models, with some of the more relevant aspects being representativeness of the observations for climatological states, and whether a given model-measurement difference addresses a sampling or a model error.

Redemann, J.↗

Time-Varying Modification of Reference Model for Adaptive Control with Performance Optimization

This paper presents a new adaptive control ap- proach that involves a performance optimization objective. The control synthesis involves the design of a performance optimizing adaptive controller from a subset of control inputs. The resulting effect of the performance optimizing adaptive controller is to modify the initial reference model into a time-varying reference model which satisfies the performance optimization requirement obtained from an optimal control problem. The time-varying reference model modification is accomplished by the real-time solutions of the time-varying Riccati and Sylvester equations coupled with the least-squares parameter estimation of the sen- sitivities of the performance metric. The effectiveness of the proposed method is demonstrated by an application of maneuver load alleviation control for a flexible aircraft.

Adaptive Control↗

Experimental Combustion and Flame Characterization of a Chemical Looping-Based Oxidative Dehydrogenation Byproduct Fuel Mixture Containing High CO2 Dilution

Abstract This study investigates the combustion performance of a CO2-rich fuel mixture containing ethane and methane as active species using a constant volume combustion chamber. This fuel is obtained as a byproduct of a chemical looping-based oxidative dehydrogenation (Cl-ODH) process ethylene production. The byproduct gas mixture has 40.79% CO2, 39.49% ethane, and 4.88% methane by weight with other minor compounds. Using this fuel for energy extraction would improve the process efficiency of the ethane to ethylene conversion. After initial combustion modeling, the gas fuel mixture was reduced to just the major species: CO2, ethane, and methane. The mixture was then tested for flammability limits and combustion performance under spark-ignition conditions. Effects of ambient conditions like temperatures between 300 and 400 K with initial pressures from 1 to 10 bar were tested. The effects of stoichiometry were tested to understand flame velocities and heat release. The fuel mixture showed an overall reduced flame velocity compared to gasoline. Instability in combustion was believed to be caused by the dissociation of ethane under elevated conditions. At higher pressures, the flame produces lower cumulative heat release. Simulations were also performed using a model tuned to replicate the operations of the combustion chamber used in the experiments. Heat release and unburnt fuel mass data were calculated to identify the discrepancies in the combustion completeness at elevated pressures. The effects of CO2 quenching the flame coupled with the increased dissociation of the fuel species can lead to up to more than 75% of the fuel mixture being unburnt. Data from this study were used to modify a small-scale spark-ignition engine to use this fuel and produce usable energy.

Energy & Fuels↗

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

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

42 ENGINEERING↗