Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 955 records · Page 53

Photoluminescence spectra of point defects in semiconductors: Validation of first-principles calculations

Optically and magnetically active point defects in semiconductors are interesting platforms for the development of solid state quantum technologies. Their optical properties are usually probed by measuring photoluminescence spectra, which provide information on excitation energies and on the interaction of electrons with lattice vibrations. We present a combined computational and experimental study of photoluminescence spectra of defects in diamond and SiC, aimed at assessing the validity of theoretical and numerical approximations used in first-principles calculations, including the use of the Franck-Condon principle and the displaced harmonic oscillator approximation. We focus on prototypical examples of solid state qubits, the divacancy centers in SiC and the nitrogen-vacancy in diamond, and we report computed photoluminescence spectra as a function of temperature that are in very good agreement with the measured ones. As expected we find that the use of hybrid functionals leads to more accurate results than semilocal functionals. Interestingly our calculations show that constrained density functional theory (CDFT) and time-dependent hybrid DFT perform equally well in describing the excited state potential energy surface of triplet states; our findings indicate that CDFT, a relatively cheap computational approach, is sufficiently accurate for the calculations of photoluminescence spectra of the defects studied here. Finally, we find that only by correcting for finite-size effects and extrapolating to the dilute limit can one obtain a good agreement between theory and experiment. Our results provide a detailed validation protocol of first-principles calculations of photoluminescence spectra, necessary both for the interpretation of experiments and for robust predictions of the electronic properties of point defects in semiconductors.

36 MATERIALS SCIENCE↗

Massively parallel computation of RCS with finite elements

One of the promising combinations of finite element approaches for scattering problems uses Whitney edge elements, spherical vector wave-absorbing boundary conditions, and bi-conjugate gradient solution for the frequency-domain near field. Each of these approaches may be criticized. Low-order elements require high mesh density, but also result in fast, reliable iterative convergence. Spherical wave-absorbing boundary conditions require additional space to be meshed beyond the most minimal near-space region, but result in fully sparse, symmetric matrices which keep storage and solution times low. Iterative solution is somewhat unpredictable and unfriendly to multiple right-hand sides, yet we find it to be uniformly fast on large problems to date, given the other two approaches. Implementation of these approaches on a distributed memory, message passing machine yields huge dividends, as full scalability to the largest machines appears assured and iterative solution times are well-behaved for large problems. We present times and solutions for computed RCS for a conducting cube and composite permeability/conducting sphere on the Intel ipsc860 with up to 16 processors solving over 200,000 unknowns. We estimate problems of approximately 10 million unknowns, encompassing 1000 cubic wavelengths, may be attempted on a currently available 512 processor machine, but would be exceedingly tedious to prepare. The most severe bottlenecks are due to the slow rate of mesh generation on non-parallel machines and the large transfer time from such a machine to the parallel processor. One solution, in progress, is to create and then distribute a coarse mesh among the processors, followed by systematic refinement within each processor. Elimination of redundant node definitions at the mesh-partition surfaces, snap-to-surface post processing of the resulting mesh for good modelling of curved surfaces, and load-balancing redistribution of new elements after the refinement are auxiliary steps expected to result in a robust low i/o system for very large finite element problems.

Parker, Jay↗

Tuning nonequilibrium phase transitions with inertia

In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dispersion relation and instability for an anisotropic nonuniform flowing plasma

A generalized formula for wave instability is developed for an anisotropic nonuniform plasma with finite flows and temperatures. Six-moment fluid equations are solved to give the analytic expression for wave instability in arbitrarily nonuniform plasmas. The analytic formula explicitly states the dependence of wave instability on the nonuniformities of number density, flow velocity, and anisotropic or isotropic pressure. The accuracy of the formalism is verified by a numerical calculation of implicit dispersion relations in complex Fourier space. The analysis shows that nonuniformity plays a critical role in plasma instability, while the flow velocity and anisotropic pressures determine the growth rate of the instability. Lastly, the instability diagram and associated instability criterion for anisotropy-driven instability are introduced as applications of the formalism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating profile stiffness and critical gradients in shaped TCV discharges using local gyrokinetic simulations of turbulent transport

The experimental observation made on the TCV tokamak of a significant confinement improvement in plasmas with negative triangularity ( δ < 0) compared to those with standard positive triangularity has been interpreted in terms of different degrees of profile stiffness (Sauter et al. Phys. Plasmas 21 055906, 2014) and/or different critical gradients. Employing the Eulerian gyrokinetic code GENE (Jenko et al. Phys. Plasmas 7 1904, 2000), profile stiffness and critical gradients are studied under TCV relevant conditions. For the considered experimental discharges, trapped electron modes (TEMs) and electron temperature gradient (ETG) modes are the dominant microinstabilities, with the latter providing a significant contribution to the non-linear electron heat fluxes near the plasma edge. Two series of simulations with different levels of realism are performed, addressing the question of profile stiffness at various radial locations. Retaining finite collisionality, impurities and electromagnetic effects, as well as the physical electron-to-ion mass ratio are all necessary in order to approach the experimental flux measurements. However, flux-tube simulations are unable to fully reproduce the TCV results, pointing towards the need to carry out radially nonlocal (global) simulations, i.e. retaining finite machine size effects, in a future study. Some conclusions about the effect of triangularity can nevertheless be drawn based on the flux-tube results. In particular, the importance of considering the sensitivity to both temperature and density gradient is shown. The flux tube results show an increase of the critical gradients towards the edge, further enhanced when δ < 0, and they also appear to indicate a reduction of profile stiffness towards plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A comparison of two types of structural optimization procedures for satisfying flutter requirements

A comparison is made of results obtained from mathematical programming and optimality criteria procedures for the minimum-mass design of typical aircraft wing structures to satisfy prescribed flutter requirements. The mathematical programming method is based on an interior penalty function approach. A rigorous optimality criterion and an intuitive optimality criterion based on uniform strain energy density are considered. An intuitive resizing procedure is used for both optimality criterion solutions. All results are calculated using the same computer program, changing only the optimization procedure. Both high- and low-aspect-ratio wings are examined. Finite elements are used for structural modeling, and the generalized coordinates for the flutter solution are based on the natural vibration modes of the structure.

Haftka, R. T.↗

Enhanced erythrocyte suspension layer stability achieved by surface tension lowering additives

In connection with a fractionation procedure involving the separation of particles, a dilute suspension of these particles in a liquid is carefully layered on a dense liquid. Under ideal conditions, the suspension forms a zone of finite thickness with a 'sharp' interface between the suspension layer and the supporting liquid. Under an applied field, e.g., gravitational or electrical, the particles in the suspension layer migrate to form different layers according to their size and/or density or according to their electrophoretic mobilities. However, in many cases the ideal conditions necessary for the fractionation process are not obtained. Many studies have been conducted to explore the reasons for suspension layer 'instability'. The present investigation represents an extension of a study conducted by Omenyi et al. (1981). An electrostatic repulsion-van der Waals mechanism was used to study the stability of fixed erythrocyte suspensions layered on a D2O cushion.

Omenyi, S. N.↗

Thick-target bremsstrahlung interpretation of short time-scale solar hard X-ray features

Steady-state analyses of bremsstrahlung hard X-ray production in solar flares are appropriate only if the lifetime of the high energy electrons in the X-ray source is much shorter than the duration of the observed X-ray burst. For a thick-target nonthermal model, this implies that a full time-dependent analysis is required when the duration of the burst is comparable to the collisional lifetime of the injected electrons, in turn set by the lengths and densities of the flaring region. In this paper we present the results of such a time-dependent analysis, and we point out that the intrinsic temporal signature of the thick-target production mechanism, caused by the finite travel time of the electrons through the target, may indeed rule out such a mechanism for extremely short duration hard X-ray events.

Emslie, A. G.↗

Electrostatic forces in planetary rings

The average charge on a particle in a particle-plasma cloud, the plasma potential inside the cloud, and the Coulomb force acting on the particle are calculated. The net repulsive electrostatic force on a particle depends on the plasma density, temperature, density of particles, particle size, and the gradient of the particle density. In a uniformly dense ring the electrostatic repulsion is zero. It is also shown that the electrostatic force acts like a pressure force, that even a collisionless ring can be stable against gravitational collapse, and that a finite ring thickness does not necessarily imply a finite velocity dispersion. A simple criterion for the importance of electrostatic forces in planetary rings is derived which involves the calculation of the vertical ring thickness which would result if only electrostatic repulsion were responsible for the finite ring thickness. Electrostatic forces are entirely negligible in the main rings of Saturn and the E and G rings. They may also be negligible in the F ring. However, the Uranian rings and Jupiter's ring seem to be very much influenced by electrostatic repulsion. In fact, electrostatic forces could support a Jovian ring which is an order of magnitude more dense than observed.

Goertz, C. K.↗

Coupling of magnetopause-boundary layer to the polar ionosphere

The plasma dynamics in the low-latitude boundary layer and its coupling to the polar ionosphere under boundary conditions at the magnetopause are investigated. In the presence of a driven plasma flow along the magnetopause, the Kelvin-Helmholtz instability can develop, leading to the formation and growth of plasma vortices in the boundary layer. The finite ionospheric conductivity leads to the decay of these vortices. The competing effect of the formation and decay of vortices leads to the formation of strong vortices only in a limited region. Several enhanced field-aligned power density regions associated with the boundary layer vortices and the upward field-aligned current (FAC) filaments can be found along the postnoon auroral oval. These enhanced field-aligned power density regions may account for the observed auroral bright spots.

Wei, C. Q.↗

Adaptive Mesh Refinement for Microelectronic Device Design

Finite element and finite volume methods are used in a variety of design simulations when it is necessary to compute fields throughout regions that contain varying materials or geometry. Convergence of the simulation can be assessed by uniformly increasing the mesh density until an observable quantity stabilizes. Depending on the electrical size of the problem, uniform refinement of the mesh may be computationally infeasible due to memory limitations. Similarly, depending on the geometric complexity of the object being modeled, uniform refinement can be inefficient since regions that do not need refinement add to the computational expense. In either case, convergence to the correct (measured) solution is not guaranteed. Adaptive mesh refinement methods attempt to selectively refine the region of the mesh that is estimated to contain proportionally higher solution errors. The refinement may be obtained by decreasing the element size (h-refinement), by increasing the order of the element (p-refinement) or by a combination of the two (h-p refinement). A successful adaptive strategy refines the mesh to produce an accurate solution measured against the correct fields without undue computational expense. This is accomplished by the use of a) reliable a posteriori error estimates, b) hierarchal elements, and c) automatic adaptive mesh generation. Adaptive methods are also useful when problems with multi-scale field variations are encountered. These occur in active electronic devices that have thin doped layers and also when mixed physics is used in the calculation. The mesh needs to be fine at and near the thin layer to capture rapid field or charge variations, but can coarsen away from these layers where field variations smoothen and charge densities are uniform. This poster will present an adaptive mesh refinement package that runs on parallel computers and is applied to specific microelectronic device simulations. Passive sensors that operate in the infrared portion of the spectrum as well as active device simulations that model charge transport and Maxwell's equations will be presented.

Cwik, Tom↗

Electronic Properties of Carbon Nanotubes and Junctions

Metallic and semiconducting Single Wall Carbon Nanotubes (CNT) have recently been characterized using scanning tunneling microscopy (STM) and the manipulation of individual CNT has been demonstrated. These developments make the prospect of using CNT as molecular wires and possibly as electronic devices an even more interesting one. We have been modeling various electronic properties such as the density of states and the transmission coefficient of CNT wires and junctions. These studies involve first calculating the stability of junctions using molecular dynamics simulations and then calculating the electronic properties using a pi-electron tight binding Hamiltonian. We have developed the expertise to calculate the electronic properties of both finite-sized CNT and CNT systems with semi-infinite boundary conditions. In this poster, we will present an overview of some of our results. The electronic application of CNT that is most promising at this time is their use as molecular wires. The conductance can however be greatly reduced because of reflection due to defects and contacts. We have modeled the transmission through CNT in the presence of two types of defects: weak uniform disorder and strong isolated scatterers. We find that the conductance is affected in significantly different manners due to these defects Junctions of CNT have also been imaged using STM. This makes it essential to derive rules for the formation of junctions between tubes of different chirality, study their relative energies and electronic properties. We have generalized the rules for connecting two different CNT and have calculated the transmission and density of states through CNT junctions. Metallic and semiconducting CNT can be joined to form a stable junction and their current versus voltage characteristics are asymmetric. CNT are deformed by the application of external forces including interactions with a substrate or other CNT. In many experiments, these deformation are expected to occur naturally. We will present some preliminary results of our calculations of the modification of CNT electronic properties as a result of deformations.

Anantram, M. P.↗

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

Mechanistic verification of empirical UO 2 fuel fracture models

Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Alloying effects on deformation induced microstructure evolution in copper

In this work, we investigated the effects of alloying elements on plastic deformation and microstructure evolution in polycrystalline copper (Cu) and Cu alloyed with 1 wt. % lead (Cu-1 %Pb). These materials were selected due to the size mismatch between Cu and Pb, with the latter forming precipitates at grain boundaries. Multi-modal characterization techniques, including neutron diffraction, electron backscatter diffraction (EBSD), and transmission electron microscopy (TEM), along with finite element simulations were employed to study the deformation behavior across multiple length scales. While both Cu and Cu-1%Pb exhibited similar macroscale response and final deformation textures, both dislocation line profile analysis and TEM revealed increased dislocation density in deformed Cu-1%Pb specimens. The presence of lead precipitates also significantly affected local plastic deformation during compression, with their influence diminishing with increasing strain. These results demonstrate the complex relationships between alloying elements, plastic deformation, microstructural evolution, and material behavior under load. The insights gained from this multi-scale and multi-technique approach contribute to the fundamental understanding of microstructural evolution in immiscible alloys and are valuable for tailoring the properties of structural materials for specific engineering applications.

36 MATERIALS SCIENCE↗

Can local white dwarfs constrain dark matter interactions?

We investigate whether nearby white dwarfs (WDs) can constrain dark matter (DM) interactions with ordinary matter. As experimental sensitivity improves, driven by the Gaia mission, the sample volume of nearby WDs has been increasing over recent years. We carefully select a sample of ten cold, isolated, non-magnetic WDs within 13 pc of the Sun. We model their carbon-oxygen interior using a finite temperature relativistic equation of state and model atmospheres to infer their core temperatures. This enables us to perform a thorough estimation of the DM capture rate and evaporation mass using actual astrophysical observations. Given the low local DM density, we focus on DM that annihilates into long-lived mediators, which escape the WD and later decay into photons. While Fermi-LAT data shows no significant gamma-ray excess, future telescopes, CTA North & South, LHAASO, SWGO, could probe DM-nucleon cross sections down to ∼ 10 -41 cm 2 for DM masses above the TeV scale. Our results are competitive with current direct detection bounds (e.g., LZ) in the multi-TeV regime. This work underscores the importance of systematic WD studies in the broader landscape of DM detection and demonstrates the synergy between astrophysical and terrestrial searches in exploring DM interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

Structural design aspects of reusable surface insulation thermal protection systems.

Low density fiber ceramic materials coated with refractory ceramics meet the requirements of reusable low weight thermal protection systems. The structural characteristics of this class of material impose unique design and analysis requirements on the application to spacecraft structural elements. Finite element type stress analysis techniques are required to adequately predict the structural response of the system. Parametric analyses have been performed to determine the response of the system to variations in geometry, and to thermal and structural load conditions. Sensitivity of coating, insulation and attachment stresses are presented and critical failure modes are identified.

Michalak, R. J.↗