New module for dualstep thermomechanical ablation simulation

Author: ShiWeiHuH

src/HyperFSI.jl

@@ -31,7 +31,7 @@ export BBTMaterial, BBTMMaterial, BBTTMaterial, BBTTAMaterial, BBTMAMaterial
 export hsource_bc!, hsource_databc!, temperature_ic!, temperature_bc!, temperature_databc!, second_bcs!, find_sec_bcs_points 
 
 # Running simulations
-export Thermstep, Thermomechstep, Dualstep, Thermstep_ablation, Thermomechstep_ablation, Flowstep, FSI_job, FSI_submit, IBM2D, Bcstruct
+export Thermstep, Thermomechstep, Dualstep, Thermstep_ablation, Thermomechstep_ablation, Dualstep_ablation, Flowstep, FSI_job, FSI_submit, IBM2D, Bcstruct
 
 
 include("IBM/boundary_counter.jl")
@@ -49,13 +49,15 @@ include("structure/time_solvers/thermomechstep.jl")
 include("structure/time_solvers/dual_timesteps.jl")
 include("structure/time_solvers/thermstep_ablation.jl")
 include("structure/time_solvers/thermomechstep_ablation.jl")
+include("structure/time_solvers/dual_timesteps_ablation.jl")
 
 include("structure/physics/bond_based_thermal_diffusion.jl")
 include("structure/physics/bond_based_thermomechanics.jl")
 include("structure/physics/bond_based_dualstep_thermomechanics.jl")
 include("structure/physics/bond_based_thermal_diffusion_temp_dependent.jl")
 include("structure/physics/bond_based_thermal_diffusion_ablation.jl")
 include("structure/physics/bond_based_thermomechanics_ablation.jl")
+include("structure/physics/bond_based_dualstep_thermomechanics_ablation.jl")
 
 include("fluid/FlowTimesolver.jl")
 include("fluid/Evolution.jl")

src/coupling_solvers/FSI_submit.jl

@@ -316,6 +316,13 @@ function FSI_submit(job::Job, mode::String, geo::AbstractPDGeometry ; kwargs...)
             ret = TMA_submit_threads(job, nthreads(), geo)
         end
         return ret        
+    elseif mode == "DTMA" # dualstep_thermomech ablation in structure
+        if Peridynamics.mpi_run()
+            ret = DTMA_submit_mpi(job)
+        else 
+            ret = DTMA_submit_threads(job, nthreads(), geo)
+        end
+        return ret           
     end
 end
 
@@ -409,6 +416,24 @@ function TMA_submit_threads(job::Job, n_chunks::Int, geo::AbstractPDGeometry)
     return dh
 end
 
+function DTMA_submit_threads(job::Job, n_chunks::Int, geo::AbstractPDGeometry)
+
+    simulation_duration = @elapsed begin
+        logo_init_logs(job.options)
+        Peridynamics.log_spatial_setup(job.options, job.spatial_setup)
+        Peridynamics.log_create_data_handler_start()
+        dh = Peridynamics.threads_data_handler(job.spatial_setup, job.time_solver, n_chunks)
+        Peridynamics.init_time_solver!(job.time_solver, dh)
+        Peridynamics.initialize!(dh, job.time_solver)
+        Peridynamics.log_create_data_handler_end()
+        Peridynamics.log_data_handler(job.options, dh)
+        Peridynamics.log_timesolver(job.options, job.time_solver)
+        solve_dual_thermomech_struct_ablation!(dh, job, geo)
+    end
+    Peridynamics.log_simulation_duration(job.options, simulation_duration)
+    return dh
+end
+
 function merge_chunk(C::Peridynamics.AbstractThreadsBodyDataHandler)
     n = C.chunks[end].system.chunk_handler.loc_points[end]
     pd_out = fill(0.0, 8, n)

src/structure/physics/ablation.jl

@@ -74,30 +74,49 @@ function find_new_bcs_idx(dh::Peridynamics.AbstractThreadsBodyDataHandler, index
     return new_add_bcs_idx, idx_map
 end
 
-# new boundary occurs after ablation, energy from ablation is divided by neighbors!
-function update_new_bcs!(chunk::Peridynamics.AbstractBodyChunk, idx_map::Dict{Int, Vector{Int}}, new_ablation_point_idx::Vector{Int})
+# new boundary occurs after ablation, only for constant flux q_const!
+function update_new_bcs_q!(chunk::Peridynamics.AbstractBodyChunk, idx_map::Dict{Int, Vector{Int}}, new_ablation_point_idx::Vector{Int})
 
     for id in new_ablation_point_idx 
         ids_new_bcs = idx_map[id]
-        energy_add_aver = 0.0
+        q_const = zeros(1,1)
         if length(ids_new_bcs) > 0
-            energy_add_aver = chunk.paramsetup.energy 
+            q_const = chunk.storage.hsource[1, id] 
+            chunk.storage.hsource[1, id] .= 0.0
 
             for j in intersect(ids_new_bcs, chunk.system.chunk_handler.loc_points)
                 j_ch = chunk.system.chunk_handler.localizer[j]
-                chunk.storage.hsource[1, j_ch] = energy_add_aver
+                chunk.storage.hsource[1, j_ch] = q_const
             end
         end 
     end
 end
 
+# ablation energy source
 function update_new_bcs!(chunk::Peridynamics.AbstractBodyChunk, new_add_bcs_idx::Vector{Int})
 
     for id in intersect(new_add_bcs_idx, chunk.system.chunk_handler.loc_points) 
         i_ch = chunk.system.chunk_handler.localizer[id]
         chunk.storage.hsource[1, i_ch] = chunk.paramsetup.energy
     end
 end
+
+# new boundary occurs after ablation, only for constant pressure p_const!
+function update_new_bcs_f!(chunk::Peridynamics.AbstractBodyChunk, idx_map::Dict{Int, Vector{Int}}, new_ablation_point_idx::Vector{Int})
+    for id in new_ablation_point_idx 
+        ids_new_bcs = idx_map[id]
+        p_const = zeros(3, 1)
+        if length(ids_new_bcs) > 0
+            p_const = chunk.storage.b_ext[:, id]
+            chunk.storage.b_ext[:, id] .= 0.0
+
+            for j in intersect(ids_new_bcs, chunk.system.chunk_handler.loc_points)
+                j_ch = chunk.system.chunk_handler.localizer[j]
+                chunk.storage.b_ext[:, j_ch] = p_const
+            end
+        end 
+    end
+end
 #=
 function find_neighbors_indices(position::Matrix{Float64}, new_ablation_point_idx::Vector{Int}, radius::Float64)
     nhs = GridNeighborhoodSearch{3}(search_radius=radius, n_points=size(position, 2))

src/structure/physics/bond_based_dualstep_thermomechanics_ablation.jl

@@ -0,0 +1,80 @@
+"""
+Bond-Based dualstep_thermomechanics Module with Temperature-Dependent Properties and Ablation
+
+This module extends the bond-based dualstep_thermomechanics model 
+by incorporating Ablation Behavior.
+
+··This module provides functionality for ablation, primarily for material point removal.
+··Uses an 'exist' state to indicate ablation status, which relates to ablation amount.
+··Specific chemical processes are configured in the ablation file, with flexible complexity.
+··Core functionality identifies new boundary points, 
+  and more critically, establishes new watertight boundaries for subsequent fluid-structure interaction.
+··Ablated points can be modeled to describe subsequent effects.
+··The variation of properties with temperature is under consideration.
+
+··Basically, Materials model are as same as BBTMAMaterial, only times itergral is different.
+"""
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:position})
+    return copy(system.position)
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:displacement})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:velocity})
+    return  zeros(3,size(system.position, 2))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:velocity_half})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:velocity_half_old})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:acceleration})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:b_int})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:b_int_old})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:b_ext})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:density_matrix})
+    return zeros(3, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:temperature})
+    return zeros(1, Peridynamics.size(system.position, 2))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:pflux})
+    return zeros(1, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:hsource})
+    return zeros(1, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:ablation_deep})
+    return zeros(1, Peridynamics.get_n_loc_points(system))
+end
+
+function Peridynamics.init_field_solver(::Dualstep_ablation, system::Peridynamics.AbstractSystem, ::Val{:exist})
+    return ones(Int, 1, size(system.position, 2))
+end
+
+
+
+
+