Skip to content

add option to use dynamic heat source#32

Open
puente-tamu wants to merge 1 commit into
ORNL-MDF:mainfrom
puente-tamu:dynamic_heat_source
Open

add option to use dynamic heat source#32
puente-tamu wants to merge 1 commit into
ORNL-MDF:mainfrom
puente-tamu:dynamic_heat_source

Conversation

@puente-tamu
Copy link
Copy Markdown

Added option to use dynamic heat source from Coleman et al. Added dynamic heat source in Finch_Solver to have m,k options to change heat source shape. Backwards compatible such that if no "shape" is specified in user input for the heat "source" the default is the original gaussian distribution. Changes were made in Finch_Inputs to account for this.

@puente-tamu puente-tamu force-pushed the dynamic_heat_source branch from 6bf7dd6 to 5603402 Compare May 28, 2026 15:27
Comment thread src/Finch_Solver.hpp

// heat source parameter constants
for ( std::size_t d = 0; d < 3; ++d )
dynamic_beam_ = ( db.source.shape == "dynamic" );
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggestion: move this logic into a new HeatSource class that gets called here in the solver. I think this new complexity warrants some separation

Comment thread src/Finch_Solver.hpp
return ( dist_to_beam[0] * dist_to_beam[0] * A_inv_[0] ) +
( dist_to_beam[1] * dist_to_beam[1] * A_inv_[1] ) +
( dist_to_beam[2] * dist_to_beam[2] * A_inv_[2] );
if ( !dynamic_beam_ )
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would expect this to look something like:

return heat_source.weight( dx, dy, dz ... );

Comment thread src/Finch_Solver.hpp
double depth_; // heat source depth (= two_sigma[2])

public:
Solver( Inputs db, LocalMeshType local_mesh )
Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Suggested change
Solver( Inputs db, LocalMeshType local_mesh )
Solver( Inputs db, LocalMeshType local_mesh, HeatSource heat_source )

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants