Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large Crystal”

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 577 records · Page 32

System and method for liquid crystal display system incorporating wire grid polarizers for large scale and large volume stereolithography

The present disclosure relates to a liquid crystal display (LCD) system. The system in one example has a light source for generating unpolarized light, and an LCD screen arranged in a path of transmittance of the unpolarized light. First and second wire grid polarizers are arranged adjacent to the LCD screen and each have a plurality of nano-scale wires, with the first and second wire grid polarizers have differing polarizations. A pitch of each of the nano-scale wires is no larger than one-third a wavelength of the unpolarized light from the light source. The wire grid polarizers create, in connection with operation of the LCD screen, a 2D light mask suitable for initiating the polymerization of an optically curable material.

Duoss, Eric B.↗

Potential productivity benefits of float-zone versus Czochralski crystal growth

Efficient mass production of single-crystal silicon is necessary for the efficient silicon solar arrays needed in the coming decade. However, it is anticipated that there will be difficulty growing such volumes of crystals using conventional Czochralski (Cz) methods. While the productivity of single crystals might increase with a crystal diameter increase, there are two obstacles to the mass production of large diameter Czochralski crystals, the long production cycle due to slow growth rate and the high heat requirements of the furnaces. Also counterproductive would be the large resistivity gradient along the growth direction of the crystals due to impurity concentration. Comparison between Float zone (FZ) and Cz crystal growth on the basis of a crystal 150 mm in diameter is on an order of two to four times in favor of the FZ method. This advantage results from high growth rates and steady-state growth while maintaining a dislocation-free condition and impurity segregation.

Abe, T.↗

Process for large-scale ammonothermal manufacturing of semipolar gallium nitride boules

Methods for large-scale manufacturing of semipolar gallium nitride boules are disclosed. The disclosed methods comprise suspending large-area single crystal seed plates in a rack, placing the rack in a large diameter autoclave or internally-heated high pressure apparatus along with ammonia and a mineralizer, and growing crystals ammonothermally. A bi-faceted growth morphology may be maintained to facilitate fabrication of large area semipolar wafers without growing thick boules.

D'Evelyn, Mark P.↗

Process for large-scale ammonothermal manufacturing of semipolar gallium nitride boules

Methods for large-scale manufacturing of semipolar gallium nitride boules are disclosed. The disclosed methods comprise suspending large-area single crystal seed plates in a rack, placing the rack in a large diameter autoclave or internally-heated high pressure apparatus along with ammonia and a mineralizer, and growing crystals ammonothermally. A bi-faceted growth morphology may be maintained to facilitate fabrication of large area semipolar wafers without growing thick boules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal Structure of Moganite and Its Anisotropic Atomic Displacement Parameters Determined by Synchrotron X-ray Diffraction and X-ray/Neutron Pair Distribution Function Analyses

The crystal structure of moganite from the Mogán formation on Gran Canaria has been re-investigated using high-resolution synchrotron X-ray diffraction (XRD) and X-ray/neutron pair distribution function (PDF) analyses. Our study for the first time reports the anisotropic atomic displacement parameters (ADPs) of a natural moganite. Rietveld analysis of synchrotron XRD data determined the crystal structure of moganite with the space group I2/a. The refined unit-cell parameters are a = 8.7363(8), b = 4.8688(5), c = 10.7203(9) Å, and β = 90.212(4)°. The ADPs of Si and O in moganite were obtained from X-ray and neutron PDF analyses. The shapes and orientations of the anisotropic ellipsoids determined from X-ray and neutron measurements are similar. The anisotropic ellipsoids for O extend along planes perpendicular to the Si-Si axis of corner-sharing SiO4 tetrahedra, suggesting precession-like movement. Neutron PDF result confirms the occurrence of OH over some of the tetrahedral sites. We postulate that moganite nanomineral is stable with respect to quartz in hypersaline water. The ADPs of moganite show a similar trend as those of quartz determined by single-crystal XRD. In short, the combined methods can provide high-quality structural parameters of moganite nanomineral, including its ADPs and extra OH position at the surface. This approach can be used as an alternative means for solving the structures of crystals that are not large enough for single-crystal XRD measurements, such as fine-grained and nanocrystalline minerals formed in various geological environments.

58 GEOSCIENCES↗

Large grained perovskite solar cells derived from single-crystal perovskite powders with enhanced ambient stability

In this study, we demonstrate the large grained perovskite solar cells prepared from precursor solution comprising single-crystal perovskite powders for the first time. Here, the resultant large grained perovskite thin film possesses negligible physical (structural) gap between each large grain and are highly crystalline as evidenced by its fan-shaped birefringence observed under polarized light, which is very different to the thin film prepared from the typical precursor route (MAI + PbI 2 ).

14 SOLAR ENERGY↗

Palladium/Beta Zeolite Passive NOx Adsorbers (PNA): Clarification of PNA Chemistry and the Effects of CO and Zeolite Crystallite Size on PNA Performance

Model 0.3 wt% and 1 wt% Pd/Beta was synthesized with nano-sized Beta crystals (average size <50 nm) and defect-free large micrometer sized crystals with similar Si/Al ratios ~13 and 15. With the aid of FTIR we reveal the formation of at least two Pd(II)-NO and one Pd(I)-NO complexes upon NO adsorption on those materials under conditions relevant to PNA. Under higher NO pressures at room temperature, we detect the formation of the first palladium (II) dinitrosyl Pd(NO)2 complex, not observed previously on the solid supports or in organometallic compounds. O2 in dry streams promotes NO storage via additional formation of nitrosyl (NO+) ions; it also leads to stabilization of significant amounts of N2O3 in BEA pores. In the presence of water, NO storage capacity is dramatically lowered. However, in the presence of CO (which is always present in exhaust gas) the performance improves. With the aid of spectroscopy and PNA measurements we show that in the presence of CO stable mixed palladium(II) carbonyl nitrosyl complex [Pd(II)(CO)(NO)] is formed and responsible for enhanced NOx storage. This points to a general conclusion that the species responsible for NOx storage on various zeolites in the presence of water and CO, are in fact, Pd(II)(CO)(NO). Decreasing the Pd loading to 0.3 wt% improves Pd dispersion and NO/Pd storage ratio. We also demonstrate that varying crystal sizes of BEA leads to significant changes in the NOx release temperature as well as NO/Pd ratio. Larger, defect free crystals store more NOx and release it at higher temperature. Furthermore, larger crystals stay significantly more active for PNA after hydrothermal aging, losing very little activity compared with defective nanocrystals. High-field 27Al NMR results show that larger BEA crystals are more resistant to de-alumination upon hydrothermal aging (HTA), and thus are more attractive for PNA storage with Pd.

Khivantsev, Konstantin↗

Large Scale DD Simulation Results for Crystal Plasticity Parameters in Fe-Cr And Fe-Ni Systems

The development of viable nuclear energy source depends on ensuring structural materials integrity. Structural materials in nuclear reactors will operate in harsh radiation conditions coupled with high level hydrogen and helium production, as well as formation of high density of point defects and defect clusters, and thus will experience severe degradation of mechanical properties. Therefore, the main objective of this work is to develop a capability that predicts aging behavior and in-service lifetime of nuclear reactor components and, thus provide an instrumental tool for tailoring materials design and development for application in future nuclear reactor technologies. Towards this end goal, the long term effort is to develop a physically based multiscale modeling hierarchy, validated and verified, to address outstanding questions regarding the effects of irradiation on materials microstructure and mechanical properties during extended service in the fission and fusion environments. The focus of the current investigation is on modern steels for use in nuclear reactors including high strength ferritic-martensitic steels (Fe-Cr-Ni alloys). The effort is to develop a predicative capability for the influence of irradiation on mechanical behavior. Irradiation hardening is related to structural information crossing different length scales, such as composition, dislocation, and crystal orientation distribution. To predict effective hardening, the influence factors along different length scales should be considered. Therefore, a hierarchical upscaling methodology is implemented in this work in which relevant information is passed between models at three scales, namely, from molecular dynamics to dislocation dynamics to dislocation-based crystal plasticity. The molecular dynamics (MD) was used to predict the dislocation mobility in body centered cubic (bcc) Fe and its Ni and Cr alloys. The results are then passed on to dislocation dynamics to predict the critical resolved shear stress (CRSS) from the evolution of local dislocation and defects. In this report the focus is on the results obtained from large scale dislocation dynamics simulations. The effect of defect density, materials structure was investigated, and evolution laws are obtained. These results will form the bases for the development of evolution and hardening laws for a dislocation-based crystal plasticity framework. The hierarchical upscaling method being developed in this project can provide a guidance tool to evaluate performance of structural materials for next-generation nuclear reactors. Combined with other tools developed in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the models developed will have more impact in improving the reliability of current reactors and affordability of new reactors.

36 MATERIALS SCIENCE↗

Sculpting 2D Crystals via Membrane Contractions before and during Solidification

When phospholipids crystallize within the otherwise fluid membranes of giant unilamellar vesicles, the resulting molecularly thin “2D” solids exhibit great variety in their morphology evolution. For example, within membranes containing moderate amounts of the crystallizing component, crystals grow with a fixed morphology depending on vesicle size. Conversely for membranes containing large amounts of the crystallizing species, we find small compact crystals on vesicles of all sizes. However, on large vesicles, growing crystals sprout flower petals that lengthen progressively. These behaviors result from two combined mechanisms: first, like other 2D solids, the shear rigidity of phospholipid crystals renders them intolerant to morphologies with nonzero Gaussian curvature. As a result and especially at elevated membrane tension, the cost of bending elasticity is reduced at the expense of line energy by the formation of flowers as opposed to compact crystals. Second, the composition-dependent tension rise during cooling relaxes via water permeation of the membrane with a time constant scaling as R2. The amount of crystal formed for a small decrease in temperature determines this composition-dependent increase in stress from thermal contractions versus solidification. Surface Evolver computations were motivated using the predicted tension evolution to develop a processing space that maps to experimental observations for initial and growing crystal morphology. Important variable groups are identified, including a scaled ratio of bending to line energy, a vesicle-size-independent group for membrane contractions, and a time constant for stress relaxation. Though processing stresses ultimately relax, the crystal morphology persists well beyond the processing window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nucleation and crystallization behaviors of nano-crystalline lithium–mica glass–ceramic prepared via sol–gel method

Graphical abstract: The effects of nucleation and crystallization treatments on nano-crystalline lithium–mica glass–ceramic, synthesized by sol–gel technique, were investigated. It was found that MgF{sub 2} crystals act as nuclei centers for the mica crystallization so that a large quantity of mica crystallites was obtained following nucleation process. The crystallization activation energy for both the un-nucleated and nucleated samples was measured as 400.2 and 229.6 kJ mol{sup −1}, respectively. The calculated Avrami exponents demonstrated that the growth mechanism of mica crystallites changes from the needle-like to three-dimensional growth with applying the appropriate nucleation treatment ▪. Highlights: ► Crystallization temperature shifts to 625 from 680 °C following nucleation process. ► Activation energy of crystallization for the nucleated specimen is 229.6 kJ mol{sup −1}. ► Crystallization activation energy for the un-nucleated specimen is 400.2 kJ mol{sup −1}. ► Needle-like growth is predominant growth mechanism for un-nucleated sample. ► Three-dimensional growth is predominant growth mechanism for nucleated sample. -- Abstract: The paper investigates the effects of nucleation and crystallization treatments on nano-crystalline lithium–mica glass–ceramics, taking the composition LiMg{sub 3}AlSi{sub 3(1+x)}O{sub 10+6x}F{sub 2} (x = 0.5) and 8 mass% MgF{sub 2} synthesized by sol–gel technique. Here, X-ray diffraction, thermal analysis and transmission electron microscopy were used to assess the structural evolutions of as-synthesized nano-crystalline lithium–mica glass–ceramics. It was found that MgF{sub 2} crystals perform as nuclei centers for the mica crystallization hence; a large quantity of mica crystallites obtained following the nucleation process at 400 °C for 12 h. For both the un-nucleated and nucleated samples, the crystallization activation energy was measured as 400.2 and 229.6 kJ mol{sup −1}, respectively. Consequently, the calculated Avrami exponents demonstrated that the growth mechanism of mica crystallites, while applying appropriate nucleation treatment, changes from needle-like to three-dimensional growth.

36 MATERIALS SCIENCE↗

Flux growth in a horizontal configuration: An analog to vapor transport growth

Flux growth of single crystals is normally performed in a vertical configuration with an upright refractory container holding the flux melt. At high temperatures, flux dissolves the charge, forming a homogeneous solution before nucleation and growth of crystals takes place under proper supersaturation generated by cooling or evaporating the flux. In this paper, we report flux growth in a horizontal configuration with a temperature gradient along the horizontal axis: a liquid transport growth analogous to the vapor transport technique. In a typical liquid transport growth, the charge is kept at the hot end of the refractory container and the flux melt dissolves the charge and transfers it to the cold end. Once the concentration of charge is above the solubility limit at the cold end, the thermodynamically stable phase nucleates and grows. Compared to the vertical flux growth, the liquid transport growth can provide a large quantity of crystals in a single growth since the charge/flux ratio is not limited by the solubility limit at the growth temperature. This technique is complementary to the vertical flux growth and can be considered when a large amount of crystals is needed but the yield from the conventional vertical flux growth is limited. Finally, we applied this technique to the growth of IrSb 3 , Mo 3 Sb 7 , and MnBi from self-flux, and the growth of FeSe, CrTe 3 , NiPSe 3 , FePSe 3 , CuInP 2 S 6 , RuCl 3 , and OsCl 4 from a halide flux.

36 MATERIALS SCIENCE↗

A low cost and large-scale synthesis of 3D photonic crystal with SP2 lattice symmetry

In this work, a novel lithographic method is proposed to prepare three-dimensional (3D) photonic crystal (PC) that is different from conventional top-down and bottom-up approaches. The method combines a 2D optical mask and off-the-axis double optical exposures to create a desirable 3D PC structure. Since the method uses only two optical exposures of a photo-resist layer, it is inherently a low-cost, high throughput and wafer-scale lithographic method. The method is implemented to make a slanted post 3D PC having the SP2 lattice symmetry. Three types of SP2 3D PC structures were successfully fabricated with a minimum feature size of d=1.5 μm over a large scale of 8x10 mm 2 , without any observable fabrication defects. The SP2 PCs are: (i) SU8 posts in air background, (ii) air pores in CdS background and (iii) Pt coated on SU8 SP2 templates. A spectroscopic study of the SP2 PCs shows select spectral regions of high reflectance, indicating the existence of a photonic stop band. This low-cost and large-scale method could enable broader technological impacts of 3D PC materials in areas such as thermo-photovoltaics and above room-temperature Bose-Einstein Condensation. Furthermore, this off-axis method could lead to the creation of an entirely new class of slanted-rod based photonic crystals, such as topological photonic crystal in 3D.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shock response of He bubbles in single crystal Cu

With large-scale molecular dynamics simulations, we investigate shock response of He nanobubbles in single crystal Cu. For sufficient bubble size or internal pressure, a prismatic dislocation loop may form around a bubble in unshocked Cu. The internal He pressure helps to stabilize the bubble against plastic deformation. However, the prismatic dislocation loops may partially heal but facilitate nucleation of new shear and prismatic dislocation loops. For strong shocks, the internal pressure also impedes internal jetting, while a bubble assists local melting; a high speed jet breaks a He bubble into pieces dispersed among Cu. Near-surface He bubbles may burst and form high velocity ejecta containing atoms and small fragments, while the ejecta velocities do not follow the three-dimensional Maxwell-Boltzmann distributions expected for thermal equilibrium. The biggest fragment size deceases with increasing shock strength. With a decrease in ligament thickness or an increase in He bubble size, the critical shock strength required for bubble bursting decreases, while the velocity range, space extension and average velocity component along the shock direction, increase. Small bubbles are more efficient in mass ejecting. Compared to voids and perfect single crystal Cu, He bubbles have pronounced effects on shock response including bubble/void collapse, Hugoniot elastic limit (HEL), deformation mechanisms, and surface jetting. HEL is the highest for perfect single crystal Cu with the same orientations, followed by He bubbles without pre-existing prismatic dislocation loops, and then voids. Complete void collapse and shear dislocations occur for embedded voids, as opposed to partial collapse, and shear and possibly prismatic dislocations for He bubbles. He bubbles lower the threshhold shock strength for ejecta formation, and increase ejecta velocity and ejected mass.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lead(II) coordination polymers based on rigid-flexible 3,5-bis-oxyacetate-benzoic acid: Structural transition driven by temperature control

Three Pb(II) complexes ([Pb{sub 3}(BOABA){sub 2}(H{sub 2}O)]·H{sub 2}O){sub n} (1), ([Pb{sub 4}(BOABA){sub 2}(µ{sub 4}-O)(H{sub 2}O){sub 2}]·H{sub 2}O){sub n} (2), and [Pb{sub 3}(BOABA){sub 2}(H{sub 2}O)]{sub n} (3) (H{sub 3}BOABA=3,5-bis-oxyacetate-benzoic acid) were obtained under the same reaction systems with different temperatures. Complexes 1 and 2 are two dimensional (2D) networks based on Pb-BOABA chains and Pb{sub 4}(µ{sub 4}-O)(COO){sub 6} SBUs, respectively. Complex 3 presents an interesting three dimensional (3D) framework, was obtained by increasing the reaction temperature. Structural transition of the crystallization products is largely dependent on the reaction temperature. Moreover, the fluorescence properties of complexes 1–3 have been investigated. - Graphical abstract: Three Pb(II) coordination polymers were obtained under the same reaction systems with different temperatures. Both of complexes 1 and 2 are 2D network. 3 presents a 3D framework based on Pb–O–C rods SBUs. The 2D to 3D structures transition between three complexes was achieved successfully by temperature control. - Highlights: • Three Pb(II) complexes were obtained under the same reaction systems with different temperatures. • Structural transition of the crystallization products is largely dependent on the reaction temperature. • The luminescence properties studies reveal that three complexes exhibit yellow fluorescence emission behavior, which might be good candidates for obtaining photoluminescent materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal Growth Scale-Up and Stability of RbSr 2 X 5 :Eu Scintillators

The discovery and development of new scintillation materials support national security needs. In these applications, large volumes of scintillator crystals are needed to achieve efficient screening for contraband. Therefore, one important step in the discovery of new scintillator compositions is testing their feasibility for scale-up and their stability. In this work, high-quality Ø22 mm crystals of two new scintillators RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu were grown via the Vertical Bridgman method, and their scintillation properties were characterized. Here, the Ø22 mm RbSr 2 I 5 :Eu crystals could be grown with fast translation rates up to 3.5 mm/h. Both RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu had high scintillation performance, including light yields of 46,000 and 61,000 ph/MeV, respectively, for Ø22 × 35 mm crystals. Additionally, properties related to physical stability were investigated, including the coefficients of thermal expansion via high-temperature X-ray diffraction (HTXRD) as well as moisture sensitivity. HTXRD confirmed the absence of solid–solid phase transitions and showed that RbSr 2 Br 5 had minimal thermal expansion anisotropy compared to RbSr 2 I 5 and some other inorganic metal halide scintillators, which favors the growth of large-sized crystals.

36 MATERIALS SCIENCE↗

A new parameterization of ice heterogeneous nucleation coupled to aerosol chemistry in WRF-Chem model version 3.5.1: evaluation through ISDAC measurements

In the Arctic, during polar night and early spring, ice clouds are separated into two leading types of ice clouds (TICs): (1) TIC1 clouds characterized by a large concentration of very small crystals and TIC2 clouds characterized by a low concentration of large ice crystals. Using a suitable parameterization of heterogeneous ice nucleation is essential for properly representing ice clouds in meteorological and climate models and subsequently understanding their interactions with aerosols and radiation. Here, we describe a new parameterization for ice crystal formation by heterogeneous nucleation in water-subsaturated conditions coupled to aerosol chemistry in the Weather Research and Forecasting model coupled with chemistry (WRF-Chem). The parameterization is implemented in the Milbrandt and Yau (2005a, b) two-moment cloud microphysics scheme, and we assess how the WRF-Chem model responds to the run-time interaction between chemistry and the new parameterization. Well-documented reference cases provided us with in situ data from the spring 2008 Indirect and Semi-Direct Aerosol Campaign (ISDAC) over Alaska. Our analysis reveals that the new parameterization clearly improves the representation of the ice water content (IWC) in polluted or unpolluted air masses and shows the poor performance of the reference parameterization in representing ice clouds with low IWC. The new parameterization is able to represent TIC1 and TIC2 microphysical characteristics at the top of the clouds, where heterogenous ice nucleation is most likely occurring, even with the known bias of simulated aerosols by WRF-Chem over the Arctic.

54 ENVIRONMENTAL SCIENCES↗

Colloids and Nucleation

The objectives of the work funded under this grant were to develop a microphotographic technique and use it to monitor the nucleation and growth of crystals of hard colloidal spheres. Special attention is given to the possible need for microgravity studies in future experiments. A number of persons have been involved in this work. A masters student, Keith Davis, began the project and developed a sheet illumination apparatus and an image processing system for detection and analysis. His work on a segmentation program for image processing was sufficient for his master's research and has been published. A post doctoral student Bernie Olivier and a graduate student Yueming He, who originally suggested the sheet illumination, were funded by another source but along with Keith made photographic series of several samples (that had been made by Keith Davis). Data extraction has been done by Keith, Bernie, Yueming and two undergraduates employed on the grant. Results are published in Langmuir. These results describe the sheet lighting technique as one which illuminates not only the Bragg scattering crystal, but all the crystals. Thus, accurate crystal counts can be made for nucleation rate measurements. The strange crystal length scale reduction, observed in small angle light scattering (SALS) studies, following the initial nucleation and growth period, has been observed directly. The Bragg scattering (and dark) crystal size decreases in the crossover region. This could be an effect due to gravitational forces or due to over- compression of the crystal during growth. Direct observations indicate a complex morphology for the resulting hard sphere crystals. The crystal edges are fairly sharp but the crystals have a large degree of internal structure. This structure is a result of (unstable) growth and not aggregation. As yet unpublished work compares growth exponents data with data obtained by SALS. The nucleation rate density is determined over a broad volume fraction range (0.505 to 0.575) and compared to classical nucleation rate calculations. The experimental nucleation rate density varies with volume fraction more slowly than expected theoretically. Gravity may influence nucleation because crystal nuclei are more dense than the surrounding fluid and may be gravitationally convected. The largest observed crystal size was near volume fraction 0.52. Near the freezing point, crystals should grow to large size because nuclei are sparse and there is plenty of room to grow without contacting other crystals. Near melting the opposite is true, with crystals growing into one another shortly after nucleation. The maximum crystal size seen near 0.52 volume fraction suggests that gravitational stresses are limiting crystal growth near freezing.

Ackerson, Bruce↗

A methodology to study crystal plasticity inside a compression test sample based on image correlation and EBSD

Modified compression tests in a coarse-grained austenitic stainless steel have been carried out in order to examine the crystal plasticity behavior for large plastic deformations. The measurements of local in-plane strains provide deeper insight into the deformation process on the local scale. These measurements are performed by digital image correlation technique (DIC) in combination with local crystal orientation measurements by using the electron backscatter diffraction technique (EBSD). Split test samples are prepared to track the strong microstructural changes during deformation, which is done incrementally in 10% steps up to 60% total macroscopic strain. The clear correlation of local strains with crystal orientation changes - e.g. in the case of mechanical twinning - permits to identify the acting deformation mechanisms. Such, experimentally determined local strain maps can be used for verification of crystal plasticity finite element method simulations (CPFEM). - Research Highlights: {yields} Method to study large strain crystal plasticity inside an austenitic FeCrNi-alloy. {yields} Correlation of local strain analyses with crystal orientation measurements. {yields} Deformation mechanism changes locally from dislocation glide to mechanical twinning. {yields} Suitable to study grain-grain interactions, slip system activation and grain boundary effects. {yields} Provide essential data for crystal plasticity FEM studies.

36 MATERIALS SCIENCE↗