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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Common Envelope Light Curves. I. Grid-code Module Calibration

The common envelope (CE) binary interaction occurs when a star transfers mass onto a companion that cannot fully accrete it. The interaction can lead to a merger of the two objects or to a close binary. The CE interaction is the gateway of all evolved compact binaries, all stellar mergers, and likely many of the stellar transients witnessed to date. CE simulations are needed to understand this interaction and to interpret stars and binaries thought to be the byproduct of this stage. At this time, simulations are unable to reproduce the few observational data available and several ideas have been put forward to address their shortcomings. The need for more definitive simulation validation is pressing and is already being fulfilled by observations from time-domain surveys. In this article, we present an initial method and its implementation for post-processing grid-based CE simulations to produce the light curve so as to compare simulations with upcoming observations. Here we implemented a zeroth order method to calculate the light emitted from CE hydrodynamic simulations carried out with the 3D hydrodynamic code Enzo used in unigrid mode. The code implements an approach for the computation of luminosity in both optically thick and optically thin regimes and is tested using the first 135 days of the CE simulation of Passy et al., where a 0.8 M {sub ⊙} red giant branch star interacts with a 0.6 M {sub ⊙} companion. This code is used to highlight two large obstacles that need to be overcome before realistic light curves can be calculated. We explain the nature of these problems and the attempted solutions and approximations in full detail to enable the next step to be identified and implemented. We also discuss our simulation in relation to recent data of transients identified as CE interactions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

RF Cavity Characterization with VORPAL

When designing a radio frequency (RF) accelerating cavity structure various figures of merit are considered before coming to a final cavity design. These figures of merit include specific field and geometry based quantities such as the ratio of the shunt impedance to the quality factor (R/Q) or the normalized peak fields in the cavity. Other important measures of cavity performance include the peak surface fields as well as possible multipacting resonances in the cavity. High fidelity simulations of these structures can provide a good estimate of these important quantities before any cavity prototypes are built. We will present VORPAL simulations of a simple pillbox structure where these quantities can be calculated analytically and compare them to the results from the VORPAL simulations. We will then use VORPAL to calculate these figures of merit and potential multipacting resonances for two cavity designs under development at Jefferson National Lab for Project X.

43 PARTICLE ACCELERATORS↗
DOE OSTI2012

Agni: Coupling Model Analysis Tools and High-Performance Subsurface Flow and Transport Simulators for Risk and Performance Assessments

Agni is open source /community /multi platform framework for model-based analyses (SA, UQ, PE, RA, DS). Agni can be applied independently from other Advanced Simulation Capability for Environmental Management (ASCEM) modules. Agni will be available as Mercurial repository (hg clone https://akuna.labworks.org/hg/Platform) including source code, manual, test and verification examples.

COMPUTER CODES↗
DOE OSTI2016

Optimization design of helium circulator diffuser in HTR-PM - 18488

Circulator is very important in HTR-PM. It needs to supply enough pressure rise to circulate cooled helium in the primary loop of HTR-PM to bring heat out. Its pressure rise requirement is the highest among the worldwide circulators. Diffuser is one main component to produce pressure rise for circulator. And the most particular for HTR-PM is that the helium must flow out from two round tubes after flowing through the diffuser, which is very different from other circulators. This special export structure increases flow loss and decreases diffuser efficiency. So it is necessary to optimize diffuser to satisfy high pressure rise requirement of HTR-PM circulator with the special export structure. Four geometrical variables of diffuser have been studied and optimized. Based on the flow and performance analysis, the optimum structure was chosen. Then the detailed simulation was done by combining it with impeller to obtain performance characteristics of the whole circulator. Finally, the simulation results were compared with experimental results. It was obtained that simulation pressure rise was 201 kPa and experimental result was 210 kPa, which verified that optimization design of diffuser was effective. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

A Neutron Streaming Problem to Test Rattlesnake Methods for TREAT

This work demonstrates the progress made at Idaho National Laboratory in the neutronic modeling of neutron streaming in the TREAT reactor core. The streaming in the slotted elements is important because of its effect on the axial distribution and energy deposition in the experiment vehicle. Two methods implemented in Rattlesnake: 1) tensor diffusion coefficients and, 2) super-homogenization, yield much improved power distributions using a diffusion solution with homogenized fuel elements. The RMS difference in the power distribution decreases below 0.661% compared to a reference Monte Carlo simulation. The maximum differences are within 1.5% of this reference. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Experimental and CFD Analysis of Condensation with Non-Condensable Gas in Multi-Tubes

After the Fukushima accident, the importance of safety system that relies on passive cooling for decay heat removal has been reemphasized. The passive containment cooling system (PCCS) is the system for cooling down a containment and lowering a containment pressure even under the station blackout. In investigating a PCCS, steam condensation in the presence of non-condensable gas is the dominant heat and mass transfer phenomenon. The non-condensable gases degrade the heat removal performance through the PCCS heat exchangers (HXs) because they are accumulated near the cooling surface and deteriorate the approach of steam by diffusion to the surface. Therefore, in order to understand steam condensation phenomenon in the PCCS, we conducted multi-tube condensation experiment with air. And we also analyzed the results by CFD simulation for the CFD-based approach to condensation phenomena prediction. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Using the MCNP6.2 Correlated Fission Multiplicity Models, CGMF and FREYA

The culmination of a multi-year, collaborative effort by researchers at Los Alamos National Laboratory (LANL), Lawrence Berkeley National Laboratory (LBNL), Lawrence Livermore National Laboratory (LLNL) and the University of Michigan, led to the delivery of two new correlated fission multiplicity models available within the latest release of MCNP{sup R} version 6.2. The needs within the nuclear nonproliferation community, particularly in simulating and predicting the neutron and gamma-ray signatures of special nuclear materials (SNM), provided the motivation for this major collaborative effort. From LANL, the CGMF code has been implemented as an upgrade to the original Cascading Gamma-ray Multiplicity (CGM) code. Additionally, the LLNL Fission Library has been upgraded from the original implementation released in prior versions of MCNP6/X to include the LBNL/LLNL Fission Reaction Event Yield Algorithm (FREYA) code. Both new codes are capable of simulating fission events with correlated neutron and gamma-ray emissions. Due to the release of these new models within MCNP6.2, the purpose of this paper is to describe the differences between these new models and existing capabilities, demonstrate their use in the context of simulation inputs and outputs, and to discuss some of the ongoing work directed toward the next generation of these features.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Radiolytic Gas Bubbles' Model and its Application to AHR's Transient Neutronics Analysis

Aqueous Homogeneous Reactor (AHR) receives much attention recently because its advantages in medical isotopes production. Fissile material and moderator is mixed in the fuel solution, leads to big negative feedback and generation of radiolytic gas bubbles; the fluidity of fuel and movement of bubbles bring out big challenge of the neutronics analysis methods. This paper aims to establish a reliable and effective gas bubbles' model, firstly a coupled Neutronics-CFD method is introduced to recognize gas bubbles' behavior characteristics, then it is simplified into mathematical expression and applied to Transient Calculation Code of Aqueous Homogeneous Reactor (TCCAHR), finally a typical AHR's transient process is simulated and compared with the experiment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Effect of Low Temperature on a 4 W/60 K Pulse-Tube Cryocooler for Cooling HgCdTe Detector

Temperature is an extremely important parameter for the material of the space-borne infrared detector. To cool an HgCdTe-infrared detector, a Stirling-type pulse-tube cryocooler (PTC) has been developed based on a great deal of numerical simulations, which are performed to investigate the thermodynamic behaviors of the PTC. The effects of different low temperatures are presented to analyze different energy flows, losses, phase shifts, and impedance matching of the PTC at a temperature range of 40–120 K, where woven wire screens are used. Finally, a high-efficiency coaxial PTC has been designed, built, and tested, operating around 60 K after a number of theoretical and experimental studies. The PTC can offer a no-load refrigeration temperature of 40 K with an input electric power of 150 W, and a cooling power of 4 W at 60 K is obtained with Carnot efficiency of 12%. In addition, a comparative study of simulation and experiment has been carried out, and some studies on reject temperatures have been presented for a thorough understanding of the PTC system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Capturing Hidden Geochemical Anomalies in Scarce Data by Fractal Analysis and Stochastic Modeling

Fractal/multifractal modeling is a widely used geomathematical approach to capturing different populations in geochemical mapping. The rationale of this methodology is based on empirical frequency density functions attained from global or local distributions. This approach is quite popular because of its simplicity and versatility; it accounts for the frequency and spatial distribution of geochemical data considering self-similarity across a range of scales. Using this technique for detection of geochemical anomalies in scarce data, however, is problematic and can lead to systematic bias in the characterization of the underlying populations. In this paper, an innovative technique is presented that provides good results without a priori assumptions. A simulation approach is adopted for fractal analysis by generating different possible distribution scenarios for the variable under study to reveal the underlying populations that are frequently hidden due to lack of data. The proposed technique is called the global simulated size–number method, and it is validated in a case study with two synthetic datasets and another case study with real dataset from the Ushtagan gold deposit in northeast Kazakhstan.

58 GEOSCIENCES↗
DOE OSTI2017

Parameter Calibration of GTN Damage Model and Formability Analysis of 22MnB5 in Hot Forming Process

Hot forming of high strength steel at elevated temperatures is an attractive technology to achieve the lightweight of vehicle body. The mechanical behavior of boron steel 22MnB5 strongly depends on the variation of temperature which makes the process design more difficult. In this paper, the Gurson–Tvergaard–Needleman (GTN) model is used to study the formability of 22MnB5 sheet at different temperatures. Firstly, the rheological behavior of 22MnB5 is analyzed through a series of hot tensile tests at a temperature range of 600-800 °C. Then, a detailed process to calibrate the damage parameters is given based on the response surface methodology and genetic algorithm method. The GTN model together with the damage parameters calibrated is then implemented to simulate the deformation and damage evolution of 22MnB5 in the process of high-temperature Nakazima test. The capability of the GTN model as a suitable tool to evaluate the sheet formability is confirmed by comparing experimental and calculated results. Finally, as a practical application, the forming limit diagram of 22MnB5 at 700 °C is constructed using the Nakazima simulation and Marciniak–Kuczynski (M–K) model, respectively. And the simulation integrated GTN model shows a higher reliability by comparing the predicted results of these two approaches with the experimental ones.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Implementing phase-covariant cloning in circuit quantum electrodynamics

An efficient scheme is proposed to implement phase-covariant quantum cloning by using a superconducting transmon qubit coupled to a microwave cavity resonator in the strong dispersive limit of circuit quantum electrodynamics (QED). By solving the master equation numerically, we plot the Wigner function and Poisson distribution of the cavity mode after each operation in the cloning transformation sequence according to two logic circuits proposed. The visualizations of the quasi-probability distribution in phase-space for the cavity mode and the occupation probability distribution in the Fock basis enable us to penetrate the evolution process of cavity mode during the phase-covariant cloning (PCC) transformation. With the help of numerical simulation method, we find out that the present cloning machine is not the isotropic model because its output fidelity depends on the polar angle and the azimuthal angle of the initial input state on the Bloch sphere. The fidelity for the actual output clone of the present scheme is slightly smaller than one in the theoretical case. The simulation results are consistent with the theoretical ones. This further corroborates our scheme based on circuit QED can implement efficiently PCC transformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Semiclassical approach to dynamics of interacting fermions

Understanding the behaviour of interacting fermions is of fundamental interest in many fields ranging from condensed matter to high energy physics. Developing numerically efficient and accurate simulation methods is an indispensable part of this. Already in equilibrium, fermions are notoriously hard to handle due to the sign problem. Out of equilibrium, an important outstanding problem is the efficient numerical simulation of the dynamics of these systems. In this work we develop a new semiclassical phase-space approach (a.k.a. the truncated Wigner approximation) for simulating the dynamics of interacting fermions in arbitrary dimensions. As fermions are essentially non-classical objects, a phase-space is constructed out of all fermionic bilinears. Classical phase-space is thus comprised of highly non-local (hidden) variables representing these bilinears, and the cost of the method is that it scales quadratic rather than linear with system size. We demonstrate the strength of the method by comparing the results to the exact quantum dynamics of fermion expansion in the Hubbard model and quantum thermalization in the Sachdev–Ye–Kitaev (SYK) model for small systems, where the semiclassics nearly perfectly reproduces correct results. We furthermore analyse fermion expansion in a larger, intractable by exact methods, 2D Hubbard model, which is directly relevant to recent cold atom experiments. - Highlights: • We propose a semiclassical method to simulate dynamics of interacting fermions. • A phase-space is constructed out of all fermionic bilinears. • We study quantum thermalization in the Sachdev–Ye–Kitaev model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Noise propagation in hybrid models of nonlinear systems: The Ginzburg–Landau equation

Every physical phenomenon can be described by multiple models with varying degrees of fidelity. The computational cost of higher fidelity models (e.g., molecular dynamics simulations) is invariably higher than that of their lower fidelity counterparts (e.g., a continuum model based on differential equations). While the former might not be suitable for large-scale simulations, the latter are not universally valid. Hybrid algorithms provide a compromise between the computational efficiency of a coarse-scale model and the representational accuracy of a fine-scale description. This is achieved by conducting a fine-scale computation in subdomains where it is absolutely required (e.g., due to a local breakdown of a continuum model) and coupling it with a coarse-scale computation in the rest of a computational domain. We analyze the effects of random fluctuations generated by the fine-scale component of a nonlinear hybrid on the hybrid's overall accuracy and stability. Two variants of the time-dependent Ginzburg–Landau equation (GLE) and their discrete representations provided by a nearest-neighbor Ising model serve as a computational testbed. Our analysis shows that coupling these descriptions in a one-dimensional simulation leads to erroneous results. Adding a random source term to the GLE provides accurate prediction of the mean behavior of the quantity of interest (magnetization). It also allows the two GLE variants to correctly capture the strength of the microscale fluctuations. Our work demonstrates the importance of fine-scale noise in hybrid simulations, and suggests the need for replacing an otherwise deterministic coarse-scale component of the hybrid with its stochastic counterpart.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Generalized non-local surface susceptibility model and Fresnel coefficients for the characterization of periodic metafilms with bianisotropic scatterers

Highlights: •Formulas for E/M fields radiated by continuous surface polarization distributions. •Non-local effective surface susceptibility model for periodic metafilms. •Generalized reflection and transmission coefficients for an arbitrary metafilm. •Successful treatment of non-planar scatterer arrays and spatial dispersion effects. -- Abstract: A non-local surface susceptibility model for the consistent description of periodic metafilms formed by arbitrarily-shaped, electrically-small, bianisotropic scatterers is developed in this paper. The rigorous scheme is based on the point-dipole approximation technique and is valid for any polarization and propagation direction of an electromagnetic wave impinging upon the metafilm, unlike existing approaches whose applicability is practically confined to very specific cases of incidence. Next, the universal form of the resulting surface susceptibility matrix is employed for the derivation of the generalized Fresnel coefficients for such surfaces, which enable the comprehensive interpretation of several significant, yet relatively unexamined, physical interactions. Essentially, these coefficients include eight distinct terms, corresponding to the co-polarized and cross-polarized reflection and transmission coefficients for the two orthogonal eigenpolarizations of a linearly-polarized incident plane wave. The above formulas are, then, utilized for the prediction of the scattering properties of metafilms with different planar and non-planar resonators, which are characterized via the featured model and two previously reported local ones. Their comparison with numerical simulation outcomes substantiates the merits of the proposed method, reveals important aspects of the underlying physics, and highlights the differences between the various modeling procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Dynamic earthquake rupture simulations on nonplanar faults embedded in 3D geometrically complex, heterogeneous elastic solids

Dynamic propagation of shear ruptures on a frictional interface in an elastic solid is a useful idealization of natural earthquakes. The conditions relating discontinuities in particle velocities across fault zones and tractions acting on the fault are often expressed as nonlinear friction laws. The corresponding initial boundary value problems are both numerically and computationally challenging. In addition, seismic waves generated by earthquake ruptures must be propagated for many wavelengths away from the fault. Therefore, reliable and efficient numerical simulations require both provably stable and high order accurate numerical methods. We present a high order accurate finite difference method for: a) enforcing nonlinear friction laws, in a consistent and provably stable manner, suitable for efficient explicit time integration; b) dynamic propagation of earthquake ruptures along nonplanar faults; and c) accurate propagation of seismic waves in heterogeneous media with free surface topography. We solve the first order form of the 3D elastic wave equation on a boundary-conforming curvilinear mesh, in terms of particle velocities and stresses that are collocated in space and time, using summation-by-parts (SBP) finite difference operators in space. Boundary and interface conditions are imposed weakly using penalties. By deriving semi-discrete energy estimates analogous to the continuous energy estimates we prove numerical stability. The finite difference stencils used in this paper are sixth order accurate in the interior and third order accurate close to the boundaries. However, the method is applicable to any spatial operator with a diagonal norm satisfying the SBP property. Time stepping is performed with a 4th order accurate explicit low storage Runge–Kutta scheme, thus yielding a globally fourth order accurate method in both space and time. We show numerical simulations on band limited self-similar fractal faults revealing the complexity of rupture dynamics on rough faults.

58 GEOSCIENCES↗
DOE OSTI2017

Multi-fidelity stochastic collocation method for computation of statistical moments

We present an efficient numerical algorithm to approximate the statistical moments of stochastic problems, in the presence of models with different fidelities. The method extends the multi-fidelity approximation method developed in . By combining the efficiency of low-fidelity models and the accuracy of high-fidelity models, our method exhibits fast convergence with a limited number of high-fidelity simulations. We establish an error bound of the method and present several numerical examples to demonstrate the efficiency and applicability of the multi-fidelity algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Interrupted in situ EBSD study of texture evolution and mechanism of surface grains in electroformed Ni after annealing with an initially duplex <100> + <111> fiber texture during uniaxial tensile deformation

Highlights: • The texture evolution of the electroformed Ni was investigated by using an interrupted in situ electron backscattered diffraction. • The initial <100> fiber component evolved into the metastable cube component with a weaker copper component. • Only a few grains of the initial <111> fiber component evolved into the copper orientation. - Abstract: The texture evolution process and mechanism of surface grains in electroformed Ni after annealing with an initially duplex <100> + <111> fiber texture during uniaxial tensile deformation were investigated by using interrupted in situ electron backscattered diffraction. The results showed that the initially duplex <100> + <111> fiber texture evolved into the metastable Cube orientation with a weak Cu orientation during the uniaxial tensile deformation. The volume fractions of <100> and <111> fiber texture decreased from 54.0% to 29.3% and 16.6% to 6.6%, respectively. The volume of the Cube orientation remained almost stable at 12%, whereas that of the Cu orientation increased from 1% to 2.9%. A simulation of the final preferred orientation based on the Taylor model was a cubic texture, which qualitatively agreed well with the experimental result. The simulated results indicated that the formed Cube and Cu textures evolved from the initial <100> and <111> fiber textures, respectively.

36 MATERIALS SCIENCE↗