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

Influence of Welding Strength Matching Coefficient and Cold Stretching on Welding Residual Stress in Austenitic Stainless Steel

Austenitic stainless steel is widely used in pressure vessels for the storage and transportation of liquid gases such as liquid nitrogen, liquid oxygen, and liquid hydrogen. Cryogenic pressure vessel manufacturing uses cold stretching technology, which relies heavily on welding joint performance, to construct lightweight and thin-walled vessels. Residual stress from welding is a primary factor in cases of austenitic stainless steel pressure vessel failure. In this paper, on the basis of Visual Environment 10.0 finite element simulation technology, the residual stress resulting from different welding strength matching coefficients (0.8, 1, 1.2, 1.4) for two S30408 plates welded with three-pass butt welds is calculated according to thermal elastoplastic theory. In addition, the stress field was calculated under a loading as high as 410 MPa and after the load was released. Path 1 was set to analyze stress along the welding line, and path 2 was set to analyze stress normal to the welding line. The welding strength matching coefficient strongly affected both the longitudinal residual stress (center of path 1) and the transverse residual stress (both ends of path 1) after the welding was completed. However, the coefficient had little effect on the longitudinal and transverse residual stress of path 2. Under the loading of 410 MPa, the longitudinal and transverse stress decreased and the stress distribution, with different welding strength matching coefficients, was less diverse. After the load was released, longitudinal and transverse stress distribution for both path 1 and path 2 decreased to a low level. Cold stretching could reduce the effect of residual stress to various degrees. Transverse strain along the stretching direction was also taken into consideration. The experimental results validated the reliability of the partial simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A novel partitioning method for block-structured adaptive meshes

We propose a novel partitioning method for block-structured adaptive meshes utilizing the meshless Lagrangian particle concept. With the observation that an optimum partitioning has high analogy to the relaxation of a multi-phase fluid to steady state, physically motivated model equations are developed to characterize the background mesh topology and are solved by multi-phase smoothed-particle hydrodynamics. In contrast to well established partitioning approaches, all optimization objectives are implicitly incorporated and achieved during the particle relaxation to stationary state. Distinct partitioning sub-domains are represented by colored particles and separated by a sharp interface with a surface tension model. In order to obtain the particle relaxation, special viscous and skin friction models, coupled with a tailored time integration algorithm are proposed. Numerical experiments show that the present method has several important properties: generation of approximately equal-sized partitions without dependence on the mesh-element type, optimized interface communication between distinct partitioning sub-domains, continuous domain decomposition which is physically localized and implicitly incremental. Therefore it is particularly suitable for load-balancing of high-performance CFD simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Numerical approach of the injection molding process of fiber-reinforced composite with considering fiber orientation

One of the most important challenges in the injection molding process of the short-glass fiber/thermoplastic composite parts is being able to predict the fiber orientation, since it controls the mechanical and the physical properties of the final parts. Folgar and Tucker included into the Jeffery equation a diffusive type of term, which introduces a phenomenological coefficient for modeling the randomizing effect of the mechanical interactions between the fibers, to predict the fiber orientation in concentrated suspensions. Their experiments indicated that this coefficient depends on the fiber volume fraction and aspect ratio. However, a definition of the fiber interaction coefficient, which is very necessary in the fiber orientation simulations, hasn't still been proven yet. Consequently, this study proposed a developed fiber interaction model that has been introduced a fiber dynamics simulation in order to obtain a global fiber interaction coefficient. This supposed that the coefficient is a sum function of the fiber concentration, aspect ratio, and angular velocity. The proposed model was incorporated into a computer aided engineering simulation package C-Mold. Short-glass fiber/polyamide-6 composites were produced in the injection molding with the fiber weight concentration of 30 wt.%, 50 wt.%, and 70 wt.%. The physical properties of these composites were examined, and their fiber orientation distributions were measured by micro-computed-tomography equipment μ-CT. The simulation results showed a good agreement with experiment results.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Design of 4-electrode optical device for application of vector electric fields to self-assembled quantum dot complexes

Self-assembled InAs quantum dots (QDs) are of great interest as components of optoelectronic devices that can operate at the quantum limit. The charge configuration, interdot coupling, and symmetry of complexes containing multiple QDs can all be tuned with applied electric fields, but the magnitude and angle of the electric field required to control each of these parameters depend on the orientation of the QD complex. We present a 4-electrode device compatible with optical excitation and emission that allows application of electric fields with arbitrary magnitudes and angles relative to isolated QD complexes. We demonstrate the electric field tunability of this device with numerical simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Folding of small knotted proteins: Insights from a mean field coarse-grained model

A small but relevant number of proteins whose native structure is known features nontrivial topology, i.e., they are knotted. Understanding the process of folding from a swollen unknotted state to the biologically relevant native conformation is, for these proteins, particularly difficult, due to their rate-limiting topological entanglement. To shed some light into this conundrum, we introduced a structure-based coarse-grained model of the protein, where the information about the folded conformation is encoded in bonded angular interactions only, which do not favor the formation of native contacts. A stochastic search scheme in parameter space is employed to identify a set of interactions that maximizes the probability to attain the knotted state. The optimal knotting pathways of the two smallest knotted proteins, obtained through this approach, are consistent with the results derived by means of coarse-grained as well as full atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

HELIOS-K: An ultrafast, open-source opacity calculator for radiative transfer

We present an ultrafast opacity calculator that we name HELIOS-K. It takes a line list as an input, computes the shape of each spectral line, and provides an option for grouping an enormous number of lines into a manageable number of bins. We implement a combination of Algorithm 916 and Gauss–Hermite quadrature to compute the Voigt profile, write the code in CUDA, and optimize the computation for graphics processing units (GPUs). We restate the theory of the k-distribution method and use it to reduce ∼10{sup 5}–10{sup 8} lines to ∼10–10{sup 4} wavenumber bins, which may then be used for radiative transfer, atmospheric retrieval and general circulation models. The choice of line-wing cutoff for the Voigt profile is a significant source of error and affects the value of the computed flux by ∼10%. This is an outstanding physical (rather than computational) problem, due to our incomplete knowledge of pressure broadening of spectral lines in the far line wings. We emphasize that this problem remains regardless of whether one performs line-by-line calculations or uses the k-distribution method and affects all calculations of exoplanetary atmospheres requiring the use of wavelength-dependent opacities. We elucidate the correlated-k approximation and demonstrate that it applies equally to inhomogeneous atmospheres with a single atomic/molecular species or homogeneous atmospheres with multiple species. Using a NVIDIA K20 GPU, HELIOS-K is capable of computing an opacity function with ∼10{sup 5} spectral lines in ∼1 s and is publicly available as part of the Exoclimes Simulation Platform (www.exoclime.org).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Testing ΛCDM cosmology at turnaround: where to look for violations of the bound?

In ΛCDM cosmology, structure formation is halted shortly after dark energy dominates the mass/energy budget of the Universe. A manifestation of this effect is that in such a cosmology the turnaround radius—the non-expanding mass shell furthest away from the center of a structure— has an upper bound. Recently, a new, local, test for the existence of dark energy in the form of a cosmological constant was proposed based on this turnaround bound. Before designing an experiment that, through high-precision determination of masses and —independently— turnaround radii, will challenge ΛCDM cosmology, we have to answer two important questions: first, when turnaround-scale structures are predicted to be close enough to their maximum size, so that a possible violation of the bound may be observable. Second, which is the best mass scale to target for possible violations of the bound. These are the questions we address in the present work. Using the Press-Schechter formalism, we find that turnaround structures have in practice already stopped forming, and consequently, the turnaround radius of structures must be very close to the maximum value today. We also find that the mass scale of ∼ 10{sup 13} M{sub ⊙} characterizes the turnaround structures that start to form in a statistically important number density today —and even at an infinite time in the future, since structure formation has almost stopped. This mass scale also separates turnaround structures with qualitatively different cosmological evolution: smaller structures are no longer readjusting their mass distribution inside the turnaround scale, they asymptotically approach their ultimate abundance from higher values, and they are common enough to have, at some epoch, experienced major mergers with structures of comparable mass; larger structures exhibit the opposite behavior. We call this mass scale the transitional mass scale and we argue that it is the optimal for the purpose outlined above. As a corollary, we explain the different accretion behavior of small and larger structures observed in already conducted numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

THE INNER DYNAMICAL MASS ACROSS GALAXY MORPHOLOGY: A WEAK SCALING WITH TOTAL STELLAR MASS

We compile a sample of dispersion-supported and rotation-supported galaxies for which the dynamical mass enclosed at 500 pc can be measured. We find that this dynamical quantity increases only slowly with stellar mass (M{sub total}(<500 pc)∝M{sub stars}{sup 0.16}) over ∼9 decades in baryonic mass and ∼3 in length scale, with a sudden upturn at the highest masses (M {sub stars} ≳ 10{sup 10} M {sub ☉}). This upturn occurs for the earliest types (S0 and Sa) of disk galaxies, and is consistent with the ''extra'' mass within 500 pc being predominantly stellar. This sudden change may be indicative of different bulge formation mechanisms between early- and late-type disks. We also compare the data to the results of recent galaxy formation simulations. These broadly track the observed trend, with some tendency to put too much mass in the inner 500 pc, depending on the details of each simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Application of Grover’s Quantum Algorithm for SDES Key Searching

The problem of finding the key of Simplified-DES (SDES)—a model of a block cipher DES—by Grover’s quantum algorithm is considered. Examples of application of Grover’s algorithm are presented. A quantum system with the minimum number of qubits is constructed that implements SDES key searching by a single pair of plaintext and ciphertext and requires only 19 qubits. This quantum circuit is simulated by using a Quipper quantum simulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Comparison of Hydrocode Simulations with Measured Shock Wave Velocities

We have conducted detailed 1- and 2-dimensional hydrodynamics calculations to assess the quality of simulations commonly made to understand various shock processes in a sample and to design shock experiments. We began with relatively simple shock experiments, where we examined the effects of the equation of state and the viscoplastic strength models. Eventually we included spallation in copper and iron and a solid-solid phase transformation in iron to assess the quality of the damage and phase transformation simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2016

X-ray simulation algorithms used in ISP

ISP is a simulation code which is sometimes used in the USNDS program. ISP is maintained by Sandia National Lab. However, the X-ray simulation algorithm used by ISP was written by scientists at LANL – mainly by Ed Fenimore with some contributions from John Sullivan and George Neuschaefer and probably others. In email to John Sullivan on July 25, 2016, Jill Rivera, ISP project lead, said “ISP uses the function xdosemeters_sim from the xgen library.” The is a fortran subroutine which is also used to simulate the X-ray response in consim (a descendant of xgen). Therefore, no separate documentation of the X-ray simulation algorithms in ISP have been written – the documentation for the consim simulation can be used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

SNe Ia AND THEIR ENVIRONMENT: THEORY AND APPLICATIONS TO SN 2014J

We present theoretical semi-analytic models for the interaction of stellar winds with the interstellar medium (ISM) or prior mass loss implemented in our code SPICE, assuming spherical symmetry and power-law ambient density profiles and using the Π-theorem. This allows us to test a wide variety of configurations, their functional dependencies, and to find classes of solutions for given observations. Here, we study Type Ia Supernova (SN Ia) surroundings of single and double degenerate systems, and their observational signatures. Winds may originate from the progenitor prior to the white dwarf (WD) stage, the WD, a donor star, or an accretion disk (AD). For M{sub Ch} explosions, the AD wind dominates and produces a low-density void several light years across, surrounded by a dense shell. The bubble explains the lack of observed interaction in late time SN light curves for, at least, several years. The shell produces narrow ISM lines Doppler shifted by 10–100 km s{sup −1}, and equivalent widths of ≈100 mÅ and ≈1 mÅ in cases of ambient environments with constant density and produced by prior mass loss, respectively. For SN2014J, both mergers and M{sub Ch} mass explosions have been suggested based on radio and narrow lines. As a consistent and most likely solution, we find an AD wind running into an environment produced by the red giant wind of the progenitor during the pre-WD stage, and a short delay, 0.013–1.4 Myr, between the WD formation and the explosion. Our framework may be applied more generally to stellar winds and star formation feedback in large scale galactic evolution simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A BARYONIC EFFECT ON THE MERGER TIMESCALE OF GALAXY CLUSTERS

Accurate estimation of the merger timescales of galaxy clusters is important for understanding the cluster merger process and further understanding the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merger process. The merger timescale decreases upon increasing the gas fraction of galaxy clusters. For example, the merger timescale is shortened by a factor of up to 3 for merging clusters with gas fractions of 0.15, compared with the timescale obtained with 0 gas fractions. The baryonic effect is significant for a wide range of merger parameters and is particularly more significant for nearly head-on mergers and high merging velocities. The baryonic effect on the merger timescale of galaxy clusters is expected to have an impact on the structure formation in the universe, such as the cluster mass function and massive substructures in galaxy clusters, and a bias of “no-gas” may exist in the results obtained from the dark matter-only cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Time delay analysis of the lensed quasar SDSS J1001+5027

We modify the algorithm we proposed in Aghamousa and Shafieloo for the time delay estimation of strongly lensed systems incorporating the weighted cross-correlation and weighted summation of correlation coefficients. We show the high performance of this algorithm by applying it to Time Delay Challenge (TDC1) simulated data. We apply then our proposed method to the light curves of the lensed quasar SDSS J1001+5027 since this system has been well studied by other groups, to compare our results with their findings. In this work we propose a new estimator, the “mirror” estimator, along with a list of criteria for reliability testing of the estimation. Our mirror estimator results are −117.1{sub −3.7}{sup +7.1} and −117.1{sub −8.8}{sup +7.2} using simple Monte Carlo simulations and simulated light curves provided by Rathna Kumar et al., respectively. Although the TDC1 simulations do not reflect the properties of the SDSS J1001+5027 light curves, using these simulations results in a smaller uncertainty, which shows that the higher quality observations can lead to a substantially more precise time delay estimation. Our time delay estimation is in agreement with the findings of the other groups for this strongly lensed system, and the difference in the size of the error bars reflects the importance of appropriate light curve simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Tuning linear response dynamics near the Dirac points in the bosonic Mott insulator

Optical lattice systems offer the possibility of creating and tuning Dirac points which are present in the tight-binding lattice dispersions. For example, such a behavior can be achieved in the staggered flux lattice or honeycomb type of lattices. Here we focus on the strongly correlated bosonic dynamics in the vicinity of Dirac points. In particular, we investigate bosonic Mott insulator phase in which quasiparticle excitations have a simple particle–hole interpretation. We show that linear response dynamics around Dirac points, can be significantly engineered at least in two ways: by the type of external perturbation or by changing the lattice properties. The key role is played by the interband transitions. Moreover, we explain that the behavior of these transitions is directly connected to different energy scales of the effective hopping amplitudes for particles and holes. Presented in this work theoretical study about tunability of linear response dynamics near the Dirac points can be directly simulated in the optical lattice systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A generalized transport-velocity formulation for smoothed particle hydrodynamics

The standard smoothed particle hydrodynamics (SPH) method suffers from tensile instability. In fluid-dynamics simulations this instability leads to particle clumping and void regions when negative pressure occurs. In solid-dynamics simulations, it results in unphysical structure fragmentation. In this work the transport-velocity formulation of Adami et al. (2013) is generalized for providing a solution of this long-standing problem. Other than imposing a global background pressure, a variable background pressure is used to modify the particle transport velocity and eliminate the tensile instability completely. Furthermore, such a modification is localized by defining a shortened smoothing length. The generalized formulation is suitable for fluid and solid materials with and without free surfaces. The results of extensive numerical tests on both fluid and solid dynamics problems indicate that the new method provides a unified approach for multi-physics SPH simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Nonlinear physics of shear Alfvén waves

Shear Alfvén waves (SAW) play fundamental roles in thermonuclear plasmas of fusion interest, since they are readily excited by energetic particles in the MeV range as well as by the thermal plasma components. Thus, understanding fluctuation induced transport in burning plasmas requires understanding nonlinear SAW physics. There exist two possible routes to nonlinear SAW physics: (i) wave-wave interactions and the resultant spectral energy transfer; (ii) nonlinear wave-particle interactions of SAW instabilities with energetic particles. Within the first route, it is advantageous to understand and describe nonlinear processes in term of proximity of the system to the Alfvénic state, where wave-wave interactions are minimized due to the cancellation of Reynolds and Maxwell stresses. Here, various wave-wave nonlinear dynamics are elucidated in terms of how they break the Alfvénic state. In particular, we discuss the qualitative and quantitative modification of the SAW parametric decay process due to finite ion compressibility and finite ion Larmor radius. We also show that toroidal geometry plays a crucial role in the nonlinear excitation of zonal structures by Alfvén eigenmodes. Within the second route, the coherent nonlinear dynamics of structures in the energetic particle phase space, by which secular resonant particle transport can occur on meso- and macro-scales, must be addressed and understood. These 'nonlinear equilibria' or 'phase-space zonal structures' dynamically evolve on characteristic (fluctuation induced) turbulent transport time scales, which are generally of the same order of the nonlinear time scale of the underlying fluctuations. In this work, we introduce the general structure of nonlinear Schrödinger equations with complex integro-differential nonlinear terms, which govern these physical processes. To elucidate all these aspects, theoretical analyses are presented together with numerical simulation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗