Calculating Local Interface Curvature for Fluid Dynamics Simulation

AI Thread Summary
The discussion focuses on calculating local interface curvature for fluid dynamics simulations using the boundary element method to solve the Laplace equation. The mean interface curvature is defined as 2H=div n, where n is the unit normal to the surface. A challenge arises when differentiating the scalar values of the unit normal, leading to confusion about the divergence calculation. The solution involves deriving a parabolic equation from three points along the interface, allowing for the calculation of curvature using the second derivative. The participant expresses appreciation for the forum as a valuable resource despite initially not receiving responses.
sam_the_man
Messages
2
Reaction score
0

Homework Statement


I am using the boundary element method to solve unknowns to the Laplace equation from classic potential flow theory for the time evolution of a fluid air interface. At each time step, I need to solve a material derivative equation numerically at every node along an interface to find the new velocity potential.

In order to calculate the material derivative, I need to calculate the local interface curvature (et al.).

Homework Equations



From text, the local (mean) interface curvature can be calculated as 2H=div n. Where H is the mean interface curvature and n is the unit normal to the surface.

The Attempt at a Solution


The divergence of a vector field is a somewhat trivial calculation, e.g.:
div F = (dF1/dx+dF2/dy+dF3/z) where each value of F is some function that can be differentiated (pde) (like x*y^2). So here is the problem/question, in the case of the unit normal, the values are scalar values (such as [1 0 1]'), therefore if I differentiate each of these values with respect to the independent variable the entire equation equals zero. No doubt I am missing something fundamental here, any advise would be greatly appreciated.
I've attached a very simple sketch of a local discretized interface with nodes, and a unit normal just to help visualize what I am working with.
 

Attachments

  • Untitled.jpg
    Untitled.jpg
    2.3 KB · Views: 394
Last edited:
Physics news on Phys.org
I've figured this out and I'm going to answer my own question.

Rather than find the unit normal for the point, it is easier to find the equation for the "curve" of interest. In the diagram I attached, I could take 3 points two either side of where I have annotated the normal. From the location of these points I can formulate an equation for the curve from the parabolic equation y=ax^2+bx+c. Since I have three points, I can find the values for a b and c for example;

y_1=ax_1^2+bx_1+c
y_2=ax_2^2+bx_2+c
y_3=ax_3^2+bx_3+c

which can be solved numerous ways for a, b and c which in turn are plugged back into the original equation.

From this the curvature can be calculated as H(x)=abs(y")/[1+(y')^2]^(3/2) at any given point and for any curve equation.

Side note: I spend quite a bit of time on here as a "spectator" generally I stumble across this site while looking for little tidbits of info. Even though I didn't receive any response for my first try here, I would say keep up the good work. This is a great resource.
 
Last edited:
Thread 'Voltmeter readings for this circuit with switches'
TL;DR Summary: I would like to know the voltmeter readings on the two resistors separately in the picture in the following cases , When one of the keys is closed When both of them are opened (Knowing that the battery has negligible internal resistance) My thoughts for the first case , one of them must be 12 volt while the other is 0 The second case we'll I think both voltmeter readings should be 12 volt since they are both parallel to the battery and they involve the key within what the...
Thread 'Correct statement about a reservoir with an outlet pipe'
The answer to this question is statements (ii) and (iv) are correct. (i) This is FALSE because the speed of water in the tap is greater than speed at the water surface (ii) I don't even understand this statement. What does the "seal" part have to do with water flowing out? Won't the water still flow out through the tap until the tank is empty whether the reservoir is sealed or not? (iii) In my opinion, this statement would be correct. Increasing the gravitational potential energy of the...
Back
Top