Wärmeübergang bei erzwungener turbulenter Strömung in Ringspalten und örtlich beliebig veränderlichen Randbedingungen zweiter Art

Zusammenfassung Das Problem des Wärmeübergangs bei turbulenter Strömung in konzentrischen Ringspalten wird für den dreidimensionalen Fall theoretisch gelöst, wobei die Wandwärmestromdichte sowohl in azimutaler als auch in axialer Richtung beliebig variiert. Die Lösung der Energiegleichung erfolgt mit der klassischen Methode der Superposition und Trennung der Variablen, wobei das dabei auftretende Sturm-Liouvillesche Eigenwertproblem numerisch gelöst wird. Zur Lösung werden Verteilungen für die Geschwindigkeit und anisotropen turbulenten Austauschgrößen verwendet, die mit dem phänomenlogischen Turbulenzmodell von Ramm berechnet wurden. Ergebnisse werden über einen weiten Bereich der Reynolds-Zahl (10⁴ Re 10⁶), der Prandtl-Zahl (0 Pr 100) und für verschiedene Radienverhältnisse diskutiert.

---

Turbulent forced convection heat transfer in annuli with arbitrarily varying boundary conditions of second kind

The problem of turbulent flow heat transfer in concentric annuli is analysed for the general threedimensional case in which the wall heat flux varies arbitrarily in both the circumferential and axial directions. The energy equation is solved using the classical method of superposition and separating variables, where the resulting Sturm-Liouville problem are evaluated numerically. The solution is based on velocity profiles and anisotropic thermal turbulent transport properties evaluated by Ramm's phenomenological turbulence model. Results are discussed over a wide range of Reynolds number (10⁴ Re 10⁶), Prandtl number (0 Pr 100) and radius ratio.

---

Bezeichnungen a,b Fourierkoeffizienten - B geometrische Funktion, [s(1-r) + r]/(1–s) - C Koeffizienten - D hydraulischer Durchmesser, 2(r₂ – r₁) - E Energietransportfunktion - f axiale Wärmestromdichteverteilung - F azimutale Wärmestromdichteverteilung - g radiale Temperaturfunktion - l Kanallänge - L dimensionslose Kanallänge, 1/D - M axialer Temperaturgradient im thermisch ausgebildeten Bereich - n harmonischer Parameter - Nu Nusselt-Zahl - Pe Péclet-Zahl - Pr Prandtl-Zahl - q Wärmestromdichte - Q dimensionslose Wärmestromdichte, q/q₀ - r dimensionslose radiale Koordinate, (R-r₁)/(r₂-r₁) - r₁,r₂ innerer und äußerer Ringspaltradius - R radiale Koordinate - Re Reynolds-Zahl - s Ringspaltverhältnis, r₁/r₂ - T dimensionslose Temperatur, 2· · (-_E/(D· q₀ - u dimensionslose Geschwindigkeit, U/U_m - U Geschwindigkeit - x dimensionslose axiale Koordinate, X/D - X axiale Koordinate - Wärmeübergangskoeffizient - _un modifizierter Eigenwert - halber Segmentwinkel - turbulente Austauschgröe - Temperatur - dimensionslose Temperaturdifferenz, T - T_m - Wärmeleitfähigkeit - _un Eigenwerte - kinematische Viskosität - azimutale Koordinate - Eigenfunktionen Indizes e thermischer Einlauf - E Eintritt bei x=0 - H Wärme - i Bedingung an der i-ten benetzten Oberfläche (i=1 – Innenrohr, i=2 - Außenrohr) - j Bedingung, wenn nur an der j-ten Oberfläche des Ringspaltes die Wärme übertragen wird (j=1,2) - ij Bedingung an der i-ten Oberfläche, wenn nur an der j-ten Oberfläche des Ringspaltes die Wärme übertragen wird (ij=11, 12, 22, 21) - m mittel - n Ordnung der Harmonischen - r radiale Richtung - u Ordnung des Eigenwertproblems - azimutale Richtung - 0 umfangskonstant - thermisch ausgebildet     

Thermal entrance region heat transfer for MHD laminar flow in parallel-plate channels with unequal wall temperatures

The problem of thermal entry heat transfer for Hartmann flow in parallel-plate channels with uniform but unequal wall temperatures considering viscous dissipation, Joule heating and axial conduction effects is approached by the eigenfunction expansion method. The series expansion coefficients for the nonorthogonal eigenfunctions are obtained by using a method for nonorthogonal series described by Kantorovich and Krylov [21]. Numerical results are obtained for the case with entrance condition parameter _o=1 and open circuit condition K=1. The parametric values of Ha=0, 2, 6, 10 and Br=0, –1 are considered for Hartmann and Brinkman numbers, respectively.

---

Zusammenfassung Das Problem der Wärmeübertragung im thermischen Einlauf einer Hartmannströmung im ebenen Spalt mit einheitlichen, aber ungleichen Wandtemperaturen wurde unter Berücksichtigung viskoser Dissipation, Joulescher Heizung und axialer Wärmeleitung mit Hilfe einer Entwicklung nach Eigenfunktionen behandelt. Die Koeffizienten der Entwicklung für nichtorthogonale Eigenfunctionen wurde nach einer Methode für nichtorthogonale Reihen nach Kantorovicz und Krylow [21] berechnet. Numerische Ergebnisse werden für den Eintrittsparameter _o=1 und die Bedingung für den offenen Stromkreis K=1 erhalten. Die Parameterwerte Ha=0, 2, 6, 10 und Br=0, –1 werden für die jeweiligen Werte der Hartmann- und der Brinckman-Zahl betrachtet.

---

Nomenclature a one-half of channel height - ¯B,B₀ magnetic field Induction vector and magnitude of applied magnetic field - Br Brinkman number, _f U_m ²/(k_c) - C_n,D_n coefficients in the series expansion of _e, see eq. 16 - c_p specific heat at constant pressure - ,E₀ electric field intensity vector and component - E_n,O_n even and odd eigenfunctions - Ha Hartmann number, (/_f)^1/2 B_o a - h₁,h₂ local heat transfer coefficients at lower and upper plates - ¯J,J_y electric current density vector and component - K external loading parameter, E_o/(B_o U_m) - k thermal conductivity - Nu₁, Nu₂ local Nusselt numbers, h₁,a/k and h₂a/k, respectively - P fluid pressure - Pe Peclet number, PrRe - Pr Prandtl number, C_{p f}/k - q₁,q₂ rates of heat transfer per unit area,^–k(T/Z)Z=–a'^–k(T/Z) Z=a respectively - Re Reynolds number, U_ma/u_f - T,T₀,T₁,T₂ fluid temperature, uniform entrance temperature, uniform but different lower and upper plate temperatures, respectively - T_b,T_m bulk temperature and (T₁+T₂)/2 - U,U_m,u axial, mean and dimensionless velocities, respectively - ¯V velocity vector - X,Z axial and transverse coordinates - x,z dimensionless coordinates - _n,_n even and odd eigenvalues - ,₀,_b dimensionless fluid, entrance and bulk temperatures, respectively - _c,_e,_f characteristic temperature difference (T₂-T_m), and dimensionless fluid temperatures, defined by eq. (10) - _e,_f magnetic permeability and viscosity of fluid - fluid density - electric conductivity - viscous dissipation function - (1-)/2     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

Fixed point theorem of nonexpansive mappings in convex metric spaces

Let X be a convex metric space with the property that every decreasing sequence of nonenply dosed subsets of X with diameters tending to has menemptyintersection. This paper proved that if T is a mapping of a elosed conver nonempty subset K of X into itself satisfying the inequality:d(T_x,T_y)≤ad(x,t)+b{d(x,T_x)+d(y,T_y)}+c{d(x,T_x)+d(y,T_y)}for all x,y in K,where 0≤a<1,b≥0,c≥0,a+c≠0 and a+2b+3c≤1, then T has a unique fixed point in K.     

Viscosity of concentrated protein solutions

The viscosity of solutions of four proteins (Bovine Serum Albumin, Ovalbumin, _s-1 Casein, Lysozyme), brought to the random coil conformation, has been measured over a large concentration range extending into the entanglement region. A master curve is obtained in the dilute and semi-dilute regions with the reduced variables and of Simha and Utracki.By using Graessley's expression for the polymer coil expansion at a given concentration in the semi-dilute region (c ^* c c ^**), a simple equation is established giving the relative viscosity _r as a function of concentrationc: forc ^* c c ^**, ln _r = 2a[]c ^*(c/c ^*)^1/2a – (2a - 1)[]c ^*; wherec ^* is the incipient overlap concentration, [] the intrinsic viscosity, anda the Mark-Houwink exponent for the polymer-solvent considered.This equation fits well the experimental results. The adjustment yields for the parametera values which are comprised between 0.6 and 0.7, as expected, for Bovine Serum Albumin and Ovalbumin, but very close to 0.5 for _s-1 Casein and Lysozyme. This can be explained by the fact that the molecular weights of the two latter proteins are lower than, or very close, the critical molecular weight; the critical molecular weight is estimated to be about 20000.     

Dynamic crack resistance of acrylic plastic

The dependence of the stressintensity factor near the tip of a growing crack in an SO120 acrylic plastic on the crackpropagation velocity K_I(:v) within the range of 10^–5 –300 m/sec is determined. Specific features of crack propagation associated with the shape of the curve K_I(v), which has discontinuities and nonuniqueness intervals, are discussed.     

Recovery-Faktor des frontal angeströmten zylindrischen Mantelthermoelementes mit ebener Stirnfläche

Zusammenfassung Zur Messung schnell veränderlicher Fluidtemperaturen mit Hilfe von Berührungsthermometern sind im allgemeinen hohe Anströmgeschwindigkeiten der Fühler notwendig, um gute Wärmeübergangsverhältnisse sicherzustellen. Dadurch erhält der Recovery-Faktor dieser Fühler besondere Bedeutung. Haupteinflußparameter sind hier die Reynolds- und die Machzahl.Es wird der Recovery-Faktor eines Mantelthermoelementes von 0,5 mm Durchmesser in Wasserdampfströmung in Abhängigkeit von der Anströmmachzahl und der mit dem Manteldurchmesser gebildeten Reynoldszahl auf experimentellem Weg ermittelt. Untersucht werden die Bereiche 0,4 M 0,8 und 3000 Re 9000. Bei niedrigen Machzahlen und hohen Reynoldszahlen sind die Recovery-Werte klein. Sie steigen bis zu einer Grenz-Machzahl mit wachsender Anströmgeschwindigkeit und bei sinkender Reynoldszahl. Oberhalb dieser Grenze besteht keine Abhängigkeit von der Machzahl. Das Recovery-Verbalten des Fühlers wird entscheidend durch die Ablösezone bestimmt, die an der Fühlerspitze auftritt.

---

Recovery factor of a cylindrical thermocouple with plane frontal area in normal flow

Generally, high flow speeds are required to assure adequate heat transfer conditions in determining rapidly changing fluid temperatures by means of contact thermometers. Consequently, the recovery factor of these temperature gages is of predominant importance with Reynolds- and flow Mach numbers as main influence parameters.In the present study, the recovery factor of an enclosed thermocouple of 0,5 mm diameter is determined in steam as a function of free stream Mach number and Reynolds number based on the tip diameter. The Mach number range 0,4 M 0,8 and Reynolds numbers 3000 Re 9000 were considered respectively. It was found that the recovery factor is relatively small with low Mach numbers and high Reynolds numbers. Increasing the Mach number and lowering the Reynolds number yields higher recovery factors. After reaching a limiting value, no effect of the flow speed is observed and the recovery factor remains constant with increasing Mach number. Of major influence is the separation zone at the tip of the thermocouple.

---

Bezeichnungen A effektiver Strömungsquerschnitt (Bild 1) - A^* effektiver engster Strömungsquerschnitt (Bild 1) - a örtliche Schallgeschwindigkeit - c spezifische Wärme - d Durchmesser des Mantelthermoelementes - d₁ Durchmesser der Thermoelement-Halterung in der Ebene A^* - d₂ Innendurchmesser der Düse, zylindrischer Teil - L charakteristische Länge, Lauflänge - M örtliche Machzahl - M^* kritische Machzahl - Pr c_p/ - Re Reynoldszahl (Gl.(7)) - Re_o a_o l/_o (Index o: gebildet mit Total werten) - r Recovery-Faktor (Gl. (2)) - T_A Referenz-Temperatur: 293,1 K - T_s statische Temperatur - T_e Eigentemperatur - T_t Totaltemperatur - u Anströmgeschwindigkeit - x Längenkoordinate, axiale Richtung - _t thermischer Ausdehnungskoeffizient des Düsenmaterials - Verdrängungsdicke der Düsengrenzschicht im Querschnitt A - Verdrängungsdicke der Düsengrenzschicht im Querschnitt A^* - J Temperaturdifferenz (Gl. (1)) - J_i i =1, 2, 3, 4 gemessene Temperaturen - Adiabatenexponent - kinematische Zähigkeit     

Flow of a viscous fluid in the initial section of channels with intense injection

Problems of flows in the initial sections of plane, circular and annular channels are solved by numerical integration of the Navier-Stokes equations on the interval 10 R_b, 3000. The initial section is considered to be the part of the flow where the local Reynolds number does not exceed the critical value.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 187–190, March–April, 1987.     

Rotating bunsen flames as solutions of the Kuramoto-Sivashinsky equation

Steadily rotating solutions of the Kuramoto-Sivashinsky equationu _t + ² u++¦u¦ ² =c ² are studied. These solutions bifurcate from the steady radial solution of the above equation. For large values ofc and angular velocities such thatN¦¦<2c<(N+1)¦¦, we show that there exists a 2N-1 family of bifurcating solutions. The proof is based on a certain generic transversality assumption. A computer-assisted proof of this assumption is given for 1N10.     

Thermal stability of composite superconducting tape under the effect of a two-dimensional dual-phase-lag heat conduction model

Thermal stability of composite superconducting tape subjected to a thermal disturbance is numerically investigated under the effect of a two-dimensional dual-phase-lag heat conduction model. It is found that the dual-phase-lag model predicts a wider stable region as compared to the predictions of the parabolic and the hyperbolic heat conduction models. The effects of different design, geometrical and operating conditions on superconducting tape thermal stability were also studied.a conductor width, (m) - A conductor cross sectional area of, (m²) - As conductor aspect ratio, (a/b) - b conductor thickness, (m) - Bi Biot number - B dimensionless disturbance Intensity - C heat capacity, (J m^–3 K^–1) - D disturbance energy density, (W m^–3) - f volume fraction of the stabilizer in the conductor - g(T) steady capacity of the Ohmic heat source, (W m^–3) - g_max Ohmic heat generation with the whole current in the stabilizer, (W m^–3) - G_max dimensionless maximum Joule heating - h convective heat transfer coefficient, (W m^–2 K^–1) - J current density, (A m^–2) - k thermal conductivity of conductor, (W m^–1 K^–1) - q conduction heat flux vector, (W m^–2) - Q dimensionless Joule heating - R relaxation times ratio (_T/2_q) - t rime, (s) - T temperature, (K) - T_c critical temperature, (K) - T_c1 current sharing temperature, (K) - T_i initial temperature, (K) - T_o ambient temperature, (K) - x, y co-ordinate defined in Fig. 1, (m) - thermal diffusivity (m² s^–1) - dimensionless time - _i dimensionless duration time - dimensionless y-variable - _o superconductor dimensionless thickness - dimensionless temperature - _c1 dimensionless current sharing temperature - 1 dimensionless maximum temperature - dimensionless disturbance energy - numerical tolerance - x width of conductor subjected to heat disturbances, (m) - y thickness of conductor subjected to heat disturbances, (m) - dimensionless x-variable - _o superconductor dimensionless width - stabilizer electrical resistivity, () - _q relaxation time of heat flux, (s) - _T relaxation time of temperature gradient, (s) - i initial - sc current sharing - max maximum - o ambient     

Untersuchung des lokalen Stoffüberganges am laminaren und turbulenten Rieselfilm

Zusammenfassung Der lokale Stoffübergang wurde in Abhängigkeit von der Meßlänge, dem Startort und der Zulaufhöhe gemessen. Der Gültigkeitsbereich der Theorie von Nusselt wird ermittelt. Die Reynolds-Zahl nahm Werte zwischen 3,86 und 2496 an. Die örtlich wirkende Hydrodynamik ist entscheidend für das Anwachsen der örtlichen Sherwood-Zahl. Die Genauigkeit aller Versuchsergebnisse kann auf ± 5% abgeschätzt werden.

---

Investigation of the local mass transfer of a laminar and turbulent falling liquid film

The local mass transfer was measured as a function of the measuring length, the starting point and the liquid height above the ring-slot. The range of the Reynolds number was 3,86 Re 2496. The validity of the Nusselt theory and the range of it is shown. The local hydrodynamic is the most important factor of the increase of the local Sherwood number. The accuracy of the measurements is ± 5%.

---

Bezeichnungen a Temperaturleitfähigkeit m²/s=/(c_p) - c Konzentration, c=¯c + c kmol/m³ - c_i0 Konzentration im Flüssigkeitskern kmol/m³ - D Diffusionskoeffizient m²/s - EL-NR Elektrodennummer - Fa Faraday-Konstante A s/kgäq=96,5·10⁶ - g Erdbeschleunigung m/s² - i_G Grenzstromdichte A/m² - u Geschwindigkeit in x-Richtung, u= + u - U Umfang des Rohres m - v Geschwindigkeit in y-Rich- m/stung, v=¯v + v - V^* Volumenstrom m³/s - x Lauflänge, Koordinate in m Strömungsrichtung - x_M Meßlänge für den Stoff-Übergang m - x_ST Startort für den Stoff-Übergang m - y Wegkoordinate senkrecht zur Rohroberfläche m - z Wertigkeit der Elektro-denreaktion kgäq/kmol - ZH Zulaufhöhe m - Wärmeübergangskoeffizient W/m²C - Stoffübergangskoeffizient m/s - Filmdicke m - Wärmeleitfähigkeit W/(mC) - kinematische Viskosität m²/s - Re=u/=V^*/U Reynolds-Zahl - Pr=/a=c_p/ Prandtl-Zahl - Sc=/D Schmidt-Zahl - Nu= / Nusselt-Zahl - Sh= /D Sherwood-Zahl - SHL lokale Sherwood-Zahl - SHM mittlere Sherwood-Zahl - - zeitlich gemittelt - örtlich gemittelt Die Durchführung der Arbeit am Institut für Verfahrens — und Kältetechnik der ETH Zürich bei Prof. Dr. P. Grassmann wurde ermöglicht durch Zuschüsse der Kommission zur Förderung der wissenschaftlichen Forschung und meiner Eltern.     

Measurement of the Coefficient of Transverse Dispersion in Flow Through Packed Beds for a Wide Range of Values of the Schmidt Number

Experimental values of the coefficient of transverse dispersion (D _T) were measured with the system 2-naphthol/water, over a range of temperatures between 293K and 373K, which corresponds to a range of values of viscosity () between 2.83×10^–4 Ns/m² and 1.01×10^–3 Ns/m² and of molecular diffusion coefficient (D _m) between 1.03×10^–9 m²/s and 5.49×10^–9 m²/s. Since the density () of water is close to 10³ kg/m³, the corresponding variation of the Schmidt number (Sc=/D _m) was in the range 1000 – 50. More than 200 experimental values of the transverse dispersion coefficient were obtained using beds of silica sand with average particle sizes (d) of 0.297 and 0.496mm, operated over a range of interstitial liquid velocities (u) between 0.1mm/s and 14mm/s and this gave a variation of the Reynolds number (Re=du/) between 0.01 and 3.5.Plots of the dimensionless coefficient of transverse dispersion (D _T/D _m) vs. the Peclet number (Pe_m=ud/D _m) based on molecular diffusion bring into evidence the influence of Sc on transverse dispersion. As the temperature is increased, the value of Sc decreases and the values of D _T/D _m gradually approach the line corresponding to gas behaviour (i.e. Sc 1), which is known to be well approximated by the equation D _T/D _m=1/+ud/12D _m, where is the tortuosity with regard to diffusion.     

Untersuchungen des Blasensiedens einiger verflüssigter Gase und ihrer binÄren Mischungen

Zusammenfassung Die Arbeit behandelt das Blasensieden an einer ebenen, horizontalen, 64·10^–4m² gro\en KupferoberflÄche von Stickstoff, Methan, Aethan und Gemischen aus Stickstoff/Methan und Methan/Aethan bei verschiedenen Drücken. Die Messwerte und empirischen Ausgleichskurven sind angegeben. Die Daten für die binÄren Gemische bestÄtigen die Gleichung von Happel und Stephan.

---

Heat transfer in nucleate boiling of liquified gases and their binary mixtures

The paper deals with pool boiling of nitrogen, methane, ethane and mixtures of nitrogen/methane and methane/ethane at different pressures. The boiling surface was a horizontal, plan copper disk of 64·10^–4m²· The measured data points and their fit by an empirical correlation are given. For mixtures the correlation of Happel and Stephan provides a good agreement with the results.

---

Formelzeichen a₁, a₂, a₃ Konstanten - b₁,b₂ Konstanten - K Konstante - n,n₁,n₂ Konstanten - p Druck Pa, M Pa - p_red reduzierter Druck (p/p_kritisch) - q^* WÄrmestromdichte W/m² - R_p Rauhtiefe m - x Molanteil der leichter flüchtigen Komponente in der Flüssigkeit - y Molanteil der leichter flüchtigen Komponente im Dampf - WÄrmeübergangskoeffi-W m^–2 K^–1zient - _i idealer WÄrmeüber-W m^–2K^–1 - _i gangskoeffizient - realer WÄrmeübergangs-W m^–2K^–1 - _r koeffizient - T Temperaturdifferenz K     

Classical limits for the coefficient of variation for the normal distribution

The exact classical limits for the coefficient of variation c for the normal distribution are derived. The hand-calculating approximated classical limits for c having high accuracy are given to meet practical engineering needs. Using Odeh and Owen's computational method and Brent's algorithm, the tables for the r-upper exact classical limits of coefficient of variation for normal distribution are calculated for the different confidence coefficient , the sample size n=1 (1)30,40,60,120, the sample coefficient of variation =0.01(0.01)0.20. It is shown that if n8, 0.20, then the -upper exact classical limits c_u for c are slightly higher than the exact fiducial limits c_u,F for c. if n>8, 0.20, then c_u–c_u,F<5×10^–6.     

Comparison of network theory predictions with stress/time data in shear and elongation for a low-density polyethylene melt

Summary Experimental data are presented which show the variation with time of the shear stress and primary normal stress difference during shear flow with a stepfunction shear rate; the material (Melt I at 150 °C) is a low-density polyethylene melt for which stress-growth and elastic recovery data in elongational flow experiments have been previously reported. A method of comparing the data with the predictions of the rubberlike-liquid constitutive equations is given, based on the use of a specially-chosen rate-of-strain invariantI _IV, defined in [4.1]. From this comparison, it is shown that the disagreement between theory and experiment is about the same for shear flow and for elongational flow, and that the extent of disagreement does not exceed 10% for short-duration flows such thatI _IV t 3.

---

Zusammenfassung Es werden Meßergebnisse über die zeitliche Änderung der Schubspannung und der ersten Normalspannungsdifferenz bei Scherfließen einer LDPE-Schmelze (Schmelze I bei 150 °C) vorgelegt. Der zeitliche Spannungsverlauf bei und die elastische Erholung nach Dehnfließen sind für dieses Material bereits früher mitgeteilt worden. Hier wird das Verhalten der Schmelze bei Scherung und bei Dehnung mit den Voraussagen der rubberlike liquid-Zustandsgleichung verglichen, wobei eine speziell gewählte InvarianteI _IV der Deformationsgeschwindigkeit verwendet wird (definiert in [4.1]). Der Vergleich zeigt Abweichungen von Theorie und Experiment, die für Scher- und Dehnfließen etwa gleich groß sind. Die Abweichungen liegen unter 10%, wenn für das Produkt ausI _IV und der Deformationszeitt der WertI _IV t = 3 nicht überschritten wird.

---

---

With 7 figures and 2 tables     

Flow birefringence of polymer melts: Application to the investigation of time dependent rheological properties

Summary Transient stresses including normal stresses, which are developed in a polymer melt by a suddenly imposed constant rate of shear, are investigated by mechanical measurement and, indirectly, with the aid of the flow birefringence technique. For the latter purpose use is made of the so-called stress-optical law, which is carefully checked.It appears that the essentially linear model of the rubberlike liquid, as proposed byLodge, is capable of describing the behaviour of polymer melts rather well, if the applied total shear does not exceed unity. In order to describe also steady state values of the stresses successfully, one should extend measurements to extremely low shear rates.These statements are verified with the aid of a method which was originally designed bySchwarzl andStruik for the practical calculation of interrelations between linear viscoelastic functions. In the present paper dynamic shear moduli are used as reference functions.

---

Zusammenfassung Mit der Zeit anwachsende Spannungen, darunter auch Normalspannungen, wie sie sich nach dem plötzlichen Anlegen einer konstanten Schergeschwindigkeit in einer Polymerschmelze entwickeln, werden mit Hilfe mechanischer Messungen und indirekt mit Hilfe der Strömungsdoppelbrechung untersucht. Für den letzteren Zweck wird das sogenannte spannungsoptische Gesetz herangezogen, dessen Gültigkeit sorgfältig überprüft wird.Es ergibt sich, daß das im Wesen lineare Modell der gummiartigen Flüssigkeit, wie es vonLodge vorgeschlagen wurde, sich recht gut zur Beschreibung des Verhaltens von Polymerschmelzen eignet, solange der im ganzen angelegte Schub den Wert Eins nicht überschreitet. Um auch stationäre Werte der Spannungen in die Beschreibung erfolgreich einzubeziehen, sollte man die Messungen bis zu extrem niedrigen Schergeschwindigkeiten ausdehnen.Die gemachten Feststellungen werden mit Hilfe einer Methode verifiziert, die vonSchwarzl undStruik ursprünglich für die praktische Berechnung von Beziehungen zwischen Zustandsfunktionen entwickelt wurde, die dem linear viskoelastischen Verhalten entsprechen. In der vorliegenden Veröffentlichung dienen die dynamischen Schubmoduln als Bezugsfunktionen.

---

a _T shift factor - B _ij Finger deformation tensor - C stress-optical coefficient, (m²/N) - f (p _jl ) undetermined scalar function - G shear modulus, (N/m²) - G(t) time dependent shear modulus, (N/m²) - G() shear storage modulus, (N/m²) - G() shear loss modulus, (N/m²) - G _r reduced shear storage modulus, (N/m²) - G _r reduced shear loss modulus, (N/m²) - H() shear relaxation time spectrum, (N/m²) - k Boltzmann constant, (Nm/°K) - n _ik refractive index tensor - p undetermined hydrostatic pressure, (N/m²) - p _ij ,p _ik stress tensor, (N/m²) - p ₂₁ shear stress, (N/m²) - p ₁₁ –p ₂₂ first normal stress difference, (N/m²) - p ₂₂ –p ₃₃ second normal stress difference, (N/m²) - q shear rate, (s^–1) - t, t time, (s) - T absolute temperature, (°K) - T ₀ reference temperature, (°K) - x the ratiot/ - x position vector of a material point after deformation, (m) - x position vector of a material point before deformation, (m) - ₀, ₁ constants in eq. [37] - ₀, ₁ constants in eq. [37] - shear deformation - (t, t) time dependent shear deformation - _ij unity tensor - n flow birefringence in the 1–2 plane - (q) non-Newtonian shear viscosity, (N s/m²) - ^* () complex dynamic viscosity, (N s/m²) - | ^* ()| absolute value of complex dynamic viscosity, (N s/m²) - () real part of complex dynamic viscosity, (N s/m²) - () imaginary part of complex dynamic viscosity, (N s/m²) - (t — t) memory function, (N/m² · s) - v number of effective chains per unit of volume, (m^–3) - temperature dependent density, (kg/m³) - ₀ density at reference temperatureT ₀, (kg/m³) - relaxation time, (s) - integration variable, (s) - (x) approximate intensity function - ₁ (x) error function - extinction angle - _m orientation angle of the stress ellipsoid - circular frequency, (s^–1) - 1 direction of flow - 2 direction of the velocity gradient - 3 indifferent direction - t time dependence The present investigation has been carried out under the auspices of the Netherlands Organization for the Advancement of Pure Research (Z. W. O.).North Atlantic Treaty Organization Science Post Doctoral Fellow.Research Fellow, Delft University of Technology.With 11 figures and 2 tables     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Temperament Development Modeled as a Nonlinear Complex Adaptive System

Using approach-withdrawal (AW) as a specific instance of temperament, a theoretical model of temperament as a complex dynamic system is proposed. Developmental contextualism (Lerner, 1998) serves as a guiding theory in determining the structural components of the system and Kauffman's (1993) Boolean models of self-organization are adapted to estimate the parameter functions. In this model P(AW) = f(, ) where P(AW) is the probability density function of an approach or a withdrawal response, ( is a standardized parameter estimate of the biological sensitivity to stimulation, and is a standardized parameter estimate of the contextual response to an approach or withdrawal response. It is theorized that the functions of ( and follow a Hill function of the forms: d /dt = (²/c² + ²) – K₁ d /dt = ( ²/c² + ²) – K₂, where K₁, K₂, and c are system constants. This results in a double sigmoid function in which at extreme values of and the system stabilizes on a steady state of either approach or withdrawal response patterns. At intermediate parameter values the probability density functions of approach and withdrawal responses are wider. Thus, AW can be modeled as representing two basins of attraction. In addition, considerations are given to the systems sensitivity to initial conditions.     

Mixed convection flow over a vertical plate with localized heating (cooling), magnetic field and suction (injection)

An analysis is carried out to study the effects of localized heating (cooling), suction (injection), buoyancy forces and magnetic field for the mixed convection flow on a heated vertical plate. The localized heating or cooling introduces a finite discontinuity in the mathematical formulation of the problem and increases its complexity. In order to overcome this difficulty, a non-uniform distribution of wall temperature is taken at finite sections of the plate. The nonlinear coupled parabolic partial differential equations governing the flow have been solved by using an implicit finite-difference scheme. The effect of the localized heating or cooling is found to be very significant on the heat transfer, but its effect on the skin friction is comparatively small. The buoyancy, magnetic and suction parameters increase the skin friction and heat transfer. The positive buoyancy force (beyond a certain value) causes an overshoot in the velocity profiles.A mass transfer constant - B magnetic field - C_fx skin friction coefficient in the x-direction - C_p specific heat at constant pressure, kJ.kg^–1.K - C_v specific heat at constant volume, kJ.kg^–1.K^–1 - E electric field - g acceleration due to gravity, 9.81 m.s^–2 - Gr Grashof number - h heat transfer coefficient, W.m².K^–1 - Ha Hartmann number - k thermal conductivity, W.m^–1.K - L characteristic length, m - M magnetic parameter - Nu_x local Nusselt number - p pressure, Pa, N.m^–2 - Pr Prandtl number - q heat flux, W.m^–2 - Re Reynolds number - Re_m magnetic Reynolds number - T temperature, K - T_o constant plate temperature, K - u,v velocity components, m.s^–1 - V characteristic velocity, m.s^–1 - x,y Cartesian coordinates - thermal diffusivity, m².s^–1 - coefficient of thermal expansion, K^–1 - , transformed similarity variables - dynamic viscosity, kg.m^–1.s^–1 - ₀ magnetic permeability - kinematic viscosity, m².s^–1 - density, kg.m^–3 - buoyancy parameter - electrical conductivity - stream function, m².s^–1 - dimensionless constant - dimensionless temperature, K - w, conditions at the wall and at infinity     

K+ uptake by root systems grown in soil under salinity: II. Sensitivity analysis

A numerical algorithm for the solution of multicomponent transport of Ca²⁺, Mg²⁺, Na⁺, K⁺, Cl^– in soil and their uptake by plant roots has been developed. The model emphasizes adsorptiondesorption due to cation exchange mechanism, dissolution-precipitation of CaCO₃, and pH changes at the root surface controlled by the anion-cation influx balance. A fully implicit finite difference scheme is used for numerical implementation. Sensitivity analysis was conducted to evaluate the effect of each parameter on nutrient uptake. Each parameter (independent of all others) was varied between 0.25 to 4 times its speculated average level. Predicted K⁺ uptake was found to be more sensitive to changes of root radius and the parameter indicating maximal influx of K⁺. Effective diffusion coefficient and soil moisture are less influential. The influence of C_aCO₃ dissolution and different kinds of boundary conditions were also considered.Nomenclature A, B, E matrices of coefficients for finite difference scheme - b _i coefficients of equation for H⁺ concentration - C ^I concentration of theI-component in water - C ₀ ^I initial concentration of theI-component - _I ^r0 concentration of theI-component at the inner side of the root surface - C _I ^r1 concentration of theI-component at the external boundary - C _Na ^cr critical concentration for Na⁺ influx into root - CEC cation exchange capacity - D ^*I effective diffusion coefficient of theI-component in soil - F ^I concentration of theI-component on the exchange complex - G vector of coefficients in finite difference scheme - h Hill's cooperativity index for K⁺ influx - h ⁰ value ofh whenC ^Na=0 - J ^I uptake of theI-component by a root length unit - J _I ^r0 influx at the root surface of theI-component - J _max maximal influx of K⁺ - J _max ⁰ value ofJ ^max whenC ^Na=0 - K _a apparent Michaelis-Menten coefficient for K⁺ influx - K _a ⁰ value ofK _a whenC ^Na=0 - K _i selectivity coefficient of the exchange reaction - P _m ^I permeability of root surface for theI-component - P_CO2 partial pressure of CO₂ - r radial distance - r ₀ root radius - r ₁ half the distance between adjacent root (external boundary) - R ^I retardation factor of the-component in mass balance equation for theI-component - t time - t ₀ initial time - T simulation time - v ₀ water radial velocity at the root surface - x _i coordinate of nodes of finite difference mesh Greek coefficient of linear change of K⁺ influx - coefficient of linear change of Na⁺ influx - _s mass density of soil solid phase - soil porosity - volumetric content of liquid in soil - _i coefficients in formulae for parameters of K⁺ influx - parameter of perturbation in finite difference scheme - gg ^I activity coefficient of theI-component - accuracy of iteration convergence - time step for finite difference scheme - steps of finite difference mesh Special Symbols [...] activity symbol