Message: Re: Stainless steel Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Stainless steel 

Forum: Geometry
Re: Question Stainless steel (Anil)
Date: 01 Jul, 2008
From: Joseph Perl <Joseph Perl>

Hi Anil,

Since there is no one definition for Stainless Steel, you have to define it yourself from an appropriate mixture of elements. Here's a method from one of my applications that defined my stainless steel among other things:

void ElectronBenchmarkDetector::DefineMaterials(){ 
	// Use NIST database for elements and materials whereever possible.
	G4NistManager* man = G4NistManager::Instance();
	man->SetVerbose(1);

	// Define elements from NIST 
	man->FindOrBuildElement("H");
	man->FindOrBuildElement("Be");
	G4Element* C  = man->FindOrBuildElement("C");
	man->FindOrBuildElement("N");
	man->FindOrBuildElement("O");
	man->FindOrBuildElement("Al");
	G4Element* Si = man->FindOrBuildElement("Si");
	man->FindOrBuildElement("Ti");
	G4Element* Cr = man->FindOrBuildElement("Cr");
	G4Element* Mn = man->FindOrBuildElement("Mn");
	G4Element* Fe = man->FindOrBuildElement("Fe");
	G4Element* Ni = man->FindOrBuildElement("Ni");
	man->FindOrBuildElement("W");
	man->FindOrBuildElement("Au");
	man->FindOrBuildElement("Pb");

	// Define pure NIST materials
	man->FindOrBuildMaterial("G4_Al");
	man->FindOrBuildMaterial("G4_Ti");
	man->FindOrBuildMaterial("G4_W");
	man->FindOrBuildMaterial("G4_Au");

	// Define other NIST materials
	man->FindOrBuildMaterial("G4_WATER");
	man->FindOrBuildMaterial("G4_KAPTON");
	G4Material* Air = man->FindOrBuildMaterial("G4_AIR");

	// Define materials not in NIST
	G4double density;
	G4int ncomponents;
	G4double fractionmass;
	G4Material* StainlessSteel = new G4Material("StainlessSteel", density= 8.06*g/cm3, ncomponents=6);
	StainlessSteel->AddElement(C, fractionmass=0.001);
	StainlessSteel->AddElement(Si, fractionmass=0.007);
	StainlessSteel->AddElement(Cr, fractionmass=0.18);
	StainlessSteel->AddElement(Mn, fractionmass=0.01);
	StainlessSteel->AddElement(Fe, fractionmass=0.712);
	StainlessSteel->AddElement(Ni, fractionmass=0.09);

	G4Material* beam = new G4Material("Beam", density= 1.e-5*g/cm3, ncomponents=1,
									  kStateGas, STP_Temperature, 2.e-2*bar);
	beam->AddMaterial(Air, fractionmass=1.);

	G4cout << *(G4Material::GetMaterialTable()) << G4endl;
}

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Stainless steel   (Anil - 02 Jul, 2008)
 Add Message Add Message
to: "Re: Stainless steel"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]