👨?🎓個人主頁：研學社的博客

💥💥💞💞歡迎來到本博客????💥💥


🏆博主優(yōu)勢：🌞🌞🌞博客內容盡量做到思維縝密，邏輯清晰，為了方便讀者。


?座右銘：行百里者，半于九十。

📋📋📋本文目錄如下：🎁🎁🎁
目錄
💥1 概述
📚2 運行結果
🎉3 參考文獻
🌈4 Matlab代碼實現(xiàn)

💥1 概述

膨脹波是流體力學的基本概念之一，指流體中擾動區(qū)與未擾動區(qū)的分界面，流體通過此界面會壓力降低。膨脹波問題中一個典型例子是超音速氣流繞凸角的定常流動，這種流動稱為普朗特-邁耶爾流動。在這種流動中，氣流逐漸從一個方向轉到另一個方向，同時不斷加速。

📚2 運行結果

部分代碼：

%%

clear

%% inputs

freq=10^3;%*(2.^([-20:20]/3)); %frequency

%% setup domain and mesh

domain.dim=[0.05 0.025]; %dimensions of domain

domain.off=[0 0]; %offset of domain

domain.fce=[1 2 3 4]; %face ids (for boundary conditions) on [-x +x -y +y]

nnd=domain.dim./0.0025+1; %number of nodes in each direction

mesh=blockmesh(nnd,domain.dim,domain.off,domain.fce);

%% loop for frequency

for f=1:length(freq)

fprintf(['computing frequency ',num2str(f),' of ',num2str(length(freq)),' ... \n']);

%% setup boundary conditions

addpath('MESH')

bcs=initbcs(mesh.nodes,mesh.elems,mesh.faces,8);

bcs=addbcface(bcs,1,1,5);

bcs=addbcface(bcs,1,2,6);

bcs=addbcface(bcs,1,3,7);

bcs=addbcface(bcs,1,4,8);

bcs=addbcface(bcs,1,1,9);

% bcs=addbcface(bcs,1,2,10);

bcs=addbcface(bcs,1,3,11);

bcs=addbcface(bcs,1,4,12);

bcs.type{1}(:)=2;

bcs.coef{1}(bcs.face{1}==1,5)=1; %face 1

bcs.coef{1}(bcs.face{1}==2,1)=1; %face 2

bcs.coef{1}(bcs.face{1}==3,2)=1; %face 3

bcs.coef{1}(bcs.face{1}==4,2)=1; %face 4

bcs.coef{1}(bcs.face{1}==5,6)=1; %face 5

bcs.coef{1}(bcs.face{1}==6,3)=1; %face 6

bcs.coef{1}(bcs.face{1}==7,4)=1; %face 7

bcs.coef{1}(bcs.face{1}==8,4)=1; %face 8

bcs.coef{1}(bcs.face{1}==9,8)=1; bcs.rhds{1}(bcs.face{1}==9)=1; %face 9

% bcs.coef{1}(bcs.face{1}==10,1)=1; %face 10

bcs.coef{1}(bcs.face{1}==11,6)=1; %face 11

bcs.coef{1}(bcs.face{1}==12,6)=1; %face 12

%% setup physics

addpath('PLANES')

air=air_properties_generic;

medium=4003;

PEM.name_mat=['Mat_porous_' num2str(medium-1000*floor(medium/1000))];

PEM.typ_mat=floor(medium/1000);

eval(['PEM=Mat_porous_' num2str(medium-1000*floor(medium/1000)),'(PEM);'])

PEM=properties_JCA(PEM,air,freq(f));

PEM=properties_PEM(PEM,air,freq(f));

physics=PEM;

%% assemble

[stiff,force,err,tol,condH]=assemble(mesh,bcs,physics,freq(f));

%% solve

fprintf('solving ... \n');

operationtime=cputime;

fprintf(' solving directly ... ');

q=full(stiff\force);

fprintf('done\n');

%% postprocess

nnds=mesh.nnds;

soln.usx(:,f)=q((1:nnds)+nnds*0);

soln.usy(:,f)=q((1:nnds)+nnds*1);

soln.utx(:,f)=q((1:nnds)+nnds*2);

soln.uty(:,f)=q((1:nnds)+nnds*3);

soln.sxx(:,f)=q((1:nnds)+nnds*4);

soln.sxy(:,f)=q((1:nnds)+nnds*5);

soln.syy(:,f)=q((1:nnds)+nnds*6);

soln.p(:,f)=q((1:nnds)+nnds*7);

% soln.cond(:,f)=condest(stiff);

us=mean(reshape(soln.usx(:,f),nnd(1),nnd(2)),2);

ut=mean(reshape(soln.utx(:,f),nnd(1),nnd(2)),2);

s=mean(reshape(soln.sxx(:,f),nnd(1),nnd(2)),2);

p=mean(reshape(soln.p(:,f),nnd(1),nnd(2)),2);

Z=p(1)/ut(1);

R=(Z-air.rho*air.c)/(Z+air.rho*air.c);

alpha(:,f)=1-abs(R)^2;

🎉3 參考文獻

部分理論來源于網(wǎng)絡，如有侵權請聯(lián)系刪除。

[1]劉占芳,嚴波,唐錄成.飽和多孔彈性材料中加速度波的傳播[J].重慶大學學報(自然科學版),1998(02):9-17.

🌈4 Matlab代碼實現(xiàn)