This website requires certain cookies to work and uses other cookies to help you have the best experience. By visiting this website, certain cookies have already been set, which you may delete and block. By closing this message or continuing to use our site, you agree to the use of cookies. Visit our updated privacy and cookie policy to learn more.
This Website Uses Cookies
By closing this message or continuing to use our site, you agree to our cookie policy. Learn More
This website requires certain cookies to work and uses other cookies to help you have the best experience. By visiting this website, certain cookies have already been set, which you may delete and block. By closing this message or continuing to use our site, you agree to the use of cookies. Visit our updated privacy and cookie policy to learn more.
Industrial Heating logo
search
cart
facebook twitter linkedin youtube
  • Sign In
  • Create Account
  • Sign Out
  • My Account
Industrial Heating logo
  • Home
  • Magazine
    • Current Issue
    • Digital Edition
    • Archives
  • News
  • Featured
    • IH Daily
    • IH MagEzine
    • Web Exclusives
    • IH Economic Indicators
    • The History of Industrial Heating
    • Heat Treatment Processes
    • Top 10 Heat-Treated Holiday Gifts
  • Topics
    • Additive Manufacturing / 3D Printing
    • Ceramics & Refractories / Insulation
    • Combustion & Burners
    • Heat Treating
    • Heat & Corrosion Resistant Materials / Composites
    • Induction Heat Treating
    • Industrial Gases & Atmospheres
    • Materials Characterization & Testing
    • Melting / Forming / Joining
    • Process Control & Instrumentation
    • Sintering / Powder Metallurgy
    • Vacuum / Surface Treatments
  • Columns
    • Editorial
    • The Heat Treat Doctor
    • Federal Triangle
    • MTI Profile
    • Academic Pulse
    • Heat Treat 5.0
    • International – Brazil
    • Next-Gen Leaders
  • Directories
    • Equipment Buyers Guide
    • Commercial Heat Treat Capabilities Directory
    • Aftermarket Parts & Services Directory
    • Materials Characterization & Testing Equipment Directory
    • Take a Tour
  • More
    • Classifieds
    • White Papers
    • Industrial Heating Bookstore
    • Organizations
    • Market Research
    • Custom Content & Marketing Services
    • FORGE Magazine
  • Multimedia
    • Podcasts
    • Videos
    • Webinars
    • Image Gallery
    • Mobile App
    • eBooks
  • Events
    • Meetings & Trade Shows
    • FNA
    • Heat Treat Show
  • Blog
    • Dan Herring - Heat Treatment
    • David Pye - Metallurgy
    • Dan Kay - Brazing
    • Debbie Aliya - Failure Analysis
    • Thomas Joseph - Intellectual Property
  • Contact
  • Advertise
  • Subscribe
    • Print & Digital Edition Subscriptions
    • eNewsletter
    • Online Registration
    • Customer Service
Home » Engineering Concepts
Commentary/Columns

Engineering Concepts

May 1, 2000
Dave VanAken
Reprints
No Comments
Diffusion by Numerical Analysis

Fig. 1 A schematic concentration profile, where the flux of atoms, J, is down the concentration gradient.
Ever been overwhelmed by the mathematics of diffusion? Stop looking for your error function table and grab your personal computer. All you need is a spreadsheet program, a little knowledge of numerical analysis and you can solve all sorts of complicated diffusion problems and that includes problems that do not have analytical solutions.

Most diffusion problems can be solved using a finite difference method of numerical analysis and Fick's first law of diffusion. Consider the solute profile shown in figure 1 and the change in composition at C1 produced by a flux into ( Jin) and a flux out of (Jout) a region of thickness h. The flux equations would look like the following:

Jin= [qin/delta(t)] = - D [(C1 - Co)/h]

Jout= [qout/delta(t)] = - D [(C2 - C1)/h]

where delta(t) is a finite time interval, D is the diffusivity of the solute, q represents the amount of solute per unit area entering or leaving the volume of interest and Ci is the solute per unit volume. The net amount of solute collected in time delta(t) will be

C1new - C1 = [(qin - qout)/h] = [(D x delta(t))/h^2 x (Co + C2 - 2C1)

This equation can be greatly simplified by setting [D x delta(t)/h^2] equal to 1/2 to give

C1new = (Co + C2/2

The only constraint is that the time interval and distance between composition measurements are fixed. This requirement is called the stability criterion for the finite difference method.



Fig. 2 An illustration of how a spreadsheet is set up to perform a numerical analysis of a carburization process.
Consider a carburization problem where the surface of a 0.22 weight percent (wt.%) carbon steel is held in a carburizing atmosphere. The carbon potential of the atmosphere is 1.1 wt.% at 925°C (1700°F) and the carbon profile is required after 8 hours. To solve this problem a spreadsheet would be set up as shown in figure 2. Only a portion of the spreadsheet is shown to indicate the simplicity of the technique and how each successive column is calculated. Each time interval (gt) requires a new column to be added and the time interval is set by the choice of distance interval (h) and the diffusivity (D). For a diffusivity of 1.59x10^-5 mm^2/sec, the gt time interval would be 314 seconds.

Fig. 3 A comparison of a spreadsheet calculation and the same problem solved using an analytical solution.
Thus, 92 columns or iterations of Æt would be required to calculate the profile for an 8 hour carburization cycle. To simplify the spreadsheet, the depth was run to 6 mm and the composition at 6 mm was kept constant throughout the 92 iterations at a value of 0.22 wt.% C. This is equivalent to a 6 mm sheet of steel that is held between two different atmospheres: one side of the steel sheet at a carbon potential of 1.1 wt.% and the other side at a potential of 0.22 wt.%. This constrains the carburization cycle to less than 10 hours. Longer times can be run by increasing the spreadsheet to incorporate greater depths, but the last row will still have to be held constant at the base carbon composition. A plot of the resulting carbon profile for the 8 hour cycle is shown in figure 3 along with the analytical solution using a mathematical approximation for the error function.

With a little effort, the spreadsheet can be turned into a powerful tool where carbon potential, base carbon content and diffusivity can be entered as variables to test new carburization cycles. Carbon leveling operations can also be modeled by lowering the surface carbon potential after the case depth is obtained. In the problem shown here, the time to attain a uniform carbon concentration, less than 0.9 wt.%, can be determined by changing the surface carbon concentration to 0.9 in column 93 and then inspecting the subsequent columns to determine the number of intervals required to level the carbon. Many spreadsheets are equipped with graphics packages and the effect of changing the variables can be immediately observed in these graphs at fixed number of iterations or times.

ih-subscribe

Recent Articles by Dave VanAken

ENGINEERING CONCEPTS: Hydrogen Embrittlement of Copper

ENGINEERING CONCEPTS: PID Controllers

ENGINEERING CONCEPTS: Transformation-Induced Residual Stresses

ENGINEERING CONCEPTS: Grains That Grow, and Secondary Recrystallization

Engineering Concepts: Segregation and Banding in Steel

Dave VanAken is a Professor, Dept. of Metallurgical Engineering,at University of Missouri-Rolla

Related Articles

ENGINEERING CONCEPTS: Relationship Between Hardness and Strength

ENGINEERING CONCEPTS: Grains That Grow, and Secondary Recrystallization

Engineering Concepts: Strength of Tempered Martensite

Engineering Concepts: Laser Shock Processing

You must login or register in order to post a comment.

Report Abusive Comment

Subscribe For Free!
  • Print & Digital Edition Subscriptions
  • eNewsletters
  • Online Registration
  • Subscription Customer Service

More Videos

Popular Stories

Editorial 2019: Reed Miller

Noel Nuggets

Industrial Heating Web Exclusives

Steel Mill Powered by Wind

1219IH-Fabrisonic-slide5

Metal Additive Manufacturing without Melting

Industrial Heating Industry News

Novelis to Expand, Upgrade Georgia Facility

Industrial Heating Industry News

Steel Dynamics Expansion Includes Induction Furnace, Billet Welder

IH Ipsen 360x184customcontent

Events

January 1, 2030

Webinar Sponsorship Information

For webinar sponsorship information, visit www.bnpevents.com/webinars or email webinars@bnpmedia.com.
View All Submit An Event

Poll

Additive Manufacturing

Has additive manufacturing had any impact on your business?
View Results Poll Archive

Products

Vacuum Heat Treatment Volume I

Vacuum Heat Treatment Volume I

See More Products

The History of Industrial Heating 1000 BC - Present Day


Industrial Heating Employment Marketplace

Industrial Heating

1219IH-cover144x192

2019 December

Check out the December 2019 issue of Industrial Heating, featuring "Metal Additive Manufacturing without Melting", "Furnaces with Tungsten Heating Elements Make High Product Quality Possible", and much more.

View More Create Account
  • Resources
    • List Rental
    • eNewsletter
    • Manufacturing Group
    • News
    • Want More?
    • Featured
    • Product / Event
    • Industry Links
    • Connect
    • Privacy Policy
    • Survey And Sample

Copyright ©2019. All Rights Reserved BNP Media.

Design, CMS, Hosting & Web Development :: ePublishing