BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Dierk Raabe

Dierk Raabe appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.softphysics.comp-phcond-mat.othercond-mat.stat-mechphysics.app-phphysics.flu-dynBiological Physicscond-mat.dis-nncond-mat.mes-hallMachine Learningphysics.med-ph

Catalog footprint

What is connected

35works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2025arXiv

Hydrogen-based direct reduction of multicomponent oxides: Insights from powder and pre-sintered precursors toward sustainable alloy design

The co-reduction of metal oxide mixtures using hydrogen as a reductant in conjunction with compaction and sintering of the evolving metallic blends offers a promising alternative toward sustainable alloy production through a single, integrated, and synergistic process. Herein, we provide fundamental insights into hydrogen-based direct reduction (HyDR) of distinct oxide precursors that differ by phase composition and morphology. Specifically, we investigate the co-reduction of multicomponent metal oxides targeting a 25Co-25Fe-25Mn-25Ni (at.%) alloy, by using either a compacted powder (mechanically mixed oxides) comprising Co3O4-Fe2O3-Mn2O3-NiO or a pre-sintered compound (chemically mixed oxides) comprising a Co,Ni-rich halite and a Fe,Mn-rich spinel. Thermogravimetric analysis (TGA) at a heating rate of 10 °C/min reveals that the reduction onset temperature for the compacted powder was ~175 °C, whereas it was significantly delayed to ~525 °C for the pre-sintered sample. Nevertheless, both sample types attained a similar reduction degree (~80%) after isothermal holding for 1 h at 700 °C. Phase analysis and microstructural characterization of reduced samples confirmed the presence of metallic Co, Fe, and Ni alongside MnO. A minor fraction of Fe remains unreduced, stabilized in the (Fe,Mn)O halite phase, in accord with thermodynamic calculations. Furthermore, ~1 wt.% of BCC phase was found only in the reduced pre-sintered sample, owing to the different reduction pathways. The kinetics and thermodynamics effects were decoupled by performing HyDR experiments on pulverized pre-sintered samples. These findings demonstrate that initial precursor states influence both the reduction behavior and the microstructural evolution, providing critical insights for the sustainable production of multicomponent alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The fate of water in hydrogen-based iron oxide reduction

Gas-solid reactions are cornerstones of many catalytic and redox processes that will underpin the energy and sustainability transition. The specific case of hydrogen-based iron oxide reduction is the foundation to render the global steel industry fossil-free, an essential target as iron production is the largest single industrial emitter of carbon dioxide. Our perception of gas-solid reactions has not only been limited by the availability of state-of-the-art techniques which can delve into the reacted solids in great structural and chemical detail, but we continue to miss an important reaction partner that defines the thermodynamics and kinetics of gas phase reactions: the gas molecules. In this investigation, we use the latest development in cryogenic atom probe tomography to study the quasi in-situ evolution of gas phase heavy water at iron-iron oxide interfaces resulting from the direct reduction of iron oxide by deuterium gas at 700°C. The findings provide new insights into the formation kinetics and location of water formed during hydrogen-based reduction of FeO, an its interaction with the ongoing redox reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A mechanically strong and ductile soft magnet with extremely low coercivity

Soft magnetic materials (SMMs) serve in electrical applications and sustainable energy supply, allowing magnetic flux variation in response to changes in applied magnetic field, at low energy loss1. The electrification of transport, households and manufacturing leads to an increase in energy consumption due to hysteresis losses2. Therefore, minimizing coercivity, which scales these losses, is crucial3. Yet, meeting this target alone is not enough: SMMs in electrical engines must withstand severe mechanical loads, i.e., the alloys need high strength and ductility4. This is a fundamental design challenge, as most methods that enhance strength introduce stress fields that can pin magnetic domains, thus increasing coercivity and hysteretic losses5. Here, we introduce an approach to overcome this dilemma. We have designed a Fe-Co-Ni-Ta-Al multicomponent alloy with ferromagnetic matrix and paramagnetic coherent nanoparticles (~91 nm size, ~55% volume fraction). They impede dislocation motion, enhancing strength and ductility. Their small size, low coherency and small magnetostatic energy create an interaction volume below the magnetic domain wall width, leading to minimal domain wall pinning, thus maintaining the soft magnetic properties. The alloy has a tensile strength of 1336 MPa at 54% tensile elongation, extremely low coercivity of 78 A/m (<1 Oe), moderate saturation magnetization of 100 Am2/kg, and high electrical resistivity of 103 μΩ u Ohm cm.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Accelerating phase-field-based simulation via machine learning

Phase-field-based models have become common in material science, mechanics, physics, biology, chemistry, and engineering for the simulation of microstructure evolution. Yet, they suffer from the drawback of being computationally very costly when applied to large, complex systems. To reduce such computational costs, a Unet-based artificial neural network is developed as a surrogate model in the current work. Training input for this network is obtained from the results of the numerical solution of initial-boundary-value problems (IBVPs) based on the Fan-Chen model for grain microstructure evolution. In particular, about 250 different simulations with varying initial order parameters are carried out and 200 frames of the time evolution of the phase fields are stored for each simulation. The network is trained with 90% of this data, taking the $i$-th frame of a simulation, i.e. order parameter field, as input, and producing the $(i+1)$-th frame as the output. Evaluation of the network is carried out with a test dataset consisting of 2200 microstructures based on different configurations than originally used for training. The trained network is applied recursively on initial order parameters to calculate the time evolution of the phase fields. The results are compared to the ones obtained from the conventional numerical solution in terms of the errors in order parameters and the system's free energy. The resulting order parameter error averaged over all points and all simulation cases is 0.005 and the relative error in the total free energy in all simulation boxes does not exceed 1%.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An artificial neural network for surrogate modeling of stress fields in viscoplastic polycrystalline materials

The purpose of this work is the development of an artificial neural network (ANN) for surrogate modeling of the mechanical response of viscoplastic grain microstructures. To this end, a U-Net-based convolutional neural network (CNN) is trained to account for the history dependence of the material behavior. The training data take the form of numerical simulation results for the von Mises stress field under quasi-static tensile loading. The trained CNN (tCNN) can accurately reproduce both the average response as well as the local von Mises stress field. The tCNN calculates the von Mises stress field of grain microstructures not included in the training dataset about 500 times faster than its calculation based on the numerical solution with a spectral solver of the corresponding initial-boundary-value problem. The tCNN is also successfully applied to other types of microstructure morphologies (e.g., matrix-inclusion type topologies) and loading levels not contained in the training dataset.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Combined modeling and experimental characterization of Mn segregation and spinodal decomposition along dislocation lines in Fe-Mn alloys

Mn enrichment at dislocations in Fe-Mn alloys due to segregation and spinodal decomposition along the dislocation line is studied via modeling and experimental characterization. To model these phenomena, both finite-deformation microscopic phase-field chemomechanics (MPFCM) and Monte Carlo molecular dynamics (MCMD) are employed. MPFCM calibration is carried out with the same Fe-Mn MEAM-based potential used in MCMD, as well as CALPHAD data. Simulation results for Mn segregation to, and spinodal decomposition along, straight screw and edge dislocations as well as dislocation loops, are compared with characterization results from atom probe tomography (APT) for two Fe-Mn alloy compositions. In contrast to classical Volterra dislocation theory, both MPFCM and MCMD predict a non-zero hydrostatic stress field in screw cores. Being of much smaller magnitude than the hydrostatic stress in straight edge cores, much less solute segregates to screw than to edge cores. In addition, the segregated amount in screw cores is below the critical concentration of 0.157 for the onset of spinodal decomposition along the line. Results from MPFCM-based modeling imply that the concentration dependence of the solute misfit distortion and resulting dependence of the elastic energy density on concentration have the strongest effect. The maximum amount of Mn segregating to straight edge dislocations predicted by MPFCM agrees well with APT results. The MPFCM model for Fe-Mn predicts little or no variation in Mn concentration along a straight dislocation line, in contrast to the APT results. As shown by the example of a dislocation loop in the current work, a change in the hydrostatic stress along the line due to changing character of dislocation does lead to a corresponding variation in Mn concentration. Such a variation in Mn concentration can also then be expected along a dislocation line with kinks or jogs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Green steel at its crossroads: hybrid hydrogen-based reduction of iron ores

Iron- and steelmaking cause ~7% of the global CO2 emissions, due to the use of carbon for the reduction of iron ores. Replacing carbon by hydrogen as the reductant offers a pathway to reduce emissions. However, production of hydrogen using renewable energy will remain a bottlenecks, because making the annual crude steel production of 1.8 billion tons sustainable requires a minimum amount of ~97 million tons of green hydrogen per year. Another fundamental aspect to make ironmaking sector more sustainable lies in an optimal utilization of green hydrogen and energy, thus reducing efforts for costly in-process hydrogen recycling. We therefore demonstrate here how the efficiency in hydrogen and energy consumption during iron ore reduction can be dramatically improved by the knowledge-based combination of two technologies: partially reducing the ore at low temperature via solid-state direct reduction (DR) to a kinetically defined degree, and subsequently melting and completely transforming it to iron under a reducing plasma (i.e. via hydrogen plasma reduction, HPR). Results suggest that an optimal transition point between these two technologies occurs where their efficiency in hydrogen utilization is equal. We found that the reduction of hematite through magnetite into wustite via DR is clean and efficient, but it gets sluggish and inefficient when iron forms at the outermost layers of the iron ore pellets. Conversely, HPR starts violent and unstable with arc delocalization, but proceeds smoothly and efficiently when processing semi-reduced oxides. We performed hybrid reduction experiments by partially reducing hematite pellets via DR at 700°C to 38% global reduction (using a standard thermogravimetry system) and subsequently transferring them to HPR, conducted with a gas mixture of Ar-10%H2 in an arc-melting furnace, to achieve conversion into liquid iron.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hierarchical nature of hydrogen-based direct reduction of iron oxides

Fossil-free ironmaking is indispensable for reducing massive anthropogenic CO2 emissions in the steel industry. Hydrogen-based direct reduction (HyDR) is among the most attractive solutions for green ironmaking, with high technology readiness. The underlying mechanisms governing this process are characterized by a complex interaction of several chemical (phase transformations), physical (transport), and mechanical (stresses) phenomena. Their interplay leads to rich microstructures, characterized by a hierarchy of defects ranging across several orders of magnitude in length, including vacancies, dislocations, internal interfaces, and free surfaces in the form of cracks and pores. These defects can all act as reaction, nucleation, and diffusion sites, shaping the overall reduction kinetics. A clear understanding of the roles and interactions of these dynamically-evolving nano-/microstructure features is missing. Gaining better insights in these effects could enable improved access to the microstructure-based design of more efficient HyDR methods, with potentially high impact on the urgently needed decarbonization in the steel industry.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hydrogen trapping and embrittlement in high-strength Al-alloys

Ever more stringent regulations on greenhouse gas emissions from transportation motivate efforts to revisit materials used for vehicles. High-strength Al-alloys often used in aircrafts could help reduce the weight of automobiles, but are susceptible to environmental degradation. Hydrogen (H) "embrittlement" is often pointed as the main culprit, however, the mechanisms underpinning failure are elusive: atomic-scale analysis of H inside an alloy remains a challenge, and this prevents deploying alloy design strategies to enhance the materials' durability. Here we successfully performed near-atomic scale analysis of H trapped in second-phase particles and at grain boundaries in a high-strength 7xxx Al-alloy. We used these observations to guide atomistic ab-initio calculations which show that the co-segregation of alloying elements and H favours grain boundary decohesion, while the strong partitioning of H into the second-phases removes solute H from the matrix, hence preventing H-embrittlement. Our insights further advance the mechanistic understanding of H-assisted embrittlement in Al-alloys, emphasizing the role of H-traps in retarding cracking and guiding new alloy design.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine learning-enabled high-entropy alloy discovery

High-entropy alloys are solid solutions of multiple principal elements, capable of reaching composition and feature regimes inaccessible for dilute materials. Discovering those with valuable properties, however, relies on serendipity, as thermodynamic alloy design rules alone often fail in high-dimensional composition spaces. Here, we propose an active-learning strategy to accelerate the design of novel high-entropy Invar alloys in a practically infinite compositional space, based on very sparse data. Our approach works as a closed-loop, integrating machine learning with density-functional theory, thermodynamic calculations, and experiments. After processing and characterizing 17 new alloys (out of millions of possible compositions), we identified 2 high-entropy Invar alloys with extremely low thermal expansion coefficients around 2*10-6 K-1 at 300 K. Our study thus opens a new pathway for the fast and automated discovery of high-entropy alloys with optimal thermal, magnetic and electrical properties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Phase-Field Modeling of Coupled Brittle-Ductile Fracture in Aluminum Alloys

Fracture in aluminum alloys with precipitates involves at least two mechanisms, namely, ductile fracture of the aluminum-rich matrix and brittle fracture of the precipitates. In this work, a coupled crystal plasticity-phase field model for mixed ductile-brittle failure modes is formulated and used to investigate the effect of precipitate morphology and size distribution on damage evolution in aluminum alloys. A thermodynamically consistent framework for elastic and plastic work dissipation in the fracture process zone is used to formulate the coupled constitutive behavior for brittle and ductile damage, respectively. Representative Volume Elements (RVE) with varying particle morphology and orientation were generated and their uni-axial loading was simulated to assess the damage resistance of the different model microstructures. For critical energy release rate $G_c$ values of the aluminum matrix ranging from 4 to 8~Jm$^{-2}$ for single crystals and from 10 to 16~Jm$^{-2}$ for polycrystals, the model predicts a change in failure modes from particle debonding to cracking followed by ductile matrix failure. The change of particle failure mechanism as a result of increased $G_c$ is observed for both, single and polycrystalline model microstructures. For the case of particle debonding, microstructures with circular or ellipsoidal particles (with the major axis perpendicular to the loading direction) show a higher ductility and fracture work compared to the other studied cases. In the case of particle cracking, microstructures with ellipsoidal particles aligned parallel to the loading axis show a higher ductility and fracture work among the investigated cases. The simulations qualitatively reproduce the experimentally observed particle failure mechanisms (cracking and debonding) for two different matrix alloy classes (commercially pure and 2xxx alloys) reinforced with ceramic particles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Understanding the degradation of a model Si-anode in Li-ion battery at the atomic-scale

Si-anodes have long been candidates thanks to an expected ten-fold increase in capacity compared to graphite. However, details of the mechanisms governing their degradation remain elusive, hindering science-guided development of long-lived Si-based anodes. Here we demonstrate how the latest developments in cryo-atom probe tomography enable the in-depth analysis of the electrode and electrolyte, and their interface at atomic-level.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Influence of microstructure and atomic-scale chemistry on iron ore reduction with hydrogen at 700°C

With 1.85 billion tons produced per year, steel is the most important material class in terms of volume and environmental impact. While steel is a sustainability enabler, for instance through lightweight design, magnetic devices, and efficient turbines, its primary production is not. For 3000 years, iron has been reduced from ores using carbon. Today 2.1 tons CO2 are produced per ton of steel, causing 30% of the global CO2 emissions in the manufacturing sector, which translates to 6.5% of the global CO2 emissions. These numbers qualify iron- and steel-making as the largest single industrial greenhouse gas emission source. The envisaged future industrial route to mitigate these CO2 emissions targets green hydrogen as a reductant. Although this reaction has been studied for decades, its kinetics is not well understood, particularly during the wustite reduction step which is dramatically slower than the hematite reduction. Many rate-limiting factors of this reaction are set by the micro- and nanostructure as well as the local chemistry of the ores. Their quantification allows knowledge-driven ore preparation and process optimization to make the hydrogen-based reduction of iron ores commercially viable, enabling the required massive CO2 mitigation to ease global warming. Here, we report on a multi-scale structure and composition analysis of iron reduced from hematite with pure H2, reaching down to near-atomic scale.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Laser-equipped gas reaction chamber for probing environmentally sensitive materials at near atomic scale

Numerous metallurgical and materials science applications depend on quantitative atomic-scale characterizations of environmentally-sensitive materials and their transient states. Studying the effect upon materials subjected to thermochemical treatments in specific gaseous atmospheres is of central importance for specifically studying a material's resistance to certain oxidative or hydrogen environments. It is also important for investigating catalytic materials, direct reduction of an oxide, particular surface science reactions or nanoparticle fabrication routes. This manuscript realizes such experiments upon a thermochemical reaction chamber called the Reacthub and allows for transferring treated materials under cryogenic & ultrahigh vacuum workflow conditions for characterisation by either atom probe or scanning Xe+/electron microscopies. Two examples are discussed in the present study. One protocol was in the deuterium gas charging (25 kPa D2 at 200 °C) of a high-manganese twinning-induced-plasticity steel and characterization of the ingress and trapping of hydrogen at various features (grain boundaries in particular) in efforts to relate this to the steel's hydrogen embrittlement susceptibility. Deuterium was successfully detected after gas charging but most contrast originated from the complex ion FeOD+ signal and the feature may be an artefact. The second example considered the direct deuterium reduction (5 kPa D2 at 700 °C) of a single crystal wuestite sample, demonstrating that under a standard thermochemical treatment causes rapid reduction upon the nanoscale. Further studies are required for complete confidence about these phenomena, but these experiments successfully demonstrate that how an ex-situ thermochemical treatment can be realised that captures environmentally-sensitive transient states that can be analysed by atomic-scale by atom probe microscope.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Magneto-electric Tuning of Pinning-Type Permanent Magnets through Atomic-Scale Engineering of Grain Boundaries

Pinning-type magnets maintaining high coercivity, i.e. the ability to sustain magnetization, at high temperature are at the core of thriving clean-energy technologies. Among these, Sm2Co17-based magnets are excellent candidates owing to their high-temperature stability. However, despite decades of efforts to optimize the intragranular microstructure, the coercivity currently only reaches 20~30% of the theoretical limits. Here, the roles of the grain-interior nanostructure and the grain boundaries in controlling coercivity are disentangled by an emerging magneto-electric approach. Through hydrogen charging/discharging by applying voltages of only ~ 1 V, the coercivity is reversibly tuned by an unprecedented value of ~ 1.3 T. In situ magneto-structural measurements and atomic-scale tracking of hydrogen atoms reveal that the segregation of hydrogen atoms at the grain boundaries, rather than the change of the crystal structure, dominates the reversible and substantial change of coercivity. Hydrogen lowers the local magnetocrystalline anisotropy and facilitates the magnetization reversal starting from the grain boundaries. Our study reveals the previously neglected critical role of grain boundaries in the conventional magnetisation-switching paradigm, suggesting a critical reconsideration of strategies to overcome the coercivity limits in permanent magnets, via for instance atomic-scale grain boundary engineering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparative study of hydrogen embrittlement resistance between additively and conventionally manufactured austenitic stainless steels

Hydrogen embrittlement in 304L (18wt.% Cr, 8-10wt.% Ni) austenitic stainless steel (ASS) fabricated by laser powder-bed-fusion (LPBF) was investigated by tensile testing after electrochemical hydrogen pre-charging and compared to conventionally available 304L ASSs with two different processing histories, (i) casting plus annealing (CA) and (ii) CA plus thermomechanical treatment (TMT). It was revealed that hydrogen-charging led to a significant reduction in ductility for the CA sample, but only a small effect of hydrogen was observed for the LPBF and CA-TMT samples. Hydrogen-assisted cracking behavior was found to be strongly linked to strain-induced martensitic transformation. In addition, the amount of alpha' martensite was much higher in the CA sample than in other samples, suggesting that severe hydrogen embrittlement can be correlated with the low mechanical stability of austenite. Detailed microstructural characterization showed that low austenite stability of the CA sample was mainly attributed to the retained content of delta ferrite and the chemical inhomogeneity inside the gamma matrix (gamma close to delta has ~2 wt.% higher Cr but ~2 wt.% lower Ni), but TMT enhanced the chemical homogeneity and thus austenite stability. By contrast, the LPBF process led directly, i.e. without any thermomechanical treatment, to a fully austenitic structure with homogeneous elemental distribution in the ASS. These results confirmed that the presence of delta and the chemical inhomogeneity inside gamma matrix, which promoted the deformation-induced martensitic transformation and the associated H enrichment at the gamma-alpha' interface, was the primary reason for the severe H-assisted failure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solute hydrogen and deuterium observed at the near atomic scale in high-strength steel

Observing solute hydrogen (H) in matter is a formidable challenge, yet, enabling quantitative imaging of H at the atomic-scale is critical to understand its deleterious influence on the mechanical strength of many metallic alloys that has resulted in many catastrophic failures of engineering parts and structures. Here, we report on the APT analysis of hydrogen (H) and deuterium (D) within the nanostructure of an ultra-high strength steel with high resistance to hydrogen embrittlement. Cold drawn, severely deformed pearlitic steel wires (Fe-0.98C-0.31Mn-0.20Si-0.20Cr-0.01Cu-0.006P-0.007S wt.%, ε=3.1) contains cementite decomposed during the pre-deformation of the alloy and ferrite. We find H and D within the decomposed cementite, and at some interfaces with the surrounding ferrite. To ascertain the origin of the H/D signal obtained in APT, we explored a series of experimental workflows including cryogenic specimen preparation and cryogenic-vacuum transfer from the preparation into a state-of-the-art atom probe. Our study points to the critical role of the preparation, i.e. the possible saturation of H-trapping sites during electrochemical polishing, how these can be alleviated by the use of an outgassing treatment, cryogenic preparation and transfer prior to charging. Accommodation of large amounts of H in the under-stoichiometric carbide likely explains the resistance of pearlite against hydrogen embrittlement.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Multi-component chemo-mechanics based on transport relations for the chemical potential

A chemo-mechanical model for a finite-strain elasto-viscoplastic material containing multiple chemical components is formulated and an efficient numerical implementation is developed to solve the resulting transport relations. The numerical solution relies on inverting the constitutive model for the chemical potential. In this work, a semi-analytical inversion for a general family of multi-component regular-solution chemical free energy models is derived. This is based on splitting the chemical free energy into a convex contribution, treated implicitly, and a non-convex contribution, treated explicitly. This results in a reformulation of the system transport equations in terms of the chemical potential rather than the composition as the independent field variable. The numerical conditioning of the reformulated system, discretised by finite elements, is shown to be significantly improved, and convergence to the Cahn-Hilliard solution is demonstrated for the case of binary spinodal decomposition. Chemo-mechanically coupled binary and ternary spinodal decomposition systems are then investigated to illustrate the effect of anisotropic elastic deformation and plastic relaxation of the resulting spinodal morphologies in more complex material systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dislocation Activities at the Martensite Phase Transformation Interface in Metastable Austenitic Stainless Steel: An In-situ TEM Study

Understanding the mechanism of martensitic transformation is of great importance in developing advanced high strength steels, especially TRansformation-Induced Plasticity (TRIP) steels. The TRIP effect leads to enhanced work-hardening rate, postponed onset of necking and excellent formability. In-situ transmission electron microscopy has been performed to systematically investigate the dynamic interactions between dislocations and alpha martensite at microscale. Local stress concentrations, e.g. from notches or dislocation pile-ups, render free edges and grain boundaries favorable nucleation sites for alpha martensite. Its growth leads to partial dislocation emission on two independent slip planes from the hetero-interface when the austenite matrix is initially free of dislocations. The kinematic analysis reveals that activating slip systems on two independent {111} planes of austenite are necessary in accommodating the interfacial mismatch strain. Full dislocation emission is generally observed inside of austenite regions that contain high density of dislocations. In both situations, phase boundary propagation generates large amounts of dislocations entering into the matrix, which renders the total deformation compatible and provide substantial strain hardening of the host phase. These moving dislocation sources enable plastic relaxation and prevent local damage accumulation by intense slipping on the softer side of the interfacial region. Thus, finely dispersed martensite distribution renders plastic deformation more uniform throughout the austenitic matrix, which explains the exceptional combination of strength and ductility of TRIP steels.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Aging in Structural Changes of Amorphous Solids: A Study of First Passage Time and Persistence Time Distribution

The time distribution of relaxation events in an aging system is investigated via molecular dynamics simulations. The focus is on the distribution functions of the first passage time, $p_1(Δt)$, and the persistence time, $p(τ)$. In contrast to previous reports, both $p_1$ and $p$ are found to evolve with time upon aging. The age dependence of the persistence time distribution is shown to be sensitive to the details of the algorithm used to extract it from particle trajectories. By updating the reference point in event detection algorithm and accounting for the event specific aging time, we uncover age dependence of $p(τ)$, hidden to previous studies. Moreover, the apparent age-dependence of $p_1$ in continuous time random walk with an age independent $p(τ)$ is shown to result from an implicit synchronization of all the random walkers at the starting time.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Unraveling the temperature dependence of the yield strength in single-crystal tungsten using atomistically-informed crystal plasticity calculations

We use a physically-based crystal plasticity model to predict the yield strength of body-centered cubic (bcc) tungsten single crystals subjected to uniaxial loading. Our model captures the thermally-activated character of screw dislocation motion and full non-Schmid effects, both of which are known to play a critical role in bcc plasticity. The model uses atomistic calculations as the sole source of constitutive information, with no parameter fitting of any kind to experimental data. Our results are in excellent agreement with experimental measurements of the yield stress as a function of temperature for a number of loading orientations. The validated methodology is then employed to calculate the temperature and strain-rate dependence of the yield strength for 231 crystallographic orientations within the standard stereographic triangle. We extract the strain-rate sensitivity of W crystals at different temperatures, and finish with the calculation of yield surfaces under biaxial loading conditions that can be used to define effective yield criteria for engineering design models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chapman-Enskog Analysis of Finite Volume Lattice Boltzmann Schemes

In this paper, we provide a systematic analysis of some finite volume lattice Boltzmann schemes in two dimensions. A complete iteration cycle in time evolution of discretized distribution functions is formally divided into collision and propagation (streaming) steps. Considering mass and momentum conserving properties of the collision step, it becomes obvious that changes in the momentum of finite volume cells is just due to the propagation step. Details of the propagation step are discussed for different approximate schemes for the evaluation of fluxes at the boundaries of the finite volume cells. Moreover, a full Chapman-Enskog analysis is conducted allowing to recover the Navier-Stokes equation. As an important result of this analysis, the relation between the lattice Boltzmann relaxation time and the kinematic viscosity of the fluid is derived for each approximate flux evaluation scheme. In particular, it is found that the constant upwind scheme leads to a positive numerical viscosity while the central scheme as well as the linear upwind scheme are free of this artifact.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Correlations of plasticity in sheared glasses

In a recent paper [S. Mandal et al., Phys. Rev. E 88, 022129 (2013)] the nature of spatial correlations of plasticity in hard sphere glasses was addressed both via computer simulations and in experiments. It was found that the experimentally obtained correlations obey a power law whereas the correlations from simulations are better fitted by an exponential decay. We here provide direct evidence--- via simulations of a hard sphere glass in 2D---that this discrepancy is a consequence of the finite system size in the 3D simulations. By extending the study to a 2D soft disk model at zero temperature, the robustness of the power-law decay in sheared amorphous solids is underlined. Deviations from a power law occur when either reducing the packing fraction towards the supercooled regime in the case of hard spheres or changing the dissipation mechanism from contact dissipation to a mean-field type drag for the case of soft disks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the spheroidized carbide dissolution and elemental partitioning in a high carbon bearing steel 100Cr6

We report on the characterization of high carbon bearing steel 100Cr6 using electron microscopy and atom probe tomography in combination with multi-component diffusion simulations (DICTRA). Scanning electron micrographs show that around 14 vol.% spheroidized carbides are formed during soft annealing and only 3 vol.% remain after dissolution into the austenitic matrix by austenitization at 1123 K (850 °C) for 300 s. The spheroidized particles are identified as (Fe, Cr)3C by transmission electron microscopy. Atom probe analyses reveal the redistribution and partitioning behaviors of elements, i.e. C, Si, Mn, Cr, Fe in both, the spheroidized carbides and the bainitic matrix in the sample isothermally heat-treated at 773 K (500 °C) after austenitization. A homogeneous distribution of C and gradual gradient of Cr was detected within the spheroidized carbides. Due to its limited diffusivity in (Fe, Cr)3C, Cr exhibits a maximum concentration at the surface of spheroidized carbides (16 at.%) and decreases gradually from surface towards the core down to a level of about 2 at.%. The atom probe results also indicate that the partially dissolved spheroidized carbides during austenitization may serve as nucleation sites for intermediate temperature cementite within bainite, which results in a relatively softer surface and harder core in spheroidized particles. This microstructure may contribute to the good wear resistance and fatigue properties

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Crossover from tumbling to tank-treading-like motion in dense simulated suspensions of red blood cells

Via computer simulations, we provide evidence that the shear rate induced red blood cell tumbling-to-tank-treading transition also occurs at quite high volume fractions, where collective effects are important. The transition takes place as the ratio of suspension stress to the characteristic cell membrane stress exceeds a certain value, independent of volume fraction and cell deformability. This value coincides with that for a transition from an orientationally less ordered to a highly ordered phase. The average cell deformation does not show any signature of the transition, but rather follows a simple scaling law independent of volume fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Maintaining the equipartition theorem in small heterogeneous molecular dynamics ensembles

It has been reported recently that the equipartition theorem is violated in molecular dynamics simulations with periodic boundary condition [Shirts et al, J. Chem. Phys. 125, 164102 (2006)]. This effect is associated with the conservation of the center of mass momentum. Here, we propose a fluctuating center of mass molecular dynamics approach (FCMMD) to solve this problem. Using the analogy to a system exchanging momentum with its surroundings, we work out --and validate via simulations-- an expression for the rate at which fluctuations shall be added to the system. The restoration of equipartition within the FCMMD is then shown both at equilibrium as well as beyond equilibrium in the linear response regime.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Polycrystal model of the mechanical behavior of a Mo-TiC30vol.% metal-ceramic composite using a 3D microstructure map obtained by a dual beam FIB-SEM

The mechanical behavior of a Mo-TiC30 vol.% ceramic-metal composite was investigated over a large temperature range (25^{\circ}C to 700^{\circ}C). High-energy X-ray tomography was used to reveal the percolation of the hard titanium carbide phase through the composite. Using a polycrystal approach for a two-phase material, finite element simulations were performed on a real 3D aggregate of the material. The 3D microstructure, used as starting configuration for the predictions, was obtained by serial-sectioning in a dual beam Focused Ion Beam (FIB)-Scanning Electron Microscope (SEM) coupled to an Electron Back Scattering Diffraction system (3D EBSD, EBSD tomography). The 3D aggregate consists of a molybdenum matrix and a percolating TiC skeleton. As most BCC metals, the molybdenum matrix phase is characterized by a change in the plasticity mechanisms with temperature. We used a polycrystal model for the BCC material, which was extended to two phases (TiC and Mo). The model parameters of the matrix were determined from experiments on pure molydenum. For all temperatures investigated, the TiC particles were considered as brittle. Gradual damage of the TiC particles was treated, based on an accumulative failure law that is approximated by an evolution of the apparent particle elastic stiffness. The model enabled us to determine the evolution of the local mechanical fields with deformation and temperature. We showed that a 3D aggregate representing the actual microstructure of the composite is required to understand the local and global mechanical properties of the studied composite.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Single particle fluctuations and directional correlations in driven hard sphere glasses

Via event driven molecular dynamics simulations and experiments, we study the packing fraction and shear-rate dependence of single particle fluctuations and dynamic correlations in hard sphere glasses under shear. At packing fractions above the glass transition, correlations increase as shear rate decreases: the exponential tail in the distribution of single particle jumps broadens and dynamic four-point correlations increase. Interestingly, however, upon decreasing the packing fraction, a broadening of the exponential tail is also observed, while dynamic heterogeneity is shown to decrease. An explanation for this behavior is proposed in terms of a competition between shear and thermal fluctuations. Building upon our previous studies [Chikkadi et al, Europhys. Lett. (2012)], we further address the issue of anisotropy of the dynamic correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Heterogeneous shear in hard sphere glasses

There is growing evidence that the flow of driven amorphous solids is not homogeneous, even if the macroscopic stress is constant across the system. Via event driven molecular dynamics simulations of a hard sphere glass, we provide the first direct evidence for a correlation between the fluctuations of the local volume-fraction and the fluctuations of the local shear rate. Higher shear rates do preferentially occur at regions of lower density and vice versa. The temporal behavior of fluctuations is governed by a characteristic time scale, which, when measured in units of strain, is independent of shear rate in the investigated range. Interestingly, the correlation volume is also roughly constant for the same range of shear rates. A possible connection between these two observations is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method

The deformation of an initially spherical capsule, freely suspended in simple shear flow, can be computed analytically in the limit of small deformations [D. Barthes-Biesel, J. M. Rallison, The Time-Dependent Deformation of a Capsule Freely Suspended in a Linear Shear Flow, J. Fluid Mech. 113 (1981) 251-267]. Those analytic approximations are used to study the influence of the mesh tessellation method, the spatial resolution, and the discrete delta function of the immersed boundary method on the numerical results obtained by a coupled immersed boundary lattice Boltzmann finite element method. For the description of the capsule membrane, a finite element method and the Skalak constitutive model [R. Skalak et al., Strain Energy Function of Red Blood Cell Membranes, Biophys. J. 13 (1973) 245-264] have been employed. Our primary goal is the investigation of the presented model for small resolutions to provide a sound basis for efficient but accurate simulations of multiple deformable particles immersed in a fluid. We come to the conclusion that details of the membrane mesh, as tessellation method and resolution, play only a minor role. The hydrodynamic resolution, i.e., the width of the discrete delta function, can significantly influence the accuracy of the simulations. The discretization of the delta function introduces an artificial length scale, which effectively changes the radius and the deformability of the capsule. We discuss possibilities of reducing the computing time of simulations of deformable objects immersed in a fluid while maintaining high accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fall and rise of small droplets on rough hydrophobic substrates

Liquid droplets on patterned hydrophobic substrates are typically observed either in the Wenzel or the Cassie state. Here we show that for droplets of comparable size to the roughness scale an additional local equilibrium state exists, where the droplet is immersed in the texture, but not yet contacts the bottom grooves. Upon evaporation, a droplet in this state enters the Cassie state, leading to a qualitatively new self-cleaning mechanism. The effect is of generic character and is expected to occur in any hydrophobic capillary wetting situation where a spherical liquid reservoir is involved.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Inertia effects and stress accumulation in a constricted duct: A combined experimental and lattice Boltzmann study

We experimentally and numerically investigate the flow of a Newtonian fluid through a constricted geometry for Reynolds numbers in the range $0.1 - 100$. The major aim is to study non-linear inertia effects at larger Reynolds numbers (>10) on the shear stress evolution in the fluid. This is of particular importance for blood flow as some biophysical processes in blood are sensitive to shear stresses, e.g., the initialization of blood clotting. We employ the lattice Boltzmann method for the simulations. The conclusion of the predictions is that the peak value of shear stress in the constriction grows disproportionally fast with the Reynolds number which leads to a non-linear shear stress accumulation. As a consequence, the combination of constricted blood vessel geometries and large Reynolds numbers may increase the risk of undesired blood clotting.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Shear stress in lattice Boltzmann simulations

A thorough study of shear stress within the lattice Boltzmann method is provided. Via standard multiscale Chapman-Enskog expansion we investigate the dependence of the error in shear stress on grid resolution showing that the shear stress obtained by the lattice Boltzmann method is second order accurate. This convergence, however, is usually spoiled by the boundary conditions. It is also investigated which value of the relaxation parameter minimizes the error. Furthermore, for simulations using velocity boundary conditions, an artificial mass increase is often observed. This is a consequence of the compressibility of the lattice Boltzmann fluid. We investigate this issue and derive an analytic expression for the time-dependence of the fluid density in terms of the Reynolds number, Mach number and a geometric factor for the case of a Poiseuille flow through a rectangular channel in three dimensions. Comparison of the analytic expression with results of lattice Boltzmann simulations shows excellent agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Profile blunting and flow blockage in a yield stress fluid: A molecular dynamics study

The flow of a simple glass forming system (a 80:20 binary Lennard-Jones mixture) through a planar channel is studied via molecular dynamics simulations. The flow is driven by an external body force similar to gravity. Previous studies show that the model exhibits both a static [Varnik et al. J. Chem. Phys. 120, 2788 (2004)] and a dynamic [F. Varnik and O. Henrich Phys. Rev. B 73, 174209 (2006)] yield stress in the glassy phase. \blue{These observations are corroborated by the present work, where we investigate how the presence of a yield stress may affect the system behavior in a Poiseuille-type flow geometry.} In particular, we observe a blunted velocity profile across the channel: A relatively wide region in the channel center flows with a constant velocity (zero shear rate) followed by a non linear change of the shear rate as the walls are approached. The observed velocity gradients are compared to those obtained from the knowledge of the shear stress across the channel and the flow-curves (stress versus shear rate), the latter being determined in our previous simulations of homogeneous shear flow. Furthermore, using the value of the (dynamic) yield stress known from previous simulations, we estimate the threshold body force for a complete arrest of the flow. Indeed, a blockage is observed as the imposed force falls below this threshold value. Small but finite shear rates are observed at stresses above the dynamic but below the static yield stress. We discuss the possible role of the \blue{stick-slip like motion} for this observation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Wetting gradient induced separation of emulsions: A combined experimental and lattice Boltzmann computer simulation study

Guided motion of emulsions is studied via combined experimental and theoretical investigations. The focus of the work is on basic issues related to driving forces generated via a step-wise (abrupt) change in wetting properties of the substrate along a given spatial direction. Experiments on binary emulsions unambiguously show that selective wettability of the one of the fluid components (water in our experiments) with respect to the two different parts of the substrate is sufficient in order to drive the separation process. These studies are accompanied by approximate analytic arguments as well as lattice Boltzmann computer simulations, focusing on effects of a wetting gradient on internal droplet dynamics as well as its relative strength compared to volumetric forces driving the fluid flow. These theoretical investigations show qualitatively different dependence of wetting gradient induced forces on contact angle and liquid volume in the case of an open substrate as opposed to a planar channel. In particular, for the parameter range of our experiments, slit geometry is found to give rise to considerably higher separation forces as compared to open substrate.

0 citations0 reviewsNo code linkedOpen access