You are now following this Submission

You will see updates in your followed content feed
You may receive emails, depending on your communication preferences

OpenBowShock 2D

Version 1.1.2.0 (16 MB) by Stanislav Karpuk

CFD solver for blunt bodies at high speed

4.0

(3)

613 Downloads

Updated 26 Feb 2016

View License

License

Share
Open in MATLAB Online
Download

Share 'OpenBowShock 2D'

Open in File Exchange

Open in MATLAB Online

Overview
Functions
Version History
Reviews (3)
Discussions (1)

Description: OpenBowShock2d calculates air flow properties around blunt bodies at supersonic and hypersonic speeds at steady state.
The code solves Euler equations using Steger-Warming and Van-Leer flux vector splitting.
for the tutorial and other information, visit http://openbowshock2dcom.com/index.html

Cite As

Stanislav Karpuk (2025). OpenBowShock 2D (https://se.mathworks.com/matlabcentral/fileexchange/54810-openbowshock-2d), MATLAB Central File Exchange. Retrieved October 31, 2025.

Requires

- C ++ -compiler. Make Sure the compiler is set up. Use http://www.mathworks.com/help/matlab/matlab_external/choose-c-or-c-compilers.html for more information - Adobe Reader

MATLAB Release Compatibility

Created with R2013a

Compatible with any release

Platform Compatibility

Windows macOS Linux

Tags Add Tags

aerodynamics blunt body cfd euler high speed steger warming supersonic

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Discover Live Editor

Create scripts with code, output, and formatted text in a single executable document.

Learn About Live Editor

Jacobian(imax,jmax,X,Y,dxi,deta)
RHS_metrics(X,Y,imax,jmax,p,q,dxi,deta,Xxi,Xeta,Yxi,Yeta)
Steger_Warming_implicit(AOA,k,dQbars,Jshock,conv,Tinf,U,V,p,ro,et,T,Yfin,Xix,Xiy,gamma,Ubar,Vbar,pbar,robar,imax,jmax,Etax,Etay,dQbar,dxi,deta,dt,alpha,J,R,theta,etbar,choice,alpha1,Aplus,Aminus,Bplus
Van_Leer(Xix,Xiy,gamma,Ubar,Vbar,pbar,robar,imax,jmax,aa,J,E,Eplus,Eminus)
Xixmetricesjacobian(imax,jmax,X,Y,dxi,deta,ym)
Xixmetricesjacobian_hyperbolic(imax,jmax,X,Y,dxi,deta,ym)
arc_function(bound,X1u,X1l,Y1u,Y1l,X2,choice,condb,split_top,split_bottom,point,coef,X0,l,handles)
blunt_body_draw(handles,bound,Y2,l)
boundary(bound,l,n,imax,betta1,Y,X,dxi,XI,choice)
boundary_conditions(jmax,etbar,Tbar,gamma,theta,alpha,robar,Ubar,Vbar,pbar)
check_shock_bound(imax,jmax,pnew,Jshock)
clustering(imax,jmax)
coefficients(Aplus,Bplus,Aminus,Bminus,C,dt,deta,dxi,imax,jmax,alpha,CAM,CA,CAP,CBM,CB,CBP)
computegrid(X,Y,imax,jmax,betta1,choice,a_grid_up,a_grid_down,b_grid,delta)
determine_surface_point(theta_s,theta,imax,R,n,jmax,choice,X,Y,Imax)
dimension(Qbarnew,jmax,J,imax,gamma,R,robar,Ubar,Vbar,pbar,etbar,Tbar,Jshock)
domain(choice,coef,invcoef,bound,l)
elliptic_function(bound,X1u,X1l,Y1u,Y1l,X2,choice,condb,split_top,split_bottom,point,coef,smin,l)
elliptic_grid(X1,Y1,imax,jmax,dxi,deta,N,top)
equilibrium_flow(imax,jmax,e,ro,p,T,gamma,A)
exponential_function(bound,X1u,X1l,Y1u,Y1l,X2,choice,condb,split_top,split_bottom,point,power,l,c,coef)
flux_shit(Xix,Xiy,Etax,Etay,robar,pbar,Ubar,Vbar,etbar,jmax,imax,J)
initial_conditions(imax,jmax,Minf,Tinf,roinf,pinf,AoA)
initial_position(dr,alpha,n,c,X,Y,imax)
linear_function(bound,X1u,X1l,Y1u,Y1l,X2,Y2,choice,condb,split_top,split_bottom,point,coef,l)
local_dt(pbar,robar,Ubar,Vbar,Xix,Xiy,Etax,Etay,gamma,imax,jmax,dxi,deta,A,dt)
logariphmic_function(bound,X1u,X1l,Y1u,Y1l,X2,choice,condb,split_top,split_bottom,point,l,e,coef)
matrix_A(Xix,Xiy,gamma,Ubar,Vbar,pbar,robar,imax,jmax,Aplus,Aminus,aa)
metrics_for_boundaries(dxi,deta,imax,jmax,X,Y,Xxi,Xeta,Yxi,Yeta,Xxi2,Yxi2,Xeta2,Yeta2,Xxieta,Yxieta,J,p,q,r,s,k)
mixed_grid(X1,Y1,imax,jmax,choice,ym,ds)
new_flux(p,ro,et,imax,jmax,U,V,J,choice,Yfin,Aplus,Aminus,Bplus,Bminus,gamma,Eplus,Eminus,Fplus,Fminus,Q,C,H,Hbar,Qbar)
newtons_poly(x,y)
newtons_poly_curve(condb,bound,l,coef,split_bottom,split_top,X2,X1lc,X1uc,Y1lc,Y1uc)
nondimension(ro,T,p,U,V,et,Tinf,jmax,imax,AS,l,Xfin,Yfin,Uinf)
openbowshock2D_main(varargin)
output_data(U,V,p,ro,et,T,Xfin,Yfin,jmax,imax,M,roinf,pinf,Uinf)
recangular_grid(X,Y,a_grid_up,a_grid_down,X2,i,imaxpr,bound,imax,betta1,jmax,N,alpha)
redo(l,condb,bound,split_top,split_bottom,s_bot,s_top,coef,X1l,X1u,Y1l,Y1u)
reverse(Ubar,Vbar,pbar,robar,etbar,Tbar,jmax,ro,U,Tinf,imax,l,Xfin,Yfin,ASbar,M,Uinf,roinf)
surfacepoint(bound,l,imax,Y,X,choice,betta1,jmax,coef,Imax,hclust)
tangencyangle(X,Y,imax,choice)
tridiagonal_solver1(CAM,CAP,CA,CBM,CB,CBP,dQbar,Eplus,Eminus,Fplus,Fminus,alpha,Hbar,dt,dxi,deta,imax,jmax,dQbars,choice,alpha2,betta2,A,B,C,RHS,Y,J)
user_defined(bound,X1u,X1l,Y1u,Y1l,X2,condb,split_top,split_bottom,coef,l,func,point,choice)
vertical_clustering(betta1,jmax,imax,R,r)
write_output(p,U,V,et,T,ro,imax,jmax,Minf,roinf,AoA,Tinf,Xfin,Yfin)

Version	Published	Release Notes
1.1.2.0	26 Feb 2016		Download

Select a Web Site

Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .

(English)
(Deutsch)
(Français)

（简体中文）
(English)

You can also select a web site from the following list

How to Get Best Site Performance

Select the China site (in Chinese or English) for best site performance. Other MathWorks country sites are not optimized for visits from your location.

Americas

América Latina (Español)
Canada (English)
United States (English)

Europe

Belgium (English)
Denmark (English)
Deutschland (Deutsch)
España (Español)
Finland (English)
France (Français)
Ireland (English)
Italia (Italiano)
Luxembourg (English)

Netherlands (English)
Norway (English)
Österreich (Deutsch)
Portugal (English)
Sweden (English)
Switzerland
United Kingdom(English)

Asia Pacific

Australia (English)
India (English)
New Zealand (English)
中国
- 简体中文Chinese
- English
日本Japanese (日本語)
한국Korean (한국어)

Contact your local office