I am working in matlab..this code is regarding
atomization of water which actually used in
evaporation using spray technology...My droplets are
of order 500 microns which reduce upto 10
microns......I am using ODE45 to solve an differential
equation...where the integration has to stop when
diameter reaches to zero are less then zero..I am
using event option 4 that purpose but its not
working.....I am getting an error msg saying 'odezero:
an event disappeared (internal error)'..I am sending
over view of my code......can anybody help me in
solving my problem..I will be very thankful if anyone
can help me....
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%hollow_cone.m%%%%%%%%%%%%%%%%%%%%%%%%
%%%% input parameters
%%% constants...
options =
odeset('RelTol',1e-3,'AbsTol',10e-8,'Events',@events_all,'Refine',4);
 [t,y,TE,YE,IE] = ode45(@atomization_temp, [0 2.4e-3],
[0.2e-3 0], options, nuv_w,...);

%%%%%%%%%%%%%%%%%%%%%%% atomization_temp.m %%%%%%%%%%
function [adot] = atomization_temp(t,y,nuv_w,...)
    
%%%calculations.......
if r <= y(1)
ad1 =-(y(1) - r)/tuv; %differential equation for
droplet dia
                           % this drop should not go
-ve
else
    ad1 = (y(1) - r)/tuv;
end
%%%%%%%%calculations....

W_f %%%velocity
adot =[ad1; W_f];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%
events_all%%%%%%%%%%%%%%%%%%%%
function [value, isterminal, direction] =
events_all(t,y,...

 value = abs(y(1)) + y(1) ; %Detect diameter = 0;
                                  
isterminal = 1; % Stop the integration

direction = -1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Thank you
hari...

iitm