Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

TTM-units #113

Open
wants to merge 2 commits into
base: develop
Choose a base branch
from
Open

TTM-units #113

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
bb73eb4
TTM unit fix
Jul 16, 2024
6784452
TTM units the other way
Jul 16, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 6 additions & 6 deletions src/program/temperature_pulse.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -106,17 +106,17 @@ double temperature_pulse_function(double function_time){
///
/// I(t) = 2/(tp delta sqrt(pi/ln 2)) exp[-(4 ln 2)(t/tp)**2] W/m^3
///
/// Assuming a nominal penetration depth of 10 nm and conversion from J/m^2 -> mJ/cm^2,
/// I(t) is now in units of s^1 m^-1 and laser power P is in units of mJ/cm^2
/// Assuming a nominal penetration depth of 10 nm and conversion from mJ/cm^2 to J/m^2
/// I(t) is now in units of s^1 m^-1 and laser power P is in units of J/m^2
///
///-----------------------------------------------------------------------------------------
double two_temperature_function(double ftime){

const double i_pump_time = 1.0/sim::pump_time;
const double reduced_time = (ftime-3.*sim::pump_time)*i_pump_time;
const double four_ln_2 = 2.77258872224; // 4 ln 2
// 2/(delta sqrt(pi/ln 2))*0.1, delta = 10 nm, J/m^2 -> mJ/cm^2 (factor 0.1)
const double two_delta_sqrt_pi_ln_2 = 9394372.787;
// 2/(delta sqrt(pi/ln 2))*0.1, delta = 10 nm, mJ/cm^2 -> (10 factor)
const double two_delta_sqrt_pi_ln_2 = 939437278.7;
const double pump=sim::pump_power*two_delta_sqrt_pi_ln_2*
exp(-four_ln_2*reduced_time*reduced_time)*i_pump_time;

Expand Down Expand Up @@ -148,8 +148,8 @@ double two_temperature_function(double ftime){
double double_pump_two_temperature_function(double ftime){

const double four_ln_2 = 2.77258872224; // 4 ln 2
// 2/(delta sqrt(pi/ln 2))*0.1, delta = 10 nm, J/m^2 -> mJ/cm^2 (factor 0.1)
const double two_delta_sqrt_pi_ln_2 = 9394372.787;
// 2/(delta sqrt(pi/ln 2))*0.1, delta = 10 nm, mJ/cm^2 -> J/m^2 (10 factor)
const double two_delta_sqrt_pi_ln_2 = 939437278.7;

const double i_pump_time1 = 1.0/sim::pump_time;
const double reduced_time1 = (ftime-3.*sim::pump_time)*i_pump_time1;
Expand Down