Single sheet

[alexruiz95]

Hello,

This is an awesome code great work. I have a question. I am going through the code and trying to see how I can modify it to have a single sheet. Is it easily possible?

[KRproject-tech]

It is time-consuming to correct this code to a single-sheet case.
Therefore, I will upload the program for a single sheet to GitHub later.

[alexruiz95]

Hmm.. Looking at new_X_func_FAST.m. I think it would be this line out_0 = out1.A1_A2_Xn + A_mat1\( d_t*[ zero_mat(:,1); Q_global]);

of this section?

`X2_matrix = [ zero_mat eye_mat;
zero_mat zero_mat];
A_mat1 = D_matrix - alpha_vd_tX2_matrix;
A_mat2 = D_matrix + (1 - alpha_v)d_tX2_matrix;

%%[2-0] 固定部のノードでは dt_q = 0, dtt_q = 0 [2-0] dt_q = 0, dtt_q = 0 at the nodes of the fixed part
not_i_vec = (1:N_q_all);
not_i_vec(i_vec) = [];

%% [3] 状態ベクトル更新 [3] State vector update

if stage == 0 %% 1回目の計算値を使いまわす．%% Reuse the first calculated value.
out1.A1_A2_Xn = A_mat1( A_mat2X_vec([ not_i_vec N_q_all+not_i_vec]) );
end
out_0 = out1.A1_A2_Xn + A_mat1( d_t[ zero_mat(:,1);
Q_global]);

out = X_vec;
out([ not_i_vec N_q_all+not_i_vec]) = out_0;
out = blkdiag( Fmat, Fmat)*out;

end`

Im looking at the variables and values the sim produces and what new_X_func_FAST.m produces.

[alexruiz95]

Additionally i have two questions,

1. what is the typical compute time on your end? It takes a long time on 16 cores.
2. Do you happen to have essentially the same code with no structure? like just the UVLM? Since I only care for the flow field and aerodynamic coefficients induced by the bending and twisting a one way FSI.

I think this code will go along way with my research if I can adapt it.

Thanks for all the help and feedback it is really appreciated.

[KRproject-tech]

1. what is the typical compute time on your end? It takes a long time on 16 cores.
   Currently checking with Ryzen 5900X.

2. Do you happen to have essentially the same code with no structure? like just the UVLM? Since I only care for the flow field and aerodynamic coefficients induced by the bending and twisting a one way FSI.
   Added mass term in the pressure acting on a sheet is modified from the force calculation given by Katz & Plotkin to achieve strong two-way coupling between the fluid and a sheet.

[alexruiz95]

Hi,

Yeah I am curious on the compute time on your end. Sometimes things work fast on one computer but not the other. Is it supposed to be super fast? For a 2sec run for a single sheet it takes 1009 seconds.

Still struggling on the time dependent displacement. Perhaps a similar question that could help me work it out is how would I apply a displacement from the start? Im looking at how we apply the displacement of zero, how would I apply at other nodes a small displacement?

Thanks for all the feedback

[KRproject-tech]

The times employed in the numerical analysis is non-dimensionalised by length and flow velocity as

The computing time under Ryzen 5900X is 2.5 hour for t^* = 30 [-].