Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modeling simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,603 records · Page 89

Dynamic modelling and simulation study for the Galileo spacecraft pulsed-mode spinup/400 N main engine burn/spindown maneuvers

Two Galileo dynamic models were developed to simulate the spinup, the main engine burn and the spindown maneuvers for the critical events of the Jupiter orbit insertion and perijove raise maneuver. The dynamic interaction between the spin thruster pulsing frequency, the science magnetometer boom flexible modes and the propellant slosh modes were studied. The prediction of safe orbit insertion was validated in flight. For the perijove raise maneuver, the simulation results indicated that propellant unporting would occur if the original spin thruster duty cycle was not updated. Based on further simulation results, a new duty cycle was selected to prevent propellant unporting. Subsequent Galileo flight data indicated that the perijove raise maneuver was successfully executed.

Ih, Che-Hang Charles↗

Dynamic Modeling and Simulation Study for the Galileo Spacecraft Pulsed-Mode Spinup/400 N Main Engine Burn/Spindown Maneuvers

Two Galileo dynamic models were developed to simulate the spinup/400-N main engine burn/spindown maneuvers for the critical events of Jupiter Orbit Insertion (JOI) and Perijove Raise Maneuver (PRM). The dynamic interaction among the spin thruster pulsing frequency science/magnetometer (SCI/MAG) boom flexible modes, and the propellant slosh modes were studied.

Galileo orbit insertion planetary spacecraft↗

Intrusive analysis for NEK5000: development of intrusive uncertainty quantification for high-dimensional, high-fidelity codes.

We study the use of lower-fidelity training data for uncertainty quantification of complex simulation models. In our approach, computationally expensive full-model outputs are approximated applying proper orthogonal-decomposition-based dimensionality reduction to the full model. A Gaussian-processes-based machine learning approach is then used to model the difference (or, rather, the correspondence) between the higher-fidelity and the lower-fidelity data. This stochastic model can be constructed by using few additional full code evaluations; it is used to calibrate arbitrarily many lower-quality outputs in order to create additional training data for sampling-based analysis. In effect, an adequate approximation of the simulation model’s response to uncertainty can be constructed at a modest computational cost. In fact, we aim at the number of full-model evaluations comparable to that required for a linear approximation (as opposed to numbers traditionally associated with sampling in high-dimensional spaces). In this report, we explain the basic algorithm, suggest some performance tuning options, and give a first characterization for the class of simulation models of nuclear engineering for which the suggested report is effective.

Roderick, O.↗

Next Generation Co-Molded One-Piece Automotive Parts

The purpose of this project was to determine the feasibility of co-molding Class A Sheet Molding Compound (SMC), structural SMC, and continuous fiber prepreg materials to produce a single piece co-molded automotive part, such as a hood. The combination of the three molding materials and the incorporation of selective design features (ribs, flanges, corrugations) was expected to eliminate the need for inner reinforcement panels, significantly reducing tooling costs and simplifying the manufacturing process. A multi-material solution would also result in significant weight savings. The scope of this project included the material characterization of the three materials, development of cure and flow simulation models, and validation of the models against parts molded with different combinations of the materials on an 11”x11” plaque tool with rib features. The intent of this project was to apply the learnings obtained from co-molding a part with simplified geometry to a Phase 2 project that would produce a single piece hood, co-molded with the same materials. Resins from INEOS Composites were provided to IDI to produce SMC and continuous fiber prepregs. Purdue University and INEOS Composites characterized the rheological and curing behavior of the different materials as well as the mechanical properties of the co-molded parts. This data was used by Purdue University to create simulation models. Models were created to predict both flow patterns and predict mechanical properties of different laminate constructions in and around the ribs. Michigan State University - Corktown validated Purdue’s models by co-molding the three materials in different combinations on a tool containing rib features provided by Century Tool. The co-molded parts were evaluated against the simulation models by Purdue, Corktown and INEOS. The project team demonstrated the following: • The ability to obtain a cohesive co-molded structure made up of Class A SMC, Structural SMC, and continuous fiber prepreg. • The ability to obtain a co-molded part with a Class A surface • Modeling of flow and fiber orientation • Modeling to predict mechanical properties of multi-material co-molded structures The predicted flow behavior and fiber orientation from simulations were then compared with the molded samples. Exact local orientation state was difficult to compare between microscopy and flow simulation, but captured general trends. Multi-material flow behavior was generally modeled well with SMC materials, but the introduction of woven material sheets requires further model development. Purdue compared the predicted versus actual mechanical properties, and INEOS Composites evaluated the surface appearance of unpainted and painted co-molded parts. The models developed by Purdue demonstrated that the mechanical properties can be predicted for parts with varying material configurations. The resulting models can be applied to future co-molded part design, tooling design, and molding conditions.

36 MATERIALS SCIENCE↗

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗

Evaluation of a High-Resolution Numerical Weather Prediction Model’s Simulated Clouds Using Observations from CloudSat, GOES-13 and in Situ Aircraft

This study aimed to assess tropical cloud properties predicted by Environment and Climate Change Canada’s Global Environmental Multiscale (GEM)modelwhen run with the Milbrandt–Yau double-moment cloud microphysical scheme and one-way nesting that culminated at a (∼300 km)2 inner domain with 0.25 km horizontal grid spacing. The assessment utilized satellite and in situ data collected during the High Ice Water Content (HIWC) and High Altitude Ice Crystals (HAIC) projects for a mesoscale convective system on 16May 2015 over French Guiana.Data fromCloud- Sat’s cloud-profiling radar and GOES-13’s imager were compared to data either simulated directly by GEM or produced by operating on GEM’s cloud data with both the CFMIP (Cloud Feedback Model Intercomparison Project) Observation Simulator Package (COSP) instrument simulator and a three-dimensional Monte Carlo solar radiative transfer model. In situ observations were made from research aircraft – Canada’s National Research Council Convair-580 and the French SAFIRE Falcon-20 – whose flight paths were aligned with CloudSat’s ground-track. Spatial and temporal shifts of clouds simulated by GEM compared well to GOES-13 imagery. There are, however, differences between simulated and observed amounts of high and low cloud. While GEM did well at predicting ranges of ice-water content (IWC) near 11 km altitude (Falcon-20), it produces too much graupel and snow near 7 km (Convair-580). This produced large differences between CloudSat’s and COSP-generated radar reflectivities and two-way attenuations. On the other hand, CloudSat’s inferred values of IWC agree well with in situ samples at both altitudes. Generally, GEM’s visible reflectances exceeded GOES-13’s on account of having produced too much low-level liquid cloud. It is expected that GEM’s disproportioning of cloud hydrometeors will improve once it includes a better representation of secondary ice production.

Qu, Zhipeng↗

Propulsion System Modeling and Simulation

The Aerospace Systems Design Laboratory at the School of Aerospace Engineering in Georgia Institute of Technology has developed a core competency that enables propulsion technology managers to make technology investment decisions substantiated by propulsion and airframe technology system studies. This method assists the designer/manager in selecting appropriate technology concepts while accounting for the presence of risk and uncertainty as well as interactions between disciplines. This capability is incorporated into a single design simulation system that is described in this paper. This propulsion system design environment is created with a commercially available software called iSIGHT, which is a generic computational framework, and with analysis programs for engine cycle, engine flowpath, mission, and economic analyses. iSIGHT is used to integrate these analysis tools within a single computer platform and facilitate information transfer amongst the various codes. The resulting modeling and simulation (M&S) environment in conjunction with the response surface method provides the designer/decision-maker an analytical means to examine the entire design space from either a subsystem and/or system perspective. The results of this paper will enable managers to analytically play what-if games to gain insight in to the benefits (and/or degradation) of changing engine cycle design parameters. Furthermore, the propulsion design space will be explored probabilistically to show the feasibility and viability of the propulsion system integrated with a vehicle.

Tai, Jimmy C. M.↗

Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery Process

This paper introduces a first-order physics-based model that accounts for the fundamental heat and mass transfer between a humid-air vapor stream on feed side to another flow stream on permeate side. The model comprises a few optional submodels for membrane mass transport; and it adopts a segment-by-segment method for discretizing heat and mass transfer governing equations for flow streams on feed and permeate sides. The model is able to simulate both dehumidifiers and energy recovery ventilators in parallel-flow, cross-flow, and counter-flow configurations. The predicted tresults are compared reasonably well with the measurements. The open-source codes are written in C++. The model and open-source codes are expected to become a fundament tool for the analysis of membrane-based dehumidification in the future.

Gao, Zhiming↗

Simulations of Dynamics and Transport during the September 2002 Antarctic Major Warming

A mechanistic model simulation initialized on 14 September 2002, forced by 100-hPa geopotential heights from Met Office analyses, reproduced the dynamical features of the 2002 Antarctic major warming. The vortex split on approx.25 September; recovery after the warming, westward and equatorward tilting vortices, and strong baroclinic zones in temperature associated with a dipole pattern of upward and downward vertical velocities were all captured in the simulation. Model results and analyses show a pattern of strong upward wave propagation throughout the warming, with zonal wind deceleration throughout the stratosphere at high latitudes before the vortex split, continuing in the middle and upper stratosphere and spreading to lower latitudes after the split. Three-dimensional Eliassen-Palm fluxes show the largest upward and poleward wave propagation in the 0(deg)-90(deg)E sector prior to the vortex split (coincident with the location of strongest cyclogenesis at the model's lower boundary), with an additional region of strong upward propagation developing near 180(deg)-270(deg)E. These characteristics are similar to those of Arctic wave-2 major warmings, except that during this warming, the vortex did not split below approx.600 K. The effects of poleward transport and mixing dominate modeled trace gas evolution through most of the mid- to high-latitude stratosphere, with a core region in the lower-stratospheric vortex where enhanced descent dominates and the vortex remains isolated. Strongly tilted vortices led to low-latitude air overlying vortex air, resulting in highly unusual trace gas profiles. Simulations driven with several meteorological datasets reproduced the major warming, but in others, stronger latitudinal gradients at high latitudes at the model boundary resulted in simulations without a complete vortex split in the midstratosphere. Numerous tests indicate very high sensitivity to the boundary fields, especially the wave-2 amplitude. Major warmings occurred for initial fields with stronger winds and larger vortices, but not smaller vortices, consistent with the initiation of wind-deceleration by upward-propagating waves near the poleward edge of the region where wave 2 can propagate above the jet core. Thus, given the observed 100-hPa boundary forcing, stratospheric preconditioning is not needed to reproduce a major warming similar to that observed. The anomalously strong forcing in the lower stratosphere can be viewed as the primary direct cause of the major warming.

winds↗

INL Report: Multiphysics Modeling and Simulation of Deimos, an Advanced Reactor Experiment

When designing a novel reactor, nuclear experiments are essential to validate the predictive capability of modeling/simulation tools and to justify the investment of a full-scale prototype. Being able to accurately predict temperature reactivity coefficients is of high importance to reactor designers as it impacts the safety and performance of the system. Deimos is a proposed graphite-moderated, beryllium-reflected, high assay low-enriched uranium (HALEU) tri-structural isotropic (TRISO) fueled experiment for NCERC. Deimos is a valuable experiment for validating the predictive performance of modeling and simulation tools since it comprises a relatively thermal energy spectrum with 84% thermal neutrons causing fission, 15% epithermal neutrons causing fission, and 1% fast neutrons causing fission. Additionally, by electrically heating the experiment, reactivity can be validated at various temperatures. Deimos will serve as a testbed for future advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Low Frequency Modulation of Extreme Temperature Regimes in a Changing Climate

The project examines long-term changes in extreme temperature episodes (ETE) associated with planetary climate modes (PCMs) in both the real atmospheric and climate model simulations. The focus is on cold air outbreaks (CAOs) and warm waves (WWs) occurring over the continental US during the past 60 winters. No significant long-term trends in either WWs or CAOs are observed over the US. The annual frequency of CAOs is affected by the (i) North Atlantic Oscillation (NAO) over the Southeast US and (ii) Pacific–North American (PNA) pattern over the Northwest US. WW frequency is influenced by the (i) NAO over the eastern US and (ii) combined influence of PNA, Pacific decadal oscillation (PDO), and ENSO over the southern US. The collective influence of PCMs accounts for as much as 50% of the regional variability in ETE frequency. During CAO (WW) events occurring over the southeast US, there are low (high) pressure anomalies at higher atmospheric levels over the southeast US with oppositely-signed pressure anomalies in the lower atmosphere over the central US. These patterns lead to anomalous northerly (for CAOs) or southerly (for WWs) flow into the southeast leading to cold or warm surface air temperature anomalies, respectively. One distinction is that CAOs involve substantial air mass transport while WW formation is more local in nature. The primary differences among event categories are in the origin and nature of the pressure anomaly features linked to ETE onset. In some cases, PCMs help to provide a favorable environment for event onset. Heat budget analyses indicate that latitudinal transport in the lower atmosphere is the main contributor to regional cooling during CAO onset. This is partly offset by adiabatic warming associated with subsiding air. Additional diagnoses reveal that this latitudinal transport is partly due to the remote physical influence of a shallow cold pool of air trapped along the east side of the Rocky Mountains. ETE and PCM behavior is also studied in (CMIP5) climate model simulations. Although the climate models considered are able to represent the overall behavior of ETEs, the frequency of WWs (CAOs) is too high (low) in many models. While all models qualitatively replicate the overall structure of the PNA pattern, a small minority of models fails to properly simulate the NAO pattern. Model shortcomings in representing the NAO and PNA patterns have important consequences for simulating associated regional variability in surface air temperature and storm track behavior. The influence of PCMs on ETEs is underestimated in most CMIP5 models. In particular, none of the models are able to accurately simulate observed linkages between ETEs and the PDO, due to a gross misrepresentation of the PDO pattern in most models. Our results indicate that predictions of future CAO and WW behavior are currently limited by the ability of climate models to accurately represent PCM characteristics. Our study also considers the behavior of PCMs known as annular modes. It is determined that north-south movements in the stratospheric jet stream (related to the Polar Annular Mode) result in long-lasting impacts upon surface weather conditions including regional air temperature anomalies. The structure and dynamics of the stratospheric northern annular mode (or SNAM, related to changes in the strength of the stratospheric jet stream) was studied in CMIP5 models. In models with poorly-resolved stratospheres, the amplitude of SNAM at stratospheric altitudes is typically too weak, consistent with weaker stratospheric jet variability. However, this distinction does not carry over to the associated tropospheric signature of SNAM. A regional analysis illustrates that most CMIP5 models (regardless of whether the stratosphere is well-resolved) have anomalously weak and eastward shifted (compared to observed SNAM events) storm track and sea level pressure anomaly patterns during SNAM events. Analyses of stratosphere–troposphere coupling reveal that large-scale wave activity in the stratosphere is anomalously weak in CMIP5 models having poorly-resolved stratospheres, due to correspondingly weak upward propagation of wave activity from the troposphere. This suggests that such models have anomalously weak stratosphere-troposphere coupling in association with SNAM events. Paradoxically, however, there is little apparent impact upon the attendant tropospheric variability.

54 ENVIRONMENTAL SCIENCES↗

Assessment of Modeling and Simulation Technical Gaps in Safety Analysis of High Burnup Accident-Tolerant Fuels

The United States nuclear industry is facing a strong challenge to maintain regulatory-required levels of safety while ensuring economic competitiveness to stay in business. Safety remains a key parameter for all aspects related to the operation of light water reactor (LWR) nuclear power plants (NPPs), and it can be achieved more economically by using a risk-informed ecosystem, such as that being developed by the Risk-Informed Systems Analysis (RISA) Pathway under the U.S. Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program. The LWRS Program is promoting a wide range of research and development activities to maximize both the safety and economically efficient performance of NPPs through improved scientific understanding, especially given that many plants are considering second license renewal. The RISA Pathway has two main goals: The deployment of methodologies and technologies that enable better representation of safety margins and the factors that contribute to cost and safety, and; The development of advanced applications that enable cost-effective plant operation. As part of the RISA Pathway, the Enhanced Resilient Plant (ERP) project refers to an NPP where safety is improved by implementing various measures, such as accident-tolerant fuels (ATF), diverse and flexible coping strategy (FLEX), enhancements to plant components and systems, incorporation of augmented or new passive cooling systems, and utilization of advanced battery technologies. The objective of the ERP research is to use novel methods and computational tools to enhance existing reactors’ safety while reducing operational costs. Many U.S. utilities are targeting implementation of ATFs instead of traditional fuel in the near future since ATFs offer benefits in terms of improved performance and cost savings. The robust properties of ATF make it possible to extend the refueling cycle from 18 to 24 months in addition to the opportunity to use less of fuel. Extensive safety assessments are required to support regulatory requirements and obtain the approvals to use ATFs and the ERP project support the industry by developing novel effective methodologies for safety evaluations. In this project, the technical gaps in the modeling and simulation of the high burnup (HBU) ATF were assessed in terms of the fuel cladding behavior during the postulated accident events. The issues were identified in modeling the cladding deformation, the hydrodynamic change due to cladding deformation and the critical heat flux (CHF). The RELAP5-3D cladding deformation model was assessed by multiple verification tests and validation with the instrumented fuel assembly (IFA) experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigation of observational error sources in multi-Doppler-radar three-dimensional variational vertical air motion retrievals

Abstract. Multi-Doppler-radar network observations have been used in different configurations over the last several decades to conduct three-dimensional wind retrievals in mesoscale convective systems. Here, the impacts of the selected radar volume coverage pattern (VCP), the sampling time for the VCP, the number of radars used, and the added value of advection correction on the retrieval of the vertical air motion in the upper part of convective clouds are examined using the Weather Research and Forecasting (WRF) model simulation, the Cloud Resolving Model Radar SIMulator (CR-SIM), and a three-dimensional variational multi-Doppler-radar retrieval technique. Comparisons between the model truth (i.e., WRF kinematic fields) and updraft properties (updraft fraction, updraft magnitude, and mass flux) retrieved from the CR-SIM-generated multi-Doppler-radar field are used to investigate these impacts. The findings are that (1) the VCP elevation strategy and sampling time have a significant effect on the retrieved updraft properties above 6  km in altitude; (2) 2  min or shorter VCPs have small impacts on the retrievals, and the errors are comparable to retrievals using a snapshot cloud field; (3) increasing the density of elevation angles in the VCP appears to be more effective to reduce the uncertainty than an addition of data from one more radar, if the VCP is performed in 2  min ; and (4) the use of dense elevation angles combined with an advection correction applied to the 2  min VCPs can effectively improve the updraft retrievals, but for longer VCP sampling periods (5  min ) the value of advection correction is challenging. This study highlights several limiting factors in the retrieval of upper-level vertical velocity from multi-Doppler-radar networks and suggests that the use of rapid-scan radars can substantially improve the quality of wind retrievals if conducted in a limited spatial domain.

54 ENVIRONMENTAL SCIENCES↗

Space-based Doppler lidar sampling strategies: Algorithm development and simulated observation experiments

Lidar Atmospheric Wind Sounder (LAWS) Simulation Models (LSM) were developed to evaluate the potential impact of global wind observations on the basic understanding of the Earth's atmosphere and on the predictive skills of current forecast models (GCM and regional scale). Fully integrated top to bottom LAWS Simulation Models for global and regional scale simulations were developed. The algorithm development incorporated the effects of aerosols, water vapor, clouds, terrain, and atmospheric turbulence into the models. Other additions include a new satellite orbiter, signal processor, line of sight uncertainty model, new Multi-Paired Algorithm and wind error analysis code. An atmospheric wind field library containing control fields, meteorological fields, phenomena fields, and new European Center for Medium Range Weather Forecasting (ECMWF) data was also added. The LSM was used to address some key LAWS issues and trades such as accuracy and interpretation of LAWS information, data density, signal strength, cloud obscuration, and temporal data resolution.

Emmitt, G. D.↗

Physics guided machine learning for multi-material decomposition of tissues from dual-energy CT scans of simulated breast models with calcifications

We introduce a physics guided data-driven method for image-based multi-material decomposition for dual-energy computed tomography (CT) scans. The method is demonstrated for CT scans of virtual human phantoms containing more than two types of tissues. The method is a physics-driven supervised learning technique. We take advantage of the mass attenuation coefficient of dense materials compared to that of muscle tissues to perform a preliminary extraction of the dense material from the images using unsupervised methods. We then perform supervised deep learning on the images processed by the extracted dense material to obtain the final multi-material tissue map. The method is demonstrated on simulated breast models with calcifications as the dense material placed amongst the muscle tissues. The physics-guided machine learning method accurately decomposes the various tissues from input images, achieving a normalized root-mean-squared error of 2.75%.

Gopalakrishnan Meena, Murali↗

COMSOL Models and Simulation Results for Closed-Loop Geothermal Systems with Thermal Enhancements

The Closed-Loop Geothermal Working Group and the National Renewable Energy Laboratory (NREL) have studied the performance of disc-shaped and linear thermal enhancements (conductive material introduced around a wellbore) in closed-loop geothermal systems. COMSOL models were developed to investigate the increase in heat production in systems with thermal enhancements compared to those without, across a range of scenarios. This dataset includes two example COMSOL models, a summary of the simulation results, and a link to an accompanying paper. Studied are idealized versions of two thermal enhancement designs: (1) existing or newly created fractures surrounding the wellbore, filled with highly conductive material, and (2) a 'fishbone structure,' where short wellbores are drilled radially outward from the main wellbore and filled with highly conductive material.

15 GEOTHERMAL ENERGY↗