Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized 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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Multi-Modular Coordination Control of HTR-PM600 Plant - 190

HTR-PM600 is a high temperature gas-cooled reactor plant based on the multi-modular scheme, i.e. six modular high temperature gas-cooled reactor (MHTGR) driving a common secondary loop system. Six MHTGR-based nuclear steam supplying system (NSSS) modules are tightly coupled by the secondary fluid flow network (FFN) which contains the secondary side of steam generators, feedwater system, main steam system, turbine generator system. The secondary FFN control, which focuses on the feedwater flowrates and main steam pressure, plays a central role in multi-modular coordination of HTR-PM600 plant. In this paper, a differential algebraic model of this secondary FFN is first proposed. Then, the corresponding flowrate-pressure control strategy is proposed, which takes the simple form of proportional-integral feedback laws. Numerical simulation results show both the feasibility and satisfactory transient performance of this controller. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

High Performance-Multidisciplinary Simulation for Regional Scale Earthquake Hazard and Risk Assessments

Moderate-to-large magnitude earthquakes represent a significant risk around the world and a major issue across the Department of Energy (DOE) mission interests, ranging from DOE's one-of-a-kind mission critical facilities to national energy systems (nuclear power plants, electric/gas distribution systems, and renewable energy production facilities). As a result, for the DOE there is a strong motivation in promoting and developing a risk-informed and performance-based approach for both the assessment of existing facilities and the design of new structures. DOE National Labs have historically made multidisciplinary contributions to the solution of challenging seismic hazard and risk problems and have been at the forefront for the development of design standards for major facilities. The current approach for estimating the ground motion hazard for critical facilities is based upon a Probabilistic Seismic Hazard Analysis (PSHA), whereby measured ground motion data are collected and homogenized into a probabilistic framework to estimate the hazard for future events. However, recent studies have highlighted the site-specific nature of earthquakes and shown the limitations of a purely empirical approach. Moreover, the design of infrastructure systems has been moving toward a performance-based approach, in which specific performance goals are specified at different levels of excitation and for different limit states. Such an approach requires the availability of robust and reliable computational tools for the nonlinear dynamic analysis of complex structural systems. Ongoing research at Lawrence Berkeley National Laboratory (LBNL) and Lawrence Livermore National Laboratory (LLNL) aims at developing an advanced simulation environment at the exascale for a coupled assessment of earthquake hazard (ground motions) and earthquake risk (structural demands)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Time Scale Effects on the Temperature Coefficient of Reactivity

For a nuclear reactor, the temperature coefficient of reactivity (TCR) is of significance for both power control and safety analysis, but it is not always a constant. For example, the TCR may be dependent upon temperature due to Doppler broadening, and may be also dependent upon time scale due to heat transfer and thermal expansion. A fast burst reactor is a favorable nuclear facility to simulate the time scale effects on the TCR. Fast burst reactors are generally unmoderated research reactors of enriched uranium metal or uranium alloys. The neutrons are fast and the Doppler broadening effects can be neglected, and the TCR is mainly contributed by fuel thermal expansion. Fast burst reactors can operate on both steady state power mode and burst power mode. In steady state power mode, the systems are assumed to be in thermal equilibrium. In burst mode, the full width at half maximum (FWHM) of the power is from tens of microseconds to a few milliseconds, and the heat transfer is generally neglected. If FWHM is short enough, thermal expansion will have a time delay due to inertial effect. This paper is aimed to study the exact time scale effects on the TCR due to thermal expansion. The TCR contributed by the Doppler broadening effects was neglected, because it is relatively small (The k{sub eff} as a function of temperature was also calculated using the Monte Carlo particle transport code JMCT by our group, and no definite difference due to the temperature change was obtained near room temperature), and there is no method to evaluate the time scale dependence of the Doppler broadening effects yet. In the time scale from 1 x 10{sup -6} s to 1 x 10{sup 5} s, the TCR of a Godiva-like spherical fast burst reactor was calculated with perturbation theory. A Godiva-like spherical fast burst reactor was studied, the calculation of the time scale effects on TCR involved thermal conduction, thermoelasticity and reactivity variation analysis. Both the thermal and the thermoelastic problems were solved using ANSYS software. ANSYS is software implementation of finite element simulations, which is widely used in engineering field. The reactivity variation due to thermal expansion was calculated with perturbation theory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Monte Carlo Study of the Core-Shell Phase Transitions in Nanostructure Systems

We present in this work the effects of the external magnetic field, the exchange anisotropy, the coupling exchange interaction and the temperature on the magnetic properties of the nanostructure systems. More exactly, the effect of the core-shell spin Ising systems has been also examined. A number of favorable phenomena such as the wetting magnetic and the layering transitions have been investigated. The Monte Carlo simulations have been used to investigate the behavior of these nanostructure systems and to obtain the interesting results presented in the form of the phase diagrams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Crystal Field Effects on Magnetic Properties of a Mixed Ferrimagnetic Ising System: a Monte Carlo Study

The effects of two different single-ion anisotropies on the magnetization of the mixed (7/2, 5/2) Ising ferrimagnetic system on a square lattice are investigated by the use of Monte Carlo simulation. The ground-state phase diagram is obtained. Compensation points are detected for certain values of the crystal fields. Different types of magnetization curves are available.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

A reaction–diffusion model for phenotypic evolution

We present a reaction–diffusion mathematical model for the evolutionary dynamics of phenotypic evolution. A detailed deduction of the equations is presented for the one-dimensional version, from which a more general model is proposed. Particular cases are studied using analytical approximations and numerical simulations. Results indicate that the approach proposed produces results that are coherent with mainstream models in evolutionary dynamics, suggesting that the reaction–diffusion model could be an alternative tool in the analysis of evolutionary dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Long-range force and moment calculations in multiresolution simulations of molecular systems

Multiresolution simulations of molecular systems such as DNAs, RNAs, and proteins are implemented using models with different resolutions ranging from a fully atomistic model to coarse-grained molecules, or even to continuum level system descriptions. For such simulations, pairwise force calculation is a serious bottleneck which can impose a prohibitive amount of computational load on the simulation if not performed wisely. Herein, we approximate the resultant force due to long-range particle-body and body-body interactions applicable to multiresolution simulations. Since the resultant force does not necessarily act through the center of mass of the body, it creates a moment about the mass center. Although this potentially important torque is neglected in many coarse-grained models which only use particle dynamics to formulate the dynamics of the system, it should be calculated and used when coarse-grained simulations are performed in a multibody scheme. Herein, the approximation for this moment due to far-field particle-body and body-body interactions is also provided.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2013

Assimilating irregularly spaced sparsely observed turbulent signals with hierarchical Bayesian reduced stochastic filters

In this paper, we consider a practical filtering approach for assimilating irregularly spaced, sparsely observed turbulent signals through a hierarchical Bayesian reduced stochastic filtering framework. The proposed hierarchical Bayesian approach consists of two steps, blending a data-driven interpolation scheme and the Mean Stochastic Model (MSM) filter. We examine the potential of using the deterministic piecewise linear interpolation scheme and the ordinary kriging scheme in interpolating irregularly spaced raw data to regularly spaced processed data and the importance of dynamical constraint (through MSM) in filtering the processed data on a numerically stiff state estimation problem. In particular, we test this approach on a two-layer quasi-geostrophic model in a two-dimensional domain with a small radius of deformation to mimic ocean turbulence. Our numerical results suggest that the dynamical constraint becomes important when the observation noise variance is large. Second, we find that the filtered estimates with ordinary kriging are superior to those with linear interpolation when observation networks are not too sparse; such robust results are found from numerical simulations with many randomly simulated irregularly spaced observation networks, various observation time intervals, and observation error variances. Third, when the observation network is very sparse, we find that both the kriging and linear interpolations are comparable.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Resolved-particle simulation by the Physalis method: Enhancements and new capabilities

We present enhancements and new capabilities of the Physalis method for simulating disperse multiphase flows using particle-resolved simulation. The current work enhances the previous method by incorporating a new type of pressure-Poisson solver that couples with a new Physalis particle pressure boundary condition scheme and a new particle interior treatment to significantly improve overall numerical efficiency. Further, we implement a more efficient method of calculating the Physalis scalar products and incorporate short-range particle interaction models. We provide validation and benchmarking for the Physalis method against experiments of a sedimenting particle and of normal wall collisions. We conclude with an illustrative simulation of 2048 particles sedimenting in a duct. In the appendix, we present a complete and self-consistent description of the analytical development and numerical methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

An unstructured-mesh finite-volume MPDATA for compressible atmospheric dynamics

An advancement of the unstructured-mesh finite-volume MPDATA (Multidimensional Positive Definite Advection Transport Algorithm) is presented that formulates the error-compensative pseudo-velocity of the scheme to rely only on face-normal advective fluxes to the dual cells, in contrast to the full vector employed in previous implementations. This is essentially achieved by expressing the temporal truncation error underlying the pseudo-velocity in a form consistent with the flux-divergence of the governing conservation law. The development is especially important for integrating fluid dynamics equations on non-rectilinear meshes whenever face-normal advective mass fluxes are employed for transport compatible with mass continuity—the latter being essential for flux-form schemes. In particular, the proposed formulation enables large-time-step semi-implicit finite-volume integration of the compressible Euler equations using MPDATA on arbitrary hybrid computational meshes. Furthermore, it facilitates multiple error-compensative iterations of the finite-volume MPDATA and improved overall accuracy. The advancement combines straightforwardly with earlier developments, such as the nonoscillatory option, the infinite-gauge variant, and moving curvilinear meshes. A comprehensive description of the scheme is provided for a hybrid horizontally-unstructured vertically-structured computational mesh for efficient global atmospheric flow modelling. The proposed finite-volume MPDATA is verified using selected 3D global atmospheric benchmark simulations, representative of hydrostatic and non-hydrostatic flow regimes. Besides the added capabilities, the scheme retains fully the efficacy of established finite-volume MPDATA formulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Lipid droplets fusion in adipocyte differentiated 3T3-L1 cells: A Monte Carlo simulation

Several human worldwide diseases like obesity, type 2 diabetes, hepatic steatosis, atherosclerosis and other metabolic pathologies are related to the excessive accumulation of lipids in cells. Lipids accumulate in spherical cellular inclusions called lipid droplets (LDs) whose sizes range from fraction to one hundred of micrometers in adipocytes. It has been suggested that LDs can grow in size due to a fusion process by which a larger LD is obtained with spherical shape and volume equal to the sum of the progenitors’ ones. In this study, the size distribution of two populations of LDs was analyzed in immature and mature (5-days differentiated) 3T3-L1 adipocytes (first and second populations, respectively) after Oil Red O staining. A Monte Carlo simulation of interaction between LDs has been developed in order to quantify the size distribution and the number of fusion events needed to obtain the distribution of the second population size starting from the first one. Four models are presented here based on different kinds of interaction: a surface weighted interaction (R2 Model), a volume weighted interaction (R3 Model), a random interaction (Random model) and an interaction related to the place where the LDs are born (Nearest Model). The last two models mimic quite well the behavior found in the experimental data. This work represents a first step in developing numerical simulations of the LDs growth process. Due to the complex phenomena involving LDs (absorption, growth through additional neutral lipid deposition in existing droplets, de novo formation and catabolism) the study focuses on the fusion process. The results suggest that, to obtain the observed size distribution, a number of fusion events comparable with the number of LDs themselves is needed. Moreover the MC approach results a powerful tool for investigating the LDs growth process. Highlights: • We evaluated the role of the fusion process in the synthesis of the lipid droplets. • We compared the size distribution of the lipid droplets in immature and mature cells. • We used the Monte Carlo simulation approach, simulating 10 thousand of fusion events. • Four different interaction models between the lipid droplets were tested. • The best model which mimics the experimental measures was selected.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Ion velocity distributions in the sheath and presheath of a biased object in plasma

Ion velocity distributions in the vicinity of a spherical object with a negative potential with respect to collisionless, source-free plasma are studied with three-dimensional numerical simulations. The ion dynamics around the object leads to distorted radial velocity distributions in the presheath and the sheath edge region. Far in the sheath, an increase in the thermal velocity in the radial direction is observed. Different potentials of the object, ion temperatures, and ion masses are considered, as well as the role of spatial and temporal resolutions in laboratory measurements of ion velocity distributions. The simulations are carried out with the DiP3D, a three-dimensional particle-in-cell numerical code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Femtosecond laser-plasma interaction with prepulse-generated liquid metal microjets

Ultrashort laser pulse interaction with a microstructured surface of a melted metal is a promising source of hard x-ray radiation. Microstructuring is achieved by a weak prepulse that produces narrow high-density microjets. As an x-ray source, the interaction of the main laser pulse with such jets is shown to be nearly two orders of magnitude more efficient than the interaction with ordinary metal targets. This paper presents the results of optical and x-ray studies of laser-plasma interaction physics under such conditions supported by numerical simulations of microjet formation and fast-electron generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Implications for the electron distribution from the stationary hydrodynamic model of a one-dimensional plasma expansion into vacuum

It is shown that the hydrodynamic model of a one-dimensional collisionless plasma expansion is contained in the kinetic description as a special case. This belongs to a specific choice for the electron distribution function. Moreover, the consequences of the use of the hydrodynamic approach regarding the temporal evolution of the electron phase space density are investigated. It turns out that only the case of a hydrodynamic description with the adiabatic constant {kappa}=3 is physically self-consistent. Numerical simulations confirm this argumentation. The analysis for the case {kappa}=3 is extended to the kinetics of a relativistic electron gas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Polymer-loaded propagating modes on a one-dimensional photonic crystal

We numerically and experimentally demonstrate that a polymer film-coated one-dimensional photonic crystal (1DPC) can sustain transverse electric (TE) polarized modes without the limit of guided layer's thickness. Our results indicate that two propagating modes are existing inside the polymer film, the first one is the TE polarized Bloch surface wave, and the second one is the TE polarized guided mode. Here in, the evolution of these two modes with change in the polymer film thickness is presented. Our numerical simulation results are in well-agreement with the experimental data obtained using back focal plane imaging.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

All-dielectric three-dimensional broadband Eaton lens with large refractive index range

We proposed a method to realize three-dimensional (3D) gradient index (GRIN) devices requiring large refractive index (RI) range with broadband performance. By combining non-resonant GRIN woodpile photonic crystals structure in the metamaterial regime with a compound liquid medium, a wide RI range (1–6.32) was fulfilled flexibly. As a proof-of-principle for the low-loss and non-dispersive method, a 3D Eaton lens was designed and fabricated based on 3D printing process. Full-wave simulation and experiment validated its omnidirectional wave bending effects in a broad bandwidth covering Ku band (12 GHz–18 GHz)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Numerical solution of the time dependent Ginzburg-Landau equations for mixed (d + s)-wave superconductors

The time-dependent Ginzburg-Landau formalism for (d + s)-wave superconductors and their representation using auxiliary fields is investigated. By using the link variable method, we then develop suitable discretization of these equations. Numerical simulations are carried out for a mesoscopic superconductor in a homogeneous perpendicular magnetic field which revealed peculiar vortex states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗