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,315 records · Page 73

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Relation between energetic and standard geodesic acoustic modes

Geodesic Acoustic Modes (GAMs) are electrostatic, axisymmetric modes which are non-linearly excited by turbulence. They can also be excited linearly by fast-particles; they are then called Energetic-particle-driven GAMs (EGAMs). Do GAMs and EGAMs belong to the same mode branch? Through a linear, analytical model, in which the fast particles are represented by a Maxwellian bump-on-tail distribution function, we find that the answer depends on several parameters. For low values of the safety factor q and for high values of the fast ion energy, the EGAM originates from the GAM. On the contrary, for high values of q and for low values of the fast ion energy, the GAM is not the mode which becomes unstable when fast particles are added: the EGAM then originates from a distinct mode, which is strongly damped in the absence of fast particles. The impact of other parameters is further explored: ratio of the ion temperature to the electron temperature, width of the fast particle distribution, mass and charge of the fast ions. The ratio between the EGAM and the GAM frequencies was found in experiments (DIII-D) and in non-linear numerical simulations (code GYSELA) to be close to 1/2: the present analytical study allows one to recover this ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron attenuation in free, neutral ethane clusters

The electron effective attenuation length (EAL) in free, neutral ethane clusters has been determined at 40 eV kinetic energy by combining carbon 1s x-ray photoelectron spectroscopy and theoretical lineshape modeling. More specifically, theory is employed to form model spectra on a grid in cluster size (N) and EAL (λ), allowing N and λ to be determined by optimizing the goodness-of-fit χ{sup 2}(N, λ) between model and observed spectra. Experimentally, the clusters were produced in an adiabatic-expansion setup using helium as the driving gas, spanning a range of 100–600 molecules in mean cluster size. The effective attenuation length was determined to be 8.4 ± 1.9 Å, in good agreement with an independent estimate of 10 Å formed on the basis of molecular electron-scattering data and Monte Carlo simulations. The aggregation state of the clusters as well as the cluster temperature and its importance to the derived EAL value are discussed in some depth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A spacecraft's own ambient environment: The role of simulation-based research

Spacecraft contamination has long been a subject of study in the rarefied gas dynamics community. Professor Mikhail Ivanov coined the term a spacecraft's 'own ambient environment' to describe the effects of natural and satellite driven processes on the conditions encountered by a spacecraft in orbit. Outgassing, thruster firings, and gas and liquid dumps all contribute to the spacecraft's contamination environment. Rarefied gas dynamic modeling techniques, such as Direct Simulation Monte Carlo, are well suited to investigate these spacebased environments. However, many advances were necessary to fully characterize the extent of this problem. A better understanding of modeling flows over large pressure ranges, for example hybrid continuum and rarefied numerical schemes, were required. Two-phase flow modeling under rarefied conditions was necessary. And the ability to model plasma flows for a new era of propulsion systems was also required. Through the work of Professor Ivanov and his team, we now have a better understanding of processes that create a spacecraft's own ambient environment and are able to better characterize these environments. Advances in numerical simulation have also spurred on the development of experimental facilities to study these effects. The relationship between numerical results and experimental advances will be explored in this manuscript.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modulation-mediated unlocking of a parametrically phase-locked spin torque oscillator

Modulation of an oscillator is crucial for its application in communication devices. While the output power and linewidth of single magnetic tunnel junction-based spin-torque oscillators (MTJ-STO) are not yet adequate for practical uses, the synchronization of such devices can overcome these limitations. Here, we investigate the modulation behavior of a parametrically synchronized MTJ-STO and show experimentally that modulation of the synchronized state preserves synchronization as long as the modulation frequency, f{sub mod}, is above a characteristic frequency, f{sub unlock}. The unlocking frequency increases with the modulation amplitude in agreement with analytical estimates and numerical simulations. These phenomena are described as a non-resonant unlocking mechanism, whose characteristics are directly related to inherent parameters of the oscillator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Research on anisotropy of fusion-produced protons and neutrons emission from high-current plasma-focus discharges

The paper concerns fast protons and neutrons from D-D fusion reactions in a Plasma-Focus-1000U facility. Measurements were performed with nuclear-track detectors arranged in “sandwiches” of an Al-foil and two PM-355 detectors separated by a polyethylene-plate. The Al-foil eliminated all primary deuterons, but was penetrable for fast fusion protons. The foil and first PM-355 detector were penetrable for fast neutrons, which were converted into recoil-protons in the polyethylene and recorded in the second PM-355 detector. The “sandwiches” were irradiated by discharges of comparable neutron-yields. Analyses of etched tracks and computer simulations of the fusion-products behavior in the detectors were performed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Combining gait optimization with passive system to increase the energy efficiency of a humanoid robot walking movement

There are several approaches to create the Humanoid robot gait planning. This problem presents a large number of unknown parameters that should be found to make the humanoid robot to walk. Optimization in simulation models can be used to find the gait based on several criteria such as energy minimization, acceleration, step length among the others. The energy consumption can also be reduced with elastic elements coupled to each joint. The presented paper addresses an optimization method, the Stretched Simulated Annealing, that runs in an accurate and stable simulation model to find the optimal gait combined with elastic elements. Final results demonstrate that optimization is a valid gait planning technique.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the influence of wall roughness in particle-laden flows

The distribution of inertial particles in turbulent flows is highly nonuniform and is governed by the local dynamics of the turbulent structures of the underlying carrier flow field. In wall-bounded flows, wall roughness strongly affects the turbulent flow field, nevertheless its effects on the particle transport in two-phase turbulent flows has been still poorly investigated. The issue is discussed here by addressing direct numerical simulations of a dilute dispersion of heavy particles in a turbulent channel flow, bounded by irregular two-dimensional rough surfaces, in the one-way coupling regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗