Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Numerical Simulation and Experimental Validation of Nondendritic Structure Formation in Magnesium Alloy Under Oscillation and Ultrasonic Vibration

In this study, the formation of nondendritic structures in the primary phase of magnesium alloy solidified under oscillation and ultrasonic vibration was investigated by numerical simulation and experimentally. The growth and motion of a dendrite during solidification was simulated by a combination of the lattice Boltzmann method and the phase-field method. The simulation and experimental results indicated that higher oscillation amplitudes and acoustic streaming made the microstructures change from dendritic to nondendritic in the α-Mg primary phases. A sufficient shear stress and an appropriate flow time in the barrel should be satisfied when a given inclined angle is selected. The effects of the flow and the thermal field on the nucleation behavior, the constitutional undercooling, and the growth morphology are also discussed. It was found that a high shear rate and high turbulence can help homogenize the temperature and the concentration fields and collide and rotate the α-Mg primary phases.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Numerical simulation of effective mechanical properties of stochastic composites with consideration for structural evolution under intensive dynamic loading

Mechanical behavior of stochastic metal-ceramic composites with the aluminum matrix under high-rate deformation at shock-wave loading is numerically simulated with consideration for structural evolution. Effective values of mechanical parameters of metal-ceramic composites Al B{sub 4}C, Al SiC, and Al Al{sub 2}O{sub 3} are evaluated depending on different concentration of ceramic inclusions.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Numerical simulation of shrouded propellers

A numerical model was developed for the evaluation of the performance characteristics of a shrouded propeller. Using this model, a computational study was carried out to investigate the feasibility of improving the aerodynamic performance of a propeller by encasing it in a shroud. The propeller blade was modeled by a segmented bound vortex positioned along the span of the blade at its quarter-chord-line. The shroud was modeled by a number of discrete vortex rings. Due to the mutual dependence of shroud and propeller vortex strengths and the propeller vortex wake an iterative scheme was employed. Three shroud configurations were considered: a cylindrical and two conical shrouds. The computed performance of the shrouded propeller was compared with that of a free propeller of identical propeller geometry. The numerical results indicated that the cylindrical shroud outperformed the conical shroud configurations for the cases considered. Furthermore, when compared to the free propeller performance, the cylindrical shroud showed a considerable performance enhancement over the free propeller. However, the improvements were found to decrease with an increase in the advance ratio and to virtually diminish at advance ratios of about 2.5.

Afjeh, Abdollah A.↗
NASA NTRS1991

Numerical simulations of microflare evolution in the solar transition region and corona

Several observers report transient ultraviolet brightenings, often referred to as microflares, in the solar atmosphere. In this paper, the results are presented of a series of one-dimensional numerical simulations examining possible relationships between microflares and the generation of dynamical chromospheric and transition region features. Low-energy and medium-energy microflares eject long-lived cool, dense gas plugs into the corona, with the gas plug traversing the loop apex in the medium energy case. In the case of high-energy microflares, the gas plug is rapidly heated to the temperature of the surrounding corona, and the results resemble the dynamics occurring in standard solar flare thick-target electron beam models.

Sterling, Alphonse C.↗
DOE OSTI2013

Numerical simulation of the main characteristics of a high-pressure DF – CO{sub 2} laser for amplification of picosecond laser pulses

The gain characteristics of the medium of a pulsed DF – CO{sub 2} laser in the ten-micron region at the working gas pressures from 1 to 2.5 atm, which were experimentally determined in [4], are numerically simulated using a scheme that includes the main chemical and relaxation processes. It is shown that the chosen scheme of processes makes it possible (i) to numerically describe the experimental data on the temporal behaviour of gains; (ii) to explain the reason for early degradation of gains (at the degree of D{sub 2} transformation of about 20%); (iii) from comparison of experimental and calculated temporal gain profiles, to determine the degree of photodissociation of F{sub 2} molecules, which is an important parameter determining the operation of pulsed chemical DF – CO{sub 2} lasers; and (iv) to predict the gain characteristics of working mixtures depending on their composition and pressure and on the initiation parameters. The predicted gains in the mixtures of the optimal composition at a pressure of 2.5 atm, a degree of dissociation of F{sub 2} molecules per flash ∼1%, and a flash duration at half-width of ∼3 μs are ∼7 m{sup -1}. (lasers)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Numerical simulation of optical vortex propagation and reflection by the methods of scalar diffraction theory

Using the equations of scalar diffraction theory we consider the formation of an optical vortex on a diffractive optical element. The algorithms are proposed for simulating the processes of propagation of spiral wavefronts in free space and their reflections from surfaces with different roughness parameters. The given approach is illustrated by the results of numerical simulations. (propagation of wave fronts)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Numerical simulation of adiabatons in electromagnetically induced transparency under quasi-resonance conditions

The evolution of adiabatons in electromagnetically induced transparency in the Λ scheme of degenerate quantum transitions J = 0 → J = 1 → J = 2 with Doppler broadening of spectral lines has been numerically simulated taking into account the effect of resonance detunings. It is shown that, in the case of linearly polarised fields, an increase in the probe-field resonance detuning (under exact-resonance conditions for the control radiation) leads to a transformation of electromagnetically induced transparency into electromagnetically induced absorption at certain stages. When the control-field resonance detuning is varied, the transparency of the medium for the probe (exactly resonant) radiation monotonically decreases with increasing detuning because of the rising role of single-photon absorption. In the case of circularly polarised control radiation and linearly polarised input probe field, a probe pulse propagating in the medium splits into two pulses with oppositely directed circular polarisations. An increase in the probe pulse resonance detuning (under exact-resonance conditions for the control radiation) leads primarily to an increase in the absorption by the medium of the probe pulse, the direction of circular polarisation for which coincides with the circular-polarisation direction for the control radiation. (nonlinear optical phenomena)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1995

Testing higher-order Lagrangian perturbation theory against numerical simulations. 2: Hierarchical models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann-Lemaitre cosmogonies is compared with numerical simulations. We study the dynamics of hierarchical models as a second step. In the first step we analyzed the performance of the Lagrangian schemes for pancake models, the difference being that in the latter models the initial power spectrum is truncated. This work probed the quasi-linear and weakly non-linear regimes. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non-linear regime. We smooth the initial data by using a variety of filter types and filter scales in order to determine the optimal performance of the analytical models, as has been done for the 'Zel'dovich-approximation' - hereafter TZA - in previous work. We find that for spectra with negative power-index the second-order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low-order statistics like the power-spectrum. However, in contrast to the results found for pancake models, where the higher-order schemes get worse than TZA at late non-linear stages and on small scales, we here find that the second-order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these results we expect that the second-order truncated Lagrangian model is especially useful for the modelling of standard dark matter models such as Hot-, Cold-, and Mixed-Dark-Matter.

Melott, A. L.↗
NASA NTRS2012

Numerical Simulation of Noise from Supersonic Jets Passing Through a Rigid Duct

The generation, propagation and radiation of sound from a perfectly expanded Mach 2.5 cold supersonic jet flowing through an enclosed rigid-walled duct with an upstream J-deflector have been numerically simulated with the aid of OVERFLOW Navier-Stokes CFD code. A one-equation turbulence model is considered. While the near-field sound sources are computed by the CFD code, the far-field sound is evaluated by Kirchhoff surface integral formulation. Predictions of the farfield directivity of the OASPL (Overall Sound Pressure Level) agree satisfactorily with the experimental data previously reported by the author. Calculations also suggest that there is significant entrainment of air into the duct, with the mass flow rate of entrained air being about three times the jet exit mass flow rate.

Kandula, Max↗
DOE OSTI2015

Contribution of the numerical simulations in the development of the French HLLW vitrification Processes - 15126

For many years, the CEA (Commissariat a l'energie Atomique et aux energies Alternatives) Marcoule France has developed various processes dedicated to radioactive waste confinement, especially vitrification processes for HLLW. For 10 years now, the numerical simulation has become an important tool for research and development projects held in the CEA-AREVA Joint Vitrification Laboratory (LCV). Induction heating, fluid mechanics and thermal simulations take part of all new R and D projects. The contributions of such simulations are, first, the enhancement of the working knowledge of existing process. Those data are very useful to define optimisation choices, for example upgrades made on the hot metallic melter used since the 90's at La Hague facility. Second, the simulations are, of course, also used at the conception stage of new processes as a tool allowing wide ranges parametric tests. This has been extensively used in the design of the cold crucible inductive melter (CCIM) commissioned in 2010 at La Hague plant. Finally, it is a powerful and relatively cheap tool for prospective studies for processes of the future. Whatever the purpose, the potential benefits are gains on the reliability, the output capacity and the life time. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1987

Contributions of numerical simulation data bases to the physics, modeling and measurement of turbulence

The use of simulation data bases for the examination of turbulent flows is an effective research tool. Studies of the structure of turbulence have been hampered by the limited number of probes and the impossibility of measuring all desired quantities. Also, flow visualization is confined to the observation of passive markers with limited field of view and contamination caused by time-history effects. Computer flow fields are a new resource for turbulence research, providing all the instantaneous flow variables in three-dimensional space. Simulation data bases also provide much-needed information for phenomenological turbulence modeling. Three dimensional velocity and pressure fields from direct simulations can be used to compute all the terms in the transport equations for the Reynolds stresses and the dissipation rate. However, only a few, geometrically simple flows have been computed by direct numerical simulation, and the inventory of simulation does not fully address the current modeling needs in complex turbulent flows. The availability of three-dimensional flow fields also poses challenges in developing new techniques for their analysis, techniques based on experimental methods, some of which are used here for the analysis of direct-simulation data bases in studies of the mechanics of turbulent flows.

Moin, Parviz↗
DOE OSTI2026

A unified neural-network framework for nucleon imaging from numerical simulations of QCD

Parton distributions encode the momentum-space structure and, in their generalizations, the spatial tomography of quarks and gluons inside hadrons, the building blocks of visible matter. We present a unified neural-network approach that learns these distributions directly from matrix elements calculated via numerical simulations of quantum chromodynamics (QCD) on the lattice by fitting two complementary inputs simultaneously: data matched to physical quantities via known momentum-space and coordinate-space formalisms. Utilizing data from both methods stabilizes the extraction and mitigates biases that can arise when either is used alone. We validate the method on controlled mock data and apply it to lattice-QCD matrix elements to extract parton distribution functions (PDFs). We show benefits of such an approach for determining the physical quantities. We further extend the framework to zero-skewness generalized parton distributions and demonstrate nucleon tomography within the same neural-network parameterization. Our results provide an adaptable and systematically improvable approach for extracting partonic distributions from Euclidean correlators. It can incorporate polarization, additional channels, and future experimental constraints from current and future facilities, such as the Electron-Ion Collider.

Hadronic Spectroscopy↗
NASA NTRS1987

Numerical simulation of unsteady viscous flows

Most unsteady viscous flows may be grouped into two categories, i.e., forced and self-sustained oscillations. Examples of forced oscillations occur in turbomachinery and in internal combustion engines while self-sustained oscillations prevail in vortex shedding, inlet buzz, and wing flutter. Numerical simulation of these phenomena was achieved due to the advancement of vector processor computers. Recent progress in the simulation of unsteady viscous flows is addressed.

Hankey, Wilbur L.↗
DOE OSTI2016

On a silicon-based photonic-crystal cavity for the near-IR region: Numerical simulation and formation technology

The production technology of a photonic-crystal cavity formed as a group of holes in a silicon strip waveguide by ion-beam etching is described. The parasitic effect associated with hole conicity which develops upon hole formation by the given technology is studied. Numerical simulation shows that the hole-conicity induced decrease in the cavity quality factor can be compensated with consideration for the hole volume. The influence of the waveguide thickness on the resonance wavelength and quality factor of the photonic-crystal cavity is analyzed.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.

Yang, Y. S.↗
DOE OSTI2016

Numerical simulations of turbulent molecular clouds regulated by radiation feedback forces. I. Star formation rate and efficiency

Radiation feedback from stellar clusters is expected to play a key role in setting the rate and efficiency of star formation in giant molecular clouds. To investigate how radiation forces influence realistic turbulent systems, we have conducted a series of numerical simulations employing the Hyperion radiation hydrodynamics solver, considering the regime that is optically thick to ultraviolet and optically thin to infrared radiation. Our model clouds cover initial surface densities between Σ{sub cl,0}∼10--300 M{sub ⊙} pc{sup −2}, with varying initial turbulence. We follow them through turbulent, self-gravitating collapse, star cluster formation, and cloud dispersal by stellar radiation. All our models display a log-normal distribution of gas surface density Σ; for an initial virial parameter α{sub vir,0}=2, the log-normal standard deviation is σ{sub lnΣ}=1--1.5 and the star formation rate coefficient ε{sub ff, ρ-bar }=0.3--0.5, both of which are sensitive to turbulence but not radiation feedback. The net star formation efficiency (SFE) ε{sub final} increases with Σ{sub cl,0} and decreases with α{sub vir,0}. We interpret these results via a simple conceptual framework, whereby steady star formation increases the radiation force, such that local gas patches at successively higher Σ become unbound. Based on this formalism (with fixed σ{sub lnΣ}), we provide an analytic upper bound on ε{sub final}, which is in good agreement with our numerical results. The final SFE depends on the distribution of Eddington ratios in the cloud and is strongly increased by the turbulent compression of gas.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Numerical simulation of boundary layers. Part 1: Weak formulation and numerical method

A numerical method designed to solve the time-dependent, three-dimensional, incompressible Navier-Stokes equations in boundary layers is presented. The fluid domain is the half-space over a flat plate, and periodic conditions are applied in the horizontal directions. The discretization is spectral. The basis functions are divergence-free and a weak formulation of the momentum equation is used, which eliminates the pressure term. An exponential mapping and Jacobi polynomials are used in the semi-infinite direction, with the irrotational component receiving special treatment. Issues related to the accuracy, stability and efficiency of the method are discussed. Very fast convergence is demonstrated on some model problems with smooth solutions. The method has also been shown to accurately resolve the fine scales of transitional and turbulent boundary layers.

Spalart, P. R.↗
DOE OSTI2023

Numerical Simulations of Flow-Induced Deflections in MITR LEU Fuel Plate Due to Channel Size Disparity

The hydromechanical stability of the fuel plates in parallel coolant channels of a Materials Testing Reactor (MTR) fuel element design is of great importance to the safety of research and test reactors. Previous analytical, experimental, and numerical efforts focused on parallel channels with the same or similar size; also, in the prior numerical simulations, the fuel plate was often assumed to be perfectly flat. This work presents the results of a fluid-structure interaction simulation performed to evaluate the flow-induced deflections of the fuel plates in the low-enriched uranium (LEU, <20 wt% 235 U) fuel element design for the conversion (from highly enriched uranium) of the Massachusetts Institute of Technology Reactor (MITR-II, also referred to as MITR). Various manufacturing and assembly tolerances of the MITR LEU elements are considered in the analysis, and the effects of channel size disparity, nonideal plate shape, and flow rate uncertainty are investigated. Results show that, for all cases analyzed, the deflection occurs toward the larger channel, and the change in any channel stripe remains small (less than 0.021 mm) compared to fabrication tolerances. In addition to simulation work, a hydraulic performance test of the MITR LEU fuel element is currently planned to support conversion to the use of LEU fuel.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗