Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulators”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Development of a Kalman filter estimator for simulation and control of particulate matter distribution of a diesel catalyzed particulate filter

The knowledge of the temperature and particulate matter mass distribution is essential for monitoring the performance and durability of a catalyzed particulate filter. A catalyzed particulate filter model was developed, and it showed capability to accurately predict temperature and particulate matter mass distribution and pressure drop across the catalyzed particulate filter. However, the high-fidelity model is computationally demanding. Therefore, a reduced order multi-zone particulate filter model was developed to reduce computational complexity with an acceptable level of accuracy. In order to develop a reduced order model, a parametric study was carried out to determine the number of zones necessary for aftertreatment control applications. The catalyzed particulate filter model was further reduced by carrying out a sensitivity study of the selected model assumptions. The reduced order multi-zone particulate filter model with 5 × 5 zones was selected to develop a catalyzed particulate filter state estimator considering its computational time and accuracy. Next, a Kalman filter–based catalyzed particulate filter estimator was developed to estimate unknown states of the catalyzed particulate filter such as temperature and particulate matter mass distribution and pressure drop (Δ P) using the sensor inputs to the engine electronic control unit and the reduced order multi-zone particulate filter model. A diesel oxidation catalyst estimator was also integrated with the catalyzed particulate filter estimator in order to provide estimates of diesel oxidation catalyst outlet concentrations of NO 2 and hydrocarbons and inlet temperature for the catalyzed particulate filter estimator. The combined diesel oxidation catalyst–catalyzed particulate filter estimator was validated for an active regeneration experiment. The validation results for catalyzed particulate filter temperature distribution showed that the root mean square temperature error by using the diesel oxidation catalyst–catalyzed particulate filter estimator is within 3.2 °C compared to the experimental data. Similarly, the Δ P estimator closely simulated the measured total Δ P and the estimated cake pressure drop error is within 0.2 kPa compared to the high-fidelity catalyzed particulate filter model.

Engineering↗
DOE OSTI2021

Novel Geochemistry Determined from High Pressure, High Temperature Simulation Experiments of Hydraulic Fracture Test Site 2

A standard method in unconventional oil and gas production is the process of hydraulic fracturing followed by a shut-in period, during which the fracture fluid remains pressurized in the reservoir for up to three weeks before production begins. Despite this widely used process, very little is known about what occurs in the reservoir during this shut-in process. In order to properly delineate potential reservoir reactions that may occur during shut-in that would lead to corrosion and scaling events, experiments were conducted in high pressure, high temperature reactors to simulate conditions of the Wolfcamp Formation in the Delaware Basin. Experimental design allowed the assessment of the effect of proppant, microbiology, and time on the mineralogy and fluid chemistry in the reservoir. Results suggest the biggest impact on fluid chemistry and shale mineralogy during shut-in is time. Analyses demonstrate dissolution of the shale material, with maximum dissolved ions occurring after 7 days of the shut-in period. After 21 days, results suggest precipitation occurs. The Delaware Basin is demonstrated to be high in sulfate content, which further increases in the fluid due to dissolution reactions during shut-in. The pressure vessel experiments suggest there was no significant contribution to reactions from microbiology during the shut-in period. Furthermore, early production samples demonstrate a significant selection of the microorganism Caminicella, a genus of which has previously been correlated to corrosion and sulfide production. Results suggest shut-in conditions may provide high sulfate concentrations that could later be utilized by a shifted microbial community to drive potential well infrastructure failure. This is the first study to incorporate microbiology, mineralogy, and fluid chemistry to investigate the fundamental geochemical reactions that occur during shut-in and early phase production of the Delaware Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/2-21URTC/D021S032R001/2477425/urtec-2021-5219-ms.pdf/1 by Carol Worster on 28 February 2022 Basin. Results from this study can complement observations from the Hydraulic Fracture Test Site 2 observations

58 GEOSCIENCES↗
DOE OSTI2017

Molecular Dynamics Simulation of Fission Fragment Damage in Nuclear Fuel and Surrogate Material

We have performed classical molecular dynamics simulations of swift heavy ion damage, typical of fission fragments, in nuclear fuel (UO2) for energy deposition per unit length of 3.9 keV/nm. We did not observe amorphization. The damage mainly consisted of isolated point defects. Only about 1% of the displacements occur on the uranium sublattice. Oxygen Frenkel pairs are an order of magnitude more numerous than uranium Frenkel pairs in the primary damage state. In contrast, previous results show that the ratio of Frenkel pairs on the two sublattices is close to the stoichiometric ratio in ceria. These differences in the primary damage state may lead to differences in radiation response of UO2 and CeO2.

Molecular dynamics, radiation damage in materials,↗
DOE OSTI2016

Simulating Fragmentation and Fluid-Induced Fracture in Disordered Media Using Random Finite-Element Meshes

Fracture and fragmentation are extremely nonlinear multiscale processes in which microscale damage mechanisms emerge at the macroscale as new fracture surfaces. Numerous numerical methods have been developed for simulating fracture initiation, propagation, and coalescence. In this paper, we present a computational approach for modeling pervasive fracture in quasi-brittle materials based on random close-packed Voronoi tessellations. Each Voronoi cell is formulated as a polyhedral finite element containing an arbitrary number of vertices and faces. Fracture surfaces are allowed to nucleate only at the intercell faces. Cohesive softening tractions are applied to new fracture surfaces in order to model the energy dissipated during fracture growth. The randomly seeded Voronoi cells provide a regularized discrete random network for representing fracture surfaces. The potential crack paths within the random network are viewed as instances of realizable crack paths within the continuum material. Mesh convergence of fracture simulations is viewed in a weak, or distributional, sense. The explicit facet representation of fractures within this approach is advantageous for modeling contact on new fracture surfaces and fluid flow within the evolving fracture network. Finally, applications of interest include fracture and fragmentation in quasi-brittle materials and geomechanical applications such as hydraulic fracturing, engineered geothermal systems, compressed-air energy storage, and carbon sequestration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

A Finite Element Model for Simulation of Carbon Dioxide Sequestration

We present a hydro-mechanical model, followed by stress, deformation, and shear-slip failure analysis for geological sequestration of carbon dioxide (CO2). The model considers the poroelastic effects by taking into account of the two-way coupling between the geomechanical response and the fluid flow process. Analytical solutions for pressure and deformation fields were derived for a typical geological sequestration scenario in our previous work. A finite element approach is introduced here for numerically solving the hydro-mechanical model with arbitrary boundary conditions. The numerical approach was built on an open-source finite element code Elmer, and results were compared to the analytical solutions. The shear-slip failure analysis was presented based on the numerical results, where the potential failure zone is identified. Information is relevant to the prediction of the maximum sustainable injection rate or pressure. The effects of caprock permeability on the fluid pressure, deformation, stress, and the shear-slip failure zone were also quantitatively studied. It was shown that a larger permeability in caprock and base rock leads to a larger uplift but a smaller shear-slip failure zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Whitepaper prepared for DOE Workshop on Integrated Simulations for Magnetic Fusion Energy Sciences: Topic E: Beyond interpretive simulations

Predictive simulations of magnetic confinement fusion and burning plasmas will require the quantification of all uncertainties and errors relating to the simulation capabilities. These include: discretization error (temporal and spatial); incomplete convergence error (nonlinear, linear, etc.); uncertainties in input data (initial conditions, boundary conditions, coefficients, thermo-physical properties, source terms, etc.); and uncertainties in the component models (the specific form and parameter values).13 Typically, there is a focus on model parameter uncertainties but uncertainties about which model to use and about the bias or discrepancy of a particular model (sometimes called model form error) often dominate parameter uncertainty. This whitepaper addresses the challenges of performing uncertainty quantification (UQ) for expensive computational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Converting DYNAMO simulations to Powersim Studio simulations

DYNAMO is a computer program for building and running 'continuous' simulation models. It was developed by the Industrial Dynamics Group at the Massachusetts Institute of Technology for simulating dynamic feedback models of business, economic, and social systems. The history of the system dynamics method since 1957 includes many classic models built in DYANMO. It was not until the late 1980s that software was built to take advantage of the rise of personal computers and graphical user interfaces that DYNAMO was supplanted. There is much learning and insight to be gained from examining the DYANMO models and their accompanying research papers. We believe that it is a worthwhile exercise to convert DYNAMO models to more recent software packages. We have made an attempt to make it easier to turn these models into a more current system dynamics software language, Powersim © Studio produced by Powersim AS 2 of Bergen, Norway. This guide shows how to convert DYNAMO syntax into Studio syntax.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

NEK5000 Assessment Milestone Report: Single-Phase Natural Circulation using Direct Numerical Simulation (DNS) & Large Eddy Simulation (LES) Methods

This study evaluates the capabilities of the NEK5000 open-source code as an alternative to the commercially available software that is currently heavily relied upon. The first three case studies initiated the author with (1) the syntax and flow of the program, (2) varying initial and boundary conditions and (3) geometries of greater complexity & dimension. Case 1 studies transient conjugate heat transfer in the entry region of a pipe. Cases 2 & 3 simulate the natural convection flow in an enclosed differentially heated cavity (DHC) of two and three dimensions, respectively. The final case study, Case 4, produces three-dimensional single-phase natural convective flow using the Boussinesq approximation within a simplified representation of an aqueous homogeneous reactor system. The Boussinesq flow in Case 4 clearly depicts both Rayleigh-Taylor and Kelvin-Helmholtz instabilities qualitatively during post-processing analysis. Further analysis of Case 4, up to the point of a quasi-steady state (i.e. thermal-equilibrium condition), is done using the LES solver method in NEK5000. The volume-averaged solution temperatures, when compared between solvers, showed no significant difference. Due to its minimal compromise in fidelity of solution resolution and inexpensive computing cost, compared to the DNS solver, the LES solver is recommended for the current application. Methodology development for the four case studies are reported and a preliminary optimization study with different polynomial resolutions for Case 4 is addressed as a resource when performing future simulations using NEK5000. Hence, this report can be used as a reference starting point to the end user of NEK5000 for natural circulation phenomena in nuclear reactor safety applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

An Approach for Simulating Transport of Nonspherical, Sublimating Particles in Rarefied Gas Flows via the Direct Simulation Monte Carlo Method

This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Simulation of Wingtip Vortex Flows with Reynolds-Averaged Navier–Stokes and Scale-Resolving Simulation Methods

The simulation of wingtip vortex flows is reported on with Reynolds-averaged Navier–Stokes (RANS) equations and scale-resolving simulation (SRS) models. These range from turbulent viscosity and Reynolds-stress model (RSM) RANS closures, to hybrid and bridging SRS methods. The aim of the study was threefold: assess the relevance of replicating the experimental flow conditions, evaluate the numerical requisites of each mathematical model, and determine their modeling accuracy in prediction of a representative wingtip vortex flow. The selected problem is the flow around a NACA 0012 wing at 10 deg of angle of attack and Reynolds number of 4.60×10 6 . The results confirm the relevance of reproducing the experimental flow conditions, in particular at the inlet boundary where the flow is not uniform. For this reason, the evaluation of modeling errors using the available measurements requires the specification of the experimental inlet conditions. On the other hand, the quantification of the modeling error indicates that solely RANS–RSM can achieve an accurate representation of the flow dynamics. Whereas the well-recognized limitations of turbulent viscosity RANS closures to deal with solid-body rotation lead to the overprediction of turbulence and consequent rapid diffusion of the wingtip vortex, the predictive use of SRS methods may reveal excessively complex and demanding.

42 ENGINEERING↗
DOE OSTI2020

Investigation of Mixing Law Efficacy for Gaseous Hydrodynamic Simulations

A computational simulation of various mixing laws for gaseous equations of state using planar traveling shocks for multiple mixtures in three dimensions is analyzed against nominal experimental data. Numerical simulations utilize the Sandia National Laboratories shock hydrodynamic code CTH and other codes including the thermochemical equilibrium code TIGER and the uncertainty qualification and sensitivity analysis code DAKOTA. The mixtures are a 1:1 and a 1:3 molar mixture of helium and sulfur hexafluoride. The mixing laws to be analyzed are the ideal gas law, Amagat's Law, Dalton's Law, the BKW EOS, the EXP6 EOS, and the JCZ3 EOS. Examination of the experimental data with TIGER revealed that the shock strength should not be strong enough to turn the mixture non-ideal as the compressibility factor z was essentially unity (z ≈ 1.02). Experimental results show that none of the equations of state are able to accurately predict the properties of the shocked mixture; similar discrepancies have been observed in previous works. Kinetic molecular theory appears to introduce a parameter which offers an explanation regarding the discrepancies. Implementation of the KMT pararneter into the EOS is left for future work.

42 ENGINEERING↗
DOE OSTI2020

Simulating the Trajectory and Biomass Growth of Free-Floating Macroalgal Cultivation Platforms along the U.S. West Coast

Trajectory tracking and macroalgal growth models were coupled to support a novel macroalgae-harvesting concept known as the Nautical Off-shore Macroalgal Autonomous Device (NOMAD). The NOMAD consists of 5 km long carbon-fiber longlines that are seeded and free float southward along the U.S. West Coast for approximately 3 months before harvesting off the California coast, taking advantage of favorable environmental conditions. The trajectory and macroalgal growth models were applied to answer planning questions pertinent to the techno-economic analysis such as identifying the preferred release location, approximate pathway, timing until harvest, and estimated growth. Trajectories were determined with the General NOAA Operational Modeling Environment (GNOME) model, using 11 years of current and wind data, determining probabilities by running nearly 40,000 Monte Carlo simulations varying the start time and location. An accompanying macroalgal growth model was used to estimate the growth of macroalgae based on the trajectory tracks and environmental forcing products, including light, temperature and nutrients. Model results show that NOMAD lines transit south in the months of April to September due to seasonal currents, taking approximately 3 months to reach Southern California. During transit, NOMAD lines are dispersed but typically avoid beaching or passing through marine sanctuaries. NOMAD lines can yield up to 30 kg wet weight per meter of cultivation line.

09 BIOMASS FUELS↗
DOE OSTI2021

As-Built Simulation of the High Flux Isotope Reactor

The Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) is an 85 MWt flux trap-type research reactor that supports key research missions, including isotope production, materials irradiation, and neutron scattering. The core consists of an inner and an outer fuel element containing 171 and 369 involute-shaped plates, respectively. The thin fuel plates consist of a U 3 O 8 -Al dispersion fuel (highly enriched), an aluminum-based filler, and aluminum cladding. The fuel meat thickness is varied across the width of the involute plate to reduce thermal flux peaks at the radial edges of the fuel elements. Some deviation from the designed fuel meat shaping is allowed during manufacturing. A homogeneity scan of each fuel plate checks for potential anomalies in the fuel distribution by scanning the surface of the plate and comparing the attenuation of the beam to calibration standards. While typical HFIR simulations use homogenized fuel regions, explicit models of the plates were developed under the Low-Enriched Uranium Conversion Program. These explicit models typically include one inner and one outer fuel plate with nominal fuel distributions, and then the plates are duplicated to fill the space of the corresponding fuel element. Therefore, data extracted from these simulations are limited to azimuthally averaged quantities. To determine the reactivity and physics impacts of an as-built outer fuel element and generate azimuthally dependent data in the element, 369 unique fuel plate models were generated and positioned. This model generates the three-dimensional (i.e., radial–axial–azimuthal) plate power profile, where the azimuthal profile is impacted by features within the adjacent control element region and beryllium reflector. For an as-built model of the outer fuel element, plate-specific homogeneity data, 235 U loading, enrichment, and channel thickness measurements were translated into the model, yielding a much more varied azimuthal power profile encompassed by uncertainty factors in analyses. These models were run with the ORNL-TN and Shift Monte Carlo tools, and they contained upwards of 500,000 cells and 100,000 unique tallies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Prediction of Metabolite Concentrations, Rate Constants and Post-Translational Regulation Using Maximum Entropy-Based Simulations with Application to Central Metabolism of Neurospora crassa

We report the application of a recently proposed approach for modeling biological systems using a maximum entropy production rate principle in lieu of having in vivo rate constants. The method is applied in four steps: (1) a new ordinary differential equation (ODE) based optimization approach based on Marcelin’s 1910 mass action equation is used to obtain the maximum entropy distribution; (2) the predicted metabolite concentrations are compared to those generally expected from experiments using a loss function from which post-translational regulation of enzymes is inferred; (3) the system is re-optimized with the inferred regulation from which rate constants are determined from the metabolite concentrations and reaction fluxes; and finally (4) a full ODE-based, mass action simulation with rate parameters and allosteric regulation is obtained. From the last step, the power characteristics and resistance of each reaction can be determined. The method is applied to the central metabolism of Neurospora crassa and the flow of material through the three competing pathways of upper glycolysis, the non-oxidative pentose phosphate pathway, and the oxidative pentose phosphate pathway are evaluated as a function of the NADP/NADPH ratio. It is predicted that regulation of phosphofructokinase (PFK) and flow through the pentose phosphate pathway are essential for preventing an extreme level of fructose 1,6-bisphophate accumulation. Such an extreme level of fructose 1,6-bisphophate would otherwise result in a glassy cytoplasm with limited diffusion, dramatically decreasing the entropy and energy production rate and, consequently, biological competitiveness.

Cannon, William R.↗
DOE OSTI2017

Precursor Evolution and Stress Corrosion Cracking Initiation of Cold-Worked Alloy 690 in Simulated Pressurized Water Reactor Primary Water

Stress corrosion crack initiation of two thermally-treated, cold-worked (CW) alloy 690 (UNS N06690) materials was investigated in 360oC simulated PWR primary water using constant load tensile (CLT) tests and blunt notch compact tension (BNCT) tests equipped with direct current potential drop (DCPD) for in-situ detection of cracking. SCC initiation was not detected by DCPD for either the 21% or 31%CW CLT specimens loaded at their yield stress after ~9,220 hours, however intergranular (IG) precursor damage and isolated surface cracks were observed on the specimens. The two 31%CW BNCT specimens loaded at moderate stress intensity after several cyclic loading ramps showed DCPD-indicated crack initiation after 10,400 hours of exposure at constant stress intensity, which was resulted from significant growth of IG cracks. The 21%CW BNCT specimens only exhibited isolated small IG surface cracks and showed no apparent DCPD change throughout the test. Post-test cross-section examinations revealed many grain boundary (GB) nano-cavities in the bulk of all the CLT and BNCT specimens particularly for the 31%CW materials. Cavities were also found along GBs extending to the surface suggesting an important role in crack nucleation. This paper provides an overview of the evolution of GB cavities and discusses their effects on crack initiation in CW alloy 690.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Impacts of aerosols on seasonal precipitation and snowpack in California based on convection-permitting WRF-Chem simulations

Here, a version of the WRF-Chem model with fully coupled aerosol–meteorology–snowpack is employed to investigate the impacts of various aerosol sources on precipitation and snowpack in California. In particular, the impacts of locally emitted anthropogenic and dust aerosols, and aerosols transported from outside California are studied. We differentiate three pathways of aerosol effects: aerosol–radiation interaction (ARI), aerosol–snow interaction (ASI), and aerosol–cloud interaction (ACI). The convection-permitting model simulations show that precipitation, snow water equivalent (SWE), and surface air temperature averaged over the whole domain (34–42° N, 117–124° W, not including ocean points) are reduced when aerosols are included, therefore reducing large biases in these variables due to the absence of aerosol effects in the model. Aerosols affect California water resources through the warming of mountaintops and the reduction of precipitation; however, different aerosol sources play different roles in changing surface temperature, precipitation, and snowpack in California by means of various weights of the three pathways. ARI by all aerosols mainly cools the surface, leading to slightly increased SWE over the mountains. Locally emitted dust aerosols warm the surface of mountaintops through ASI, in which the reduced snow albedo associated with dusty snow leads to more surface absorption of solar radiation and reduced SWE. Transported aerosols and local anthropogenic aerosols play a dominant role in increasing nonprecipitating clouds but reducing precipitation through ACI, leading to reduced SWE and runoff on the Sierra Nevada, as well as the warming of mountaintops associated with decreased SWE and hence lower surface albedo. The average changes in surface temperature from October 2012 to June 2013 are about –0.19 and 0.22 K for the whole domain and over mountaintops, respectively. Overall, the averaged reduction during October to June is about 7 % for precipitation, 3 % for SWE, and 7 % for surface runoff for the whole domain, while the corresponding numbers are 12, 10, and 10 % for the mountaintops. The reduction in SWE is more significant in a dry year, with 9 % for the whole domain and 16 % for the mountaintops. The maximum reduction of ~20 % in precipitation occurs in May and is associated with the maximum aerosol loading, leading to the largest decrease in SWE and surface runoff over that period. It is also found that dust aerosols can cause early snowmelt on the mountaintops and reduced surface runoff after April.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Simulation Speed Analysis and Improvements of Modelica Models for Building Energy Simulation

This paper presents an approach for speeding up Modelica models. Insight is provided into how Modelica models are solved and what determines the tool’s computational speed. Aspects such as algebraic loops, code efficiency and integrator choice are discussed. This is illustrated using simple building simulation examples and Dymola. The generality of the work is in some cases verified using OpenModelica. Using this approach, a medium sized office building including building envelope, heating ventilation and air conditioning (HVAC) systems and control strategy can be simulated at a speed five hundred times faster than real time.

Jorissen, Filip↗
DOE OSTI2012

Molecular simulation as a tool for studying lignin

Lignocellulosic biomass provides a sustainable source of sugars for biofuel and biomaterial production. However, biomass resistance to degradation imposes difficulties for economical conversion of plant carbohydrates to fermentable sugars. One of the key contributors to recalcitrance is lignin. Understanding the properties of lignin macromolecules in the cell wall matrix is useful for manipulating biomass structure to generate more easily degradable biomass. Along with experimental techniques such as 2D-NMR and mass spectrometry, computational techniques can be useful for characterizing the structural and energetic properties of the biomass assembly and its individual constituents. Here, we provide a brief introduction to lignin, review some of the recent, relevant scientific literature, and give our perspectives on the role of molecular simulation in understanding lignin structure.

09 BIOMASS FUELS↗