Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weighted relaxation”

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 55 records · Page 3

A Coulomb collision algorithm for weighted particle simulations

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Miller, Ronald H.↗

On relaxations of the max k -cut problem formulations

Here, a tight continuous relaxation is a crucial factor in solving mixed integer formulations of many NP-hard combinatorial optimization problems. The (weighted) max k-cut problem is a fundamental combinatorial optimization problem with multiple notorious mixed integer optimization formulations. In this paper, we explore four existing mixed integer optimization formulations of the max k-cut problem. Specifically, we show that the continuous relaxation of a binary quadratic optimization formulation of the problem is: (i) stronger than the continuous relaxation of two mixed integer linear optimization formulations and (ii) at least as strong as the continuous relaxation of a mixed integer semidefinite optimization formulation. We also conduct a set of experiments on multiple sets of instances of the max k-cut problem using state-of-the-art solvers that empirically confirm the theoretical results in item (i). Furthermore, these numerical results illustrate the advances in the efficiency of global non-convex quadratic optimization solvers and more general mixed integer nonlinear optimization solvers. As a result, these solvers provide a promising option to solve combinatorial optimization problems. Our codes and data are available on GitHub.

97 MATHEMATICS AND COMPUTING↗

ICDARTS: Improving the Stability of Cyclic DARTS

Cyclic DARTS (CDARTS) is a Differentiable Architecture Search (DARTS)-based approach to neural architecture search (NAS) that uses a cyclic feedback mechanism to train search and evaluation networks concurrently. This training protocol aims to optimize the search process and evaluate the deep evaluation network comprised of discretized candidate operations. However, this approach introduces a loss function for the evaluation network dependent on the search network. The dissimilarity between the evaluation network’s loss function used during the search and retraining phases results in a search network that is a sub-optimal proxy for the final evaluation network accessed during retraining. We present a revised approach that removes the dependency of the evaluation network weights upon those of the search network. In addition, we introduce a modified process for relaxing the search network’s zero operations that allows these operations to be retained in the final evaluation networks.

Herron, Emily↗

Warping of unsymmetric cross-ply graphite/epoxy laminates

Warping in unsymmetric graphite/epoxy laminates was studied with particular attention given to the change of residual stresses resulting from long term environmental exposure. Square, cured prepreg sheets were measured for edge deflection with a cathetometer, then quartered and remeasured. Two postcuring durations were then used, 7.5 and one hr at 177 C; varying cooldown rates after curing were used for other samples, and one set was stored in vacuum at 75 C. Maximum deflections and weight changes were measured periodically at room temperature. Average curvatures, the effect of postcure, and the effect of long-term exposure were determined. Larger panels exhibited cylindrical warping and smaller panels underwent anticlastic warping. The deflections were related to weight changes, i.e. moisture absorption, and the lower the moisture content, the higher the deflection. Relaxation of residual stresses at 75 C was neglibible after 220 days.

Hahn, H. T.↗

Adaptive array antenna for satellite cellular and direct broadcast communications

Adaptive phased-array antennas provide cost-effective implementation of large, light weight apertures with high directivity and precise beamshape control. Adaptive self-calibration allows for relaxation of all mechanical tolerances across the aperture and electrical component tolerances, providing high performance with a low-cost, lightweight array, even in the presence of large physical distortions. Beam-shape is programmable and adaptable to changes in technical and operational requirements. Adaptive digital beam-forming eliminates uplink contention by allowing a single electronically steerable antenna to service a large number of receivers with beams which adaptively focus on one source while eliminating interference from others. A large, adaptively calibrated and fully programmable aperture can also provide precise beam shape control for power-efficient direct broadcast from space. Advanced adaptive digital beamforming technologies are described for: (1) electronic compensation of aperture distortion, (2) multiple receiver adaptive space-time processing, and (3) downlink beam-shape control. Cost considerations for space-based array applications are also discussed.

Horton, Charles R.↗

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

In this study, a method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob–Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable and Memory-Efficient Algorithms for Controlling Networked Epidemic Processes Using Multiplicative Weights Update Method

We study the problem of designing scalable algorithms to find effective intervention strategies for controlling stochastic epidemic processes on networks. This is a common problem arising in agent based models for epidemic spread. Previous approaches to this problem focus on either heuristics with no guarantees or approximation algorithms that scale only to networks corresponding to county-sized populations, typically, with less than a million nodes. In particular, the mathematical-programming based approaches need to solve the Linear Program (LP) relaxation of the problem using an LP solver, which restricts the scalability of this approach. In this work, we overcome this restriction by designing an algorithm that adapts the multiplicative weights update (MWU) framework, along with the sample average approximation (SAA) technique, to approximately solve the linear program (LP) relaxation for the problem. To scale this approach further, we provide a memory-efficient algorithm that enables scaling to large networks, corresponding to country-size populations, with over 300 million nodes and 30 billion edges. Furthermore, we show that this approach provides near-optimal solutions to the LP in practice.

Sambaturu, Prathyush↗

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ultrafast bipolar conductivity driven by intense single-cycle terahertz pulses in a topological insulator Bi 2 Se 3

We present the measurement of ultrafast terahertz (THz) conductivity of a 50 nm thick topological insulator Bi 2 Se 3 at low temperature 5 K by using a THz pump and THz probe spectroscopy scheme. Here, the THz conductivity driven by a single-cycle intense THz pump exhibits a bipolar lineshape and frequency-dependent relaxation dynamics due to the separated surface and bulk charge carriers with distinctly different scattering rates and relaxation times. We also demonstrate THz pump field-dependent excitation of the surface and bulk spectral weights, which suggests one of the optimal THz excitation conditions for generating the most charge carriers from the topology-protected surface state relative to the bulk. This also allows THz control of surface and bulk transport channels in a selective way.

36 MATERIALS SCIENCE↗

Tuning Organic Semiconductor Packing and Morphology through Non-equilibrium Solution Processing (Final Report)

Organic semiconductors (OSCs) are a promising candidate to produce low-cost, large area, and flexible electronics. There has previously been successful development of organic field effect transistors (OFETs), photovoltaics (OPV), and bioelectronics using OSCs. Solution processing of these materials allows for the fabrication of large-area devices in the kinetic crystallization regime. The charge transport capabilities have been shown to depend on the morphology and molecular packing of the OSC within thin films. Our hypothesis is that solution processing conditions may significantly impact OSC morphology and as a result the charge transport. Therefore, this work has focused on obtaining a fundamental understanding of the morphology and molecular packing of OSCs for optimal device performance. Through our work, we have gained a better understanding of how these structural conditions were influenced by the solution processing conditions. Our studies have led to a more systematic understanding of the various parameters that impact OSC morphology. Solution processing leads to non-equilibrium films, thus allowing for the formation of diverse morphologies that are inaccessible by other fabrication methods. Previously, our group has focused on tuning the morphology of small-molecular OSCs [53, 55-57]. However, little work had been devoted to polymer OSCs. There was a lack of detailed studies characterizing the solution-state of polymer OSCs in terms of their conformation, degree of entanglement, polymer aggregation, and chain relaxation dynamics. The solution state properties are also likely to be influenced by the rigidity and molecular weight of the polymer OSCs investigated. Thus, we have focused our attention on gaining a better insight of polymer OSC films and solution processing methods. Through this proposal, our approach is to investigate the correlation between solution-state properties and the morphology of the resulting polymer OSC films. We worked three specific aims: investigate the effects of 1) polymer OSC solution-state properties on final film morphology; 2) molecular additives on the solution-state and final film properties; 3) controlled pre-aggregation in the solution-state in the final film morphology. More rigid and planar polymer backbones should promote interchain charge transport and more efficient interchain hopping between polymer chains. Through tailoring the polymer backbone, polymer sidechains, and molecular weight, we expected the altered solution-state properties to affect the final film morphology. In addition, reducing the entanglements and promoting chain alignment will likely prevent charge carrier trapping through conformation disorders. We thus studied different mechanisms, such as molecular additives, to reduce entanglements in solution. Devices fabricated from these solutions were expected to have improved charge transport abilities. In addition to tailoring the solution-state characteristics of the polymer OSCs, we investigated the effect of solution processing methods on the molecular packing and morphology of polymer OSC films. Two main solution processing techniques, e.g., spin-coating and solution shearing, were employed to fabricate OFETs. Spin coating was employed to prepare OSC films, which produce isotropic films. This technique creates several parameters to tune such as spin coating speed, acceleration, and time. Additionally, our research group developed the solution shearing method, which consists of the solution initially sandwiched between two plates. By sliding the top plate, the solution front is exposed and drying begins. This technique allows for the creation of aligned large crystalline domains. Solution shearing consists of different processing parameters which can affect the final morphology, such as shearing speed, substrate temperature and temperature gradient, distance between both plates, and the tilt angle of the top plate. Due to the complexity of the polymer systems investigated, numerous techniques were employed to characterize the polymer OSC solution state and films. All the materials were characterized using the DOE supported synchrotron X-ray scattering facilities at the Stanford Synchrotron Radiation Lightsource (SSRL). Using grazing incidence X-ray scattering (GIXS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) techniques, the crystalline structure and molecular orientations of the thin films were measured. Optical absorption (UV-Vis) spectroscopy was employed to determine the aggregation state in solution and films of the polymer systems. Polarized UV-Vis also allowed for the determination of the relative degree of polymer chain alignment for solution sheared films. Additionally, various other techniques were used to investigate other properties within the film, such as atomic force microscopy (AFM) and solution rheology. Finally, the device performance is quantified through the fabrication and characterization of OFETs, which will highlight the effects of morphology on charge transport.

36 MATERIALS SCIENCE↗

Neutron scattering study of fluctuating and static spin correlations in the anisotropic spin glass Fe 2 Ti O 5

The anisotropic spin-glass transition, in which spin freezing is observed only along the c axis in pseudobrookite Fe 2 TiO 5 , has long been perplexing because the Fe 3+ moments (d 5 ) are expected to be isotropic. Recently, neutron diffraction demonstrated that surfboard-shaped antiferromagnetic nanoregions coalesce above the glass transition temperature T g ≈55K, and a model was proposed in which the freezing of the surfboard magnetization fluctuations leads to the anisotropic spin-glass state. Given this model, we have carried out high-resolution inelastic neutron scattering measurements of the spin-spin correlations to understand the temperature dependence of the intrasurfboard spin dynamics on neutron (picosecond) timescales. Here, in this paper, we report on the temperature-dependence of the spin fluctuations measured from single-crystal Fe 2 TiO 5 . Strong quasi-elastic magnetic scattering, arising from intrasurfboard correlations, is observed well above T g . The spin fluctuations possess a steep energy–wave vector relation and are indicative of strong exchange interactions, consistent with the large Curie-Weiss temperature. As the temperature approaches T g from above, a shift in spectral weight from inelastic to elastic scattering is observed. At various temperatures between 4 and 300 K, a characteristic relaxation rate of the fluctuations is determined. Despite the freezing of most of the spin correlations, an inelastic contribution remains even at base temperature, signifying the presence of fluctuating intrasurfboard spin correlations to at least T/T g ≈0.1 , consistent with an energy landscape that is a hybrid between conventional and geometrically frustrated spin glasses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comprehensive Dynamics in a Polyelectrolyte Complex Coacervate

The linear viscoelastic response, LVR, of a hydrated polyelectrolyte complex coacervate, PEC, was evaluated over a range of frequencies, temperatures, and salt concentrations. The PEC was a nearly stoichiometric blend of a quaternary ammonium poly([3-(methacrylamido)propyl]trimethylammonium chloride), PMAPTAC, and poly(2-acrylamido-2-methyl-1-propanesulfonic acid sodium salt), PAMPS, an aliphatic sulfonate, selected because they remain fully charged over the conditions of use. Narrow molecular weight distribution polyelectrolytes were prepared using fractionation techniques. A partially deuterated version of PMAPTAC was incorporated to determine its coil radius of gyration, R g , within PECs using small-angle neutron scattering. Chain dimensions were determined to be Gaussian with a Kuhn length of 2.37 nm, which remained constant from 25 to 65 °C. The LVR for a series of matched molecular weight PECs, mostly above the entanglement threshold, exhibited crossovers of modulus versus frequency classically attributed to the reptation time, relaxation between entanglements, and the relaxation of a Kuhn length of units (the “monomer” time). The scaling for zero shear viscosity, η 0 , versus chain length, N, was η 0 ~ N 3.1 , in agreement with “sticky reptation” theory. The lifetime and activation energy, E p , of a pair between polyanion and polycation repeat units, Pol + Pol – , were determined from diffusion coefficients of salt ions within the PEC. The activation energy for LVR of salt-free PECs was 2E p , showing that the key mechanism limiting the dynamics of undoped PECs is pair exchange. An FTIR technique was used to distinguish whether SCN– acts as a counterion or a co-ion within PECs. Doping of PECs with NaSCN breaks Pol + Pol – pairing efficiently, which decreases effective cross-linking and decreases viscosity. Finally, an equation was derived that quantitatively predicts this effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toughening mechanism in elastomer-modified epoxy resins, part 2

The role of matrix ductility on the toughenability and toughening mechanism of elastomer-modified DGEBRA epoxies was investigated. Matrix ductility was varied by using epoxide resins of varying epoxide monomer molecular weights. These epoxide resins were cured using 4,4' diaminodiphenyl sulfone (DDS) and, in some cases, modified with 10% HYCAR(r)CTBN 1300X8. Fracture roughness values for the neat epoxies were found to be almost independent on the monomer molecular weight of the epoxide resin used. However, it was found that the fracture toughness of the elastomer-modified epoxies was very dependent upon the epoxide monomer molecular weight. Tensile dilatometry indicated that the toughening mechanism, when present, is similar to the mechanisms found for the piperidine cured epoxies in Part 1. SEM and OM corroborate this finding. Dynamic mechanical studies were conducted to shed light on the toughenability of the epoxies. The time-dependent small strain behavior of these epoxies were separated into their bulk and shear components. The bulk component is related to brittle fracture, whereas the shear component is related to yielding. It can be shown that the rates of shear and bulk strain energy buildup for a given stress are uniquely determined by the values of Poisson's ratio, nu. It was found that nu increases as the monomer molecular weight of the epoxide resin used increases. This increase in nu can be associated with the low temperature beta relaxation. The effect of increasing cross-link density is to shift the beta relaxation to higher temperatures and to decrease the magnitude of the beta relaxation. Thus, increasing cross-link density decreases nu and increases the tendency towards brittle fracture.

Yee, A. F.↗

The application of high temperature superconductors to space electrical power distribution components

Some important space based electrical power distribution systems and components are examined to determine what might be achieved with the introduction of high temperature superconductors (HTS). Components that are compared in a before and after fashion include transformers, transmission lines, and capacitors. It is concluded that HTS has its greatest effect on the weight associated with transmission lines, where the weight penalty could be reduced by as much as 130 kg/kW/km of cable. Transformers, because 28 percent of their mass is in the conductor, are reduced in weight by the same factor. Capacitors are helped the least with only negligible savings possible. Finally, because HTS can relax the requirement to use alternating current in order to reduce conductor mass, it will be possible to generate significant savings by eliminating most transformers and capacitors.

Aron, Paul R.↗

The application of high temperature superconductors to space electrical power distribution components

Some important space based electrical power distribution systems and components are examined to determine what might be achieved with the introduction of high temperature superconductors (HTS). Components that are compared in a before-and-after fashion include transformers, transmission lines, and capacitors. It is concluded that HTS has its greatest effect on the weight associated with transmission lines, where the weight penalty could be reduced by as much as 130 kg/kW/km of cable. Transformers, because 28 percent of their mass is in the conductor, are reduced in weight by the same factor. Capacitors are helped the least with only negligible savings possible. Finally, because HTS can relax the requirement to use alternating current in order to reduce conductor mass, it will be possible to generate significant savings by eliminating most transformers and capacitors.

Aron, Paul R.↗

Employing Relaxed Smoothness Constraints on Imaginary Part of Refractive Index in AERONET Aerosol Retrieval Algorithm

In the Aerosol Robotic Network (AERONET) retrieval algorithm, smoothness constraints on the imaginary part of the refractive index provide control of retrieved spectral dependence of aerosol absorption by preventing the inversion code from fitting the noise in optical measurements and thus avoiding unrealistic oscillations of retrievals with wavelength. The history of implementation of the smoothness constraints in the AERONET retrieval algorithm is discussed. It is shown that the latest version of the smoothness constraints on the imaginary part of refractive index, termed standard and employed by Version 3 of the retrieval algorithm, should be modified to account for strong variability of light absorption by brown-carbon-containing aerosols in UV through mid-visible parts of the solar spectrum. In Version 3 strong spectral constraints were imposed at high values of the Ångström exponent (440–870 nm) since black carbon was assumed to be the primary absorber, while the constraints became increasingly relaxed as aerosol exponent deceased to allow for wavelength dependence of absorption for dust aerosols. The new version of the smoothness constraints on the imaginary part of the refractive index assigns different weights to different pairs of wavelengths, which are the same for all values of the Ångström exponent. For example, in the case of four-wavelength input, the weights assigned to short-wavelength pairs (440–675, 675–870 nm) are small so that smoothness constraints do not suppress natural spectral variability of the imaginary part of the refractive index. At longer wavelengths (870–1020 nm), however, the weight is 10 times higher to provide additional constraints on the imaginary part of refractive index retrievals of aerosols with a high Ångström exponent due to low sensitivity to aerosol absorption for longer channels at relatively low aerosol optical depths. The effect of applying the new version of smoothness constraints, termed relaxed, on retrievals of single-scattering albedo is analyzed for case studies of different aerosol types: black- and brown-carbon-containing fine mode aerosols, mineral dust coarse mode aerosols, and urban industrial fine mode aerosol. It is shown that for brown-carbon-containing aerosols employing the relaxed smoothness constraints resulted in significant reduction in retrieved single-scattering albedo and spectral residual errors (compared to standard) at the short wavelengths. For example, biomass burning smoke cases showed a reduction in single-scattering albedo and spectral residual error at 380 nm of ∼ 0.033 and ∼ 17 %, respectively, for the Rexburg site and ∼ 0.04 and ∼ 12.7 % for the Rimrock site, both AERONET sites in Idaho, USA. For a site with very high levels of black-carbon-containing aerosols (Mongu, Zambia), the effect of modification in the smoothness constraints was minor. For mineral dust aerosols at small Ångström exponent values (Mezaira site, UAE), the spectral constraint on the imaginary part of the refractive index was already relaxed in Version 3; therefore the new relaxed constraint results in minimal change. In the case of weakly absorbing urban industrial aerosols at the GSFC site, there are significant changes in retrieved single-scattering albedo using relaxed assumption, especially reductions at longer wavelengths: ∼ 0.016 and ∼ 0.02 at 875 and 1020 nm, respectively, for 440 nm aerosol optical depth (AOD) ∼ 0.3. The modification of smoothness constraints on the imaginary part of the refractive index has a minor effect on retrievals of other aerosol parameters such as the real part of the refractive index and parameters of the aerosol size distribution. The implementation of the relaxed smoothness constraints on the imaginary part of the refractive index in the next version of the AERONET inversion algorithm will produce significant impacts at some sites in short wavelength channels (380 and 440 nm) for some biomass burning smoke cases with significant brown carbon content and possibly in mid-visible channels (500 and 675 nm) to near-infrared channels (870 to 1020 nm) for some urban industrial aerosol types. However, most differences in single-scattering albedo retrievals between those applying the new relaxed constraint and the standard constraint will be within the uncertainty of the single-scattering albedo retrievals, depending on the level of aerosol optical depth, Ångström exponent, brown carbon content and wavelength.

Aliaksandr Siniuk↗

Magnetic activity complexes, thermal relaxation oscillations, and the dynamics of the azimuthal magnetic field of a star

The large amount of magnetic flux cycling through an activity complex in a period of six months suggests that the azimuthal magnetic field in the convective zone of the sun is not less than 3000 Gauss and probably is considerably more. A field of 3 kiloGauss or more has the effect of blocking the convective heat transport, producing a cool shadow above the azimuthal field and causing heat to accumulate below. The incremental weight of the cool shadow is sufficient to overpower the magnetic buoyancy, pressing the field against the bottom of the convective zone. The accumulation of heat underneath produces a thermal relaxation oscillation that erupts to the surface at intervals of the order of a week, supplying the intermittent emergence of magnetic flux observed in the activity complex.

Parker, E. N.↗