![]() |
Princeton Fusion Systems MATLAB Toolboxes
MATLAB toolboxes for fusion energy analysis
|
| ▼ Coils | |
| ► Demos | Coils demos |
| ► Inductances | |
| CoAxForces.m | Forces between PCC and LTS for various axial distances Published Demo |
| FluxLoopDemos.m | Flux conservation of series vs |
| ► Mass | |
| AxialArrayHTS.m | PFRC HTS magnet array mass including cryocoolers Published Demo Compute the total mass as a function of wire temperature to see the relative impact of cryocooler mass |
| MagnetMassComparison.m | Magnet structural mass model comparison, shell vs |
| StructureMassDemo.m | Mass calculation demo Published Demo Demonstrates SolenoidStructuralMass |
| ► Pancake | |
| AxialPancake.m | Axial magnet pancake design script Published Demo The alpha and beta shape parameters define the stresses and the fields at the inner magnet wall |
| NozzlePancake.m | Nozzle magnet pancake design script Published Demo The alpha and beta shape parameters define the stresses and the fields at the inner magnet wall |
| PancakeMagnetModel.m | Draw a magnet using Tube Published Demo This consists of a set of double pancakes, each of which is drawn using Tube |
| ► SingleLoop | |
| CoilSizing.m | Coil sizing analysis Published Demo Amperium wire: Type 8502-350, copper laminated |
| CoilStudy.m | PFRC coil currents design study Published Demo Calculate the radii for an ellipsoidal solenoid and calculate the currents to maintain a desired average field |
| CompensationCoil.m | Create a compensating coil for the DFD dipole Published Demo Uses AWG3 wire |
| CurrentShapeDemo.m | Compare current profiles for PFRC magnets Published Demo Three current shapes are compared for fixed magnet radii |
| WireFieldEstimates.m | Field strength and currents for PFRC coils Published Demo |
| ► Wire | |
| NbTiWireDemo.m | NbTi model of a particular Supercon wire Published Demo The wire is Supercon 54S33 with properties: |
| ► Forces | Self and between-coil forces |
| CoilSelfForce.m | Calculate the self force for a circular coil with a given current |
| CoilStresses.m | CoilStresses>Demo Shows a example similar to Fig 4.2 in the reference |
| CoilStressFiniteSolenoid.m | First check comparison to infinite solenoid |
| CoilStressSolenoid.m | Stresses in an infinite solenoid |
| SolenoidForces.m | Total force between coaxial current loops |
| Help | |
| ► Inductances | |
| CoAxBField.m | Determine Bz |
| CoAxConserveFluxes.m | Compare effects of having multiple individually flux conserving loops (left side of plot) and having a small number of multiply-wound flux conserving loops in series (right side of plot) |
| CoAxEllipseRR.m | Filaments of a coil composed of elliptical Rigid-Rotor current profile |
| CoAxFluxField.m | Unpack into column vectors |
| CoAxFluxRadius.m | CoAxFluxRadius>Demo |
| CoAxMutualInductance.m | Determine the mutual inductance of coil 2 on coil 1 |
| CoAxPlotBZ.m | Plot magnitude of Bz (|Bz|) or Br contours |
| CoAxPlotFlux.m | Plot flux contours at sense points given coil positions and currents |
| CoAxRectangle.m | Produce filament parameters of a coil with a rectangular cross section |
| CoAxSelfInductance.m | Determine the self inductance of a coil |
| CoAxStruct.m | Coil data structure |
| CoAxTestParticle.m | Set up demo inputs |
| ► Magnetic | |
| CoAxFluxInChunks.m | Fprintf('CoAxFluxField '); tic; phi1 = CoAxFluxField(rS,zS,rC,zC,IC); toc; fprintf('Discrepancy was %.3g ',max(max(abs(phi-phi1)))/mean(mean(phi))) |
| CoAxL2.m | Evaluate some simplified models, until you get to measure an L-2 coil: |
| CoAxL2Pack.m | Produce filament parameters of an L-2 coil pack |
| ConstBFlux.m | Enclosed flux assuming constant B |
| ► Mass | Coil mass estimates |
| CoilTurns.m | Determine the number of turns and mass for a coil of round wire |
| MagnetArrayMass.m | Default HTS properties |
| MagnetDataStructure.m | Data.dims = [0.1 0.1]; |
| SolenoidStructuralMass.m | Structural mass for a discrete solenoid |
| ► Pancake | Pancake-wound magnet models |
| Bpancake.m | Magnetic field of a pancake coil |
| PancakeMagnetMass.m | Size a pancake magnet to produce a desired B given the wire properties |
| PancakeMagnetTurns.m | Rectangular wire |
| REffPancake.m | Effective radius of a thin pancake coil |
| ► SingleLoop | Single loop current models |
| ConcentricCoils.m | Field for coaxial coils (poorly named) |
| ConcentricCoilsFiniteWidth.m | Field for concentric coils with a finite width |
| CurrentForField.m | Current to produce a desired field at the center of a wire loop |
| FieldCurrentLoop.m | Field along the central axis of a single current loop |
| ForceBetweenCurrentLoops.m | Total force between concentric current loops |
| ► Solenoid | Discrete solenoid model |
| EllipsoidalSolenoid.m | Demo |
| FLinearCurrent.m | Demo |
| FSteppedCurrent.m | Demo |
| ScanMagnetArray.m | Individual magnet fields |
| SolveSolenoid.m | Solve for the maximum current that produces the target peak field |
| ► Wire | HTS and LTS wire models |
| NbTiHysteresis.m | Determine the heating due to hysteresis in an AC NbTi coil |
| NbTiModel.m | Compute the critical current and temperature for NbTi wire |
| NbTiPenetrationField.m | Determine the penetration field of an AC NbTi coil |
| SuperconductorDataStructure.m | Width and thickness in mm |
| SuperPowerHTS.m | Properties of SuperPower HTS 12 mm tape |
| docs | |
| ▼ FusionEngine | |
| ► Cryo | Cryocooler sizing |
| Cryocooler.m | Demo |
| CryocoolerBrake.m | Demo |
| HeatCapacityAluminum.m | Check that H is good enough: |
| ThermalConductivityG10.m | Thermal conductivity of G-10 Cryo at a given temperature |
| ThermalConductivityInconel.m | Give the thermal conductivity of Inconel at a given temperature |
| ► Demos | Fusion engine demos |
| ► Cryo | Cryocooler sizing demos |
| CryocoolerAnalysis.m | Cryocooler analysis Published Demo Plots for three values of lift heat from 50 to 100 K cold temperature |
| CryoModelComparison.m | Cryocooler model comparison Published Demo Compare the two models at 60 K, by computing input power |
| HTSAxialMagnets.m | HTS Axial coil with cryocooler mass analysis Published Demo Calculate the coil and cryocooler mass for a magnet with a specific inner radius and target central field, as a function of operating temperature |
| LTSCoolerModel.m | Cryocooler analysis at 4 K, LTS magnet temperatures Published Demo See also: CryocoolerBrake |
| MagnetMassWithCryo.m | Coil with cryocooler mass analysis Published Demo Calculate the coil and cryocooler mass for a magnet with a specific inner radius and target central field, as a function of operating temperature |
| ► FRC | FRC demos |
| DrawFRCs.m | Draw FRCs Published Demo See also FRCFieldLines, DrawFRCLines, NewFig |
| ► HeatTransfer | Heat exchanger and thermal demos |
| CoilHeatTransfer.m | Total flux to the magnets through a layer of MLI Published Demo Compute temperatures on a radiation shield with MLI |
| HeatExchangerDesign.m | Study the heat exchanger between the plasma and the magnets Published Demo The heat exchanger is modeled as a flat plate with a uniform flux q on one side and a cooling gas flow |
| ShellThermalAnalysis.m | Shell analysis Published Demo Design a MLI blanket to shield the solar radiation |
| SizeExchangerChannel.m | Size the channels for heat transfer with He/Xe Published Demo Create a set of circumferential channels around the fusion region |
| ThermalCapacity.m | Thermal capacity Nitrogen and CO2 Published Demo |
| Total1DHeatTransfer.m | Complete thermal analysis from fusion to magnets Published Demo |
| ► PlasmaPhysics | Physics demos |
| CrossSectionDemo.m | Compare fusion cross-section methods Published Demo |
| TripleProduct.m | Compute the triple product for plasma parameters (script) Published Demo |
| ► Reactor | Point design demos |
| BremsstrahlungAnalysis.m | Bremmstrahlung analysis Published Demo |
| BremsstrahlungComparison.m | Compare Bremsstrahlung for several fuels |
| DHe3Reactor.m | D-He3 fusion power and magnetic field from electron temperature Published Demo Demonstrates Bremsstrahlung, FusionPower, and Synchrotron, assuming ions and electrons are at the same temperature |
| FRCPointDesign.m | Generates a point design for a steady-state, D-He3 FRC reactor |
| FuelPowerModel.m | D-3He fuel mix analysis demo creating contour plots Published Demo Uses a fixed electron temperature for all computations |
| FusionLosses.m | Look at fusion losses for several fuels Published Demo Demonstrates Bremsstrahlung and Synchrotron |
| FusionPowerDemo.m | Explore D-3He fusion power model Published Demo Explore combinations of D and 3He temperature and fuel ratios |
| PowerPlantAnalysis.m | Analyze a D-T power plant |
| SigmaVPlot.m | Plot collision crosssections for D-D, D-T and D-He3 reactions |
| SynchrotronAnalysis.m | Synchrotron analysis of reflection Published Demo Compare Bingren's model with the NRL formulary for PFRC relevant parameters |
| SynchrotronDemo.m | Synchrotron radiation demo Published Demo Reproduce the chart in the reference |
| ► RFSystem | Radio-frequency system demos |
| PowerAmplifierDemo.m | Class E power amplifier design demo Published Demo |
| RFMass.m | RF subsystem mass estimation demo Published Demo |
| ► Shielding | Neutron and xray shielding demos |
| NeutronAnalysis.m | Fast and thermal neutron analysis for LiH and B4C Published Demo Want to first scatter the neutrons until they are thermalized, then capture them |
| NeutronPowerDemo.m | Ratio of neutron power to fusion power for D3He Published Demo For PFRC parameters: Te = 30 keV, Td = 70 keV, and T(he3) ranges from 70 to 140 keV |
| SCO2XRay.m | X-ray absorption for supercritial CO2 Published Demo |
| Shielding.m | Shielding mass simple calculations Published Demo Calculate simple cylinders of potential shielding materials |
| ThicknessCalc.m | Shielding thickness calculations from wall loading Published Demo Results of a simple attenuation model |
| XRayAttenuation.m | Demonstrate x-ray attenuation for a Xenon/Helium mix and Titanium Published Demo X-rays are from bremsstrahulung |
| ► Tokamak | Tokamak design demos |
| FusionReactorDesignDemo.m | Fusion Reactor Design demo Published Demo Reproduce the figures/tables from the original reference paper |
| ► Economics | Cost modeling |
| CapacityFactor.m | Demo |
| CRF.m | Demo |
| FuelCost.m | Demo |
| LCOE.m | Demo |
| LCOF.m | Demo |
| MachineCost.m | Infer the cost of a machine size given a few sizes and costs |
| OvernightCapitalCost.m | Demo |
| ► FRC | FRC relations |
| BetaBarnes.m | Compute average beta for an FRC using Barnes relation |
| BVacuum.m | Vacuum field calculated from magnet and separatrix radii |
| DrawFRCLines.m | Draw FRC field lines, Solov'ev model Published Demo Create a figure and draw the lines for the given radius and elongation |
| FRCFieldLines.m | Demo |
| FRCInertia.m | FRCInertia>>Demo |
| ► HeatTransfer | Heat tranfer relations |
| RadiationShield.m | Computes temperatures on a radiation shield with MLI |
| Help | |
| ► Magnets | Magnet fields and masses |
| MagneticFieldCurrentLoop.m | Computes fields for multiple current loops |
| MagnetMass.m | Demo |
| MagnetMassVirial.m | Demo |
| MagnetStoredEnergy.m | Demo |
| NozzleRadius.m | Demo |
| ► PlasmaPhysics | Plasma physics functions |
| AlfvenSpeed.m | Demo |
| AtomicNumber.m | Demo |
| BremsstrahlungFrequency.m | Demo |
| BremsstrahlungWavelength.m | Demo |
| BurnupTime.m | Reactant burnup time |
| CollisionTimes.m | Ne converted to /cm3, Te to eV |
| CoulombLog.m | Coulomb log for electron-ion collisions |
| CoulombLogIons.m | Coulomb log for ion-ion collisions |
| CrosssectionClassical.m | Demo |
| CrossSectionFromRest.m | Demo |
| CyclotronFrequency.m | Demo |
| DebyeLength.m | Demo |
| DHe3Reaction.m | Demo |
| EffectiveCharge.m | Demo D-He3 |
| FusionReaction.m | Demo |
| FusionReactionDD.m | Demo |
| GasPressure.m | Demo |
| GyroRadius.m | Demo |
| IonCyclotronFrequency.m | Demo |
| IonizationEnergy.m | Demo |
| IonizationEnergyPlasma.m | Demo |
| MaxwellianV.m | Demo |
| NumberOfElectrons.m | Demo |
| PhotonEnergy.m | Demo |
| PlasmaConstant.m | PSS internal file version information |
| PlasmaFrequency.m | Demo |
| PlasmaHeating.m | Demo |
| PlasmaPressure.m | Demo |
| PlasmaRadius.m | Demo |
| PlasmaVolume.m | Demo |
| PowerIonsToElectrons.m | Demo |
| ReactionString.m | Demo |
| SahaEquation.m | |
| SigmaV.m | Demo |
| SynchrotronFrequency.m | Speed of light |
| ► Reactor | Fusion reactor modeling |
| Beta.m | Permeability of free space |
| BField.m | Demo |
| Bremsstrahlung.m | Demo |
| CentralPressure.m | Demo |
| CentralTemperature.m | Demos |
| DensityFromFusionPower.m | Demo |
| DivertorCLR.m | Demo |
| DMassFromPower.m | Demo |
| FirstWallThermalPower.m | Demo |
| FusionPower.m | Demo |
| FusionPowerDHe3.m | Demo |
| He3MassFromPower.m | Demo |
| MagneticPressure.m | Demo |
| MagnetRadius.m | Demo |
| NetFusionPower.m | Demo |
| NumberDensityFromBeta.m | Demo |
| PlasmaInternalEnergy.m | Demos |
| PointFusionReaction.m | PSS internal file version information |
| PowerDD.m | Demo |
| PowerDHe3.m | Demo |
| PowerDT.m | Demo |
| RadiusSphericalTorus.m | |
| ReactorField.m | 2017 |
| Synchrotron.m | Demo |
| VolumeAveragedPressure.m | Demo |
| ► RFSystem | Inductor and amplifier modeling |
| ClassERFAmplifier.m | Demo |
| Inductor.m | Model an inductor |
| InductorDesign.m | Design an inductor (fmincon) |
| SquareWave.m | Generates a square wave with a magnitude either 0 or 1 |
| ► Shielding | Shielding attenuation |
| MacroCrossSection.m | LiH for thermal neutrons |
| NumberDensity.m | LiH |
| ShieldingThickness.m | Shielding thickness from neutron attenuation model |
| XRayMassAttenCoeff.m | Computes the X-ray mass attenuation coefficient |
| ► Tokamak | Tokamak beta and point modeling |
| EdgeLocalizedModes.m | Demo |
| FusionReactorDesign.m | Unpack parameters |
| SteadyStateTokamak.m | Inductances and resistances |
| TokamakBeta.m | Demo |
| TPVHRHS.m | PSS internal file version information |
| ► Transport | Transport loss relations |
| ClassicalConfinementTime.m | Question: how handle all the collisions in the plasma? |
| CollisionsTrubnikov.m | Electrons and hydrogen |
| IonElectronEquilibration.m | The coulomb log for electron/ion collisions |
| IonIonCollisionRate.m | The coulomb log is a scalar for a pair of ions |
| LawsonProduct.m | Use a volume of 1 m3 |
| PowerIE.m | Nu = IonElectronEquilibration( Te, nI, Ti, zI, aI ) |
| TransportPower.m | TODO: add electron component |
