Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional random fields”

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 73 records · Page 4

Effect of Mach number and volume fraction in air-shock interacting with a bed of randomly distributed spherical particles

Shock-particle interaction is a fundamental problem in many engineering applications, with the dynamics being heavily influenced by the incident-shock Mach number and the particle volume fraction. In this paper we present fully resolved inviscid simulations of an incident-shock wave traveling through a bed of randomly distributed spherical particles. We vary the strength of the incident shock along with the particle volume fraction in order to study the complex wave interaction during shock-particle interaction. In this study we are interested in the early-time behavior during which the particles do not move and hence in our simulations all the particles are fixed in space. We compute the streamwise average of flow field quantities to generate the x - t contour plots to study the unsteady oscillations inside the particle bed. We observe that the transmitted shock slows down under certain conditions and it is partly due to tortuosity and partly due to weakening caused by energy dissipation. We also present the force histories of the streamwise drag and lift forces for all the particles. The random distribution of particles leads to high variability in the drag force experienced by the particles. We compute the mean peak drag force as a function of the streamwise location to study the mean behavior of the transmitted shock. Based on our findings, we propose simple modifications to improve the current point-particle models used in Euler-Lagrange simulations of shock interacting with a bed of particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Influence of cooling rate in planar thermally assisted magnetic random access memory: Improved writeability due to spin-transfer-torque influence

This paper investigates the effect of a controlled cooling rate on magnetic field reversal assisted by spin transfer torque (STT) in thermally assisted magnetic random access memory. By using a gradual linear decrease of the voltage at the end of the write pulse, the STT decays more slowly or at least at the same rate as the temperature. This condition is necessary to make sure that the storage layer magnetization remains in the desired written direction during cooling of the cell. The influence of the write current pulse decay rate was investigated on two exchange biased synthetic ferrimagnet (SyF) electrodes. For a NiFe based electrode, a significant improvement in writing reproducibility was observed using a gradual linear voltage transition. The write error rate decreases by a factor of 10 when increasing the write pulse fall-time from ∼3 ns to 70 ns. For comparison, a second CoFe/NiFe based electrode was also reversed by magnetic field assisted by STT. In this case, no difference between sharp and linear write pulse fall shape was observed. We attribute this observation to the higher thermal stability of the CoFe/NiFe electrode during cooling. In real-time measurements of the magnetization reversal, it was found that Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling in the SyF electrode vanishes for the highest pulse voltages that were used due to the high temperature reached during write. As a result, during the cooling phase, the final state is reached through a spin-flop transition of the SyF storage layer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient and robust quantum random number generation by photon number detection

We present an efficient and robust quantum random number generator based upon high-rate room temperature photon number detection. We employ an electric field-modulated silicon avalanche photodiode, a type of device particularly suited to high-rate photon number detection with excellent photon number resolution to detect, without an applied dead-time, up to 4 photons from the optical pulses emitted by a laser. By both measuring and modeling the response of the detector to the incident photons, we are able to determine the illumination conditions that achieve an optimal bit rate that we show is robust against variation in the photon flux. We extract random bits from the detected photon numbers with an efficiency of 99% corresponding to 1.97 bits per detected photon number yielding a bit rate of 143 Mbit/s, and verify that the extracted bits pass stringent statistical tests for randomness. Our scheme is highly scalable and has the potential of multi-Gbit/s bit rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Velocity overshoot decay mechanisms in compound semiconductor field-effect transistors with a submicron characteristic length

Velocity overshoot is a critically important nonstationary effect utilized for the enhanced performance of submicron field-effect devices fabricated with high-electron-mobility compound semiconductors. However, the physical mechanisms of velocity overshoot decay dynamics in the devices are not known in detail. Therefore, a numerical analysis is conducted typically for a submicron GaAs metal-semiconductor field-effect transistor in order to elucidate the physical mechanisms. It is found that there exist three different mechanisms, depending on device bias conditions. Specifically, at large drain biases corresponding to the saturation drain current (dc) region, the velocity overshoot suddenly begins to drop very sensitively due to the onset of a rapid decrease of the momentum relaxation time, not the mobility, arising from the effect of velocity-randomizing intervalley scattering. It then continues to drop rapidly and decays completely by severe mobility reduction due to intervalley scattering. On the other hand, at small drain biases corresponding to the linear dc region, the velocity overshoot suddenly begins to drop very sensitively due to the onset of a rapid increase of thermal energy diffusion by electrons in the channel of the gate. It then continues to drop rapidly for a certain channel distance due to the increasing thermal energy diffusion effect, and later completely decays by a sharply decreasing electric field. Moreover, at drain biases close to a dc saturation voltage, the mechanism is a mixture of the above two bias conditions. It is suggested that a large secondary-valley energy separation is essential to increase the performance of submicron devices.

36 MATERIALS SCIENCE↗

Dataset for SIAM MPI23 Project "Model inversion for complex physical systems using low-dimensional surrogates"

Training and testing data for the SIAM MPI23 workshop "Model inversion for complex physical systems using low-dimensional surrogates". The data consists of ensembles of input and output pairs corresponding to queries of a 2D saturated groundwater flow model. The inputs consist of random vectors which map to discretized model parameter fields via a Kosambi-Karhunen-Loève expansion. Output corresponds to discretized pressure fields.

basis adaptation, conditional Karhunen-Loève expan↗

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗

Charged particle dynamics in multiple colliding electromagnetic waves. Survey of random walk, Lévy flights, limit circles, attractors and structurally determinate patterns

The multiple colliding laser pulse concept formulated by Bulanov et al. ( Phys. Rev. Lett ., vol. 104, 2010 b , 220404) is beneficial for achieving an extremely high amplitude of coherent electromagnetic field. Since the topology of electric and magnetic fields of multiple colliding laser pulses oscillating in time is far from trivial and the radiation friction effects are significant in the high field limit, the dynamics of charged particles interacting with the multiple colliding laser pulses demonstrates remarkable features corresponding to random walk trajectories, limit circles, attractors, regular patterns and Lévy flights. Lastly, under extremely high intensity conditions the nonlinear dissipation mechanism stabilizes the particle motion resulting in the charged particle trajectory being located within narrow regions and in the occurrence of a new class of regular patterns made by the particle ensembles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Chiral streamers

The interaction of time-varying electromagnetic fields and solid, liquid, and gaseous matter may lead to electrical breakdown phenomena through the excitation of ionization waves or streamers that control the dynamics of localized plasma propagation through the media. The streamers usually propagate along straight lines, either between random points in space or along a certain direction in a guided mode. Here, we report on a new type of plasma discharges with the regular helical propagation pattern driven by a pulsed dc voltage in nitrogen at sub-atmospheric-pressure conditions. The helical guided streamers, named chiral streamers or chi-streamers, are excited without any external magnetic fields, which commonly cause helical plasma motions. We also demonstrate a hybrid propagation mode involving the interchangeable chiral streamers and the straight-line propagating plasmas. High-speed, time-resolved optical imaging reveals that the chiral streamers and the hybrid patterns are made of spatially localized discrete plasma bullets, similar to the straight-line guided streamers. These results may enable effective control of propagation of confined plasmas and electromagnetic energy along pre-determined, potentially deterministic paths, which have important implications for the development of next-generation plasma-based radiation sources, communication devices, and medical treatments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗

A computational study of how surfaces affect slip family activity

Plastic deformation behavior is most conveniently assessed by characterization on a surface, but whether such observations are representative of bulk properties is uncertain. Motivated by reported inconsistencies in slip resistance probed at different depths, we investigated (i) whether the average slip family activity is affected by the presence of a surface and (ii) how the kinematic nature of available slip families influences a potential surface effect. The slip family activity as a function of distance from the surface was extracted from full-field crystal plasticity simulations of random polycrystalline hexagonal close-packed (HCP) and body-centered cubic (BCC) metals as examples of mixed in contrast to universally-high numbers of slip systems per family. Under certain conditions, a deviation from bulk slip activity is observed up to about two grains from the surface. For the easiest (least slip-resistant) family, a surface effect of decreasing activity with depth emerges if the number of slip systems falls below about six. For harder families, slip activity always increases with depth. Furthermore, these phenomena are explained on the basis of varying constraints with depth in connection with the kinematic properties of slip families in the material.

36 MATERIALS SCIENCE↗

The Importance of an Ensemble Approach for Modeling Aerosol‐Convection Interactions

Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.

aerosol-convection interactions↗

Fish Response to Turbulence Generated Using Multiple Randomly Actuated Synthetic Jet Arrays

Hydroelectric power stations generate turbulent flow conditions, which represent a potentially significant hydraulic stressor to fish passing through the turbine system. A test facility has been developed using two randomly actuated synthetic jet arrays (RASJAs) of 25 independent submersible pumps to generate a turbulent flow field for biological dose-response testing. The novel elements of this approach include the ability to control the exposure duration within a test volume due to low mean flow velocity as well as the capacity to scale the turbulence levels as a function of pump capacity. Juvenile Chinook salmon (Oncorhynchus tshawytscha) were subjected to the turbulent flow regime with average turbulence kinetic energy per unit mass of 0.089 m 2 /s 2 for periods of 2 min and 10 min. No significant loss of equilibrium or disorientation was observed after exposure for either duration at the level of turbulence achieved in this prototype. Further scaling of this approach is required to generate a complete dose-response relationship.

59 BASIC BIOLOGICAL SCIENCES↗

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigation of LSM-YSZ Composite Cathode Performance Degradation with a Multistep Charge Transfer Model

The long-term performance degradation of LSM/YSZ composite cathodes is investigated via calibrated multiphysics simulation with a multistep oxygen reduction reaction (ORR) model and structural coarsening data from a phase field study. Multiphysics simulations, including a multistep ORR mechanism, are first developed for better understanding of the limiting processes in solid oxide fuel cells. The multistep ORR mechanism includes two parallel pathways, surface pathway and bulk pathway, which consist of elementary reaction steps. The multiphysics model is simultaneously calibrated with experimental polarization curves and impedance behavior for various air/fuel supply conditions. To our knowledge, this is the first time that a multistep ORR model is calibrated with such datasets. Next, the calibrated simulations are utilized to simulate a 2D half-cell constructed with measured microstructural data and random heterogeneity. Finally, the long-term performance degradation of the half-cell is predicted by the calibrated multiphysics model coupled with structural coarsening trends simulated using a phase field-based coarsening model. Degradation of both polarization curves and impedance behavior is investigated. Thorough analyses, including changes of contributions from different pathways, the resistance components, and overall reaction order, are performed to provide more insights into cathode performance degradation due to grain coarsening phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution of the magnetorotational instability on initially tangled magnetic fields

The initial magnetic field of previous magnetorotational instability (MRI) simulations has always included a significant system-scale component, even if stochastic. However, it is of conceptual and practical interest to assess whether the MRI can grow when the initial field is turbulent. The ubiquitous presence of turbulent or random flows in astrophysical plasmas generically leads to a small-scale dynamo (SSD), which would provide initial seed turbulent velocity and magnetic fields in the plasma that becomes an accretion disc. Can the MRI grow from these more realistic initial conditions? To address this, we supply a standard shearing box with isotropically forced SSD generated magnetic and velocity fields as initial conditions and remove the forcing. We find that if the initially supplied fields are too weak or too incoherent, they decay from the initial turbulent cascade faster than they can grow via the MRI. When the initially supplied fields are sufficient to allow MRI growth and sustenance, the saturated stresses, large-scale fields and power spectra match those of the standard zero net flux MRI simulation with an initial large-scale vertical field.

79 ASTRONOMY AND ASTROPHYSICS↗

The Helicity Sign of Flux Transfer Event Flux Ropes and Its Relationship to the Guide Field and Hall Physics in Magnetic Reconnection at the Magnetopause

Flux Transfer Events (FTEs) are transient magnetic flux ropes typically found at the Earth's magnetopause on the dayside. While it is known that FTEs are generated by magnetic reconnection, it remains unclear how the details of magnetic reconnection controls their properties. A recent study showed that the helicity sign of FTEs positively correlates with the east-west (B y ) component of the Interplanetary Magnetic Field (IMF). With data from the Cluster and Magnetospheric Multiscale missions, we performed a statistical study of 166 quasi force-free FTEs. We focus on their helicity sign and possible association with upstream solar wind conditions and local magnetic reconnection properties. Using both in situ data and magnetic shear modeling, we find that FTEs whose helicity sign corresponds to the IMF B y are associated with moderate magnetic shears while those that do not correspond to the IMF B y are associated with higher magnetic shears. While uncertainty in IMF propagation to the magnetopause may lead to randomness in the determination of the flux rope core field and helicity, we rather propose that for small IMF B y , which corresponds to high shear and low guide field, the Hall pattern of magnetic reconnection determines the FTE core field and helicity sign. In that context we explain how the temporal sequence of multiple X-line formation and the reconnection rate are important in determining the flux rope helicity sign. This work highlights a fundamental connection between kinetic processes at work in magnetic reconnection and the macroscale structure of FTEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Molecular bond stabilization in the strong-field dissociation of O 2 +

We theoretically examine the rotational and vibrational dynamics of O 2 + molecular ions exposed to intense, short laser pulses for conditions realized in contemporary pump-probe experiments. We solve the time-dependent Schrödinger equation within the Born-Oppenheimer approximation for an initial distribution of randomly aligned molecular ions. For fixed peak intensities, our numerical results show that total, angle-integrated O 2 + → O( 3 P) + O + ( 4 S 0 ) dissociation yields do not monotonically increase with increasing infrared-probe pulse duration. We find this pulse-duration-dependent stabilization to be consistent with the transient trapping of nuclear probability density in a light-induced (bond-hardening) potential-energy surface and robust against rotational excitation. We analyze this stabilization effect and its underlying bond-hardening mechanism (i) in the time domain, by following the evolution of partial nuclear probability densities associated with the dipole-coupled O 2 + (a 4 Π u ) and O 2 + (f 4 Π g ) cationic states, and (ii) in the frequency domain, by examining rovibrational quantum-beat spectra for the evolution of the partial nuclear probability densities associated with these states. In this work, our analysis reveals the characteristic timescale for the bond-hardening mechanism in O 2 + and explains the onset of bond stabilization for sufficiently long pulse durations.

74 ATOMIC AND MOLECULAR PHYSICS↗

A non‐intrusive domain‐decomposition model reduction method for linear steady‐state partial differential equations with random coefficients

Abstract Domain decomposition methods have been proved to be an effective strategy to reduce the dimension of parametric partial differential equations (PDEs). However, existing domain decomposition methods for parametric PDEs are usually intrusive, which means domain decomposition based solvers need to be implemented from scratch for each target parametric PDE. To address this issue, we develop a new non‐intrusive domain‐decomposition model reduction method for linear steady‐state PDEs with random‐field coefficients. As a variant of our previous work by Mu and Zhang, the new method only needs access to the final linear system, that is, the global stiffness matrix and the right hand side, of a deterministic PDE solver, in order to build a domain‐decomposition‐based reduced model without intrusive implementation from scratch. The key idea is to remove the interface condition between sub‐domains and rely on the correlation between columns of the linear system to couple the sub‐domains. The non‐intrusive feature enables the applicability of the proposed method to a broader class of uncertainty quantification problems, where many legacy codes/solvers can be fully reused by our method. Two numerical examples including diffusion equations with random diffusivity and convection‐dominated transport with random velocity, are provided to demonstrate the effectiveness and efficiency of our method.

Zhang, Guannan↗