NANO- SCALE SOLIDIFICATION OF ALLOYS Justyna Czerwinska, Erwan Deriaz, Agnieszka Słowicka, Tomasz A. Kowalewski Department of Mechanics and Physics of Fluids Institute of Fundamental Technological Research, Polish Academy of Sciences MOTIVATION EU PROJECT NADIA - DESIGN AND MANUFACTURING OF LIGHT ALLOYS BY INTELLIGENT PROCESSING
MOLECULAR DYNAMICS SIMULATIONS
A CLOSE STUDY OF ALUMINUM ALLOYS
Assumptions: σ 1 2 σ 6 U ( r ) = ε − r r ε, σ - parameters of the Lennard- Jones potential r – distance between molecules
Al- 21.5%Si obs er ved cryst als
• • • •
spherical non- polar particles interactions: Lennard- Jones potential only neighbour particles interact calculations: program LAMMPS Sandia National Laboratories
by
σ (from HUT - NADIA rep ort Ap ril 20 0 8)
To p red ict alloy m ech an ical p rop ert ies an d t o avoid exp en s ive exp erim en t al t est s we em p loy n u m erical sim u lat ion s b y m olecu lar d yn am ics. Firs t , n u m er ical exp er im en t s are p erform ed t o obt ain t h e p h en om en ological m et al liq u id p r op ert ies (t ran s p ort coefficien t s ), t h e s olid ificat ion ch aract erist ics (p h as e d iagram , crys t alliz at ion fr on t s p eed ). Fin ally, t h e m ech an ical p rop er t ies of t h e res u lt in g alloy will b e evalu at ed .
MODELING CRYSTALLIZATION FCC crystal (Al) cs t - www.n rl.n avy.m il/ lat t ice
T= 2 0 0 0 K Len n a r d - Jon es cr ys t a l for m s FCC s t r u ct u r es
T= 1 0 8 9 K
T= 9 1 2 K
T= 6 8 4 K
Molecu la r d yn a m ics cr ys t a lliz a t ion s im u la t ion wit h Len n a r d - Jon es p ot en t ia l, liq u id a t om s a r e colou r ed in b lu e a n d cr ys t a lliz ed a t om s in gr een
Gr een - Ku bo es tim ation of d iffu s ion in liq u id u s ed t o valid at e t h e m od el (tem p er atu r e u n its ar e ar bit rar y)
MECHANICAL PROPERTIES Swit ch in g t o Em b ed d ed At om Mod el, we t es t t h e m ech an ical p r op ert ies of t h e s olid ified m et als we ob t ain ed .
Mea s u r ed en er gy d is s ip a t ion d u r in g s olid ifica t ion a n d m elt in g of a Al- 1 7 %Si- 4 %Cu + 3 0 %Al2 O3 a lloy [HUT NADIA r ep or t ].
En er gy d is s ip a t ion d u r in g t h e coolin g of a n h om ogen eou s Len n a r d - Jo n es m a t er ia l in a m olecu la r d yn a m ics n u m er ica l s im u la t ion .
MODELING METAL LIQUID chemical element Al Cu Mg Sn Si We m od el t h e t h e m os t com m on atomic mass (g/mol) 27 63.6 24.3 118.7 28.1 elem en t s u sed t o creat e alu m in u m alloys : cop p er (Cu ), m agn es iu m (Mg) 2390 7990 1560 7010 2100 density (kg/m3) an d silicon (Si). We also ad d t in (Sn ) Fusion temperature (K) 933 1358 923 505 1687 t o t est ou r m od el. Al Cu Mg Sn Si We d et erm in e exp erim en t ally (wit h M- D LJ coef ε 2.7480 4.7150 2.0270 0.5040 5.1720 s im u lat ion s ) t h e Len n ard - Jon es coefficien t s ε LJ coef σ 1.1336 1.0023 1.2977 1.4554 1.2105 an d σ cor resp on d in g t o t h es e m et als . Density ρ 0.7654 1.0726 0.5435 0.4679 0.6304
Mat erial s t ren gt h t es t d evice
MODELING ALLOYS Mat erial u n d ergoin g a s t r es s at t h e n an o- s cale Us in g m ixin g ru les: = 1 2
= 1 2 /2
we sim u lat ed t h e solid ificat ion of a Al- 17%Si- 4%Cu alloy – figu re on t h e righ t – an d d et ect ed t h e liq u id u s an d t h e s olid u s u s in g a geom et rical m et h od . Th e res u lt s obt ain ed are in agreem en t wit h t h e exp erim en t al resu lt s. Fou r com p on en t Len n ard - Jon es cr ys t al
Th e s im u lat ion s h ows a lin ear resp on s e t o s t rain s t res s for t es t ed m at erial Pe rs pe ctiv e s : - r ecreat e accu rat ely t h e cryst al s t ru ct u r e of alloys , - t es t for in n ovat ive m et al m ixt u res an d alloyin g p r oces ses , - fin d in t erest in g alloys recip es.
Dec 16, 2008 - â Centrale Marseille/LATP, Technopôle de Château-Gombert, 38 rue F. ... differentiation of wavelet bases that shows that the above families of ...
Feb 16, 2006 - that forms a basis (sometime orthogonal) of a functionnal .... 8. B. #йи is an orthogonal basis change (orthogonal matrix). ..... Scale separation.
the example of the merging of three vortices in dimension 2. ..... condition: δt ⤠Cδx2. This study is inspired by the approach of Temam in his book ... (3.5). = 1 δt. [. PJ. (â« (n+1)δt nδt. âtu dt. )] â (νÎun+1 â P [(un · â)un]),
Jul 26, 2013 - It presents an axial symmetry along (Ox). To be consistent ..... p=1, q=5, alpha=9/5 ... with skew-symmetric discretizations, SIAM J. Numer. Anal.
We present an original adaptive scheme using a dynamically refined grid for the ...... physical space while in the gravitational case it collapses, allowing the ...
Oct 25, 2010 - finally we show how this stability criterion extends to non-linear equations, ...... Pδx the orthogonal projector onto the space of discretization. ..... We extend the scope of application of the stability conditions (3.9) to other ca
Apr 17, 2008 - âInstitute of Mathematics, Polish Academy of Sciences. ul. ..... In the vector case, the algorithm should be applied on each component of the vector ...... For higher thresholdings, the effects of the anisotropy begin to be visible.
Engineering and Materials Science, University of Minnesota, Minneapolis, ... lens as a function of P, the load applied. This technique gives the dependence of the interfacial adhesion energy, or strain energy release rate, G, on the crack propagation
The first and the third summands of the. Lagrangian describe the free motion of electron, the last summand describes the same of photons, and the term Ïϵ. A.
This behavior depends on both the solid fraction and the shear rate. ... Keywords: Casting; Semi-solid alloys; Theory & modeling. 1. .... The driving force of.
Abstract. In this paper we investigate the use of divergence-free wavelet bases for the Coherent. Vortex Extraction (CVE) of turbulent flows. We begin with a short ...
adsorbate. The adsorption process generally involves chemisorption in which the gas ... diffusion coefficient is large when the molecules mix rapidly. The first ...
Mar 1, 2010 - copy is furnished to the author for internal non-commercial research ... First applications of the Craya decomposition to the dynamics of turbulent fluctua- .... The motivation is to dispose of an efficient algorithm for local CrayaâH
Jun 20, 2005 - (NFLV) is based on total internal reflection-fluorescence recovery after photobleaching ... ing a high surface energy ( s > 72 mJ=m2). Rms rough- .... layers could be the origin of a better momentum transfer and thus of a larger ...
These IBMs have focused on early life his- ... off the west coast of southern Africa is a highly dynamic system. Its main ...... In: Oceanography and Marine Biology.
Dec 1, 2002 - substrate, as measured through atomic-force microscopy techniques. .... the glass transition temperature Tg of the latex, while the elastic ...
Jul 15, 2002 - The oxidation kinetics of a non planar silicon surface were already of great .... kinetic effects by computing local stresses in the sphere. Our results are ...... 1 D.-B. Kao, J. McVitte, W. Nix, and K. Saraswat, IEEE Trans. Electron.
deposited material remained at the electrode surface, even following .... (EG&G, Princeton Applied Research) or a multipotentiostat .... core level spectra were referenced to the Au 4f7/2 binding energy ...... the Internet at http://pubs.acs.org.
Circassians at the service of the king of Jordan. They are Cherkessians or Chechnyans, originate from the Caucasus and have been guarding the safety of ...