Документ взят из кэша поисковой машины. Адрес оригинального документа : http://acoustics.phys.msu.ru/teachers/lecture/matmodel.pdf Дата изменения: Thu Oct 18 14:53:30 2007 Дата индексирования: Mon Oct 1 20:15:36 2012 Кодировка:

Any statement below is my private opinion. This talk contains a lot of controversial point
1st question: Does the Engineering need for Mathematics today? 2nd : If yes, what kind of Mathematics seems to be preferable? Well-known statement: "Some area of knowledge can be

named as a SCIENCE if and only if this area applies the MATHEMATICAL LANGUAGE"
Is this LANGUAGE necessary for any Engineer to communicate with his colleagues? Yes, but that is only a minor constituent. Basic statement: "Mathematics is much clever, than

Mathematician"
It means, that after the scientist or engineer developed new mathematical model and designed corresponding equations, the solution to these equations produces new information unpredictable for the author! Consider examples


To derive equations, one could say only 3 sentences: 1. Conservation of mass takes place in the volume where the motion is localized: m = d V = const or in differential form: r + div u = 0 t 2. Conservation of linear momentum takes place in the same volume: r r P = u dV = const or in differential form: ( p ik + uiuk + ik ) = 0 ui +
t x
k

Example 1: Fluid dynamics

3.

Viscosity leads to absorption coordinates and the motion is the fluid. From this sentences ik =

of energy, if velocity depends on different from steady-state rotation of we derive the viscous stress tensor: ul ul ui uk 2 + ik - ik + xl xl xk xi 3

Solutions to these equations describe thousands of new phenomena completely unknown at the derivation of this mathematical model! These equations are used to calculate operation regimes of aircrafts, ships, rockets, industrial flows, weather forecast, etc., etc...


Example 2: Maxwell's equations

1.Coulomb's law 2.Biot-Savart-Laplace law 3.Faraday law In addition, the existence of displacement current (alternating current can flow through the capacity) and the lack of magnetic charge are taken into account.

r r 4 div D = 4 , rot H = c r r 1 div B = 0, rot E = - c

r r 1 D j+ c t r B . t

Only a few results of experiments offer the possibility to derive these equations:

Interesting: Maxwell derived these equations using a primitive mechanical model: system of rotating gears with liquid between them
Strikingly! These primitive ideas led to universal equations describing everything: the operation of How one can study methods of Radio, TV, Mobile phones, properties of new design of new mathematical engineering materials, etc.etc. models and get the unique experience? I can demonstrate the simple exercise


A simplest humanitarian example
Philosophy of a poor painter: The more money, the less happiness Solution to this system:

M

dH = -aM , dt

dM = +bH dt

t
debts Solution to this system:

Periods of Comfortable life alternate with distress

Philosophy of a businessman: More money - more happiness

M

dH = + aM , dt

dM = +bH dt

t
laid-down capital

Amount of money increases exponentially!


Simple example of design of mathematical model (REAL Mathematical Modeling)
Simplest problem of demography: variation of male and female population with account of birth and death only.

dM M = +bMF - tM dt Increases with
Temporal change in number of Males increase in frequency of intimate relations between M and F Decreases because of mortality

dF M = + gMF - dt tF
The same ideas lead to the 2nd equation for Females

Notations: M ­ number of Males, F ­ number of Females, t M , t F -corresponding life intervals, b, g = (1 - b) -probabilities of birth of boys and girls. The simplest solution is the stationary one, when M and F are constant. Requirement M=F leads to simple formula: tF If life interval for women is 70 years and for men 50 years b= t F + t M (Russia), we calculate b =7/12. It means, that from among 12 of born babies, 7 must be boys and 5 must be girls


A little of my own experience in real mathematical modeling
1. Parametric underwater sonar (echo depth-sounding, fishery, mine location in marine sediments) was produced industrially after the engineering design method based on nonlinear acoustic equations was developed (See details in the book: B.K.Novikov, O.V.Rudenko, V.I.Timoshenko "Nonlinear Underwater Acoustics", American Inst.of Physics, 1987 ­ translated from Russian ed., 1981)



p p b 2 p c z - c 3 p - 2c 3 2 = 2 p


Signal transmission from air to underwater
V.E.Gusev, A.A.Karabutov. Laser optoacoustics. New York, Academic Press, 1994

Laser Beam

Linear and nonlinear (for high-intensity laser beam) optics equations Thermodynamic nonlinear equations describing light-sound transformation through thermal expansion, evaporation, etc

Light-sound transformation area
Sound wave

Acoustic equations


Modulated radiation pressure for shear wave excitation
Shear wave
Ultrasonic pulse

Palpation to detect Tumor: is that region of increased density?

No, density of tumor is the same as density of surrounding tissue Application to medical diagnostics: Shear Wave Elasticity Imaging Sarvazyan and Rudenko. US Patent 5,810,731 (Sept.22, 1998)


Mathematical model for «Shear Wave Elasticity Imaging»
Equation for nonlinear ultrasonic beam in tissue



p p 2 -3p - x c 2


0



c ( ) p( x , - )d = p 2

r2 r2 p( x = 0, r , t ) = p 0 2 ( t ) sin t + 2cd a
Averaging solution gives Radiation Pressure

Boundary condition on focusing transducer

Fx =

b p 3 c

2

This pressure is pulsating because US is modulated in kHz frequencies. This vibration excites the shear wave in tissue

2 sx 2 - (ct + ) s x = Fx 2 t t


Studies in Fluid Dynamics of water jet are connected with two main problems: 1. Flows in contracting nozzle forming high-speed jet. 2. Instability of water jet in air. Results are briefly described below.
Instability leads to breakdown of jet and to catastrophic decrease in its cutting ability.
Several physical factors were taken into account: (i) capillary forces at water free surface determined by both surface tension of pure water and polymer fibers added artificially to improve stability; (ii) aerodynamic forces appearing at the streamlining of jet irregularities.
Length of instability Nozzle

It was shown, that jet velocities around sound speed in air (330 m/s) are very undesirable, because instability increases explosively due to the wave resonance phenomenon. Next undesirable factor is the presence of air bubbles in liquid.

For instability length the solution is derived:
L-1 = ins I1 (k R) 1 U L I 0 (k R)
2 4

(a - c)2

b 1 + b2 sin arctan , a -c 2
2

U 2 J 0 J1 + Y0Y1 a = Ak , J12 + Y12 c=

b = Ak

U

2 1 , k R J12 + Y12

2

k
R
2

(k

2

R2 -1

)


To form the high-speed jet having velocities in air from 200 to 900 m/s, it is necessary to accelerate liquid in contracting nozzle by difference in pressure between input and output cross-sections. This fast (nonlinear) flow is described rr r r by equations:

(u ) u + u = -( p / 0 ), div u= 0

Exact solution depending on 3 constants was derived: 2 2 2 - a 2 + 2 + C3 1 + D a + a2 + 2 1 + uz = W , W = 2 2 2 2 2 - a + + C3 1 + z a2 + 2 2 1+
2 2 1 W 2 a - 2 1 + - C3 2 + ur = z 2 1 + 2 + C3 2
1/ 2

2

(

)

(

2

)



, =

1 0.5



1, 2

=

p0 0 (r - r
2 2

4 QL2

22 1

)

r22 - r12 L ln 1 + 4L2 z1, 2 4L2 z1, 2 1+ 2 2 r2 - r1 L

r z

This exact solution gives optimum shape of nozzle and formulas connecting main geometrical parameters of nozzle and water flow: 2

Q = r12 p00

r22 - r12 - 2L



0.2 An important equation exists connecting the total mass stream (outlet

z

discharge) and geometrical parameters of the nozzle. Strong restriction on dimensions:

2 Lr1 > r22 - r12

It means that nozzle length must be large enough. For example, if output radius is 1mm, and input radius is 1 cm, the nozzle length must be not less than 5 cm. Otherwise the stream considered above cannot be realized. This restriction was discovered experimentally by Water Jet Sweden and explained by this theory.


MATHEMATICS PHYSICS ENGINEERING Real Mathematical Modeling:
Level 1: Level 2: Designing of New equations to describe processes in Physics and Engineering

Exact Solution: the best way to extract
information from the mathematical model

Approximate calculations and asymptotic Level 3: methods ­ are used if mathematicians cannot
find the exact solution Level 4: Computer simulation (incorrect interpretation of "Mathematical Modeling") is used to solve complicated problems after main phenomena are studied at Level 2 and Level 3. Direct computer analysis of existing Models using unknown software leads often to principal mistakes