Copper Shock Simulation of Yield Strength

ISEF Category: Materials Science

Subcategory: Computation and Theory  ·  Difficulty: Advanced  ·  Setup: University Lab  ·  Time: Full Year

The Hook

Copper can get stronger when its internal structure gets smaller. That sounds backward, but atoms do not follow the same rules as a backpack or a bridge. In a shock, tiny defects can decide when the metal starts to deform. You can test that idea with molecular dynamics and see how strain rate changes yield strength.

What Is It?

This project studies how nano-twinned copper behaves when a shock wave hits it. Nano-twins are thin mirror-like layers inside a crystal. Think of them like neat stacks of cards inside the metal. Those layers can block or steer dislocations, which are the tiny defects that let metals bend instead of snap.

Molecular dynamics, or MD, is a simulation method that tracks atoms as they move. LAMMPS is a free MD code that many researchers use. In your project, you would build a copper model with twin boundaries, apply shock loading, and watch when dislocations appear. The main question is how the loading rate changes the stress needed to start plastic deformation, which is when the material begins to permanently change shape.

Why This Is a Good Topic

This is a strong science fair topic because you can change one variable, measure a clear outcome, and compare your results across several simulation runs. It connects to real problems in impact resistance, aerospace materials, and protective design. You can learn how to run atomistic simulations, define controls, read stress-strain data, and fit scaling trends, all skills that matter in materials research.

Research Questions

• How does strain rate change the yield strength of nano-twinned copper under shock loading?
• What is the effect of twin spacing on the stress needed for dislocation nucleation?
• Does the number of twin boundaries change the first location where plastic deformation begins?
• To what extent does shock direction relative to the twin planes alter yield behavior?
• Which loading rates produce the largest gap between elastic response and permanent deformation?
• What is the effect of sample size on the measured strain-rate scaling in the simulation?

Basic Materials

• Computer with enough memory to run atomistic simulations.
• LAMMPS installed on a local machine or access to a university cluster.
• Atomistic visualization tool such as OVITO or VMD.
• Python with NumPy, SciPy, and Matplotlib for data analysis.
• Text editor or notebook for recording simulation settings.
• Reference paper on nano-twinned copper and shock response.
• Spreadsheet software for tracking runs and outputs.

Advanced Materials

• High-performance computing access for longer LAMMPS runs.
• Prebuilt copper interatomic potential file validated for shock or defect studies.
• OVITO Pro or open-source OVITO for defect analysis and coordination maps.
• Python with pymatgen, ASE, NumPy, SciPy, and pandas for workflow and analysis.
• Version control software such as Git for tracking model changes.
• Cluster job scheduler access and submission scripts.
• Published benchmark datasets for copper shock or dislocation nucleation.

Software & Tools

• LAMMPS: Runs the molecular-dynamics simulation of copper under shock loading.
• OVITO: Visualizes atoms, dislocations, and twin boundaries after each run.
• Python: Processes stress, strain, and yield data, then fits scaling trends.
• ImageJ: Measures or compares visual features in exported simulation frames when needed.
• Excel or Google Sheets: Tracks simulation conditions and organizes output values.

Experiment Steps

1. Define the exact question you want the simulation to answer, then choose one primary variable, such as strain rate or twin spacing.
2. Build a simple copper model with nano-twin boundaries and decide how you will verify that the structure is correct before loading begins.
3. Plan a set of loading conditions that changes only one factor at a time, so you can compare outcomes fairly.
4. Choose the signal that marks yield, such as the first big deviation from linear stress response or the first clear defect burst.
5. Design your analysis workflow before you run the full batch, including how you will extract stress-strain curves and detect dislocation onset.
6. Set up repeat runs and controls, then decide how you will test whether any scaling trend stays consistent across sample sizes or orientations.

Common Pitfalls

• Using a copper potential that was not tested for shock or defect behavior, which can give fake yield trends.
• Comparing runs with different cell sizes without normalizing the geometry, which mixes size effects with strain-rate effects.
• Calling the first noisy stress spike “yield,” which can mistake vibration for real plastic deformation.
• Ignoring whether the twin boundaries actually stayed intact, which can hide a different deformation path.
• Running only one simulation per condition, which makes random atomic noise look like a real trend.

What Makes This Competitive

A stronger project goes beyond a single stress-strain plot. You can compare more than one twin spacing, shock direction, or sample size, then test whether the same scaling law still holds. Clear defect analysis helps, too, because judges want to see that you know why the curve changed, not just that it changed. A careful uncertainty analysis and a clean comparison to published results can push the work much higher.

Project Variations

• Study how twin spacing changes dislocation nucleation in nano-twinned copper under the same shock conditions.
• Compare nano-twinned copper with ordinary polycrystalline copper to test whether twin boundaries delay yield.
• Analyze how shock direction relative to the twin planes changes the yield-strength scaling in the same model.

Learn More

• LAMMPS documentation: Search the official LAMMPS manual for shock loading, metal units, and defect-related examples.
• OVITO tutorials: Use the official OVITO site to learn how to identify dislocations and visualize crystal defects.
• NIST Materials Data Repository: Search for copper structure and property references, plus benchmark-style materials data.
• PubMed: Search for review articles on dislocation nucleation, strain-rate sensitivity, and shock deformation in metals.
• Acta Materialia: Search recent peer-reviewed papers on nano-twinned copper, dislocation activity, and deformation mechanisms.
• MIT OpenCourseWare: Look for materials science and computational materials courses that cover crystal defects and simulation basics.