Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “contact model”

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 145 records · Page 8

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

New Approach Methodology for Assessing Inhalation Risks of a Contact Respiratory Cytotoxicant: Computational Fluid Dynamics-Based Aerosol Dosimetry Modeling for Cross-Species and In Vitro Comparisons

Regulatory agencies are considering alternative approaches to assessing inhalation toxicity that utilizes in vitro studies with human cells and in silico modeling in lieu of additional animal studies. In support of this goal, computational fluid-particle dynamics models were developed to estimate site-specific deposition of inhaled aerosols containing the fungicide, chlorothalonil, in the rat and human for comparisons to prior rat inhalation studies and new human in vitro studies. Under bioassay conditions, the deposition was predicted to be greatest at the front of the rat nose followed by the anterior transitional epithelium and larynx corresponding to regions most sensitive to local contact irritation and cytotoxicity. For humans, simulations of aerosol deposition covering potential occupational or residential exposures (1–50 µm diameter) were conducted using nasal and oral breathing. Aerosols in the 1–5 µm range readily penetrated the deep region of the human lung following both oral and nasal breathing. Under actual use conditions (aerosol formulations >10 µm), the majority of deposited doses were in the upper conducting airways. Beyond the nose or mouth, the greatest deposition in the pharynx, larynx, trachea, and bronchi was predicted for aerosols in the 10–20 µm size range. Only small amounts of aerosols >20 µm penetrated past the pharyngeal region. Using the ICRP clearance model, local retained tissue dose metrics including maximal concentrations and areas under the curve were calculated for each airway region following repeated occupational exposures. These results are directly comparable with benchmark doses from in vitro toxicity studies in human cells leading to estimated human equivalent concentrations that reduce the reliance on animals for risk assessments.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Impacts of Non-Ideal Back Contact on Capacitance Measurements in CdTe Solar Cells

CdTe solar cells suffer from a non-ideal back contact which can strongly affect the equivalent circuit model and complicate capacitance measurements. Here, five different back contact materials are deposited on identical CdTe absorbers and their influence on capacitance measurements is characterized. A five-element equivalent circuit model is shown to represent a CdTe solar cell with a non-ideal back contact, and capacitance-frequency (C-f) measurements on each sample clearly show the influence of this contact. Capacitance-voltage (C-V) measurements on each sample show the influence of frequency on measured capacitance. In the AlGaOx/Au sample, temperature dependent C-fs show the back barrier height to be 240 meV, and ongoing measurements will give a comparison of back barrier heights for all samples. This study provides insight into the efficacy of various back contacts and highlights potential errors in the interpretation of capacitance measurements due to the presence of the back contact.

cadmium compounds↗

Changes in Hydrogen Concentration and Defect State Density at the Poly-Si/SiOx/c-Si Interface Due to Firing

We determined the density of defect states of poly-Si/SiOx/c-Si junctions featuring a wet chemical interfacial oxide from lifetime measurements using the MarcoPOLO model to calculate recombination and contact resistance in poly-Si/SiOx/c-Si-junctions. In samples that did not receive any hydrogen treatment, the Dit,cSi is about 2 × 1012 cm-2 eV1 before firing and rises to 3–7 × 1012 cm2 eV1 during firing at measured peak temperatures between 620 °C and 863 °C. To address the question of why AlOx/SiNy stacks in contrast to pure SiNy layers for hydrogenation during firing provides better passivation quality, we have measured the hydrogen concentrations at the poly-Si/SiOx/c-Si interface as a function of AlOx layer thickness and compared these to J0 and calculated Dit,c-Si values. We observe an increase of the hydrogen concentration at the SiOx/c-Si interface upon firing as a function of the firing temperature that exceeds the defect concentrations at the interface several times. However, the AlOx layer thickness appears to cause an increase in hydrogen concentration at the SiOx/c-Si interface in these samples rather than exhibiting a hydrogen blocking property.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Impacts of Non-Ideal Back Contact on Capacitance Measurements in CdTe Solar Cells

CdTe solar cells suffer from a non-ideal back contact which can strongly affect the equivalent circuit model and complicate capacitance measurements. Here, five different back contact materials are deposited on identical CdTe absorbers and their influence on capacitance measurements is characterized. A five-element equivalent circuit model is shown to represent a CdTe solar cell with a non-ideal back contact, and capacitance-frequency (C-f) measurements on each sample clearly show the influence of this contact. Capacitance-voltage (C-V) measurements on each sample show the influence of frequency on measured capacitance. In the AlGaOx/Au sample, temperature dependent C-fs show the back barrier height to be 240 meV, and ongoing measurements will give a comparison of back barrier heights for all samples. This study provides insight into the efficacy of various back contacts and highlights potential errors in the interpretation of capacitance measurements due to the presence of the back contact.

14 SOLAR ENERGY↗

Development of a co-anneal process for double-side TOPCon precursor fabricated by ex-situ POCI3 and APCVD boron diffusion

The aim of this study was to develop a simple and industrially attractive co-anneal process to fabricate a high-quality DS-TOPCon precursor with textured n-TOPCon on front and planar p- TOPCon on rear by ex-situ POCI3 and APCVD boron diffusion. This requires only one high temperature anneal with no additional masking steps. Excellent iVOC of 733mV and iFF of ~86% were achieved after SiNX passivation on both sides, prior to contact firing. Our device modeling projects that this precursor in combination with a manufacturing-friendly poly-Si patterning technique on front can enable> 25% DS-TOPCon cells at low-cost.

APCVD↗

Dynamic condensation model of rolling droplets for high-performance heat transfer

Condensation is a critical process during vapor-liquid phase change in relation to heat transfer. To achieve a high heat transfer coefficient, the classical model for dropwise condensation requires a low contact angle and low contact angle hysteresis, failing to align with experimental observations on a hydrophobic and slippery quasi-liquid surface (QLS). Here, we report a dynamic condensation model that incorporates high-frequency condensate removal by emphasizing the role of timescale during droplet growth and shedding. Our model agrees well with the experimental result that a surface with high contact angle and low contact angle hysteresis promotes condensation, particularly during rolling-propelled condensate removal. Particle image velocimetry reveals that rolling droplets on a hydrophobic QLS exhibit 4-fold higher shedding speeds than the sliding droplets on a hydrophilic QLS, leading to significant heat transfer enhancement. This work deepens our theoretical understanding of condensation heat transfer and provides advanced physics-informed design rationales for water and energy systems.

42 ENGINEERING↗

Structure determination of the HgcAB complex using metagenome sequence data: insights into microbial mercury methylation

Bacteria and archaea possessing the hgcAB gene pair methylate inorganic mercury (Hg) to form highly toxic methylmercury. HgcA consists of a corrinoid binding domain and a transmembrane domain, and HgcB is a dicluster ferredoxin. However, their detailed structure and function have not been thoroughly characterized. We modeled the HgcAB complex by combining metagenome sequence data mining, coevolution analysis, and Rosetta structure calculations. In addition, we overexpressed HgcA and HgcB in Escherichia coli, confirmed spectroscopically that they bind cobalamin and [4Fe-4S] clusters, respectively, and incorporated these cofactors into the structural model. Surprisingly, the two domains of HgcA do not interact with each other, but HgcB forms extensive contacts with both domains. The model suggests that conserved cysteines in HgcB are involved in shuttling HgII, methylmercury, or both. These findings refine our understanding of the mechanism of Hg methylation and expand the known repertoire of corrinoid methyltransferases in nature.

59 BASIC BIOLOGICAL SCIENCES↗

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗

Electromagnetic signals of inelastic dark matter scattering

Light dark sectors in thermal contact with the Standard Model can naturally produce the observed relic dark matter abundance and are the targets of a broad experimental search program. A key light dark sector model is the pseudo-Dirac fermion with a dark photon mediator. The dynamics of the fermionic excited states are often neglected. We consider scenarios in which a nontrivial abundance of excited states is produced and their subsequent de-excitation yields interesting electromagnetic signals in direct detection experiments. We study three mechanisms of populating the excited state: a primordial excited fraction, a component up-scattered in the Sun, and a component up-scattered in the Earth. We find that the fractional abundance of primordial excited states is generically depleted to exponentially small fractions in the early universe. Nonetheless, this abundance can produce observable signals in current dark matter searches. MeV-scale dark matter with thermal cross sections and higher can be probed by down-scattering following excitation in the Sun. Up-scatters of GeV-scale dark matter in the Earth can give rise to signals in current and upcoming terrestrial experiments and X-ray observations. We comment on the possible relevance of these scenarios to the recent excess in XENON1T.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Peridigm Meshfree Peridynamics Code

Abstract Peridigm is a meshfree peridynamics code written in C++ for use on large-scale parallel computers. It was originally developed at Sandia National Laboratories and is currently managed as an open-source, community driven software project. Its primary features include bond-based, state-based, and non-ordinary state-based constitutive models, bond failure laws, contact, and support for explicit and implicit time integration. To date, Peridigm has been used primarily by methods developers focused on solid mechanics and material failure. Peridigm utilizes foundational software components from Sandia’s Trilinos project and was designed for extensibility. This paper provides an overview of the solution methods implemented in Peridigm , a discussion of its software infrastructure, and demonstrates the use of Peridigm for the solution of several example problems.

97 MATHEMATICS AND COMPUTING↗

The homogenous alternative to biomineralization: Zn- and Mn-rich materials enable sharp organismal “tools” that reduce force requirements

We measured hardness, modulus of elasticity, and, for the first time, loss tangent, energy of fracture, abrasion resistance, and impact resistance of zinc- and manganese-enriched materials from fangs, stings and other “tools” of an ant, spider, scorpion and nereid worm. The mechanical properties of the Zn- and Mn-materials tended to cluster together between plain and biomineralized “tool” materials, with the hardness reaching, and most abrasion resistance values exceeding, those of calcified salmon teeth and crab claws. Atom probe tomography indicated that Zn was distributed homogeneously on a nanometer scale and likely bound as individual atoms to more than ¼ of the protein residues in ant mandibular teeth. This homogeneity appears to enable sharper, more precisely sculpted “tools” than materials with biomineral inclusions do, and also eliminates interfaces with the inclusions that could be susceptible to fracture. Based on contact mechanics and simplified models, we hypothesize that, relative to plain materials, the higher elastic modulus, hardness and abrasion resistance minimize temporary or permanent tool blunting, resulting in a roughly 2/3 reduction in the force, energy, and muscle mass required to initiate puncture of stiff materials, and even greater force reductions when the cumulative effects of abrasion are considered. We suggest that the sharpness-related force reductions lead to significant energy savings, and can also enable organisms, especially smaller ones, to puncture, cut, and grasp objects that would not be accessible with plain or biomineralized “tools”.

59 BASIC BIOLOGICAL SCIENCES↗

Four-dimensional chromosome reconstruction elucidates the spatiotemporal reorganization of the mammalian X chromosome

Significance Breakthrough technological advances such as Hi-C have enabled large-scale sequencing of pairwise genomic contacts. However, translating the 2D contacts to 3D spatial models remains challenging. Here we develop the “4DHiC” modeling method to infer 3D information from 2D datasets across time and construct a 3D model of the inactive X-chromosome (Xi). Our model not only demonstrates a spatial phase-separation between A/B and S1/S2 compartments to form “hemispheres,” but also reveals persistent smaller-scale structures hidden under the Xi superstructure. We also simulate Xist RNA spreading dynamics and deduce that the RNA spreads between hemispheres via both surface and core to establish the Xi. Our 4DHiC will serve a growing demand for data-driven methodologies that extract spatial information from biochemical interaction data.

Lappala, Anna↗

Extracting forces from noisy dynamics in dusty plasmas

Extracting environmental forces from noisy data is a common yet challenging task in complex physical systems. Machine learning (ML) represents a robust approach to this problem, yet is mostly tested on simulated data with known parameters. Here we use supervised ML to extract the electrostatic, dissipative, and stochastic forces acting on micron-sized charged particles levitated in an argon plasma (dusty plasma). By tracking the sub-pixel motion of particles in subsequent images, we successfully estimated these forces from their random motion. The experiments contained important sources of non-Gaussian noise, such as drift and pixel-locking, representing a data mismatch from methods used to analyze simulated data with purely Gaussian noise. Our model was trained on simulated particle trajectories that included all of these artifacts, and used more than 100 dynamical and statistical features, resulting in a prediction with 50\% better accuracy than conventional methods. Lastly, in systems with two interacting particles, the model provided non-contact measurements of the particle charge and Debye length in the plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of Combined Analytical Methods to Predict Slip in Cylindrical Roller Bearings

This paper presents a combination of models that together calculate the cage and roller speeds of a cylindrical roller bearing. The models consider elastohydrodynamic lubrication and contact elasticity between the roller and raceway, roller centrifugal forces, hydrodynamic lubrication at the cage pocket, and frictional forces. Using these models, the predicted cage and roller speeds and the extent of slip are compared to measurements acquired on cylindrical roller bearings in a commercial gearbox in steady-state and transient operating conditions of a wind turbine. In steady-state conditions at low wind speeds and low lubricant temperatures, cage and roller slip up to 60% occur in the loaded zone of the bearing. In the unloaded zone, up to 80% roller slip occurs. Cage and roller slip then decrease as the lubricant temperature and wind speed increases. In general, the analytical model results match experimental measurements within 10% for lubricant temperatures above 40 degrees C and wind speeds over 10 meters per second. The analytical model is further evaluated during a transient start-up event and highly dynamic emergency stop event and also used to examine changes in the bearing design or lubricant properties. Furthermore, roller and cage slip in cylindrical bearings is a combined effect of the bearing design, applied load, shaft speed, and lubricant properties and temperature, and can be quickly evaluated with the combined analytical models.

17 WIND ENERGY↗

Development of a Flotation Simulator That Can Predict Grade vs. Recovery Curves from Mineral Liberation Data

Flotation is a kinetic process designed to separate mineral fines by control of hydrophobicity. In general, both the recovery and selectivity of a flotation process increase with contact angle. It is difficult, however, to model flotation kinetics using the thermodynamic parameter. In the present work, we consider flotation as heterocoagulation so that the bubble-particle interactions occurring in the pulp phase of a flotation cell can be modeled using the hydrophobic force as a kinetic parameter that can represent the role of contact angles in flotation. Here, the pulp-phase model developed in this manner has been combined with a froth-phase model developed earlier by considering the effect of particles in film thinning kinetics and bubble coarsening (Park et al, 2018). The combined model has been used to develop a flotation simulator that can be used to construct grade vs. recovery curves from feed characteristics (e.g., mineral liberation and particle size) and various operating parameters. The simulator has been tested successfully against a set of bench-scale continuous flotation test results reported in the literature (dos Santos and Galery, 2018).

58 GEOSCIENCES↗

Trap-Assisted Dopant Compensation Prevents Shunting in Poly-Si Passivating Interdigitated Back Contact Silicon Solar Cells

Using a trap-assisted compensation model, we explain why polycrystalline Si (poly-Si) passivating contacts are able to achieve low leakage current between the doped fingers of interdigitated back contact (IBC) monocrystalline Si solar cells despite mixing of boron and phosphorus dopants in the isolation region. The fill factor of IBC solar cells is strongly affected by the electrical isolation region between n- and p-type fingers, as this region is critical in minimizing shunting losses. During fabrication of monocrystalline Si solar cells with poly-Si passivating contacts, the intrinsic poly-Si isolation region inevitably gets contaminated with both p- and n-type dopants. Using dopant profiles measured with time-of-flight secondary ion mass spectrometry and scanning probe measurements of the isolation region, we demonstrate that despite the dopant spreading during cell processing, a well-compensated region between the doped fingers exists that prevents shunting. The trap-assisted dopant compensation mechanism significantly widens the compensated region to tens of microns, where the residual dopant densities are below the trap density. This enables a high-resistivity region, resulting in low shunt current. Using one-dimensional (1-D) finite element diode simulations, we identify the design parameters and experimental conditions under which a sufficiently resistive region can form. Furthermore, our measurements of 2-D local resistivity and work function maps across the isolation region using scanning spreading resistance microscopy and Kelvin probe force microscopy demonstrate the existence of a highly resistive, wide compensated region and confirm our proposed compensation mechanism. For our structures, this region is ~25 µm in width within a ~150 µm wide finger isolation region with nearly 3 orders of magnitude higher resistivity than the regions dominated by a single type of dopant.

14 SOLAR ENERGY↗

Nonlinear model of infection wavy oscillation of COVID-19 in Japan based on diffusion kinetics

The infectious propagation of SARS-CoV-2 is continuing worldwide, and specifically, Japan is facing severe circumstances. Medical resource maintenance and action limitations remain the central measures. An analysis of long-term follow-up reports in Japan shows that the infection number follows a unique wavy oscillation, increasing and decreasing over time. However, only a few studies explain the infection wavy oscillation. This study introduces a novel nonlinear mathematical model of the new infection wavy oscillation by applying the macromolecule diffusion theory. In this model, the diffusion coefficient that depends on population density gives nonlinearity in infection propagation. As a result, our model accurately simulated infection wavy oscillations, and the infection wavy oscillation frequency and amplitude were closely linked with the recovery rate of infected individuals. In conclusion, our model provides a novel nonlinear contact infection analysis framework.

60 APPLIED LIFE SCIENCES↗